Changeset 1244 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/conf_phys.F90

Timestamp:

Sep 17, 2009, 4:59:40 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Le calcul de ratqs faisait osciller le modele et donnait une solution
asymptotique bidon et dépendant fortement du pas de temps.

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/conf_phys.F90 (modified) (6 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/conf_phys.F90

@@ -11 +11 @@
                        fact_cldcon, facttemps,ok_newmicro,iflag_radia,&
                        iflag_cldcon, &
-                       iflag_ratqs,ratqsbas,ratqshaut, &
+                       iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
                        ok_ade, ok_aie, aerosol_couple, &
                        flag_aerosol, new_aod, &
@@ -65 +65 @@
   LOGICAL              :: new_aod
   REAL                 :: bl95_b0, bl95_b1
-  real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut
+  real                 :: fact_cldcon, facttemps,ratqsbas,ratqshaut,tau_ratqs
   integer              :: iflag_cldcon
   integer              :: iflag_ratqs
@@ -81 +81 @@
   real,SAVE           :: fact_cldcon_omp, facttemps_omp,ratqsbas_omp
   real,SAVE           :: ratqshaut_omp
+  real,SAVE           :: tau_ratqs_omp
   integer,SAVE        :: iflag_radia_omp
   integer,SAVE        :: iflag_rrtm_omp
@@ -729 +730 @@
   ratqshaut_omp = 0.3
   call getin('ratqshaut',ratqshaut_omp)
+
+!Config Key  = tau_ratqs
+!Config Desc =  
+!Config Def  = 1800.
+!Config Help = 
+!
+  tau_ratqs_omp = 1800.
+  call getin('tau_ratqs',tau_ratqs_omp)
 
 !
@@ -1375 +1384 @@
     ratqsbas = ratqsbas_omp
     ratqshaut = ratqshaut_omp
+    tau_ratqs = tau_ratqs_omp
+
     iflag_radia = iflag_radia_omp
     iflag_rrtm = iflag_rrtm_omp
@@ -1500 +1511 @@
   write(numout,*)' ratqsbas = ',ratqsbas 
   write(numout,*)' ratqshaut = ',ratqshaut 
+  write(numout,*)' tau_ratqs = ',tau_ratqs 
   write(numout,*)' top_height = ',top_height 
   write(numout,*)' overlap = ',overlap