Changeset 963 for LMDZ4/trunk/libf/phylmd/conf_phys.F90

Timestamp:

Jun 9, 2008, 11:54:29 AM (17 years ago)

Author:

Laurent Fairhead

Message:

On passe les parametres lus dans parameter_mix.data et supcrit.data dans le
physiq.def
LF

File:

• LMDZ4/trunk/libf/phylmd/conf_phys.F90 (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/conf_phys.F90

@@ -22 +22 @@
  include "nuage.h"
  include "YOMCST.h"
+ include "YOMCST2.h"
 !IM : on inclut/initialise les taux de CH4, N2O, CFC11 et CFC12
 include "clesphys.h"
@@ -83 +84 @@
   integer,SAVE :: iflag_coupl_omp,iflag_clos_omp,iflag_wake_omp
   integer,SAVE :: iflag_cvl_sigd_omp
+  REAL, SAVE :: supcrit1_omp, supcrit2_omp
+  INTEGER, SAVE :: iflag_mix_omp
+  real, save :: scut_omp, qqa1_omp, qqa2_omp, gammas_omp, Fmax_omp, alphas_omp
 
   REAL,SAVE :: R_ecc_omp,R_peri_omp,R_incl_omp,solaire_omp,co2_ppm_omp
@@ -991 +995 @@
   ok_slab_sicOBS_omp = .true.
   call getin('ok_slab_sicOBS', ok_slab_sicOBS_omp)
+!
+!
+! PARAMETERS FOR CONVECTIVE INHIBITION BY TROPOS. DRYNESS
+!
+!Config Key  = supcrit1
+!Config Desc =
+!Config Def  = .540
+!Config Help =
+!
+  supcrit1_omp = .540
+  call getin('supcrit1',supcrit1_omp)
+
+!
+!Config Key  = supcrit2
+!Config Desc =
+!Config Def  = .600
+!Config Help =
+!
+  supcrit2_omp = .600
+  call getin('supcrit2',supcrit2_omp)
+
+!
+! PARAMETERS FOR THE MIXING DISTRIBUTION
+!
+!
+!Config Key  = iflag_mix
+!Config Desc =
+!Config Def  = 1
+!Config Help =
+!
+  iflag_mix_omp = 1
+  call getin('iflag_mix',iflag_mix_omp)
+
+!
+!Config Key  = scut
+!Config Desc =
+!Config Def  = 0.95
+!Config Help =
+!
+  scut_omp = 0.95
+  call getin('scut',scut_omp)
+
+!
+!Config Key  = qqa1
+!Config Desc =
+!Config Def  = 1.0
+!Config Help =
+!
+  qqa1_omp = 1.0
+  call getin('qqa1',qqa1_omp)
+
+!
+!Config Key  = qqa2
+!Config Desc =
+!Config Def  = 0.0
+!Config Help =
+!
+  qqa2_omp = 0.0
+  call getin('qqa2',qqa2_omp)
+
+!
+!Config Key  = gammas
+!Config Desc =
+!Config Def  = 0.05
+!Config Help =
+!
+  gammas_omp = 0.05
+  call getin('gammas',gammas_omp)
+
+!
+!Config Key  = Fmax
+!Config Desc =
+!Config Def  = 0.65
+!Config Help =
+!
+  Fmax_omp = 0.65
+  call getin('Fmax',Fmax_omp)
+
+!
+!Config Key  = alphas  
+!Config Desc =
+!Config Def  = -5.
+!Config Help =
+!
+  alphas_omp = -5.
+  call getin('alphas',alphas_omp)
+
 
 
@@ -1097 +1188 @@
     cvl_corr = cvl_corr_omp
     ok_lic_melt = ok_lic_melt_omp
+    supcrit1 = supcrit1_omp
+    supcrit2 = supcrit2_omp
+    iflag_mix = iflag_mix_omp
+    scut = scut_omp
+    qqa1 = qqa1_omp
+    qqa2 = qqa2_omp
+    gammas = gammas_omp
+    Fmax = Fmax_omp
+    alphas = alphas_omp
 
 !$OMP MASTER
@@ -1178 +1278 @@
   write(numout,*)' inertie_ice = ', inertie_ice
   write(numout,*)' inertie_sno = ', inertie_sno
+  write(numout,*)' supcrit1 = ', supcrit1
+  write(numout,*)' supcrit2 = ', supcrit2
+  write(numout,*)' iflag_mix = ', iflag_mix
+  write(numout,*)' scut = ', scut
+  write(numout,*)' qqa1 = ', qqa1
+  write(numout,*)' qqa2 = ', qqa2
+  write(numout,*)' gammas = ', gammas
+  write(numout,*)' Fmax = ', Fmax
+  write(numout,*)' alphas = ', alphas
 
   write(numout,*)' lonmin lonmax latmin latmax bilKP_ins =',&