anchor.diff:       &         qold,maux,aux,t,dt)
         endif
   
  -c ---------------------------------------------------------------
  -c
  -c     # Assign flux registers fp_chombo, etc.  The reason we do this, rather
  -c     # then pass fp_chombo, etc directly into step2, is that fp and
  -c     # fp_chombo (fm and fm_chombo, etc.) are dimensioned differently.
  -c
  -
  -
  -      refratio = mx/mxc
  -      if (mcapa > 0) then
  -c        # with kappa
  -         do i = 1,mx+1
  -            ic = (i-1)/refratio + 1
  -            do j = 1,my
  -               jc = (j-1)/refratio + 1
  -               do m = 1,meqn
  -                  fp_chombo(i,j,m) = fp(i,j,m)/aux(i,j,mcapa)
  -                  fm_chombo(i,j,m) = fm(i,j,m)/aux(i-1,j,mcapa)
  -                  fpc_chombo(i,j,m) =
  -     &                  fp(i,j,m)/aux_coarse(ic,jc,mcapa)
  -                  fmc_chombo(i,j,m) =
  -     &                  fm(i,j,m)/aux_coarse(ic-1,jc,mcapa)
  -                  if (i .eq. 1) then
  -                     qadd_x(i,j,m) =
  -     &                     qadd_x(i,j,m)/aux_coarse(ic-1,jc,mcapa)
  -                  elseif (i .eq. mx+1) then
  -                     qadd_x(i,j,m) =
  -     &                     qadd_x(i,j,m)/aux_coarse(ic,jc,mcapa)
  -                  endif
  -               enddo
  -            enddo
  -         enddo
  -
  -         do i = 1,mx
  -            ic = (i-1)/refratio + 1
  -            do j = 1,my+1
  -               jc = (j-1)/refratio + 1
  -               do m = 1,meqn
  -                  gp_chombo(i,j,m) = gp(i,j,m)/aux(i,j,mcapa)
  -                  gm_chombo(i,j,m) = gm(i,j,m)/aux(i,j-1,mcapa)
  -                  gpc_chombo(i,j,m) = gp(i,j,m)/aux_coarse(ic,jc,mcapa)
  -                  gmc_chombo(i,j,m) =
  -     &                  gm(i,j,m)/aux_coarse(ic,jc-1,mcapa)
  -                  if (j .eq. 1) then
  -                     qadd_y(i,j,m) =
  -     &                     qadd_y(i,j,m)/aux_coarse(ic,jc-1,mcapa)
  -                  elseif (j .eq. my+1) then
  -                     qadd_y(i,j,m) =
  -     &                     qadd_y(i,j,m)/aux_coarse(ic,jc,mcapa)
  -                  endif
  -               enddo
  -            enddo
  -         enddo
  -      else
  -c        # no kappa
  -         do i = 1,mx+1
  -            do j = 1,my
  -               do m = 1,meqn
  -                  fp_chombo(i,j,m) = fp(i,j,m)
  -                  fm_chombo(i,j,m) = fm(i,j,m)
  -                  fpc_chombo(i,j,m) = fp(i,j,m)
  -                  fmc_chombo(i,j,m) = fm(i,j,m)
  -               enddo
  -            enddo
  -         enddo
  -
  -         do i = 1,mx
  -            do j = 1,my+1
  -               do m = 1,meqn
  -                  gp_chombo(i,j,m) = gp(i,j,m)
  -                  gm_chombo(i,j,m) = gm(i,j,m)
  -                  gpc_chombo(i,j,m) = gp(i,j,m)
  -                  gmc_chombo(i,j,m) = gm(i,j,m)
  -               enddo
  -            enddo
  -         enddo
  -      endif
  -
  -c ---------------------------------------------------------
  -
  -      if (debug) then
  -         write(iunit,*) 'qold...'
  -         do i = 1,mx
  -            do j = 1,my
  -               write(iunit,'(2I5,4E24.16)') i,j,qold(i,j,1),
  -     &               fm(i,j,1), fp(i,j,1)
  -            enddo
  -            write(iunit,*) ' '
  -         enddo
  -      endif
  -
  -      if (debug) then
  -         close(iunit)
  -      endif
  -
  -      if (debug) then
  -         write(6,*) 'In Clawpatch2; input a number to continue'
  -         read(5,*) some_number
  -      endif
  -
         end
   
   
  
positive.code:      &      i0auxr
        integer i,j,m, refratio, ic, jc
  
  c     # For debugging
        integer iunit, some_number
        logical debug
  
  c     Needed by Riemann solvers.
        double precision dtcom, dxcom,dycom,tcom
        integer icom, jcom
        common/comxyt/dtcom,dxcom,dycom,tcom,icom,jcom
  
  c -----------------------------------------------------------------
  c     # Fix auxtype - convert from integers to characters
  
negative.code:      &         'size of maxmaux'
           stop
        endif
  
  c ------------------------------------------------------------------
  c     # Set variables in common block, needed by step2.f and flux2.f
  