Changeset 2325 for LMDZ5/trunk/libf/phylmd/physiq.F90

Timestamp:

Jul 8, 2015, 10:17:34 PM (9 years ago)

Author:

oboucher

Message:

Enabling aerosols for the case RRTM + NSW=2
In this case we use the old aerosol properties

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -3107 +3107 @@
      !jq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
      IF (flag_aerosol .gt. 0) THEN
-        IF (iflag_rrtm .EQ. 0) THEN !--old radiation
+         IF (iflag_rrtm .EQ. 0) THEN !--old radiation
            IF (.NOT. aerosol_couple) THEN
               !
@@ -3124 +3124 @@
 !
 #ifdef CPP_RRTM
+           IF (NSW.EQ.6) THEN
+!--new aerosol properties
+!
              CALL readaerosol_optic_rrtm( debut, aerosol_couple, &
              new_aod, flag_aerosol, itap, jD_cur-jD_ref, &
@@ -3131 +3134 @@
              tausum_aero, tau3d_aero)
 
-             CALL aeropt_lw_rrtm
+           ELSE IF (NSW.EQ.2) THEN 
+!--for now we use the old aerosol properties
+!
+              CALL readaerosol_optic( &
+                   debut, new_aod, flag_aerosol, itap, jD_cur-jD_ref, &
+                   pdtphys, pplay, paprs, t_seri, rhcl, presnivs,  &
+                   mass_solu_aero, mass_solu_aero_pi,  &
+                   tau_aero, piz_aero, cg_aero,  &
+                   tausum_aero, tau3d_aero)
+!
+                   !--natural aerosols
+                   tau_aero_sw_rrtm(:,:,1,:)=tau_aero(:,:,3,:)
+                   piz_aero_sw_rrtm(:,:,1,:)=piz_aero(:,:,3,:)
+                   cg_aero_sw_rrtm (:,:,1,:)=cg_aero (:,:,3,:)
+                   !--all aerosols
+                   tau_aero_sw_rrtm(:,:,2,:)=tau_aero(:,:,2,:)
+                   piz_aero_sw_rrtm(:,:,2,:)=piz_aero(:,:,2,:)
+                   cg_aero_sw_rrtm (:,:,2,:)=cg_aero (:,:,2,:)
+           ELSE
+              abort_message='Only NSW=2 or 6 are possible with aerosols and iflag_rrtm=1'
+              call abort_physic(modname,abort_message,1)
+           ENDIF
+
+           CALL aeropt_lw_rrtm
+!
 #else
-
-              abort_message='You should compile with -rrtm if running with iflag_rrtm=1'
-              call abort_physic(modname,abort_message,1)
+           abort_message='You should compile with -rrtm if running with iflag_rrtm=1'
+           call abort_physic(modname,abort_message,1)
 #endif
               !