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

Timestamp:

Apr 14, 2017, 4:42:31 PM (7 years ago)

Author:

oboucher

Message:

Introducing dry AOD diagnostics for the total aerosols and specieswise
The calculations are only performed if the diagnostics are requested

File:

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

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -92 +92 @@
        topswcf_aero,solswcf_aero,   &
        tausum_aero,tau3d_aero,      &
+       drytausum_aero,              &
        !
        !variables CFMIP2/CMIP5
@@ -1580 +1581 @@
            !!!  xios_field_is_active("toplwai").OR.xios_field_is_active("sollwai") !))) & 
            swaero_diag=.TRUE. 
+
+       !--setting up dryaod_diag to TRUE in XIOS case 
+       DO naero = 1, naero_tot-1
+         IF (xios_field_is_active("dryod550_"//name_aero_tau(naero)) dryaod_diag=.TRUE. 
+       ENDDO
 #endif 
 
@@ -3371 +3377 @@
                         tr_seri, mass_solu_aero, mass_solu_aero_pi,  &
                         tau_aero_sw_rrtm, piz_aero_sw_rrtm, cg_aero_sw_rrtm,  &
-                        tausum_aero, tau3d_aero)
+                        tausum_aero, drytausum_aero, tau3d_aero)
 #endif
                    !
@@ -3411 +3417 @@
        ELSE   !--flag_aerosol = 0 
           tausum_aero(:,:,:) = 0.
+          drytausum_aero(:,:) = 0.
           mass_solu_aero(:,:) = 0.
           mass_solu_aero_pi(:,:) = 0.
@@ -3678 +3685 @@
 #ifndef CPP_XIOS
           !--OB 30/05/2016 modified 21/10/2016
-          !--here we return swaero_diag to FALSE
+          !--here we return swaero_diag and dryaod_diag to FALSE
           !--and histdef will switch it back to TRUE if necessary
           !--this is necessary to get the right swaero at first step
           !--but only in the case of no XIOS as XIOS is covered elsewhere 
           IF (debut) swaero_diag = .FALSE.
+          IF (debut) dryaod_diag = .FALSE.
 #endif
           !