check.f.html | |
Source file: check.f | |
Directory: /home/rjl/git/rjleveque/clawpack-4.x/amrclaw/2d/lib | |
Converted: Tue Jul 26 2011 at 12:59:08 using clawcode2html | |
This documentation file will not reflect any later changes in the source file. |
c c --------------------------------------------------------- c subroutine check(nsteps,time,nvar,naux) c c :::::::::::::::::::::: CHECK ::::::::::::::::::::::::::::::::; c check point routine - can only call at end of coarse grid cycle c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::; implicit double precision (a-h,o-z) integer tchkunit parameter (tchkunit = 12) character chkname*12 character tchkname*12 include "call.i" c c ### make the file name showing the time step c chkname = 'fort.chkxxxx' tchkname = 'fort.tckxxxx' nstp = nsteps do 20 ipos = 12, 9, -1 idigit = mod(nstp,10) chkname(ipos:ipos) = char(ichar('0') + idigit) tchkname(ipos:ipos) = char(ichar('0') + idigit) nstp = nstp / 10 20 continue open(unit=tchkunit,file=tchkname,status='unknown', . form='formatted') open(unit=chkunit,file=chkname,status='unknown', . form='unformatted') c c ### dump the data c write(tchkunit,*) 'Checkpoint file at time t = ',time write(tchkunit,*) 'alloc size memsize = ',memsize write(tchkunit,*) 'Number of steps taken = ',nsteps write(chkunit) lenmax,lendim,memsize write(chkunit) (alloc(i),i=1,lendim) write(chkunit) hxposs,hyposs,possk,icheck write(chkunit) lfree,lenf write(chkunit) rnode,node,lstart,newstl,listsp,tol, 1 ibuff,mstart,ndfree,lfine,iorder,mxnest, 2 intratx,intraty,kratio,iregsz,jregsz, 2 iregst,jregst,iregend,jregend, 3 kcheck,nsteps, 3 time,matlabu write(chkunit) evol,rvol,rvoll,lentot,tmass0,cflmax c close(chkunit) close(tchkunit) c return end