Changeset 836 for codes/icosagcm/devel/Python

Timestamp:

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

Author:

dubos

Message:

devel : Cp(T) thermodynamics (TBC)

File:

• codes/icosagcm/devel/Python/src/kernels_caldyn_hydro.jin (modified) (1 diff)

codes/icosagcm/devel/Python/src/kernels_caldyn_hydro.jin

@@ -65 +65 @@
       gv = (g*Rd)*temp_ik/p_ik   ! specific volume v = Rd*T/p
     END_GEOPOT
+  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)
+    BALANCE( rhodz(THECELL) )
+    GEOPOT('temp_ik')
+      Cp_ik = nu*( theta(CELL,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
+    END_GEOPOT
   CASE(thermo_moist)
     BALANCE( rhodz(THECELL)*(1.+theta(THECELL,2)) )