Changeset 2787 for LMDZ5/branches/testing/libf/phylmd/rrtm/aeropt_lw_rrtm.F90

Timestamp:

Jan 30, 2017, 5:54:45 PM (8 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r2727:2785 into testing branch

Location:

LMDZ5/branches/testing

Files:

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

LMDZ5/branches/testing/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 phys_state_var_mod, ONLY: tau_aero_lw_rrtm
+  USE YOERAD, ONLY: NLW
+  USE YOMCST, ONLY: RG
 
   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
+  CHARACTER*20 modname
   !
   !--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  /
+  DATA alpha_abs_CIDUST_16bands /                         &
+  0.001, 0.003, 0.005, 0.006, 0.012, 0.030, 0.148, 0.098, &
+  0.017, 0.053, 0.031, 0.008, 0.010, 0.011, 0.013, 0.015  /
+  !
+  modname='aeropt_lw_rrtm'
   !
   IF (NLW.NE.nbands_lw_rrtm) THEN
-     print *,'Erreur NLW doit etre egal a 16 pour cette routine'
-     stop
+    CALL abort_physic(modname,'Erreur NLW doit etre egal a 16 pour cette routine',1)
   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(:,:,1,inu) = MAX(1.e-15,tmp_var_pi*alpha_abs_CIDUST_16bands(inu))
+           tau_aero_lw_rrtm(i,k,1,inu) = 1.e-15  !--test
+           !
+         ENDDO
+      ENDDO
       !
     ENDDO
+    ! 
+   ENDIF
+   !
+  ELSE !--no aerosol LW effects
     !
-    !--avoid very small values
-    tau_aero_lw_rrtm = MAX(tau_aero_lw_rrtm,1.e-15)
-    ! 
-  ELSE !--not aerosol_couple
-    !
-    !--no LW effects if not coupled to INCA
     tau_aero_lw_rrtm = 1.e-15
   ENDIF