anchor.diff:          END DO
   
   !!        !set normal direction
  -         if (ixy.eq.1) then
  +         if (ixy .eq. 1) then
               mu=2
               mv=3
            else
  
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
  