anchor.diff:           open(iunit,file='clawpatch2.out')
            write(iunit,'(A,I5)') 'mx     = ', mx
            write(iunit,'(A,I5)') 'my     = ', my
  -         write(iunit,'(A,I5)') 'mxc    = ', mxc
  -         write(iunit,'(A,I5)') 'myc    = ', myc
            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 -----------------------------------------------------------------------
  -c     # Set up boundary conditions
  -      do i = 1,4
  -         if (intersectsBoundary(i) .eq. 0) then
  -            mthbc_amr(i) = -1
  -         else
  -            mthbc_amr(i) = mthbc(i)
  -         endif
  -      enddo
  -
  -      if (debug) then
  -c         write(iunit,'(A,4I4)') 'mthbc_tmp: ',(mthbc_amr(i),i=1,4)
  -      endif
  -
  -      call bc2(maxmx, maxmy, meqn, mbc, mx, my, xlower, ylower,
  -     &      dx, dy, qold, maux, aux, t, dt, mthbc_amr)
  -
  -c ----------------------------------------------------------------------
  -c     # Before qold is modified, handle non-conservative fix-up at
  -c     # coarse/fine boundaries by modifying qadd_x, y, z.
  -
  -      if (level .gt. 0) then
  -         i0wave  = 1
  -         i0s     = i0wave +  (maxm+2*mbc)*meqn*mwaves
  -         i0amdq  = i0s +     (maxm+2*mbc)*mwaves
  -         i0apdq  = i0amdq +  (maxm+2*mbc)*meqn
  -         i0ql    = i0apdq +  (maxm+2*mbc)*meqn
  -         i0qr    = i0ql +    (maxm+2*mbc)*meqn
  -         i0auxl  = i0qr +    (maxm+2*mbc)*meqn
  -         i0auxr  = i0auxl +  (maxm+2*mbc)*maux
  -         mused   = i0auxr +  (maxm+2*mbc)*maux - 1
  -c
  -         if (mused .gt. mwork) then
  -            write(6,*) 'Not enough work space in ClawPatch2'
  -            write(6,*) 'mused = ', mused, '   mwork =',mwork
  -            stop
  -         endif
  -
  -         call qad2(maxm,meqn,mwaves,mbc,mx,my,maux,qold,qold_coarse,
  -     &         aux, aux_coarse,
  -     &         intersectsBoundary, mxc,myc,dt,dx,dy,
  -     &         work(i0ql),work(i0qr), work(i0auxl),work(i0auxr),
  -     &         work(i0wave),work(i0s),work(i0amdq),work(i0apdq),
  -     &         qadd_x,qadd_y, auxtype)
  -
         endif
   
   c -----------------------------------------------------------------------
  
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:      &         qold,maux,aux,t,dt)
        endif
  
        end
  
  
  