Skip to content
Snippets Groups Projects
Commit a48250fd authored by Dave Whipp's avatar Dave Whipp
Browse files

Removed unused declared variables; Added dilatation rate stats to list of output files

parent 0ebf39ee
No related branches found
No related tags found
No related merge requests found
......@@ -20,7 +20,6 @@ include 'mpif.h'
type(parameters) params
character(*) aktion
character*40 :: string
character*2 :: cis
integer iproc,ierr,is,us
......@@ -52,6 +51,7 @@ if (params%doDoRuRe .and. iproc.eq.0) then
open (828,file='./DORURE/temp_stats.dat',status='replace') ; close(828)
open (829,file='./DORURE/e3d_stats.dat',status='replace') ; close(829)
open (830,file='./DORURE/dt_stats.dat',status='replace') ; close(830)
open (831,file='./DORURE/dilatrtemp_stats.dat',status='replace') ; close(831)
do is=1,params%ns
call int_to_char(cis,2,is)
! us=600+is
......@@ -174,7 +174,7 @@ logical doDoRuRe
integer sizediag
double precision diag(sizediag)
integer ndof,istep,iter,iter_nl
integer ierr,iproc,nproc
integer ierr,iproc
call mpi_comm_rank (mpi_comm_world,iproc,ierr)
if (doDoRuRe .and. iproc.eq.0) then
......@@ -260,6 +260,14 @@ if (doDoRuRe .and. iproc.eq.0) then
write(805,'(3es13.4)') istep+dble(iter)*.01d0+dble(iter_nl)*.0001d0,minval(q),maxval(q)
close(805)
open (831,file ='./DORURE/dilatrtemp_stats.dat', &
position='append', &
form ='formatted', &
access ='sequential', &
status ='old')
write(831,'(3es13.4)') istep+dble(iter)*.01d0+dble(iter_nl)*.0001d0,minval(dilatrtemp),maxval(dilatrtemp)
close(831)
end if
end subroutine DoRuRe_leaf_stats
......@@ -335,7 +343,7 @@ implicit none
include 'mpif.h'
logical doDoRuRe
integer istep,iter,iter_nl,nl,n,nz,infog,nb
integer istep,iter,iter_nl,n,nz,infog,nb
double precision rinfog
real time1,time2,time3
integer iproc,ierr
......@@ -522,7 +530,7 @@ include 'mpif.h'
type (octreesolve) osolve
type (parameters) params
integer istep,iter,iter_nl
integer ierr,iproc,i,ie,k,il,nlev,err
integer ierr,iproc,i,k,err
integer,dimension(:),allocatable::levs
integer,dimension(:),allocatable::levelstat
double precision phi(8),vol,voltot
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment