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

Timestamp:

Jun 12, 2015, 11:38:57 PM (10 years ago)

Author:

musat

Message:

RRTM :
Correction bug interface avec rrtm (radlwsw_m.F90)
Ajouter la possibilite d'utilsation d'anciennes proprites optique dans le cas RRTM 2bandes (readaerosol_optic_rrtm.F90)
Correction rrtm (rrtm_rtrn1a_140gp.F90)

COSP :
Modifications pour traiter les valeurs indefines :
mise a 0 en cas d'utilisation de IOIPSL
mise a la valeur recuperee dans .xml en cas d'utilisation de XIOS

Location:

LMDZ5/trunk/libf/phylmd/rrtm

Files:

• readaerosol_optic_rrtm.F90 (modified) (2 diffs)
• rrtm_rtrn1a_140gp.F90 (modified) (1 diff)

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

@@ -303 +303 @@
   END DO
 
+! AI 12 juin 2015
+! Modif Olivier pour prendre en compte des prop optiq aero ancienne pour RRTM 2bds 
+
+IF (NSW.EQ.2) THEN
+!--old aerosol properties
+
+   ! old aeropt_2bands for rrtm
+   CALL aeropt_2bands( &
+         pdel, m_allaer, pdtphys, rhcl, &
+         tau_aero, piz_aero, cg_aero,   &
+         m_allaer_pi, flag_aerosol, &
+         pplay, t_seri, presnivs)
+
+   ! aeropt_5wv only for validation and diagnostics.
+   CALL aeropt_5wv(                    &
+        pdel, m_allaer,                &
+        pdtphys, rhcl, aerindex,       &
+        flag_aerosol, pplay, t_seri,   &
+        tausum_aero, tau3d_aero, presnivs)
+
+ELSEIF (NSW.EQ.6) THEN
+!--new aerosol propertie
   ! aeropt_6bands for rrtm
   CALL aeropt_6bands_rrtm( &
@@ -316 +338 @@
        flag_aerosol, pplay, t_seri,   &
        tausum_aero, tau3d_aero )
+
+ELSE
+
+    PRINT *,'Cas NSW non prevu pour RRTM - NSW=',NSW
+    STOP
+
+ENDIF
 
   ! Diagnostics calculation for CMIP5 protocol

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

@@ -496 +496 @@
 !    Z_FACCMB2D(I_LEV-1) = Z_FACCLD1D(I_LEV-1) * Z_FACCLR2D(I_LEV) *&
 !     & (1.0_JPRB - Z_CLDFRAC(I_LEV+1))  
-    if (istcldd(i_lev).ne.1) then
+    if (istcldd(i_lev).ne.1.and.i_lev.ne.0) then
        z_faccmb1d(i_lev-1) = max(0.,min(z_cldfrac(i_lev+1)-z_cldfrac(i_lev), &
                             z_cldfrac(i_lev-1)-z_cldfrac(i_lev)))