anchor.diff:  
            cflgrid = dmax1(cflgrid,cfl1d)
   c     
  -c     # update fluxes for use in AMR:
  -c     # NOTE : We update ghost cell values for subcycling
  +c        # update fluxes for use in AMR:
  +c        # NOTE : We update ghost cell values for subcycling
            do  m=1,meqn
               do  i=2-mbc,mx+mbc
                  fm(i,j,m) = fm(i,j,m) + faddm(i,m)
  
positive.code:             end do
           end do
        end do
  c     
  c     
        return
        end
  
negative.code:          endif
  c     
  c     
  c        # Store the value of j along this slice in the common block
  c        # comxyt in case it is needed in the Riemann solver (for
  c        # variable coefficient problems)
           jcom = j
  c     
  c        # compute modifications fadd and gadd to fluxes along this slice:
           ixy = 1
           call flux2(1,maxm,meqn,maux,mbc,mx,
       &        q1d,dtdx1d,aux1,aux2,aux3,
  