Changeset 1525 for LMDZ5/trunk/libf/phylmd/newmicro.F

Timestamp:

May 25, 2011, 12:55:27 PM (13 years ago)

Author:

idelkadi

Message:

Modifications concerning the cloud scheme:

1. In newmicro, it now possible to read a min and max effective radius of ice particles from physiq.def : rei_min and rei_max which were initially set to 3.5 and 61.29 microns in the code.

concerns : conf_phys.F90, nuage.h, newmicro.F

2. In physiq.F, in case of combination of iflag_cldcon>=5 (A. Jam cloud scheme) and iflag_coupl=5

concerns : physiq.F and thermcell_main.F90

File:

• LMDZ5/trunk/libf/phylmd/newmicro.F (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/newmicro.F

@@ -155 +155 @@
       REAL diff_paprs(klon, klev), zfice1, zfice2(klon, klev)
       REAL rad_chaud_tab(klon, klev), zflwp_var, zfiwp_var
+      REAL d_rei_dt
 
 ! Abderrahmane oct 2009
       Real reliq(klon, klev), reice(klon, klev)
 
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! FH : 2011/05/24
+!
+! rei = ( rei_max - rei_min ) * T(°C) / 81.4 + rei_max
+! to be used for a temperature in celcius T(°C) < 0
+! rei=rei_min for T(°C) < -81.4
+!
+! Calcul de la pente de la relation entre rayon effective des cristaux
+! et la température.
+! Pour retrouver les résultats numériques de la version d'origine,
+! on impose 0.71 quand on est proche de 0.71
+
+      d_rei_dt=(rei_max-rei_min)/81.4
+      if (abs(d_rei_dt-0.71)<1.e-4) d_rei_dt=0.71
+!      print*,'d_rei_dT ',d_rei_dt,rei_min,rei_max
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 c
 c Calculer l'epaisseur optique et l'emmissivite des nuages
@@ -294 +311 @@
                   
                      tc = t(i,k)-273.15
-                     rei = 0.71*tc + 61.29 
-                     if (tc.le.-81.4) rei = 3.5 
+                     rei = d_rei_dt*tc + rei_max
+                     if (tc.le.-81.4) rei = rei_min
                      if (zflwp_var.eq.0.) radius = 1. 
                      if (zfiwp_var.eq.0. .or. rei.le.0.) rei = 1. 
@@ -321 +338 @@
 c     asymptotical value of 3.5 microns at T<-81.4 C added to be 
 c     consistent with observations of Heymsfield et al. 1986]:
+c  2011/05/24 : rei_min = 3.5 becomes a free parameter as well as rei_max=61.29
                   tc = t(i,k)-273.15
-                  rei = 0.71*tc + 61.29 
-                  if (tc.le.-81.4) rei = 3.5 
+                  rei = d_rei_dt*tc + rei_max
+                  if (tc.le.-81.4) rei = rei_min
 c     -- cloud optical thickness :