anchor.diff:                       uin(i,j,k,UMZ) * grav(i,j,k,3)
   
                  ! rotation
  -               SrE = SrE + omega**2*(uin(i,j,k,UMX)*x + uin(i,j,k,UMY)*y)
  +               !SrE = SrE + omega**2*(uin(i,j,k,UMX)*x + uin(i,j,k,UMY)*y)
   
                  uout(i,j,k,UEDEN) = uout(i,j,k,UEDEN) + SrE * dt
   
  
positive.code: 
        problo(:) = 0.0d0
  
        do k = lo(3),hi(3)
           z = problo(3) + dx(3)*(float(k)+0.5d0) - center(3)
  
  
negative.code:                SrW = rho * grav(i,j,k,3)
  
                 ! rotation
                 !SrU = SrU + 2.d0*uin(i,j,k,UMY)*omega + rho*omega**2*x
                 !SrV = SrV - 2.d0*uin(i,j,k,UMX)*omega + rho*omega**2*y
  
                 ! update the momenta
                 uout(i,j,k,UMX)   = uout(i,j,k,UMX) + SrU * dt
  