Changeset 2738 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Jan 2, 2017, 11:24:30 PM (7 years ago)

Author:

oboucher

Message:

Swapping the order of CSSO4 and ASSO4 aerosols (and fixing an issue on aerindex). Preparing the ground for nitrate aerosols (coarse soluble, accumulation soluble, coarse insoluble). Modifying the LW aeropt routine so that it is compatible with both INCA and climatological aerosols (for dust only). Adding a new flag ok_alw for activating aerosol direct LW effect (for dust only). This change is bit comparable for flag_aerosol=6, flag_rrtm=1, NSW=6.

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (9 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -980 +980 @@
     ! Parameters
     LOGICAL ok_ade, ok_aie    ! Apply aerosol (in)direct effects or not
+    LOGICAL ok_alw            ! Apply aerosol LW effect or not
     LOGICAL ok_cdnc ! ok cloud droplet number concentration (O. Boucher 01-2013)
     REAL bl95_b0, bl95_b1   ! Parameter in Boucher and Lohmann (1995)
-    SAVE ok_ade, ok_aie, ok_cdnc, bl95_b0, bl95_b1
-    !$OMP THREADPRIVATE(ok_ade, ok_aie, ok_cdnc, bl95_b0, bl95_b1)
+    SAVE ok_ade, ok_aie, ok_alw, ok_cdnc, bl95_b0, bl95_b1
+    !$OMP THREADPRIVATE(ok_ade, ok_aie, ok_alw, ok_cdnc, bl95_b0, bl95_b1)
     LOGICAL, SAVE :: aerosol_couple ! true  : calcul des aerosols dans INCA
     ! false : lecture des aerosol dans un fichier
@@ -1142 +1143 @@
             fact_cldcon, facttemps,ok_newmicro,iflag_radia, &
             iflag_cld_th,iflag_ratqs,ratqsbas,ratqshaut,tau_ratqs, &
-            ok_ade, ok_aie, ok_cdnc, aerosol_couple,  &
+            ok_ade, ok_aie, ok_alw, ok_cdnc, aerosol_couple,  &
             flag_aerosol, flag_aerosol_strat, new_aod, &
             flag_bc_internal_mixture, bl95_b0, bl95_b1, &
@@ -1587 +1588 @@
             jjmp1,nlevSTD,clevSTD,rlevSTD, dtime,ok_veget, &
             type_ocean,iflag_pbl,iflag_pbl_split,ok_mensuel,ok_journe, &
-            ok_hf,ok_instan,ok_LES,ok_ade,ok_aie,  &
+            ok_hf,ok_instan,ok_LES,ok_ade,ok_aie, &
             read_climoz, phys_out_filestations, &
             new_aod, aerosol_couple, &
@@ -3293 +3294 @@
        !jq - introduce the aerosol direct and first indirect radiative forcings
        !jq - Johannes Quaas, 27/11/2003 (quaas@lmd.jussieu.fr)
-       IF (flag_aerosol .gt. 0) THEN
+       IF (flag_aerosol .GT. 0) THEN
           IF (iflag_rrtm .EQ. 0) THEN !--old radiation
              IF (.NOT. aerosol_couple) THEN
@@ -3312 +3313 @@
 #ifdef CPP_RRTM
                 IF (NSW.EQ.6) THEN
-                   !--new aerosol properties
+                   !--new aerosol properties SW and LW
                    !
-                   CALL readaerosol_optic_rrtm( debut, aerosol_couple, &
+                   CALL readaerosol_optic_rrtm( debut, aerosol_couple, ok_alw, &
                         new_aod, flag_aerosol, flag_bc_internal_mixture, itap, jD_cur-jD_ref, &
                         pdtphys, pplay, paprs, t_seri, rhcl, presnivs,  &
@@ -3320 +3321 @@
                         tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm,  &
                         tausum_aero, tau3d_aero)
-
+                   !
                 ELSE IF (NSW.EQ.2) THEN 
                    !--for now we use the old aerosol properties
@@ -3339 +3340 @@
                    piz_aero_sw_rrtm(:,:,2,:)=piz_aero(:,:,2,:)
                    cg_aero_sw_rrtm (:,:,2,:)=cg_aero (:,:,2,:)
+                   !
+                   !--no LW optics
+                   tau_aero_lw_rrtm(:,:,:,:) = 1.e-15
+                   !
                 ELSE
                    abort_message='Only NSW=2 or 6 are possible with ' &
@@ -3344 +3349 @@
                    CALL abort_physic(modname,abort_message,1)
                 ENDIF
-
-                !--call LW optical properties for tropospheric aerosols 
-                !--only works for INCA aerosol (aerosol_couple = TRUE)
-                CALL aeropt_lw_rrtm(aerosol_couple,paprs,tr_seri)
-                !
 #else
                 abort_message='You should compile with -rrtm if running ' &
@@ -3357 +3357 @@
              ENDIF
           ENDIF
-       ELSE
+       ELSE   !--flag_aerosol = 0 
           tausum_aero(:,:,:) = 0.
           mass_solu_aero(:,:) = 0.