Index: trunk/ICOSA_LMDZ/src/phyvenus/interface_icosa_lmdz.f90
===================================================================
--- trunk/ICOSA_LMDZ/src/phyvenus/interface_icosa_lmdz.f90	(revision 2556)
+++ trunk/ICOSA_LMDZ/src/phyvenus/interface_icosa_lmdz.f90	(revision 2591)
@@ -199,4 +199,5 @@
 
   ! For Cp(T)
+  LOGICAL :: cpofT
   REAL(rstd) :: nu_venus
   REAL(rstd) :: t0_venus
@@ -415,10 +416,15 @@
 ! 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)