Changeset 836 for codes/icosagcm/devel/src/kernels_hex

Timestamp:

05/03/19 13:25:21 (5 years ago)

Author:

dubos

Message:

devel : Cp(T) thermodynamics (TBC)

File:

• codes/icosagcm/devel/src/kernels_hex/compute_geopot.k90 (modified) (1 diff)

codes/icosagcm/devel/src/kernels_hex/compute_geopot.k90

@@ -82 +82 @@
          END DO
       END DO
+   CASE(thermo_variable_Cp)
+      ! thermodynamics with variable Cp
+      ! Cp.dT = dh = Tds + vdp
+      ! pv = RT
+      ! => ds = (dh+v.dp)/T = Cp.dT/T - R dp/p
+      ! Cp(T) = Cp0 * (T/T0)^nu
+      ! => s(p,T) = Cp(T)/nu - R log(p/preff)
+      ! h = Cp(T).T/(nu+1)
+      !DIR$ SIMD
+      DO ij=ij_omp_begin_ext,ij_omp_end_ext
+         pk(ij,llm) = ptop + .5*g* rhodz(ij,llm)
+      END DO
+      DO l = llm-1,1,-1
+         !DIR$ SIMD
+         DO ij=ij_omp_begin_ext,ij_omp_end_ext
+            pk(ij,l) = pk(ij,l+1) + (.5*g)*( rhodz(ij,l) + rhodz(ij,l+1) )
+         END DO
+      END DO
+      IF(caldyn_eta == eta_lag) THEN
+         !DIR$ SIMD
+         DO ij=ij_omp_begin_ext,ij_omp_end_ext
+            ps(ij) = pk(ij,1) + .5*g* rhodz(ij,1)
+         END DO
+      END IF
+      DO l = 1,llm
+         !DIR$ SIMD
+         DO ij=ij_omp_begin_ext,ij_omp_end_ext
+            p_ik = pk(ij,l)
+            Cp_ik = nu*( theta(ij,l,1) + Rd*log(p_ik/preff) )
+            temp_ik = Treff* (Cp_ik/cpp)**(1./nu)
+            gv = (g*Rd)*temp_ik/p_ik ! specific volume v = Rd*T/p
+            pk(ij,l) = temp_ik
+            geopot(ij,l+1) = geopot(ij,l) + gv*rhodz(ij,l)
+         END DO
+      END DO
    CASE(thermo_moist)
       !DIR$ SIMD