anchor.diff:  
                    else
   
  -                   vx = adv(domlo(1),j,k,UMX)/adv(domlo(1),j,k,URHO)
  -                   vy = adv(domlo(1),j,k,UMY)/adv(domlo(1),j,k,URHO)
  -                   vz = adv(domlo(1),j,k,UMZ)/adv(domlo(1),j,k,URHO)
  +                   vx = adv(domhi(1),j,k,UMX)/adv(domhi(1),j,k,URHO)
  +                   vy = adv(domhi(1),j,k,UMY)/adv(domhi(1),j,k,URHO)
  +                   vz = adv(domhi(1),j,k,UMZ)/adv(domhi(1),j,k,URHO)
    
                    endif
   
  
positive.code: 
                   else
  
                     vx = adv(domlo(1),j,k,UMX)/adv(domlo(1),j,k,URHO)
                     vy = adv(domlo(1),j,k,UMY)/adv(domlo(1),j,k,URHO)
                     vz = adv(domlo(1),j,k,UMZ)/adv(domlo(1),j,k,URHO)
   
                   endif
  
  
negative.code: 
                   else
  
                     vx = adv(i,domlo(2),k,UMX)/adv(i,domlo(2),k,URHO)
                     vy = adv(i,domlo(2),k,UMY)/adv(i,domlo(2),k,URHO)
                     vz = adv(i,domlo(2),k,UMZ)/adv(i,domlo(2),k,URHO)
   
                   endif
  
  