Changeset 4171 for LMDZ6/branches/LMDZ-ECRAD/libf/phylmd/newmicro.F90

Timestamp:

Jun 17, 2022, 4:24:49 PM (2 years ago)

Author:

lguez

Message:

Sync latest trunk changes to branch LMDZ-ECRAD.

Location:

LMDZ6/branches/LMDZ-ECRAD

Files:

• . (modified) (1 prop)
• libf/phylmd/newmicro.F90 (modified) (11 diffs)

LMDZ6/branches/LMDZ-ECRAD/libf/phylmd/newmicro.F90

@@ -1 +1 @@
 ! $Id$
 
-SUBROUTINE newmicro(flag_aerosol, ok_cdnc, bl95_b0, bl95_b1, paprs, pplay, t, pqlwp, pclc, &
+SUBROUTINE newmicro(flag_aerosol, ok_cdnc, bl95_b0, bl95_b1, paprs, pplay, t, pqlwp, picefra, pclc, &
     pcltau, pclemi, pch, pcl, pcm, pct, pctlwp, xflwp, xfiwp, xflwc, xfiwc, &
-    mass_solu_aero, mass_solu_aero_pi, pcldtaupi, re, fl, reliq, reice, &
+    mass_solu_aero, mass_solu_aero_pi, pcldtaupi, latitude_deg,re, fl, reliq, reice, &
     reliq_pi, reice_pi)
 
@@ -9 +9 @@
   USE phys_local_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, &
       reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra,  & 
-      zfice, dNovrN
+      zfice, dNovrN, ptconv
   USE phys_state_var_mod, ONLY: rnebcon, clwcon
   USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14)
   USE ioipsl_getin_p_mod, ONLY : getin_p
   USE print_control_mod, ONLY: lunout
+  USE lscp_tools_mod, only: icefrac_lscp
+
 
 
@@ -29 +31 @@
 
   ! t-------input-R-temperature
-  ! pqlwp---input-R-eau liquide nuageuse dans l'atmosphere dans la partie
-  ! nuageuse (kg/kg)
+  ! pqlwp---input-R-eau liquide nuageuse dans l'atmosphere dans la maille (kg/kg)
+  ! picefra--input-R-fraction de glace dans les nuages
   ! pclc----input-R-couverture nuageuse pour le rayonnement (0 a 1)
   ! mass_solu_aero-----input-R-total mass concentration for all soluble
@@ -38 +40 @@
   ! bl95_b0-input-R-a PARAMETER, may be varied for tests (s-sea, l-land)
   ! bl95_b1-input-R-a PARAMETER, may be varied for tests (    -"-      )
+  ! latitude_deg-input latitude in degrees
 
   ! re------output-R-Cloud droplet effective radius multiplied by fl [um]
@@ -58 +61 @@
   include "radepsi.h"
   include "radopt.h"
+  include "clesphys.h"
 
   ! choix de l'hypothese de recouvrement nuageuse via radopt.h (IM, 19.07.2016)
@@ -81 +85 @@
   REAL t(klon, klev)
   REAL pclc(klon, klev)
-  REAL pqlwp(klon, klev)
+  REAL pqlwp(klon, klev), picefra(klon,klev)
   REAL pcltau(klon, klev)
   REAL pclemi(klon, klev)
   REAL pcldtaupi(klon, klev)
+  REAL latitude_deg(klon)
 
   REAL pct(klon)
@@ -124 +129 @@
   REAL, PARAMETER :: t_glace_max_old = 273.13
 
-  REAL rel, tc, rei
+  REAL rel, tc, rei, iwc, dei, deimin, deimax
   REAL k_ice0, k_ice, df
   PARAMETER (k_ice0=0.005) ! units=m2/g
@@ -148 +153 @@
   ! jq-end
   ! IM cf. CR:parametres supplementaires
+  REAL dzfice(klon,klev)
   REAL zclear(klon)
   REAL zcloud(klon)
@@ -229 +235 @@
   ELSE ! of IF (iflag_t_glace.EQ.0)
     DO k = 1, klev
-        CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:,k))
- 
-
-        ! JBM: icefrac_lsc is now contained icefrac_lsc_mod
+
+! JBM: icefrac_lsc is now contained icefrac_lsc_mod
 !       zfice(i, k) = icefrac_lsc(t(i,k), t_glace_min, &
 !                                 t_glace_max, exposant_glace)
