check.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	check (nsteps, time, nvar, naux)

Function/Subroutine Documentation

check()

subroutine check	(		nsteps,
			time,
			nvar,
			naux
	)

check point routine - can only call at end of coarse grid cycle

Definition at line 5 of file check.f.

References amr_module::alloc, amr_module::avenumgrids, amr_module::cflmax, amr_module::check_a, amr_module::checkpt_style, amr_module::chkunit, amr_module::dbugunit, amr_module::evol, amr_module::hxposs, amr_module::hyposs, amr_module::ibuff, amr_module::icheck, amr_module::intratx, amr_module::intraty, amr_module::iorder, amr_module::iregend, amr_module::iregridcount, amr_module::iregst, amr_module::iregsz, amr_module::jregend, amr_module::jregst, amr_module::jregsz, amr_module::kcheck, amr_module::kratio, amr_module::lendim, amr_module::lenf, amr_module::lenmax, amr_module::lentot, amr_module::lfine, amr_module::lfree, amr_module::listsp, amr_module::lstart, amr_module::matlabu, amr_module::memsize, amr_module::mstart, amr_module::mxnest, amr_module::ndfree, amr_module::ndfree_bnd, amr_module::newstl, amr_module::node, gauges_module::num_gauges, amr_module::numgrids, amr_module::outunit, amr_module::possk, gauges_module::print_gauges_and_reset_nextloc(), amr_module::rnode, amr_module::rvol, amr_module::rvoll, amr_module::timebound, amr_module::timeboundcpu, amr_module::timeregridding, amr_module::timeregriddingcpu, amr_module::timestepgrid, amr_module::timestepgridcpu, amr_module::timetick, amr_module::timetickcpu, amr_module::timevalout, amr_module::timevaloutcpu, amr_module::tmass0, amr_module::tol, amr_module::tvoll, and amr_module::tvollcpu.

Referenced by flagcheck(), and tick().

c
c :::::::::::::::::::::: CHECK ::::::::::::::::::::::::::::::::;
c :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::;

      use amr_module
      use gauges_module, only:  num_gauges
      use gauges_module, only: print_gauges_and_reset_nextloc

      implicit double precision (a-h,o-z)
      integer tchkunit
      parameter(tchkunit = 13)
      character  chkname*13
      character  tchkname*13

      write(6,601) time,nsteps
 601  format('Creating checkpoint file at t = ',e16.9,'  nsteps = ',i5)
c
      if (checkpt_style < 0) then

c         #  Alternate between two sets of files, overwriting the oldest
c         #  one, so that files do not accumulate with frequent checkpoints.
c
c         # Note that logical check_a is stored in amr_module, initialized
c         # in amr2 and perhaps reset properly in restrt.

          if (check_a) then
              chkname = 'fort.chkaaaaa'
              tchkname = 'fort.tckaaaaa'
            else
              chkname = 'fort.chkbbbbb'
              tchkname = 'fort.tckbbbbb'
            endif
          check_a = .not. check_a   ! to use other file next time

      else

c         # make a new checkpoint file for this time
c         ###  make the file name showing the time step
c
          chkname = 'fort.chkxxxxx'
          tchkname = 'fort.tckxxxxx'
          nstp = nsteps
          do 20 ipos = 13, 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
      endif

      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(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,ndfree_bnd,lfine,iorder,mxnest,
     2          intratx,intraty,kratio,iregsz,jregsz,
     2          iregst,jregst,iregend,jregend, 
     3          numgrids,kcheck,nsteps,
     3          time,matlabu
      write(chkunit) avenumgrids, iregridcount,
     1               evol,rvol,rvoll,lentot,tmass0,cflmax,
     2               tvoll,tvollcpu,timetick,timetickcpu,
     3               timestepgrid,timestepgridcpu,
     4               timebound,timeboundcpu,
     5               timeregridding,timeregriddingcpu,
     6               timevalout,timevaloutcpu 
c
      close(chkunit)

c     # flush open running output files fort.amr, fort.gauge, fort.debug
c     # so if code dies it will at least have output up to this checkpoint time

      flush(outunit)        ! defined in amr_module.f90
      flush(dbugunit)       ! defined in amr_module.f90
c      flush(OUTGAUGEUNIT)   ! defined in gauges_module.f90
      do ii = 1, num_gauges
         call print_gauges_and_reset_nextloc(ii)
      end do

c     # write the time stamp file last so it's not updated until data is
c     # all dumped, in case of crash mid-dump.
      write(tchkunit,*) 'Checkpoint file at time t = ',time
      write(tchkunit,*) 'alloc size memsize = ',memsize
      write(tchkunit,*) 'Number of steps taken = ',nsteps
      close(tchkunit)
c
      return

Here is the call graph for this function:

Here is the caller graph for this function: