Changeset 2542 for LMDZ5/branches/testing/libf/phylmd/rrtm

Timestamp:

Jun 6, 2016, 4:04:57 PM (8 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r2487:2541 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/rrtm/aeropt_lw_rrtm.F90 (modified) (1 diff)
• libf/phylmd/rrtm/read_rsun_rrtm.F90 (copied) (copied from LMDZ5/trunk/libf/phylmd/rrtm/read_rsun_rrtm.F90)
• libf/phylmd/rrtm/readaerosol_optic_rrtm.F90 (modified) (1 diff)
• libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90 (copied) (copied from LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato1_rrtm.F90)
• libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90 (copied) (copied from LMDZ5/trunk/libf/phylmd/rrtm/readaerosolstrato2_rrtm.F90)
• libf/phylmd/rrtm/readaerosolstrato_rrtm.F90 (deleted)
• libf/phylmd/rrtm/recmwf_aero.F90 (modified) (6 diffs)

LMDZ5/branches/testing/libf/phylmd/rrtm/aeropt_lw_rrtm.F90

@@ -1 +1 @@
 !
 ! aeropt_lw_rrtm.F90 2014-05-13 C. Kleinschmitt
+!                    2016-05-03 O. Boucher 
 !
+! This routine feeds aerosol LW properties to RRTM
+! we only consider absorption (not scattering)
 
-SUBROUTINE AEROPT_LW_RRTM
+SUBROUTINE AEROPT_LW_RRTM(aerosol_couple,paprs,tr_seri)
 
+  USE dimphy
+  USE aero_mod
+  USE infotrac_phy
   USE phys_state_var_mod, ONLY : tau_aero_lw_rrtm
+  USE YOERAD, ONLY : NLW
 
   IMPLICIT NONE
 
-  tau_aero_lw_rrtm(:,:,:,:) = 1.e-15
+  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
+  ! 
+  !
+  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
 
+  !
+  !--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  /
+  !
+  IF (NLW.NE.nbands_lw_rrtm) THEN
+     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
+    !
+    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
+      !
+    ENDDO
+    !
+    !--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
+  !
 END SUBROUTINE AEROPT_LW_RRTM

LMDZ5/branches/testing/libf/phylmd/rrtm/readaerosol_optic_rrtm.F90

@@ -19 +19 @@
        load_tmp1,load_tmp2,load_tmp3,load_tmp4,load_tmp5,load_tmp6,load_tmp7
 
-  USE infotrac
+  USE infotrac_phy
 
   IMPLICIT NONE

LMDZ5/branches/testing/libf/phylmd/rrtm/recmwf_aero.F90

@@ -161 +161 @@
 IMPLICIT NONE
 INCLUDE "clesphys.h"
-
 
 INTEGER(KIND=JPIM),INTENT(IN)    :: KPROMA 
@@ -306 +305 @@
 ! not ok_ade and not ok_aie = 1 to proceed
 ! therefore the cases have the following corresponding switches
-! 1 = not ok_ade or not ok_aie
-! 2 = ok_aie
-! 3 = ok_ade
+! 1 = not ok_ade and not ok_aie OR not ok_ade and ok_aie and swaero_diag OR ok_ade and not ok_aie and swaero_diag
+! 2 = not ok_ade and ok_aie OR ok_aie and ok_ade and swaero_diag
+! 3 = ok_ade and not ok_aie OR ok_aie and ok_ade and swaero_diag
 ! 4 = ok_ade and ok_aie
 ! 5 = no aerosol feedback wanted or no aerosol at all 
@@ -445 +444 @@
 IF (flag_aerosol .GT. 0 .OR. flag_aerosol_strat) THEN
 
-IF ( .not. ok_ade .or. .not. ok_aie ) THEN
+!--Case 1
+IF ( ( .not. ok_ade .AND. .not. ok_aie ) .OR.             & 
+   & ( .not. ok_ade .AND. ok_aie .AND. swaero_diag ) .OR. & 
+   & ( ok_ade .AND. .not. ok_aie .AND. swaero_diag ) ) THEN
 
 ! natural aerosols for direct and indirect effect 
@@ -483 +485 @@
 ENDIF 
 
-IF (ok_aie) THEN
+!--Case 2
+IF ( ( .not. ok_ade .AND. ok_aie ) .OR. & 
+   & ( ok_ade .AND. ok_aie .AND. swaero_diag ) ) THEN
 
 ! natural aerosols for direct indirect effect 
@@ -521 +525 @@
 ENDIF ! ok_aie      
 
-IF (ok_ade) THEN
+!--Case 3
+IF ( ( ok_ade .AND. .not. ok_aie ) .OR. &
+   & ( ok_ade .AND. ok_aie .AND. swaero_diag ) ) THEN
 
 ! direct effect of total aerosol activated
@@ -559 +565 @@
 ENDIF !-end ok_ade
 
+!--Case 4
 IF (ok_ade .and. ok_aie) THEN