--- trunk/Sources/phylmd/nuage.f	2017/02/27 15:44:55	213
+++ trunk/Sources/phylmd/nuage.f	2017/03/30 14:25:18	217
@@ -5,14 +5,12 @@
 contains
 
   SUBROUTINE nuage (paprs, pplay, t, pqlwp, pclc, pcltau, pclemi, pch, pcl, &
-       pcm, pct, pctlwp, ok_aie, sulfate, sulfate_pi, bl95_b0, bl95_b1, &
-       cldtaupi, re, fl)
+       pcm, pct, pctlwp)
 
     ! From LMDZ4/libf/phylmd/nuage.F, version 1.1.1.1, 2004/05/19 12:53:07
 
     use dimphy, only: klon, klev
-    use nr_util, only: pi
-    use SUPHEC_M, only: rd, rg
+    use SUPHEC_M, only: rg
 
     ! Author: Z. X. Li (LMD/CNRS)
     ! Date: 1993/09/10
@@ -30,41 +28,10 @@
     REAL, intent(inout):: pclc(klon, klev)
     ! couverture nuageuse pour le rayonnement (0 \`a 1)
 
-    REAL pcltau(klon, klev), pclemi(klon, klev)
-    ! pcltau--output-R-epaisseur optique des nuages
-    ! pclemi--output-R-emissivite des nuages (0 a 1)
-
+    REAL, intent(out):: pcltau(klon, klev) ! \'epaisseur optique des nuages
+    real pclemi(klon, klev) ! pclemi--output-R-emissivite des nuages (0 a 1)
     REAL pch(klon), pcl(klon), pcm(klon), pct(klon), pctlwp(klon)
 
-    !jq for the aerosol indirect effect
-    !jq introduced by Johannes Quaas (quaas@lmd.jussieu.fr), 27/11/2003
-    !jq
-    LOGICAL ok_aie ! Apply AIE or not?
-    ! ok_aie--input-L-apply aerosol indirect effect or not
-
-    REAL sulfate(klon, klev) ! sulfate aerosol mass concentration [ug m-3]
-    ! sulfate-input-R-sulfate aerosol mass concentration [um/m^3]
-
-    REAL, intent(in):: sulfate_pi(klon, klev) ! sulfate aerosol mass concentration [ug m-3] (pre-industrial value)
-    ! sulfate_pi-input-R-dito, pre-industrial value
-
-    REAL bl95_b0, bl95_b1 ! Parameter in B&L 95-Formula
-    ! 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 (-"-)
-
-    REAL cldtaupi(klon, klev) ! pre-industrial cloud opt thickness for diag
-    ! cldtaupi-output-R-pre-industrial value of cloud optical thickness,
-    ! needed for the diagnostics of the aerosol indirect
-    ! radiative forcing (see radlwsw)
-
-    REAL re(klon, klev) ! cloud droplet effective radius [um]
-    ! re------output-R-Cloud droplet effective radius multiplied by fl [um]
-
-    REAL fl(klon, klev) ! xliq * rneb (denominator to re ; fraction of liquid water clouds within the grid cell)
-    ! fl------output-R-Denominator to re, introduced to avoid problems in
-    ! the averaging of the output. fl is the fraction of liquid
-    ! water clouds within a grid cell
-
     ! Local:
 
     LOGICAL lo
@@ -86,9 +53,6 @@
     INTEGER nexpo ! exponentiel pour glace/eau
     PARAMETER (nexpo=6)
 
-    REAL cdnc(klon, klev) ! cloud droplet number concentration [m-3]
-    REAL cdnc_pi(klon, klev) ! cloud droplet number concentration [m-3] (pi value)
-
     !--------------------------------------------------------------------
     
     ! Calculer l'epaisseur optique et l'emmissivite des nuages
@@ -105,57 +69,6 @@
           zfice = MIN(MAX(zfice, 0.0), 1.0)
           zfice = zfice**nexpo
 
-          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. * PI * 1000. * cdnc(i, k)))**(1./3.)
-
-             ! Convert to um. CDR shall be at least 3 um.
-
-             rad_chaud = MAX(rad_chaud*1.e6, 3.)
-
-             ! 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)
-
-             ! 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 = MAX(1.1e6 * ((pqlwp(i, k)*pplay(i, k)/(RD*T(i, k))) &
-                  / (4./3.*PI*1000.*cdnc_pi(i, k)))**(1./3.), &
-                  3.) * (1.-zfice) + rad_froid * zfice
-             cldtaupi(i, k) = 3.0/2.0 * zflwp / radius
-
-          END IF ! ok_aie
-
           radius = rad_chaud * (1.-zfice) + rad_froid * zfice
           coef = coef_chau * (1.-zfice) + coef_froi * zfice
           pcltau(i, k) = 3.0/2.0 * zflwp / radius
@@ -164,8 +77,6 @@
           IF (lo) pclc(i, k) = 0.0
           IF (lo) pcltau(i, k) = 0.0
           IF (lo) pclemi(i, k) = 0.0
-
-          IF (.NOT.ok_aie) cldtaupi(i, k)=pcltau(i, k)
        END DO
     END DO