-      DO i = 1, klon
+
+      IF (ok_new_lscp) THEN
+          CALL icefrac_lscp(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:,k),dzfice(:,k))
+      ELSE
+          CALL icefrac_lsc(klon,t(:,k),pplay(:,k)/paprs(:,1),zfice(:,k))
+      ENDIF
+
+      DO i = 1, klon
+        
+        IF ((.NOT. ptconv(i,k)) .AND. ok_new_lscp .AND. ok_icefra_lscp) THEN
+        ! EV: take the ice fraction directly from the lscp code
+        ! consistent only for non convective grid points
+        ! critical for mixed phase clouds
+            zfice(i,k)=picefra(i,k)
+        ENDIF
+
         ! -layer calculation
         rhodz(i, k) = (paprs(i,k)-paprs(i,k+1))/rg ! kg/m2
@@ -321 +340 @@
             rhodz(i, k)
           zfiwp_var = 1000.*zfice(i, k)*pqlwp(i, k)/pclc(i, k)*rhodz(i, k)
-          tc = t(i, k) - 273.15
-          rei = d_rei_dt*tc + rei_max
-          IF (tc<=-81.4) rei = rei_min
+          ! Calculation of ice cloud effective radius in micron
+          IF (iflag_rei .EQ. 1) THEN
+            ! when we account for precipitation in the radiation scheme,
+            ! It is recommended to use the rei formula from Sun and Rikkus 1999 with a revision
+            ! from Sun 2001 (as in the IFS model)
+            iwc=zfice(i, k)*pqlwp(i, k)/pclc(i,k)*zrho(i,k)*1000. !in cloud ice water content in g/m3
+            dei=(1.2351+0.0105*(t(i,k)-273.15))*(45.8966*(iwc**0.2214) + &
+               & 0.7957*(iwc**0.2535)*(t(i,k)-83.15))
+            !deimax=155.0
+            !deimin=20.+40*cos(abs(latitude_deg(i))/180.*RPI)
+            !Etienne: deimax and deimin controled by rei_max and rei_min in physiq.def
+            deimax=rei_max*2.0
+            deimin=2.0*rei_min+40*cos(abs(latitude_deg(i))/180.*RPI) 
+            dei=min(dei,deimax)
+            dei=max(dei,deimin)
+            rei=3.*sqrt(3.)/8.*dei
+           ELSE
+            ! Default
+            ! for ice clouds: as a function of the ambiant temperature
+            ! [formula used by Iacobellis and Somerville (2000), with an
+            ! asymptotical value of 3.5 microns at T<-81.4 C added to be
+            ! consistent with observations of Heymsfield et al. 1986]:
+            ! 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 = d_rei_dt*tc + rei_max
+            IF (tc<=-81.4) rei = rei_min
+           ENDIF
 
           ! -- cloud optical thickness :
@@ -401 +445 @@
         rel = rad_chaud(i, k)
 
-        ! for ice clouds: as a function of the ambiant temperature
-        ! [formula used by Iacobellis and Somerville (2000), with an
-        ! asymptotical value of 3.5 microns at T<-81.4 C added to be
-        ! consistent with observations of Heymsfield et al. 1986]:
-        ! 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 = d_rei_dt*tc + rei_max
-        IF (tc<=-81.4) rei = rei_min
-
+        ! Calculation of ice cloud effective radius in micron
+
+
+        IF (iflag_rei .GT. 0) THEN
+
+            ! when we account for precipitation in the radiation scheme,
+            ! we use the rei formula from Sun and Rikkus 1999 with a revision
+            ! from Sun 2001 (as in the IFS model)
+            iwc=zfice(i, k)*pqlwp(i, k)/pclc(i,k)*zrho(i,k)*1000. !in cloud ice water content in g/m3
+            dei=(1.2351+0.0105*(t(i,k)-273.15))*(45.8966*(iwc**0.2214) + &
+               &0.7957*(iwc**0.2535)*(t(i,k)-83.15))
+            !deimax=155.0
+            !deimin=20.+40*cos(abs(latitude_deg(i))/180.*RPI)
+            !Etienne: deimax and deimin controled by rei_max and rei_min in physiq.def
+            deimax=rei_max*2.0
+            deimin=2.0*rei_min+40*cos(abs(latitude_deg(i))/180.*RPI) 
+            dei=min(dei,deimax)
+            dei=max(dei,deimin)
+            rei=3.*sqrt(3.)/8.*dei
+       
+        ELSE
+            ! Default
+            ! for ice clouds: as a function of the ambiant temperature
+            ! [formula used by Iacobellis and Somerville (2000), with an
+            ! asymptotical value of 3.5 microns at T<-81.4 C added to be
+            ! consistent with observations of Heymsfield et al. 1986]:
+            ! 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 = d_rei_dt*tc + rei_max
+            IF (tc<=-81.4) rei = rei_min
+        ENDIF
         ! -- cloud optical thickness :
         ! [for liquid clouds, traditional formula,