Fix src/Nd/restrt.f so that frame number is not advanced if output_t0…

… is True

In this case it re-writes the last frame from the original calculation
when restarting, and it should reuse that frame number.  For historical
reasons the frame number is called matlabu, stored in amr_module.f90

src/1d/restrt.f

@@ -11,6 +11,7 @@ subroutine restrt(nsteps,time,nvar,naux)
 
       logical foundFile
       dimension intrtx(maxlv),intrtt(maxlv)
+      integer matlabu0
 c
 c :::::::::::::::::::::::::::: RESTRT ::::::::::::::::::::::::::::::::
 c read back in the check point files written by subr. check.
@@ -39,6 +40,8 @@ subroutine restrt(nsteps,time,nvar,naux)
       open(rstunit,file=trim(rstfile),status='old',form='unformatted')
       rewind rstunit
 
+      matlabu0 = matlabu ! value set by amr2, 0 if output_t0, else 1
+
       read(rstunit) lenmax,lendim,isize
 
 c     # need to allocate for dynamic memory:
@@ -58,6 +61,8 @@ subroutine restrt(nsteps,time,nvar,naux)
 
       close(rstunit) 
 
+      matlabu = matlabu - 1 + matlabu0 ! redo previous frame if output_t0
+
       write(outunit,100) nsteps,time
       write(6,100) nsteps,time
  100  format(/,' RESTARTING the calculation after ',i5,' steps',

src/2d/restrt.f

@@ -11,6 +11,7 @@ subroutine restrt(nsteps,time,nvar,naux)
 
       logical foundFile
       dimension intrtx(maxlv),intrty(maxlv),intrtt(maxlv)
+      integer matlabu0
 c
 c :::::::::::::::::::::::::::: RESTRT ::::::::::::::::::::::::::::::::
 !> read back in the check point files written by subr. check.
@@ -39,6 +40,8 @@ subroutine restrt(nsteps,time,nvar,naux)
       open(rstunit,file=trim(rstfile),status='old',form='unformatted')
       rewind rstunit
 
+      matlabu0 = matlabu ! value set by amr2, 0 if output_t0, else 1
+
       !!read(rstunit) lenmax,lendim,isize
       !! new version has flexible node size, so need to read current size maxgr
       read(rstunit) lenmax,lendim,isize,maxgr
@@ -66,6 +69,8 @@ subroutine restrt(nsteps,time,nvar,naux)
 
       close(rstunit) 
 
+      matlabu = matlabu - 1 + matlabu0 ! redo previous frame if output_t0
+
       write(outunit,100) nsteps,time
       write(6,100) nsteps,time
  100  format(/,' RESTARTING the calculation after ',i5,' steps',

src/3d/restrt.f

@@ -10,6 +10,8 @@ subroutine restrt(nsteps,time,nvar,naux)
 
       logical foundFile
       dimension intrtx(maxlv),intrty(maxlv),intrtz(maxlv),intrtt(maxlv)
+      integer matlabu0
+
 c
 c :::::::::::::::::::::::::::: RESTRT ::::::::::::::::::::::::::::::::
 c read back in the check point files written by subr. check.
@@ -39,6 +41,9 @@ subroutine restrt(nsteps,time,nvar,naux)
       open(rstunit,file=trim(rstfile),status='old',form='unformatted')
       rewind rstunit
 
+      matlabu0 = matlabu ! value set in amr3, 0 if output_t0, else 1
+
+      !!read(rstunit) lenmax,lendim,isize
       !! new version has flexible node size, so need to read current size maxgr
       read(rstunit) lenmax,lendim,isize,maxgr
 
@@ -60,6 +65,8 @@ subroutine restrt(nsteps,time,nvar,naux)
 
       close(rstunit) 
 
+      matlabu = matlabu - 1 + matlabu0 ! redo previous frame if output_t0
+
       write(outunit,100) nsteps,time
       write(6,100) nsteps,time
  100  format(/,' RESTARTING the calculation after ',i5,' steps',