diff --git a/src/DoRuRe.f90 b/src/DoRuRe.f90 index 831c669da497dca780aad0acbe41951bdba65442..6a94c36319a1ceadcc8fb6d4320a5a19a797b3e5 100644 --- a/src/DoRuRe.f90 +++ b/src/DoRuRe.f90 @@ -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