Changeset 2854 for LMDZ5/trunk/libf/phylmd/rrtm

Timestamp:

Apr 14, 2017, 4:42:31 PM (7 years ago)

Author:

oboucher

Message:

Introducing dry AOD diagnostics for the total aerosols and specieswise
The calculations are only performed if the diagnostics are requested

Location:

LMDZ5/trunk/libf/phylmd/rrtm

Files:

• aeropt_5wv_rrtm.F90 (modified) (8 diffs)
• readaerosol_optic_rrtm.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/rrtm/aeropt_5wv_rrtm.F90

@@ -8 +8 @@
    flag_bc_internal_mixture, & 
    pplay, t_seri,            &
-   tausum, tau )
+   tausum, drytausum, tau )
 
   USE DIMPHY
   USE aero_mod
-  USE phys_local_var_mod, ONLY: od443aer,od550aer,od865aer,ec550aer,od550lt1aer,abs550aer
+  USE phys_local_var_mod, ONLY: od443aer,od550aer,dryod550aer,od865aer,ec550aer,od550lt1aer,abs550aer
+  USE phys_output_var_mod, ONLY: dryaod_diag
   USE YOMCST, ONLY: RD,RG
 
@@ -67 +68 @@
   REAL, DIMENSION(klon), INTENT(OUT)                      :: ai      ! POLDER aerosol index 
   REAL, DIMENSION(klon,nwave,naero_tot), INTENT(OUT)      :: tausum
+  REAL, DIMENSION(klon,naero_tot), INTENT(OUT)            :: drytausum
   REAL, DIMENSION(klon,klev,nwave,naero_tot), INTENT(OUT) :: tau
   !
@@ -342 +344 @@
   ! Initialisations
   ai(:) = 0.
+  abs550aer(:)=0.0
+  drytausum(:,:) = 0.
   tausum(:,:,:) = 0.
   tau(:,:,:,:)=0.
-
-  abs550aer(:)=0.0
 
   DO k=1, klev
@@ -476 +478 @@
     DO la=1,las
 
-    !--only 443, 550, 670 and 865 nm are used
+    !--only 443, 550, and 865 nm are used
     !--to save time 670 and AI are not computed for CMIP6
     !IF (la.NE.la443.AND.la.NE.la550.AND.la.NE.la670.AND.la.NE.la865) CYCLE
@@ -511 +513 @@
              tausum(i,la,aerindex)=tausum(i,la,aerindex)+tau(i,k,la,aerindex)
 
+             IF (la.EQ.la550.AND.dryaod_diag) THEN 
+                tau_ae5wv_int = alpha_MG_5wv(1,la,classbc)
+                drytausum(i,aerindex)=drytausum(i,aerindex)+m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*tau_ae5wv_int*fac
+             ENDIF
+
              IF (la.EQ.la550) THEN 
                 abs_ae5wv_int = abs_MG_5wv(RH_num(i,k),la,classbc)+DELTA(i,k)* &
@@ -532 +539 @@
               tausum(i,la,aerindex)=tausum(i,la,aerindex)+tau(i,k,la,aerindex)
 
+              IF (la.EQ.la550.AND.dryaod_diag) THEN 
+                 tau_ae5wv_int = alpha_aers_5wv(1,la,spsol)
+                 drytausum(i,aerindex)=drytausum(i,aerindex)+m_allaer(i,k,aerindex)/1.e6*zdh(i,k)*tau_ae5wv_int*fac
+              ENDIF
+
               IF (la.EQ.la550) THEN 
                  abs_ae5wv_int = abs_aers_5wv(RH_num(i,k),la,spsol)+DELTA(i,k)* &
@@ -554 +566 @@
             tausum(i,la,aerindex)= tausum(i,la,aerindex)+tau(i,k,la,aerindex)
 
+            IF (la.EQ.la550.AND.dryaod_diag) THEN 
+              drytausum(i,aerindex)= drytausum(i,aerindex)+tau(i,k,la,aerindex)
+            ENDIF
+
             IF (la.EQ.la550) THEN 
                abs_ae5wv_int = abs_aeri_5wv(la,spsol)
@@ -573 +589 @@
   od865aer(:)=SUM(tausum(:,la865,:),dim=2)
 
+!--dry AOD calculation for diagnostics la=la550
+  dryod550aer(:)=SUM(drytausum(:,:),dim=2)
+
 !--extinction coefficient for diagnostic
   ec550aer(:,:)=SUM(tau(:,:,la550,:),dim=3)/zdh(:,:)

LMDZ5/trunk/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

@@ -6 +6 @@
      tr_seri, mass_solu_aero, mass_solu_aero_pi, &
      tau_aero, piz_aero, cg_aero, &
-     tausum_aero, tau3d_aero )
+     tausum_aero, drytausum_aero, tau3d_aero )
 
   ! This routine will :
@@ -52 +52 @@
   REAL, DIMENSION(klon,klev,2,NSW), INTENT(OUT) :: cg_aero     ! asymmetry parameter aerosol
   REAL, DIMENSION(klon,nwave,naero_tot), INTENT(OUT)       :: tausum_aero
+  REAL, DIMENSION(klon,naero_tot), INTENT(OUT)             :: drytausum_aero
   REAL, DIMENSION(klon,klev,nwave,naero_tot), INTENT(OUT)  :: tau3d_aero
 
@@ -320 +321 @@
        flag_bc_internal_mixture,     & 
        pplay, t_seri,                &
-       tausum_aero, tau3d_aero )
+       tausum_aero, drytausum_aero, tau3d_aero )
 
   !--call LW optical properties for tropospheric aerosols