Changeset 1894 for LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

Timestamp:

Oct 22, 2013, 11:37:11 AM (11 years ago)

Author:

idelkadi

Message:

1) Patchs for idealized terra- and aqua-planets simulations.
ablsno.F90, imposing snow albedo at albsno0 if >= 0. (default val. : -1)
pbl_surface_mod.F90 : zero evaporation if either qsol0 (for terra-planets) or

evap0 (for aqua-planets) = 0.

2) New flag iflag_fisrtilp_qsat=0/1 for fisrtilp.F90, fisrtilp.h:

temperature dependency of qsat for iterative procedure
temporary flag to switch between the old (CMIP5, =0) and new versions (=1) versions :
zcond(i) = MAX(0.0,zqn(i)-zqs(i))/(1.+iflag_fisrtilp_qsat*zdqs(i))
=1 should be avoided because it creates significant supersaturation.
Both options are wrong (ongoing work ...)

clesphys.h, conf_phys_m.F90 : new parameters, iflag_fisrtilp_qsat, albsno0 and evap0

File:

• LMDZ5/trunk/libf/phylmd/conf_phys_m.F90 (modified) (8 diffs)

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -148 +148 @@
   REAL,SAVE :: inertie_sol_omp,inertie_sno_omp,inertie_ice_omp
   REAL,SAVE :: qsol0_omp
+  REAL,SAVE :: evap0_omp
+  REAL,SAVE :: albsno0_omp
   REAL      :: solarlong0
   REAL,SAVE :: solarlong0_omp
@@ -175 +177 @@
   INTEGER, SAVE :: nbapp_rad_omp, iflag_con_omp
   INTEGER, SAVE :: iflag_ener_conserv_omp
+  INTEGER, SAVE :: iflag_fisrtilp_qsat_omp
   LOGICAL,SAVE :: ok_strato_omp
   LOGICAL,SAVE :: ok_hines_omp
@@ -686 +689 @@
        CALL getin('iflag_ener_conserv',iflag_ener_conserv_omp)
 
+!Config  Key  = iflag_fisrtilp_qsat
+!Config  Desc = Flag de convection
+!Config  Def  = 1
+!Config  Help = Flag  pour la convection les options suivantes existent :
+!Config         -1 pour Kinetic energy correction
+!Config         1  conservation kinetic and enthalpy
+       iflag_fisrtilp_qsat_omp = 0
+       CALL getin('iflag_fisrtilp_qsat',iflag_fisrtilp_qsat_omp)
+
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -897 +909 @@
 !
 !-----------------------------------------------------------------------
-!  Valeur imposee de l'humidite du sol pour le modele bucket.
-!Config Key  = qsol0
-!Config Desc =  
-!Config Def  = -1.
-!Config Help = 
-!
+!  Valeur imposee pour configuration idealisees
+!Config Key  = qsol0 pour le bucket, evap0 pour aquaplanetes, albsno0
+! Default value -1 to activate the full computation
   qsol0_omp = -1.
   call getin('qsol0',qsol0_omp)
+  evap0_omp = -1.
+  call getin('evap0',evap0_omp)
+  albsno0_omp = -1.
+  call getin('albsno0',albsno0_omp)
 !
 !-----------------------------------------------------------------------
@@ -1680 +1693 @@
     iflag_con = iflag_con_omp
     iflag_ener_conserv = iflag_ener_conserv_omp
+    iflag_fisrtilp_qsat = iflag_fisrtilp_qsat_omp
 
     epmax = epmax_omp
@@ -1695 +1709 @@
     solarlong0 = solarlong0_omp
     qsol0 = qsol0_omp
+    evap0 = evap0_omp
+    albsno0 = albsno0_omp
     inertie_sol = inertie_sol_omp
     inertie_ice = inertie_ice_omp
@@ -1895 +1911 @@
   write(lunout,*)'iflag_con=',iflag_con
   write(lunout,*)'iflag_ener_conserv=',iflag_ener_conserv
+  write(lunout,*)'iflag_fisrtilp_qsat=',iflag_fisrtilp_qsat
   write(lunout,*)' epmax = ', epmax
   write(lunout,*)' ok_adj_ema = ', ok_adj_ema
@@ -1966 +1983 @@
   write(lunout,*)' solarlong0 = ', solarlong0
   write(lunout,*)' qsol0 = ', qsol0
+  write(lunout,*)' evap0 = ', evap0
+  write(lunout,*)' albsno0 = ', albsno0
   write(lunout,*)' inertie_sol = ', inertie_sol
   write(lunout,*)' inertie_ice = ', inertie_ice