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

Timestamp:

Mar 30, 2015, 11:08:45 AM (10 years ago)

Author:

jyg

Message:

1/ Introduction of two variables in the ".def" files: (i) cvl_sig2feed is
the top of the convective feeding layer in sigma coordinates (D=0.97);
(ii) cvl_comp_threshold is the threshold fraction of convective points
below which compression occurs (D=1.).
2/ Corrections of various bugs revealed by the changes in compression:

• correct bugs in cv3a_uncompress.F90 for 3 fields used for convective

scavenging.

• add a reset to zero of "sig" and "w0" for non-convective points

(cva_driver.F90).

• in cv3_routines.F90, correct bounds of a few loops in cv3_undilute2,

correct the reset of the no-convection counter in cv3_yield.

• in phys_output_write_mod.F90, correct output of wdtrainA and wdtrainM.

3/ Improve declarations in various subroutines.

Modified files:

conema3.h
cv3param.h
cv3p1_closure.F90
conf_phys_m.F90
cv3a_compress.F90
phys_output_write_mod.F90
cv3_routines.F90
concvl.F90
cva_driver.F90
cv3a_uncompress.F90

File:

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

LMDZ5/trunk/libf/phylmd/conf_phys_m.F90

@@ -191 +191 @@
     REAL,SAVE :: ecrit_LES_omp
     REAL,SAVE :: ecrit_tra_omp
+    REAL,SAVE :: cvl_comp_threshold_omp
+    REAL,SAVE :: cvl_sig2feed_omp
     REAL,SAVE :: cvl_corr_omp
     LOGICAL,SAVE :: ok_lic_melt_omp
@@ -764 +766 @@
     ! KE
     !
+
+    !Config key  = cvl_comp_threshold
+    !Config Desc = maximum fraction of convective points enabling compression
+    !Config Def  = 1.00
+    !Config Help = fields are compressed when less than a fraction cvl_comp_threshold
+    !Config Help = of the points is convective.
+    cvl_comp_threshold_omp = 1.00
+    CALL getin('cvl_comp_threshold', cvl_comp_threshold_omp)
+
+    !Config key  = cvl_sig2feed
+    !Config Desc = sigma coordinate at top of feeding layer
+    !Config Def  = 0.97
+    !Config Help = deep convection is fed by the layer extending from the surface (pressure ps)
+    !Config Help = and cvl_sig2feed*ps.
+    cvl_sig2feed_omp = 0.97
+    CALL getin('cvl_sig2feed', cvl_sig2feed_omp)
 
     !Config key  = cvl_corr
@@ -2034 +2052 @@
     ecrit_tra = ecrit_tra_omp
     ecrit_reg = ecrit_reg_omp
+    cvl_comp_threshold = cvl_comp_threshold_omp
+    cvl_sig2feed = cvl_sig2feed_omp
     cvl_corr = cvl_corr_omp
     ok_lic_melt = ok_lic_melt_omp
@@ -2135 +2155 @@
     write(lunout,*)' RN2O_per = ',RN2O_per,' RCFC11_per = ', RCFC11_per
     write(lunout,*)' RCFC12_per = ',RCFC12_per
+    write(lunout,*)' cvl_comp_threshold=', cvl_comp_threshold
+    write(lunout,*)' cvl_sig2feed=', cvl_sig2feed
     write(lunout,*)' cvl_corr=', cvl_corr
     write(lunout,*)'ok_lic_melt=', ok_lic_melt