Changeset 1667 for LMDZ5/trunk/libf/phylmd/physiq.F

Timestamp:

Oct 11, 2012, 12:42:12 PM (12 years ago)

Author:

idelkadi

Message:

Reorganisation des differents cas d'utilisation des aerosols avec les flag ok_ade, ok_aie, flag_aerosol

File:

• LMDZ5/trunk/libf/phylmd/physiq.F (modified) (8 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F

@@ -1247 +1247 @@
       LOGICAL, SAVE :: mskocean_beta
 c$OMP THREADPRIVATE(mskocean_beta)
-      REAL, dimension(klon, klev) :: beta       ! facteur sur cldtaurad et cldemirad pour evaluer les retros liees aux CRF
-      REAL, dimension(klon, klev) :: cldtaurad  ! epaisseur optique pour radlwsw,COSP
-      REAL, dimension(klon, klev) :: cldemirad  ! emissivite pour radlwsw,COSP 
+      REAL, dimension(klon, klev) :: beta         ! facteur sur cldtaurad et cldemirad pour evaluer les retros liees aux CRF
+      REAL, dimension(klon, klev) :: cldtaurad    ! epaisseur optique pour radlwsw,COSP
+      REAL, dimension(klon, klev) :: cldtaupirad  ! epaisseur optique pour radlwsw,COSP cas pre-industrial
+      REAL, dimension(klon, klev) :: cldemirad    ! emissivite pour radlwsw,COSP 
       INTEGER :: nbtr_tmp ! Number of tracer inside concvl
       REAL, dimension(klon,klev) :: sh_in ! Specific humidity entering in phytrac
@@ -3099 +3100 @@
 cjq - introduce the aerosol direct and first indirect radiative forcings
 cjq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
-      IF (ok_ade.OR.ok_aie) THEN
+      IF (flag_aerosol .gt. 0) THEN
          IF (.NOT. aerosol_couple)
      &        CALL readaerosol_optic(
@@ -3384 +3385 @@
 cIM betaCRF
 c
-      cldtaurad = cldtau
-      cldemirad = cldemi
+      cldtaurad   = cldtau
+      cldtaupirad = cldtaupi
+      cldemirad   = cldemi
 c
       if(lon1_beta.EQ.-180..AND.lon2_beta.EQ.180..AND.
@@ -3402 +3404 @@
          beta(i,k) = beta(i,k) * pctsrf(i,is_oce)
         endif
-        cldtaurad(i,k) = cldtau(i,k) * beta(i,k)
-        cldemirad(i,k) = cldemi(i,k) * beta(i,k)
+        cldtaurad(i,k)   = cldtau(i,k) * beta(i,k)
+        cldtaupirad(i,k) = cldtaupi(i,k) * beta(i,k)
+        cldemirad(i,k)   = cldemi(i,k) * beta(i,k)
        ENDDO
        ENDDO 
@@ -3424 +3427 @@
           beta(i,k) = beta(i,k) * pctsrf(i,is_oce)
          endif
-        cldtaurad(i,k) = cldtau(i,k) * beta(i,k)
-        cldemirad(i,k) = cldemi(i,k) * beta(i,k)
+        cldtaurad(i,k)   = cldtau(i,k) * beta(i,k)
+        cldtaupirad(i,k) = cldtaupi(i,k) * beta(i,k)
+        cldemirad(i,k)   = cldemi(i,k) * beta(i,k)
         endif
 c
@@ -3474 +3478 @@
      s        topsw_aero, topsw0_aero,
      s        solsw_aero, solsw0_aero,
-     e        cldtaupi,
+     e        cldtaupirad,
      s        topswai_aero, solswai_aero)
             
@@ -3497 +3501 @@
      e        t_seri,q_seri,wo,
      e        cldfra, cldemirad, cldtaurad,
-     e        ok_ade, ok_aie,
+     e        ok_ade, ok_aie, flag_aerosol,
      e        tau_aero, piz_aero, cg_aero,
-     e        cldtaupi,new_aod,
+     e        cldtaupirad,new_aod,
      e        zqsat, flwc, fiwc,
      s        heat,heat0,cool,cool0,radsol,albpla,
@@ -3538 +3542 @@
      e        t_seri,q_seri,wo,
      e        cldfra, cldemi, cldtau,
-     e        ok_ade, ok_aie,
+     e        ok_ade, ok_aie, flag_aerosol,
      e        tau_aero, piz_aero, cg_aero,
      e        cldtaupi,new_aod,