Changeset 836 for codes/icosagcm/devel/src/unstructured/data_unstructured.F90

Timestamp:

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

Author:

dubos

Message:

devel : Cp(T) thermodynamics (TBC)

File:

• codes/icosagcm/devel/src/unstructured/data_unstructured.F90 (modified) (4 diffs)

codes/icosagcm/devel/src/unstructured/data_unstructured.F90

@@ -1 +1 @@
 MODULE data_unstructured_mod
   USE ISO_C_BINDING
+  USE earth_const, ONLY : thermo_theta
   USE mpipara, ONLY : is_mpi_master
   USE grid_param, ONLY : llm, nqdyn
@@ -13 +14 @@
 
   INTEGER, PARAMETER :: eta_mass=1, eta_lag=2, &
-       thermo_theta=1, thermo_entropy=2, thermo_moist=3, thermo_boussinesq=4, &
+!       thermo_theta=1, thermo_entropy=2, thermo_moist=3, thermo_boussinesq=4, &
        caldyn_vert_cons=1, max_nb_stage=5
-  INDEX,  BIND(C) :: caldyn_thermo=thermo_theta, caldyn_eta=eta_lag, &
+  INDEX,  BIND(C) :: caldyn_eta=eta_lag, &
        caldyn_vert_variant=caldyn_vert_cons, nb_threads=0, nb_stage=0
+!  INDEX,  BIND(C) :: caldyn_thermo=thermo_theta, caldyn_eta=eta_lag, &
+!       caldyn_vert_variant=caldyn_vert_cons, nb_threads=0, nb_stage=0
   LOGICAL(C_BOOL), BIND(C) :: hydrostatic=.TRUE., debug_on=.FALSE.
   LOGICAL(C_BOOL), BIND(C, NAME='debug_hevi_solver') :: debug_hevi_solver_=.TRUE.
@@ -40 +43 @@
   TIME, PARAMETER :: print_trace_interval = 1.
   TIME, BIND(C) :: elapsed
-  NUM, BIND(C) :: g, ptop, cpp, cppv, Rd, Rv, preff, Treff, pbot, Phi_bot, rho_bot
-  NUM :: kappa
+  NUM, BIND(C) :: ptop, pbot, Phi_bot, rho_bot
+!  NUM, BIND(C) :: g, ptop, cpp, cppv, Rd, Rv, preff, Treff, pbot, Phi_bot, rho_bot
+!  NUM :: kappa
   NUM1(max_nb_stage), BIND(C)              :: tauj       ! diagonal of fast Butcher tableau
   NUM2(max_nb_stage,max_nb_stage), BIND(C) :: cslj, cflj ! slow and fast modified Butcher tableaus
@@ -252 +256 @@
   !
   SUBROUTINE init_params() BINDC(init_params)
+    USE earth_const
     kappa = Rd/cpp
     IF(is_mpi_master) THEN