Changeset 417 for codes/icosagcm/trunk/src/omega.f90

Timestamp:

06/10/16 22:28:02 (8 years ago)

Author:

dubos

Message:

Fixed computation of omega

File:

• codes/icosagcm/trunk/src/omega.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/omega.f90

@@ -146 +146 @@
                   + le(ij+u_rdown)*u(ij+u_rdown,l)*( ne(ij,rdown)*ps(ij) + ne(ij+t_rdown,lup)*ps(ij+t_rdown) )
              ugradps = .5*(bp(l)+bp(l+1)) *ugradps/(-4*Ai(ij))  ! sign convention as in Ringler et al. 2010, Eq. 22 p.3072
-             w( ij, l) =  ugradps - .5*(convm( ij,l+1)+convm(ij,l)) 
+             w( ij, l) =  ugradps - g*.5*(convm( ij,l+1)+convm(ij,l)) 
           ENDDO
         ENDDO