Changeset 1894 for LMDZ5/trunk/libf/phylmd/pbl_surface_mod.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/pbl_surface_mod.F90 (modified) (4 diffs)

LMDZ5/trunk/libf/phylmd/pbl_surface_mod.F90

@@ -269 +269 @@
     INCLUDE "YOETHF.h"
     INCLUDE "temps.h"
-!****************************************************************************************
-! Declarations specifiques pour le 1D. A reprendre 
     INCLUDE "flux_arp.h"
-! Input variables
 !****************************************************************************************
     REAL,                         INTENT(IN)        :: dtime   ! time interval (s)
@@ -512 +509 @@
 
 !****************************************************************************************
+
 ! End of declarations
 !****************************************************************************************
@@ -562 +560 @@
 ! Force soil water content to qsol0 if qsol0>0 and VEGET=F (use bucket
 ! instead of ORCHIDEE)
-    IF (qsol0>0.) THEN
+    IF (qsol0>=0.) THEN
       PRINT*,'WARNING : On impose qsol=',qsol0
       qsol(:)=qsol0
@@ -983 +981 @@
 !
 !****************************************************************************************
+
+       if (evap0>=0.) then
+          yevap(:)=evap0
+          yevap(:)=RLVTT*evap0
+       endif
+
+
        y_d_ts(1:knon)   = ytsurf_new(1:knon) - yts(1:knon)