Changeset 837 for codes/icosagcm/devel/src/kernels_unst

Timestamp:

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

Author:

dubos

Message:

devel : Cp(T) thermodynamics (TBC)

File:

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

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

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