Changeset 401 for codes/icosagcm/trunk/src/earth_const.f90

Timestamp:

06/08/16 01:51:21 (8 years ago)

Author:

dubos

Message:

Introduced entropy as prognostic variable - tested with JW06

File:

• codes/icosagcm/trunk/src/earth_const.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/earth_const.f90

@@ -9 +9 @@
   REAL(rstd),SAVE :: kappa=0.2857143
   REAL(rstd),SAVE :: cpp=1004.70885
+  REAL(rstd),SAVE :: cppv=1004.70885 ! FIXME
+  REAL(rstd),SAVE :: Rv=461.5
+  REAL(rstd),SAVE :: Treff=273.
   REAL(rstd),SAVE :: preff=101325.
   REAL(rstd),SAVE :: pa=50000.
@@ -15 +18 @@
   REAL(rstd),SAVE :: gas_constant = 8.3144621 
   REAL(rstd),SAVE :: mu                 ! molar mass of the atmosphere
+
+  INTEGER, PARAMETER,PUBLIC :: thermo_theta=1, thermo_entropy=2, thermo_moist=3
+  INTEGER, PUBLIC :: caldyn_thermo
+  !$OMP THREADPRIVATE(caldyn_thermo) 
 
   LOGICAL, SAVE :: boussinesq, hydrostatic
@@ -33 +40 @@
     CALL getin("cpp",cpp)  
     CALL getin("preff",preff)  
+    CALL getin("Treff",Treff)  
     CALL getin("scale_height",scale_height)