anchor.diff:  !!             pR = auxl(pressure_index, i)
   !!         end if
   
  -         hvL = qr(mv,i-1) 
  -         hvR = ql(mv,i)        
  +         hvL = qr(i-1,mv) 
  +         hvR = ql(i,mv)        
   
           !!check for wet/dry boundary
           IF (hR .GT. dry_tolerance) THEN
  
positive.code:     integer ii_com, jj_com
      common /common_ii/ ii_com, jj_com
  
      use_simple = .true.
  
  !!  no pressure forcing
      pL = 0
  
negative.code:             qr(i-1,3) = 0.d0
          ENDIF
  
          IF (ql(i,1) .LT. 0.d0) THEN
              ql(i,1) = 0.d0
              ql(i,2) = 0.d0
              ql(i,3) = 0.d0
  