Changeset 2738 for LMDZ5/trunk/libf/phylmd/rrtm/aeropt_lw_rrtm.F90

Timestamp:

Jan 2, 2017, 11:24:30 PM (8 years ago)

Author:

oboucher

Message:

Swapping the order of CSSO4 and ASSO4 aerosols (and fixing an issue on aerindex). Preparing the ground for nitrate aerosols (coarse soluble, accumulation soluble, coarse insoluble). Modifying the LW aeropt routine so that it is compatible with both INCA and climatological aerosols (for dust only). Adding a new flag ok_alw for activating aerosol direct LW effect (for dust only). This change is bit comparable for flag_aerosol=6, flag_rrtm=1, NSW=6.

File:

• LMDZ5/trunk/libf/phylmd/rrtm/aeropt_lw_rrtm.F90 (modified) (1 diff)

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

@@ -2 +2 @@
 ! aeropt_lw_rrtm.F90 2014-05-13 C. Kleinschmitt
 !                    2016-05-03 O. Boucher 
+!                    2016-12-17 O. Boucher
 !
 ! This routine feeds aerosol LW properties to RRTM
 ! we only consider absorption (not scattering)
+! we only consider dust for now
 
-SUBROUTINE AEROPT_LW_RRTM(aerosol_couple,paprs,tr_seri)
+SUBROUTINE AEROPT_LW_RRTM(ok_alw, pdel, zrho, flag_aerosol, m_allaer, m_allaer_pi)
 
   USE dimphy
   USE aero_mod
-  USE infotrac_phy
   USE phys_state_var_mod, ONLY : tau_aero_lw_rrtm
   USE YOERAD, ONLY : NLW
+  USE YOMCST
 
   IMPLICIT NONE
 
-  INCLUDE "YOMCST.h"
   INCLUDE "clesphys.h"
   !
   ! Input arguments:
   !
-  LOGICAL, INTENT(IN)                         :: aerosol_couple
-  REAL, DIMENSION(klon,klev+1), INTENT(in)    :: paprs
-  REAL, DIMENSION(klon,klev,nbtr), INTENT(in) :: tr_seri
-  !
-  REAL, DIMENSION(klon,klev) :: zdp, mass_temp
-  ! 
+  LOGICAL, INTENT(IN)                                 :: ok_alw
+  INTEGER, INTENT(IN)                                 :: flag_aerosol
+  REAL, DIMENSION(klon,klev),     INTENT(in)          :: pdel, zrho
+  REAL, DIMENSION(klon,klev,naero_tot),   INTENT(in)  :: m_allaer, m_allaer_pi
   !
   INTEGER inu, i, k
-  INTEGER :: id_ASBCM, id_ASPOMM, id_ASSO4M, id_ASMSAM, id_CSSO4M, id_CSMSAM, id_SSSSM
-  INTEGER :: id_CSSSM, id_ASSSM, id_CIDUSTM, id_AIBCM, id_AIPOMM, id_ASNO3M, id_CSNO3M, id_CINO3M
-
+  REAL :: zdh(klon,klev)
+  REAL :: tmp_var, tmp_var_pi
   !
   !--absorption coefficient for CIDUST
   REAL:: alpha_abs_CIDUST_16bands(nbands_lw_rrtm)   !--unit m2/g 
   DATA alpha_abs_CIDUST_16bands /              & 
-       0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, & 
-       0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0  /
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, &
+  0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0 /
   !
   IF (NLW.NE.nbands_lw_rrtm) THEN
-     print *,'Erreur NLW doit etre egal a 16 pour cette routine'
-     stop
+     PRINT *,'Erreur NLW doit etre egal a 16 pour cette routine'
+     STOP
   ENDIF
   ! 
-  IF (aerosol_couple) THEN 
-  ! 
-      do i=1,nbtr
-      select case(trim(solsym(i)))
-         case ("ASBCM")
-            id_ASBCM = i
-         case ("ASPOMM")
-            id_ASPOMM = i
-         case ("ASSO4M")
-            id_ASSO4M = i
-         case ("ASMSAM")
-            id_ASMSAM = i
-         case ("CSSO4M")
-            id_CSSO4M = i
-         case ("CSMSAM")
-            id_CSMSAM = i
-         case ("SSSSM")
-            id_SSSSM = i
-         case ("CSSSM")
-            id_CSSSM = i
-         case ("ASSSM")
-            id_ASSSM = i
-         case ("CIDUSTM")
-            id_CIDUSTM = i
-         case ("AIBCM")
-            id_AIBCM = i
-         case ("AIPOMM")
-            id_AIPOMM = i
-         case ("ASNO3M")
-            id_ASNO3M = i
-         case ("CSNO3M")
-            id_CSNO3M = i
-         case ("CINO3M")
-            id_CINO3M = i
-         end select
-    enddo
+  IF (ok_alw) THEN                                   !--aerosol LW effects
+   !
+   IF (flag_aerosol.EQ.5.OR.flag_aerosol.EQ.6) THEN  !--Dust
+    !
+    zdh(:,:)=pdel(:,:)/(RG*zrho(:,:))      !-- m
     !
     DO k=1, klev
-       zdp(:,k) = (paprs(:,k)-paprs(:,k+1))/RG !--kg/m2
-    ENDDO
-    !
-    !--for now only consider dust
-    mass_temp(:,:)=tr_seri(:,:,id_CIDUSTM)  !--kg/kg
-    !
-    DO inu=1,NLW
-      !
-      !--total aerosol
-      tau_aero_lw_rrtm(:,:,2,inu) = mass_temp(:,:)*zdp(:,:)*1000.*alpha_abs_CIDUST_16bands(inu)
-      !--no aerosol at all
-      tau_aero_lw_rrtm(:,:,1,inu) = 0.0
+      DO i=1, klon
+         !
+         tmp_var   =m_allaer(i,k,id_CIDUSTM_phy)   /1.e6*zdh(i,k) !-- g/m2
+         tmp_var_pi=m_allaer_pi(i,k,id_CIDUSTM_phy)/1.e6*zdh(i,k) !-- g/m2
+         !
+         DO inu=1, NLW
+           !
+           !--total aerosol
+           tau_aero_lw_rrtm(i,k,2,inu) = MAX(1.e-15,tmp_var*alpha_abs_CIDUST_16bands(inu))
+           !--natural aerosol 
+           tau_aero_lw_rrtm(i,k,1,inu) = MAX(1.e-15,tmp_var_pi*alpha_abs_CIDUST_16bands(inu))
+           !
+         ENDDO
+      ENDDO
       !
     ENDDO
     !
-    !--avoid very small values
-    tau_aero_lw_rrtm = MAX(tau_aero_lw_rrtm,1.e-15)
-    ! 
-  ELSE !--not aerosol_couple
+   ENDIF
+   !
+  ELSE !--no aerosol LW effects
     !
-    !--no LW effects if not coupled to INCA
     tau_aero_lw_rrtm = 1.e-15
   ENDIF