Changeset 4188 for LMDZ6/branches/LMDZ_ECRad/libf/phylmd/radlwsw_m.F90

Timestamp:

Jun 29, 2022, 9:23:14 AM (2 years ago)

Author:

idelkadi

Message:

Implementation of the Ecrad radiative transfer code in the LMD model (continued) :
Integration of aerosols (direct effect)

File:

• LMDZ6/branches/LMDZ_ECRad/libf/phylmd/radlwsw_m.F90 (modified) (5 diffs)

LMDZ6/branches/LMDZ_ECRad/libf/phylmd/radlwsw_m.F90

@@ -78 +78 @@
 ! Besoin pour ECRAD de pctsrf, zmasq, longitude, altitude
 #ifdef CPP_ECRAD
+      USE phys_local_var_mod, ONLY: rhcl, m_allaer
       USE geometry_mod, ONLY: latitude, longitude
       USE phys_state_var_mod, ONLY: pctsrf
@@ -323 +324 @@
                ZQ_SNOW(klon,klev)           ! Snow cloud mass mixing ratio (kg/kg) ?
   REAL(KIND=8) ZAEROSOL_OLD(KLON,6,KLEV), &  ! 
-               ZAEROSOL(KLON,KLEV,naero_tot) !
+               ZAEROSOL(KLON,KLEV,naero_grp) !
 ! OUTPUTS
   REAL(KIND=8) ZFLUX_DIR(klon), &           ! Direct compt of surf flux into horizontal plane
@@ -1181 +1182 @@
 !
 ! AI ATTENTION Aerosols A REVOIR
-!      DO i = 1, kdlon
-!      DO k = 1, kflev
-!      DO kk= 1, naero_tot 
+      DO i = 1, kdlon
+      DO k = 1, kflev
+      DO kk= 1, naero_grp
 !      DO kk=1, NSW
 !
@@ -1194 +1195 @@
 !      PCGA_NAT(i,kflev+1-k,kk)=cg_aero_sw_rrtm(i,k,1,kk)
 !       ZAEROSOL(i,kflev+1-k,kk)=m_allaer(i,k,kk)
-!
-!      ENDDO
-!      ENDDO
-!      ENDDO
+       ZAEROSOL(i,kflev+1-k,kk)=m_allaer(i,k,kk)
+!
+      ENDDO
+      ENDDO
+      ENDDO
 !-end OB
 !
@@ -1345 +1347 @@
  
       CALL RADIATION_SCHEME &
-      & (ist, iend, klon, klev, naero_tot, NSW, &
+      & (ist, iend, klon, klev, naero_grp, NSW, &
 ! ??? naero_tot
       & day_cur, current_time, &