Changeset 1279 for LMDZ4/trunk/libf/phylmd/newmicro.F

Timestamp:

Dec 10, 2009, 10:02:56 AM (15 years ago)

Author:

Laurent Fairhead

Message:

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

Location:

LMDZ4/trunk

Files:

• . (modified) (1 prop)
• libf/phylmd/newmicro.F (modified) (12 diffs)

LMDZ4/trunk/libf/phylmd/newmicro.F

@@ -1 @@
-!
-! $Header$
-!
+! $Id$
+!     
       SUBROUTINE newmicro (paprs, pplay,ok_newmicro,
      .                  t, pqlwp, pclc, pcltau, pclemi,
@@ -7 +6 @@
      s                  xflwp, xfiwp, xflwc, xfiwc,
      e                  ok_aie, 
-     e                  sulfate, sulfate_pi, 
+     e                  mass_solu_aero, mass_solu_aero_pi, 
      e                  bl95_b0, bl95_b1,
-     s                  cldtaupi, re, fl)
+     s                  cldtaupi, re, fl, reliq, reice)
+
       USE dimphy
       IMPLICIT none
@@ -22 +22 @@
 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 mass_solu_aero-----input-R-total mass concentration for all soluble aerosols[ug/m^3]
+c mass_solu_aero_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 (    -"-      )
@@ -94 +94 @@
       LOGICAL ok_a1lwpdep       ! a1 LWP dependent?
       
-      REAL sulfate(klon, klev)  ! sulfate aerosol mass concentration [ug m-3]
+      REAL mass_solu_aero(klon, klev)    ! total mass concentration for all soluble aerosols [ug m-3]
+      REAL mass_solu_aero_pi(klon, klev) ! - " - (pre-industrial value)
       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)
@@ -119 +119 @@
       REAL diff_paprs(klon, klev), zfice1, zfice2(klon, klev)
       REAL rad_chaud_tab(klon, klev), zflwp_var, zfiwp_var
+
+! Abderrahmane oct 2009
+      Real reliq(klon, klev), reice(klon, klev)
 
 c
@@ -157 +160 @@
                                 !             
                   cdnc(i,k) = 10.**(bl95_b0+bl95_b1*
-     &                 log(MAX(sulfate(i,k),1.e-4))/log(10.))*1.e6 !-m-3
+     &               log(MAX(mass_solu_aero(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]
@@ -165 +168 @@
                                 !
                   cdnc_pi(i,k) = 10.**(bl95_b0+bl95_b1*
-     &                 log(MAX(sulfate_pi(i,k),1.e-4))/log(10.))*1.e6 !-m-3
+     &               log(MAX(mass_solu_aero_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)))
                ENDDO
@@ -221 +224 @@
                   re(i,k) = rad_chaud_tab(i,k)*fl(i,k)
                   
+                  rel = 0.
+                  rei = 0.
                   pclc(i,k) = 0.0
                   pcltau(i,k) = 0.0
@@ -252 +257 @@
                      cldtaupi(i,k) = 3.0/2.0 * zflwp_var / radius
      &                    + zfiwp_var * (3.448e-03  + 2.431/rei)
+
                   ENDIF         ! ok_aie
                                 ! For output diagnostics
@@ -280 +286 @@
 c     for ice clouds, Ebert & Curry (1992)] 
                   
-                  if (zflwp_var.eq.0.) rel = 1. 
-                  if (zfiwp_var.eq.0. .or. rei.le.0.) rei = 1. 
-                  pcltau(i,k) = 3.0/2.0 * ( zflwp_var/rel )
+                 if (zflwp_var.eq.0.) rel = 1. 
+                 if (zfiwp_var.eq.0. .or. rei.le.0.) rei = 1. 
+                 pcltau(i,k) = 3.0/2.0 * ( zflwp_var/rel )
      &                 + zfiwp_var * (3.448e-03  + 2.431/rei)
 c     -- cloud infrared emissivity:
@@ -296 +302 @@
 
                ENDIF
-               
+              reliq(i,k)=rel
+              reice(i,k)=rei 
+!              if (i.eq.1) then
+!              print*,'Dans newmicro rel, rei :',rel, rei
+!              print*,'Dans newmicro reliq, reice :',
+!     $             reliq(i,k),reice(i,k)
+!              endif
+
             ENDDO
          ENDDO
@@ -400 +413 @@
             DO i = 1, klon
                zclear(i)=zclear(i)*(1.-MAX(pclc(i,k),zcloud(i)))
-     &              /(1.-MIN(zcloud(i),1.-ZEPSEC))
+     &              /(1.-MIN(real(zcloud(i), kind=8),1.-ZEPSEC))
                pct(i)=1.-zclear(i) 
                IF (pplay(i,k).LE.cetahb*paprs(i,1)) THEN
                   pch(i) = pch(i)*(1.-MAX(pclc(i,k),zcloud(i)))
-     &                 /(1.-MIN(zcloud(i),1.-ZEPSEC))
+     &                 /(1.-MIN(real(zcloud(i), kind=8),1.-ZEPSEC))
                ELSE IF (pplay(i,k).GT.cetahb*paprs(i,1) .AND.
      &                 pplay(i,k).LE.cetamb*paprs(i,1)) THEN
                   pcm(i) = pcm(i)*(1.-MAX(pclc(i,k),zcloud(i)))
-     &                 /(1.-MIN(zcloud(i),1.-ZEPSEC))
+     &                 /(1.-MIN(real(zcloud(i), kind=8),1.-ZEPSEC))
                ELSE IF (pplay(i,k).GT.cetamb*paprs(i,1)) THEN 
                   pcl(i) = pcl(i)*(1.-MAX(pclc(i,k),zcloud(i)))
-     &                 /(1.-MIN(zcloud(i),1.-ZEPSEC))
+     &                 /(1.-MIN(real(zcloud(i), kind=8),1.-ZEPSEC))
                endif
                zcloud(i)=pclc(i,k)