anchor.diff:       &         'size of maxmaux'
            stop
         endif
  -
  -      do m = 1,maux
  -         if (auxtype_int(m) .eq. 1) then
  -            auxtype(m) = "xleft"
  -         elseif (auxtype_int(m) .eq. 2) then
  -            auxtype(m) = "yleft"
  -         elseif (auxtype_int(m) .eq. 3) then
  -            auxtype(m) = "zleft"
  -         elseif (auxtype_int(m) .eq. 4) then
  -            auxtype(m) = "center"
  -         elseif (auxtype_int(m) .eq. 5) then
  -            auxtype(m) = "capacity"
  -         endif
  -      enddo
  -
   
   c ------------------------------------------------------------------
   c     # Set variables in common block, needed by step2.f and flux2.f
  
positive.code:       subroutine ClawPatch2(maxm, meqn, maux, mbc, method, mthlim,
       &      mcapa, mwaves, mx, my, qold, aux, dx, dy, dt, cfl,
       &      work, mwork, qold_coarse, aux_coarse, qadd_x, qadd_y,
       &      auxtype_int, xlower, ylower,
       &      intersectsBoundary, level,mthbc,t, mxc, myc, fp,fm, gp, gm,
       &      fp_chombo, fm_chombo, gp_chombo,gm_chombo,
       &      fpc_chombo, fmc_chombo, gpc_chombo,gmc_chombo)
  
        implicit none
  
        external rpn2,rpt2
  
        integer maxm,meqn,maux,mbc,mcapa,mwaves,mx,my, mwork,
       &      maxmx, maxmy, mxc, myc, intersectsBoundary(4),
       &      mthbc(4), mthbc_amr(4),level
        double precision dx,dy,dt,cfl, xlower, ylower, t,
       &      dtdx, dtdy
  
        integer method(7), mthlim(mwaves)
        integer auxtype_int(maux)
        integer maxmaux
        parameter(maxmaux = 100)
        character*10 auxtype(maxmaux)
        double precision work(mwork)
  
        double precision   qold(1-mbc:mx+mbc, 1-mbc:my+mbc, meqn)
        double precision   qold_coarse(0:mxc+1, 0:myc+1, meqn)
  
        double precision  aux(1-mbc:mx+mbc, 1-mbc:my+mbc, maux)
        double precision  aux_coarse(0:mxc+1,0:myc+1,maux)
  
  c     Conservative numerical fluxes returned here.
        double precision fp(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
        double precision fm(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
        double precision gp(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
        double precision gm(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
        double precision fp_chombo(mx+1,my,  meqn)
        double precision fm_chombo(mx+1,my,  meqn)
        double precision gp_chombo(mx,  my+1,meqn)
        double precision gm_chombo(mx,  my+1,meqn)
        double precision fpc_chombo(mx+1,my,  meqn)
        double precision fmc_chombo(mx+1,my,  meqn)
        double precision gpc_chombo(mx,  my+1,meqn)
        double precision gmc_chombo(mx,  my+1,meqn)
        double precision qadd_x(mx+1,my,  meqn)
        double precision qadd_y(mx,  my+1,meqn)
  
  c     # Local variables
        integer i0faddm, i0faddp, i0gaddm, i0gaddp,
  
negative.code:          open(iunit,file='clawpatch2.out')
           write(iunit,'(A,I5)') 'mx     = ', mx
           write(iunit,'(A,I5)') 'my     = ', my
           write(iunit,'(A,I5)') 'meqn   = ', meqn
           write(iunit,'(A,I5)') 'mwaves = ', mwaves
           write(iunit,'(A,I5)') 'mbc    = ', mbc
           write(iunit,'(A,I5)') 'maux   = ', maux
           write(iunit,'(A,I5)') 'mcapa  = ', mcapa
        endif
  
  c -----------------------------------------------------------------------
  