anchor.diff:                 gupdate = 0.5d0*dtdx1d(i) * bpasdq(m,i)
                  gaddm(m,i,2) = gaddm(m,i,2) - gupdate
                  gaddp(m,i,2) = gaddp(m,i,2) - gupdate
  -  180          continue
  +         end do
  +      end do
   c
     999 continue
         return
  
positive.code:                gupdate = 0.5d0*dtdx1d(i-1) * bpasdq(m,i)
                 gaddm(m,i-1,2) = gaddm(m,i-1,2) - gupdate
                 gaddp(m,i-1,2) = gaddp(m,i-1,2) - gupdate
    160          continue
  c
  c     # split the right-going flux difference into down-going and up-going:
        call rpt2(ixy,2,maxm,meqn,mwaves,maux,mbc,mx,
  
negative.code:       if (limit) call clawpack5_inlinelimiter(maxm,meqn,mwaves,mbc,mx,
       &      wave,s,mthlim)
  c
        do  i = 1, mx+1
  c
  c        # For correction terms below, need average of dtdx in cell
  c        # i-1 and i.  Compute these and overwrite dtdx1d:
  