Index: LMDZ6/trunk/libf/phylmd/StratAer/calcaerosolstrato_rrtm.f90 =================================================================== --- LMDZ6/trunk/libf/phylmd/StratAer/calcaerosolstrato_rrtm.f90 (revision 5557) +++ LMDZ6/trunk/libf/phylmd/StratAer/calcaerosolstrato_rrtm.f90 (revision 5559) @@ -39,29 +39,10 @@ tau_strat_wave=0.0 tau_lw_abs_rrtm=0.0 - + +!-- init tau_strat vars + tau_strat_550(:,:) =0.0 + tau_strat_1020(:,:)=0.0 + CALL miecalc_aer(tau_strat, piz_strat, cg_strat, tau_strat_wave, tau_lw_abs_rrtm, paprs, debut) - -!!--test CK: deactivate radiative effect of aerosol -! tau_strat=0.0 -! piz_strat=0.0 -! cg_strat=0.0 -! tau_strat_wave=0.0 -! tau_lw_abs_rrtm=0.0 - -!--test CK: deactivate SW radiative effect of aerosol (but leave LW) -! tau_strat=0.0 -! piz_strat=0.0 -! cg_strat=0.0 - -! DO wave=1, nwave_sw -! tau_strat_wave(:,:,wave)=0.0 -! ENDDO - -!--test CK: deactivate LW radiative effect of aerosol (but leave SW) -! tau_lw_abs_rrtm=0.0 - -! DO wave=nwave_sw+1, nwave_sw+nwave_lw -! tau_strat_wave(:,:,wave)=0.0 -! ENDDO !--total vertical aod at the 5 SW + 1 LW wavelengths Index: LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.f90 =================================================================== --- LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.f90 (revision 5557) +++ LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.f90 (revision 5559) @@ -118,5 +118,10 @@ IF (is_mpi_root.AND.is_omp_root) THEN - + + OCS_lifetime(:,:)=0.0 + SO2_lifetime(:,:)=0.0 + H2SO4_lifetime(:,:)=0.0 + O3_clim(:,:)=0.0 + !--init ncdf variables IF(flag_newclim_file) THEN @@ -332,4 +337,7 @@ !--regridding tracer concentration on the vertical + budg_3D_backgr_ocs(:,:)=0.0 + budg_3D_backgr_so2(:,:)=0.0 + DO i=1, klon DO k=1, klev Index: LMDZ6/trunk/libf/phylmd/StratAer/miecalc_aer.f90 =================================================================== --- LMDZ6/trunk/libf/phylmd/StratAer/miecalc_aer.f90 (revision 5557) +++ LMDZ6/trunk/libf/phylmd/StratAer/miecalc_aer.f90 (revision 5559) @@ -227,4 +227,9 @@ 50000.000, 0.2000, 1.49800, 1.0000E-08 /), (/nb_lambda_h2so4,4/), order=(/2,1/) ) + ! init + piz_bin(:,:)=0.0 + alpha_bin(:,:)=0.0 + cg_bin(:,:)=0.0 + !--compute particle radius for a composition of 75% H2SO4 / 25% H2O at T=293K DO bin_number=1, nbtr_bin @@ -332,5 +337,5 @@ DO bin=1, Nbin !---loop on size bins - + r_lower=exp(log(rmin)+FLOAT(bin-1)/FLOAT(Nbin)*(log(rmax)-log(rmin))) r_upper=exp(log(rmin)+FLOAT(bin)/FLOAT(Nbin)*(log(rmax)-log(rmin))) Index: LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.f90 =================================================================== --- LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.f90 (revision 5557) +++ LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.f90 (revision 5559) @@ -159,5 +159,5 @@ USE mod_phys_lmdz_para, ONLY : is_master USE infotrac_phy, ONLY: id_OCS_strat,id_SO2_strat,id_H2SO4_strat,nbtr_sulgas,nbtr_bin - USE phys_local_var_mod, ONLY : mdw + USE phys_local_var_mod, ONLY : mdw,R2SO4,R2SO4B,DENSO4,DENSO4B,f_r_wet,f_r_wetB USE aerophys, ONLY: mdwmin, V_rat USE yomcst_mod_h , ONLY : RPI @@ -205,4 +205,13 @@ nAerErupt = 1 ; nSpeciesErupt = 1 ifreqroc=2 ; flh2o=0 + + ! array init + mdw(:)=0. + R2SO4(:,:)=0. + R2SO4B(:,:,:)=0. + DENSO4(:,:)=0. + DENSO4B(:,:,:)=0. + f_r_wet(:,:)=0. + f_r_wetB(:,:,:)=0. !============= Read params ============= Index: LMDZ6/trunk/libf/phylmd/StratAer/sulfate_aer_mod.f90 =================================================================== --- LMDZ6/trunk/libf/phylmd/StratAer/sulfate_aer_mod.f90 (revision 5557) +++ LMDZ6/trunk/libf/phylmd/StratAer/sulfate_aer_mod.f90 (revision 5559) @@ -27,9 +27,10 @@ USE dimphy, ONLY : klon,klev ! nb of longitude and altitude bands USE infotrac_phy, ONLY : nbtr_bin - USE aerophys + USE aerophys, ONLY : mAIRmol,mH2Omol,dens_aer_dry,rgas USE phys_local_var_mod, ONLY: R2SO4, R2SO4B, DENSO4, DENSO4B, f_r_wet, f_r_wetB USE strataer_local_var_mod, ONLY: RRSI ! WARNING: in phys_local_var_mod R2SO4B, DENSO4B, f_r_wetB (klon,klev,nbtr_bin) ! and dens_aer_dry must be declared somewhere + USE print_control_mod, ONLY : lunout IMPLICIT NONE @@ -90,5 +91,5 @@ ! Loop on bin radius (RRSI in cm) DO IK=1,nbtr_bin - + ! *** H2SO4-H2O curved surface - Kelvin effect factor *** ! wet radius (m) (RRSI(IK) in [cm]) Index: LMDZ6/trunk/libf/phylmd/physiq_mod.F90 =================================================================== --- LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (revision 5557) +++ LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (revision 5559) @@ -114,6 +114,4 @@ ptrop, ttrop, ztrop, gravit, itroprep, Z1, Z2, fac, B, chemini_rep, chemtime_rep, coord_hyb_rep, & rtime - USE strataer_local_var_mod - USE strataer_emiss_mod, ONLY: strataer_emiss_init USE time_phylmdz_mod, ONLY: annee_ref, day_ini, day_ref, start_time USE vertical_layers_mod, ONLY: aps, bps, ap, bp @@ -126,8 +124,10 @@ USE phys_local_var_mod, ONLY: d_q_emiss - USE strataer_local_var_mod + USE strataer_local_var_mod, ONLY: strataer_init,flag_emit,flh2o,ok_qemiss,flag_verbose_strataer, & + year_emit_vol,mth_emit_vol,day_emit_vol,nErupt,nAerErupt,injdur,m_H2O_emiss_vol_daily,m_H2O_emiss_vol, & + ponde_lonlat_vol,altemiss_vol,sigma_alt_vol,xlat_min_vol,xlat_max_vol,xlon_min_vol,xlon_max_vol USE strataer_nuc_mod, ONLY: strataer_nuc_init USE strataer_emiss_mod, ONLY: strataer_emiss_init - + USE lmdz_xios, ONLY: xios_update_calendar, xios_context_finalize USE lmdz_xios, ONLY: xios_get_field_attr, xios_field_is_active, xios_context