Changeset 467 for LMDZ.3.3/branches/rel-LF/libf/phylmd/conf_phys.F90

Timestamp:

Aug 6, 2003, 4:50:49 PM (21 years ago)

Author:

lmdzadmin

Message:

Modifs sur les seuils (cdrag etc...), inclusion des diagnostics ISCCP par Ionela
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/conf_phys.F90 (modified) (7 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/conf_phys.F90

@@ -168 +168 @@
 !
 !valeur AMIP II
-  RCH4 = 1.65E-06* 16.043/28.97
+!OK  RCH4 = 1.65E-06* 16.043/28.97
 ! RCH4 = 9.137366240938903E-07
 !
 !ancienne valeur
 ! RCH4 = 1.72E-06* 16.043/28.97
-  call getin('RCH4', RCH4)
+!OK call getin('RCH4', RCH4)
+  CH4_ppb = 1650.
+  call getin('CH4_ppb', CH4_ppb)
+  RCH4 = CH4_ppb * 1.0E-09 * 16.043/28.97
 !!
 !Config Key  = RN2O
@@ -182 +185 @@
 !
 !valeur AMIP II
-  RN2O = 306.E-09* 44.013/28.97
+!OK  RN2O = 306.E-09* 44.013/28.97
 ! RN2O = 4.648939592682085E-07
 !
 !ancienne valeur
 ! RN2O = 310.E-09* 44.013/28.97
-  call getin('RN2O', RN2O)
+!OK  call getin('RN2O', RN2O)
+  N2O_ppb=306.
+  call getin('N2O_ppb', N2O_ppb)
+  RN2O = N2O_ppb * 1.0E-09 * 44.013/28.97
 !!
 !Config Key  = RCFC11
@@ -195 +201 @@
 !               
 !
-  RCFC11 = 280.E-12* 137.3686/28.97
+!OK RCFC11 = 280.E-12* 137.3686/28.97
+  CFC11_ppt = 280.
+  call getin('CFC11_ppt',CFC11_ppt)
+  RCFC11=CFC11_ppt* 1.0E-12 * 137.3686/28.97
 ! RCFC11 = 1.327690990680013E-09
-  call getin('RCFC11', RCFC11)
+!OK call getin('RCFC11', RCFC11)
 !!
 !Config Key  = RCFC12
@@ -205 +214 @@
 !               
 !
-  RCFC12 = 484.E-12* 120.9140/28.97
+!OK RCFC12 = 484.E-12* 120.9140/28.97
+  CFC12_ppt = 484.
+  call getin('CFC12_ppt',CFC12_ppt)
+  RCFC12 = CFC12_ppt * 1.0E-12 * 120.9140/28.97
 ! RCFC12 = 2.020102726958923E-09
-  call getin('RCFC12', RCFC12)
+!OK call getin('RCFC12', RCFC12)
 !!
 !! Constante solaire & Parametres orbitaux & taux gaz effet de serre END
@@ -386 +398 @@
   rad_chau2 = 9.0
   call getin('rad_chau2',rad_chau2)
+
+!
+!Config Key  = top_height
+!Config Desc =
+!Config Def  = 3
+!Config Help =
+!
+  top_height = 3
+  call getin('top_height',top_height)
+
+!
+!Config Key  = overlap
+!Config Desc =
+!Config Def  = 3
+!Config Help =
+!
+  overlap = 3
+  call getin('overlap',overlap)
+
+
+!IM
+!
+!Config Key  = cdmmax
+!Config Desc =
+!Config Def  = 1.3E-3
+!Config Help =
+!
+  cdmmax = 1.3E-3
+  call getin('cdmmax',cdmmax)
+
+!
+!Config Key  = cdhmax
+!Config Desc =
+!Config Def  = 1.1E-3
+!Config Help =
+!
+  cdhmax = 1.1E-3
+  call getin('cdhmax',cdhmax)
 
 !
@@ -415 +465 @@
   write(numout,*)' co2_ppm =',co2_ppm
   write(numout,*)' RCO2 = ',RCO2
-  write(numout,*)' RCH4 = ',RCH4
-  write(numout,*)' RN2O =  ',RN2O
-  write(numout,*)' RCFC11 =  ',RCFC11
-  write(numout,*)' RCFC12 =  ',RCFC12
+  write(numout,*)' CH4_ppb =',CH4_ppb,' RCH4 = ',RCH4
+  write(numout,*)' N2O_ppb =',N2O_ppb,' RN2O =  ',RN2O
+  write(numout,*)' CFC11_ppt=',CFC11_ppt,' RCFC11 =  ',RCFC11
+  write(numout,*)' CFC12_ppt=',CFC12_ppt,' RCFC12 =  ',RCFC12
 !IM constantes physiques END
   write(numout,*)' epmax = ', epmax
@@ -438 +488 @@
   write(numout,*)' ratqsbas = ',ratqsbas 
   write(numout,*)' ratqshaut = ',ratqshaut 
+  write(numout,*)' top_height = ',top_height 
+  write(numout,*)' overlap = ',overlap 
 
   return