Changeset 2591 for trunk/ICOSA_LMDZ

Timestamp:

Dec 6, 2021, 9:19:50 AM (3 years ago)

Author:

emillour

Message:

Venus GCM:
Add proper handling of "cpofT" case in the dynamico interface.
EM

File:

• trunk/ICOSA_LMDZ/src/phyvenus/interface_icosa_lmdz.f90 (modified) (2 diffs)

trunk/ICOSA_LMDZ/src/phyvenus/interface_icosa_lmdz.f90

@@ -199 +199 @@
 
   ! For Cp(T)
+  LOGICAL :: cpofT
   REAL(rstd) :: nu_venus
   REAL(rstd) :: t0_venus
@@ -415 +416 @@
 ! with constant Cp. 
 ! Conversion Theta/T takes this into account (see below)
-! IF CP(T)
-!    nu_venus=0.35
-!    t0_venus=460.
-! IF CP CONSTANT
-    nu_venus=0.
-    t0_venus=0.
+    cpofT=.false. ! default, assume Cp to be constant in physics
+    CALL getin('cpofT',cpofT)
+
+    IF (cpofT) THEN
+      nu_venus=0.35
+      t0_venus=460.
+    ELSE
+      nu_venus=0.
+      t0_venus=0.
+    ENDIF
+
     CALL init_cpdet_phy(cpp,nu_venus,t0_venus)