Changeset 837 for codes/icosagcm/devel/src/kernels_hex/caldyn_fast.k90

Timestamp:

05/03/19 13:40:15 (5 years ago)

Author:

dubos

Message:

devel : Cp(T) thermodynamics (TBC)

File:

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

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

@@ -27 +27 @@
          END DO
       END DO
+   CASE(thermo_variable_Cp)
+      ! thermodynamics with variable Cp
+      ! Cp(T) = Cp0 * (T/T0)^nu
+      ! => h = Cp(T).T/(nu+1)
+      DO l = ll_begin, ll_end
+         !DIR$ SIMD
+         DO ij=ij_begin, ij_end
+            berni(ij,l) = .5*(geopot(ij,l)+geopot(ij,l+1))
+            cp_ik = cpp*(pk(ij,l)/Treff)**nu
+            berni(ij,l) = berni(ij,l) + pk(ij,l)*(cp_ik/(nu+1.)-theta(ij,l,1)) ! Gibbs = h-Ts = T(Cp/(nu+1)-s)
+         END DO
+      END DO
    CASE DEFAULT
       PRINT *, 'Unsupported value of caldyn_thermo : ',caldyn_thermo ! FIXME