Changeset 836 for codes/icosagcm/devel/src/base/earth_const.f90

Timestamp:

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

Author:

dubos

Message:

devel : Cp(T) thermodynamics (TBC)

File:

• codes/icosagcm/devel/src/base/earth_const.f90 (modified) (1 diff)

codes/icosagcm/devel/src/base/earth_const.f90

@@ -7 +7 @@
 ! hence the variables set here are not THREADPRIVATE
 
-  REAL(rstd) :: g=9.80616
-  REAL(rstd) :: radius=6.37122E6
+! Some variable are BIND(C) to be accessible from Python
+
+  REAL(rstd), BIND(C) :: g=9.80616    ! acceleration of gravity assumed constant
+  REAL(rstd):: radius=6.37122E6
   REAL(rstd) :: omega=7.292E-5
   REAL(rstd) :: scale_factor=1.
   REAL(rstd),PARAMETER :: daysec=86400
 
-  REAL(rstd) :: kappa=0.2857143
-  REAL(rstd) :: cpp=1004.70885
-  REAL(rstd) :: nu=0.35 ! exponent in variable-Cp law Cp=cpp*(T/Treff)^nu
-  REAL(rstd) :: cppv=1860.
-  REAL(rstd) :: Rv=461.5
-  REAL(rstd) :: Treff=273.
-  REAL(rstd) :: preff=101325.
-  REAL(rstd) :: pa=50000. ! default value set to preff/2 by disvert_std
+  REAL(rstd), BIND(C) ::   Treff=273.,     &   ! reference temperature, used in definition of entropy
+                           preff=101325.,  &   ! reference pressure, used in definition of entropy and theta
+                           cpp=1004.70885, &   ! specific heat at constant pressure at t=Treff for dry air
+                           cppv=1860.,     &   ! specific heat at constant pressure at t=Treff for water vapor
+                           Rd,             &   ! specific perfect gas constant for dry air
+                           Rv=461.5            ! specific perfect gas constant for water vapor
+
+  REAL(rstd) :: nu=0.35     ! exponent in variable-Cp law Cp=cpp*(T/Treff)^nu
+  REAL(rstd) :: pa=50000.   ! default value set to preff/2 by disvert_std
   REAL(rstd) :: scale_height=8000. ! atmospheric scale height (m)
   REAL(rstd) :: gas_constant = 8.3144621 
-  REAL(rstd) :: Rd, mu             ! specific perfect gas constant and molar mass (?)
+  REAL(rstd) :: mu             ! molar mass (?)
 
-  INTEGER, PUBLIC, PARAMETER :: thermo_none=-99, thermo_theta=1, thermo_entropy=2, thermo_variable_Cp=3, &
+  REAL(rstd) :: kappa=0.2857143
+
+  INTEGER, PARAMETER :: thermo_none=-99, thermo_theta=1, thermo_entropy=2, thermo_variable_Cp=3, &
        thermo_moist=4, thermo_boussinesq=5, thermo_dry=10, thermo_fake_moist=11, thermo_moist_debug=100
-  LOGICAL, PUBLIC :: boussinesq
-  INTEGER, PUBLIC :: caldyn_thermo, physics_thermo
+  LOGICAL :: boussinesq
+  INTEGER, BIND(C) :: caldyn_thermo, physics_thermo
 
 CONTAINS