Changeset 3425 for LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/physiq_mod.F90

Timestamp:

Dec 13, 2018, 1:07:06 PM (6 years ago)

Author:

Laurent Fairhead

Message:

Followup to r3408 modifications for VolMIP
NL

File:

• LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/physiq_mod.F90 (modified) (3 diffs)

LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/physiq_mod.F90

@@ -1352 +1352 @@
 
        IF (type_trac == 'inca') THEN
+          ! jg : initialisation jusqu'au ces variables sont dans restart
+          ccm(:,:,:) = 0.
+          tau_aero(:,:,:,:) = 0.
+          piz_aero(:,:,:,:) = 0.
+          cg_aero(:,:,:,:) = 0.
 
           config_inca='none' ! default
@@ -3603 +3608 @@
 #else
                    !--climatologies or INCA aerosols
-                   CALL readaerosol_optic_rrtm( debut, aerosol_couple, ok_alw, &
+                   CALL readaerosol_optic_rrtm( debut, aerosol_couple, ok_alw, ok_volcan, &
                         new_aod, flag_aerosol, flag_bc_internal_mixture, itap, jD_cur-jD_ref, &
                         pdtphys, pplay, paprs, t_seri, rhcl, presnivs,  &
@@ -3693 +3698 @@
              CALL readaerosolstrato1_rrtm(debut)
             ELSEIF (flag_aerosol_strat.EQ.2) THEN
-             CALL readaerosolstrato2_rrtm(debut)
+             CALL readaerosolstrato2_rrtm(debut, ok_volcan)
             ELSE
              abort_message='flag_aerosol_strat must equal 1 or 2 for rrtm=1'