Changeset 517 for LMDZ.3.3/branches/rel-LF/libf/phylmd/newmicro.F

Timestamp:

Apr 16, 2004, 5:43:38 PM (20 years ago)

Author:

lmdzadmin

Message:

Inclusion des modifications de O. Boucher et de J. Quaas pour le calcul des
premiers effets directs et indirects dus aux aerosols
LF

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/newmicro.F (modified) (5 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/newmicro.F

@@ -1 +1 @@
       SUBROUTINE newmicro (paprs, pplay,ok_newmicro,
      .                  t, pqlwp, pclc, pcltau, pclemi,
-cIM    .                  pch, pcl, pcm, pct, pctlwp)
      .                  pch, pcl, pcm, pct, pctlwp,
-     .                  xflwp, xfiwp, xflwc, xfiwc)
-
+     s                  xflwp, xfiwp, xflwc, xfiwc,
+     e                  ok_aie, 
+     e                  sulfate, sulfate_pi, 
+     e                  bl95_b0, bl95_b1,
+     s                  cldtaupi, re, fl)
       IMPLICIT none
 c======================================================================
@@ -15 +17 @@
 c pclc----input-R-couverture nuageuse pour le rayonnement (0 a 1)
 c 
+c ok_aie--input-L-apply aerosol indirect effect or not
+c sulfate-input-R-sulfate aerosol mass concentration [um/m^3]
+c sulfate_pi-input-R-dito, pre-industrial value
+c bl95_b0-input-R-a parameter, may be varied for tests (s-sea, l-land)
+c bl95_b1-input-R-a parameter, may be varied for tests (    -"-      )
+c      
+c cldtaupi-output-R-pre-industrial value of cloud optical thickness, 
+c                   needed for the diagnostics of the aerosol indirect 
+c                   radiative forcing (see radlwsw)
+c re------output-R-Cloud droplet effective radius multiplied by fl [um]
+c fl------output-R-Denominator to re, introduced to avoid problems in
+c                  the averaging of the output. fl is the fraction of liquid
+c                  water clouds within a grid cell           
 c pcltau--output-R-epaisseur optique des nuages
 c pclemi--output-R-emissivite des nuages (0 a 1)
@@ -66 +81 @@
       parameter (DF=1.66) ! diffusivity factor
 c sb --
-
+cjq for the aerosol indirect effect
+cjq introduced by Johannes Quaas (quaas@lmd.jussieu.fr), 27/11/2003
+cjq      
+      LOGICAL ok_aie            ! Apply AIE or not?
+      LOGICAL ok_a1lwpdep       ! a1 LWP dependent?
+      
+      REAL sulfate(klon, klev)  ! sulfate aerosol mass concentration [ug m-3]
+      REAL cdnc(klon, klev)     ! cloud droplet number concentration [m-3]
+      REAL re(klon, klev)       ! cloud droplet effective radius [um]
+      REAL sulfate_pi(klon, klev)  ! sulfate aerosol mass concentration [ug m-3] (pre-industrial value)
+      REAL cdnc_pi(klon, klev)     ! cloud droplet number concentration [m-3] (pi value)
+      REAL re_pi(klon, klev)       ! cloud droplet effective radius [um] (pi value)
+      
+      REAL fl(klon, klev)       ! xliq * rneb (denominator to re; fraction of liquid water clouds within the grid cell)
+      
+      REAL bl95_b0, bl95_b1     ! Parameter in B&L 95-Formula
+      
+      REAL cldtaupi(klon, klev) ! pre-industrial cloud opt thickness for diag
+cjq-end    
 c
 c Calculer l'epaisseur optique et l'emmissivite des nuages
@@ -119 +152 @@
 
 c for liquid water clouds: 
+         IF (ok_aie) THEN
+            ! Formula "D" of Boucher and Lohmann, Tellus, 1995
+            !             
+            cdnc(i,k) = 10.**(bl95_b0+bl95_b1*
+     .           log(MAX(sulfate(i,k),1.e-4))/log(10.))*1.e6 !-m-3
+            ! Cloud droplet number concentration (CDNC) is restricted
+            ! to be within [20, 1000 cm^3]
+            ! 
+            cdnc(i,k)=MIN(1000.e6,MAX(20.e6,cdnc(i,k)))
+            !
+            !
+            cdnc_pi(i,k) = 10.**(bl95_b0+bl95_b1*
+     .           log(MAX(sulfate_pi(i,k),1.e-4))/log(10.))*1.e6 !-m-3
+            cdnc_pi(i,k)=MIN(1000.e6,MAX(20.e6,cdnc_pi(i,k)))
+            !            
+            !
+            ! air density: pplay(i,k) / (RD * zT(i,k)) 
+            ! factor 1.1: derive effective radius from volume-mean radius
+            ! factor 1000 is the water density
+            ! _chaud means that this is the CDR for liquid water clouds
+            !
+            rad_chaud = 
+     .           1.1 * ( (pqlwp(i,k) * pplay(i,k) / (RD * T(i,k)) )  
+     .               / (4./3. * RPI * 1000. * cdnc(i,k)) )**(1./3.)
+            !
+            ! Convert to um. CDR shall be at least 3 um.
+            !
+c           rad_chaud = MAX(rad_chaud*1.e6, 3.) 
+            rad_chaud = MAX(rad_chaud*1.e6, 5.) 
+            
+            ! Pre-industrial cloud opt thickness
+            !
+            ! "radius" is calculated as rad_chaud above (plus the 
+            ! ice cloud contribution) but using cdnc_pi instead of
+            ! cdnc.
+            radius = 
+     .           1.1 * ( (pqlwp(i,k) * pplay(i,k) / (RD * T(i,k)) )  
+     .               / (4./3. * RPI * 1000. * cdnc_pi(i,k)) )**(1./3.)
+            radius = MAX(radius*1.e6, 3.) 
+            
+            tc = t(i,k)-273.15
+            rei = 0.71*tc + 61.29 
+            if (tc.le.-81.4) rei = 3.5 
+            if (zflwp(i).eq.0.) radius = 1. 
+            if (zfiwp(i).eq.0. .or. rei.le.0.) rei = 1. 
+            cldtaupi(i,k) = 3.0/2.0 * zflwp(i) / radius
+     .             + zfiwp(i) * (3.448e-03  + 2.431/rei)
+         ENDIF                  ! ok_aie
+         ! For output diagnostics
+         !
+         ! Cloud droplet effective radius [um]
+         !
+         ! we multiply here with f * xl (fraction of liquid water
+         ! clouds in the grid cell) to avoid problems in the
+         ! averaging of the output.
+         ! In the output of IOIPSL, derive the real cloud droplet 
+         ! effective radius as re/fl
+         !
+         fl(i,k) = pclc(i,k)*(1.-zfice)            
+         re(i,k) = rad_chaud*fl(i,k)
+            
+c-jq end         
+         
          rel = rad_chaud
-
 c for ice clouds: as a function of the ambiant temperature
 c [formula used by Iacobellis and Somerville (2000), with an 
@@ -156 +251 @@
          IF (lo) pcltau(i,k) = 0.0
          IF (lo) pclemi(i,k) = 0.0
+         
+         IF (lo) cldtaupi(i,k) = 0.0
+         IF (.NOT.ok_aie) cldtaupi(i,k)=pcltau(i,k)            
       ENDDO
       ENDDO