Changeset 499 for codes/icosagcm/trunk/src/caldyn_kernels_hevi.f90

Timestamp:

11/27/16 00:24:42 (8 years ago)

Author:

dubos

Message:

Bugfix : NH - now passes DCMIP2012 (flat terrain), Bubble, DCMIP2016

File:

• codes/icosagcm/trunk/src/caldyn_kernels_hevi.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/caldyn_kernels_hevi.f90

@@ -363 +363 @@
        DO ij=ij_begin_ext,ij_end_ext
           pk(ij,l) = cpp*((pk(ij,l)/preff)**kappa) ! other formulae possible if exponentiation is slow
-          berni(ij) = (-.5*g*g)*(              &
+          berni(ij) = (-.25*g*g)*(              &
                  (W(ij,l)/m_il(ij,l))**2       &
                + (W(ij,l+1)/m_il(ij,l+1))**2 )
@@ -816 +816 @@
                     le_de(ij+u_ldown)*u(ij+u_ldown,l)**2 +    &
                     le_de(ij+u_rdown)*u(ij+u_rdown,l)**2 ) &
-               - .5*( gradPhi2(ij,l)  *w_il(ij,l)**2 +   &
+               - .25*( gradPhi2(ij,l)  *w_il(ij,l)**2 +   &
                       gradPhi2(ij,l+1)*w_il(ij,l+1)**2 )
        END DO