Changeset 658 for codes/icosagcm/devel/src/kernels_unst/compute_NH_geopot.k90

Timestamp:

12/30/17 02:00:38 (7 years ago)

Author:

dubos

Message:

devel/unstructured : updated kernels

File:

• codes/icosagcm/devel/src/kernels_unst/compute_NH_geopot.k90 (modified) (2 diffs)

codes/icosagcm/devel/src/kernels_unst/compute_NH_geopot.k90

@@ -4 +4 @@
    g2=g*g
    gm2 = 1./g2
+   vreff = Treff*cpp/preff*kappa
    gamma = 1./(1.-kappa)
    !$OMP BARRIER
@@ -35 +36 @@
             DO l = 1,llm
                rho_ij = (g*m_ik(l,ij))/(Phi_il(l+1,ij)-Phi_il(l,ij))
-               X_ij = Treff*exp(theta(l,ij)/cpp) ! theta = Tref.exp(s/Cp)
-               X_ij = (cpp/preff)*kappa*X_ij*rho_ij
-               p_ik(l,ij) = preff*(X_ij**gamma)
+               X_ij = log(vreff*rho_ij) + theta(l,ij)/cpp
+               p_ik(l,ij) = preff*exp(X_ij*gamma)
                c2_mik = gamma*p_ik(l,ij)/(rho_ij*m_ik(l,ij)) ! c^2 = gamma*R*T = gamma*p/rho
                A_ik(l,ij) = c2_mik*(tau/g*rho_ij)**2