Changeset 2231 for LMDZ5/trunk/libf

Timestamp:

Mar 12, 2015, 5:46:56 PM (9 years ago)

Author:

oboucher

Message:

Putting minimum values for SW and LW aerosol optical depth values
Putting default value for aerosol single scattering albedo if no aerosol
Doing this consistently in all cases covered by flag_aerosol and flag_aerosol_strat

Location:

LMDZ5/trunk/libf/phylmd

Files:

• physiq.F90 (modified) (1 diff)
• rrtm/aeropt_6bands_rrtm.F90 (modified) (4 diffs)
• rrtm/aeropt_lw_rrtm.F90 (modified) (1 diff)
• rrtm/readaerosolstrato_rrtm.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -3026 +3026 @@
         tausum_aero(:,:,:) = 0.
         IF (iflag_rrtm .EQ. 0) THEN !--old radiation
-           tau_aero(:,:,:,:) = 0.
-           piz_aero(:,:,:,:) = 0.
+           tau_aero(:,:,:,:) = 1.e-15
+           piz_aero(:,:,:,:) = 1.
            cg_aero(:,:,:,:)  = 0.
         ELSE
-           tau_aero_sw_rrtm(:,:,:,:)=0.0
-           piz_aero_sw_rrtm(:,:,:,:)=0.0
-           cg_aero_sw_rrtm(:,:,:,:)=0.0
+           tau_aero_sw_rrtm(:,:,:,:) = 1.e-15
+           tau_aero_lw_rrtm(:,:,:,:) = 1.e-15
+           piz_aero_sw_rrtm(:,:,:,:) = 1.0
+           cg_aero_sw_rrtm(:,:,:,:)  = 0.0
         ENDIF
      ENDIF

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

@@ -473 +473 @@
                                tau_ae(i,k,id_ASSSM_phy,inu)+tau_ae(i,k,id_CSSSM_phy,inu)+   &
                                tau_ae(i,k,id_SSSSM_phy,inu)+ tau_ae(i,k,id_CIDUSTM_phy,inu)
-         tau_allaer(i,k,2,inu)=MAX(tau_allaer(i,k,2,inu),1e-5)
+         tau_allaer(i,k,2,inu)=MAX(tau_allaer(i,k,2,inu),1e-15)
 
          piz_allaer(i,k,2,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ &
@@ -486 +486 @@
                                 tau_ae(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) &
                                 /tau_allaer(i,k,2,inu)
-         piz_allaer(i,k,2,inu)=MAX(piz_allaer(i,k,2,inu),0.1)
+         piz_allaer(i,k,2,inu)=MAX(piz_allaer(i,k,2,inu),0.01)
 
          cg_allaer(i,k,2,inu)=(tau_ae(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)*cg_ae(i,k,id_ASSO4M_phy,inu)+ &
@@ -506 +506 @@
                                tau_ae_pi(i,k,id_ASSSM_phy,inu)+tau_ae_pi(i,k,id_CSSSM_phy,inu)+   &
                                tau_ae_pi(i,k,id_SSSSM_phy,inu)+ tau_ae_pi(i,k,id_CIDUSTM_phy,inu)
-         tau_allaer(i,k,1,inu)=MAX(tau_allaer(i,k,1,inu),1e-5)
+         tau_allaer(i,k,1,inu)=MAX(tau_allaer(i,k,1,inu),1e-15)
 
          piz_allaer(i,k,1,inu)=(tau_ae_pi(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)+ &
@@ -519 +519 @@
                                 tau_ae_pi(i,k,id_CIDUSTM_phy,inu)*piz_ae(i,k,id_CIDUSTM_phy,inu)) &
                                 /tau_allaer(i,k,1,inu)
-         piz_allaer(i,k,1,inu)=MAX(piz_allaer(i,k,1,inu),0.1)
+         piz_allaer(i,k,1,inu)=MAX(piz_allaer(i,k,1,inu),0.01)
 
          cg_allaer(i,k,1,inu)=(tau_ae_pi(i,k,id_ASSO4M_phy,inu)*piz_ae(i,k,id_ASSO4M_phy,inu)*cg_ae(i,k,id_ASSO4M_phy,inu)+ &

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

@@ -9 +9 @@
   IMPLICIT NONE
 
-  tau_aero_lw_rrtm(:,:,:,:)=0.0
+  tau_aero_lw_rrtm(:,:,:,:) = 1.e-15
 
 END SUBROUTINE AEROPT_LW_RRTM

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

@@ -211 +211 @@
     ENDDO
 
+!--default SSA value if there is no aerosol
+!--to avoid 0 values that seems to cause some problem to RRTM
+    WHERE (tau_aero_sw_rrtm.LT.1.e-14)
+      piz_aero_sw_rrtm = 1.0
+    ENDWHERE
+
+!--in principle this should not be necessary 
+!--as these variables have min values already but just in case
+!--put 1e-15 min value to both SW and LW AOD
+    tau_aero_sw_rrtm = MAX(tau_aero_sw_rrtm,1.e-15)
+    tau_aero_lw_rrtm = MAX(tau_aero_lw_rrtm,1.e-15)
+
 end subroutine readaerosolstrato_rrtm