Changeset 5559


Ignore:
Timestamp:
Feb 24, 2025, 3:39:32 PM (7 hours ago)
Author:
lebasn
Message:

StratAer?: clean duplicated module loading, specify loading vars from module + variable initialisation for debug mode.

Location:
LMDZ6/trunk/libf/phylmd
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/StratAer/calcaerosolstrato_rrtm.f90

    r5338 r5559  
    3939  tau_strat_wave=0.0
    4040  tau_lw_abs_rrtm=0.0
    41 
     41 
     42!-- init tau_strat vars
     43  tau_strat_550(:,:) =0.0
     44  tau_strat_1020(:,:)=0.0
     45 
    4246  CALL miecalc_aer(tau_strat, piz_strat, cg_strat, tau_strat_wave, tau_lw_abs_rrtm, paprs, debut)
    43 
    44 !!--test CK: deactivate radiative effect of aerosol
    45 !  tau_strat=0.0
    46 !  piz_strat=0.0
    47 !  cg_strat=0.0
    48 !  tau_strat_wave=0.0
    49 !  tau_lw_abs_rrtm=0.0
    50 
    51 !--test CK: deactivate SW radiative effect of aerosol (but leave LW)
    52 !  tau_strat=0.0
    53 !  piz_strat=0.0
    54 !  cg_strat=0.0
    55 
    56 !  DO wave=1, nwave_sw
    57 !  tau_strat_wave(:,:,wave)=0.0
    58 !  ENDDO
    59 
    60 !--test CK: deactivate LW radiative effect of aerosol (but leave SW)
    61 !  tau_lw_abs_rrtm=0.0
    62 
    63 !  DO wave=nwave_sw+1, nwave_sw+nwave_lw
    64 !  tau_strat_wave(:,:,wave)=0.0
    65 !  ENDDO
    6647
    6748!--total vertical aod at the 5 SW + 1 LW wavelengths

  • LMDZ6/trunk/libf/phylmd/StratAer/interp_sulf_input.f90

    r5338 r5559  
    118118
    119119  IF (is_mpi_root.AND.is_omp_root) THEN
    120 
     120   
     121    OCS_lifetime(:,:)=0.0
     122    SO2_lifetime(:,:)=0.0
     123    H2SO4_lifetime(:,:)=0.0
     124    O3_clim(:,:)=0.0
     125     
    121126    !--init ncdf variables
    122127    IF(flag_newclim_file) THEN
    … …  
    332337
    333338!--regridding tracer concentration on the vertical
     339  budg_3D_backgr_ocs(:,:)=0.0
     340  budg_3D_backgr_so2(:,:)=0.0
     341 
    334342  DO i=1, klon
    335343    DO k=1, klev

  • LMDZ6/trunk/libf/phylmd/StratAer/miecalc_aer.f90

    r5268 r5559  
    227227    50000.000,    0.2000,   1.49800,   1.0000E-08 /), (/nb_lambda_h2so4,4/), order=(/2,1/) )
    228228     
     229    ! init
     230    piz_bin(:,:)=0.0
     231    alpha_bin(:,:)=0.0
     232    cg_bin(:,:)=0.0
     233   
    229234    !--compute particle radius for a composition of 75% H2SO4 / 25% H2O at T=293K
    230235    DO bin_number=1, nbtr_bin
    … …  
    332337
    333338      DO bin=1, Nbin !---loop on size bins
    334 
     339     
    335340      r_lower=exp(log(rmin)+FLOAT(bin-1)/FLOAT(Nbin)*(log(rmax)-log(rmin)))
    336341      r_upper=exp(log(rmin)+FLOAT(bin)/FLOAT(Nbin)*(log(rmax)-log(rmin)))

  • LMDZ6/trunk/libf/phylmd/StratAer/strataer_local_var_mod.f90

    r5268 r5559  
    159159    USE mod_phys_lmdz_para, ONLY : is_master
    160160    USE infotrac_phy, ONLY: id_OCS_strat,id_SO2_strat,id_H2SO4_strat,nbtr_sulgas,nbtr_bin
    161     USE phys_local_var_mod, ONLY : mdw
     161    USE phys_local_var_mod, ONLY : mdw,R2SO4,R2SO4B,DENSO4,DENSO4B,f_r_wet,f_r_wetB
    162162    USE aerophys, ONLY: mdwmin, V_rat
    163163    USE yomcst_mod_h  , ONLY : RPI
    … …  
    205205    nAerErupt = 1 ; nSpeciesErupt = 1
    206206    ifreqroc=2 ; flh2o=0
     207   
     208    ! array init
     209    mdw(:)=0.
     210    R2SO4(:,:)=0.
     211    R2SO4B(:,:,:)=0.
     212    DENSO4(:,:)=0.
     213    DENSO4B(:,:,:)=0.
     214    f_r_wet(:,:)=0.
     215    f_r_wetB(:,:,:)=0.
    207216   
    208217    !============= Read params =============

  • LMDZ6/trunk/libf/phylmd/StratAer/sulfate_aer_mod.f90

    r5268 r5559  
    2727      USE dimphy, ONLY : klon,klev ! nb of longitude and altitude bands
    2828      USE infotrac_phy, ONLY : nbtr_bin
    29       USE aerophys
     29      USE aerophys, ONLY : mAIRmol,mH2Omol,dens_aer_dry,rgas
    3030      USE phys_local_var_mod, ONLY: R2SO4, R2SO4B, DENSO4, DENSO4B, f_r_wet, f_r_wetB
    3131      USE strataer_local_var_mod, ONLY: RRSI
    3232!     WARNING: in phys_local_var_mod R2SO4B, DENSO4B, f_r_wetB (klon,klev,nbtr_bin)
    3333!          and dens_aer_dry must be declared somewhere
     34      USE print_control_mod, ONLY : lunout
    3435   
    3536      IMPLICIT NONE
    … …  
    9091!          Loop on bin radius (RRSI in cm)
    9192           DO IK=1,nbtr_bin
    92  
     93
    9394!      ***   H2SO4-H2O curved surface - Kelvin effect factor ***
    9495!            wet radius (m) (RRSI(IK) in [cm])

  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r5554 r5559  
    114114                        ptrop, ttrop, ztrop, gravit, itroprep, Z1, Z2, fac, B, chemini_rep, chemtime_rep, coord_hyb_rep, &
    115115            rtime
    116     USE strataer_local_var_mod
    117     USE strataer_emiss_mod, ONLY: strataer_emiss_init
    118116    USE time_phylmdz_mod,    ONLY: annee_ref, day_ini, day_ref, start_time
    119117    USE vertical_layers_mod, ONLY: aps, bps, ap, bp
    … …  
    126124
    127125    USE phys_local_var_mod, ONLY: d_q_emiss
    128     USE strataer_local_var_mod
     126    USE strataer_local_var_mod, ONLY: strataer_init,flag_emit,flh2o,ok_qemiss,flag_verbose_strataer, &
     127         year_emit_vol,mth_emit_vol,day_emit_vol,nErupt,nAerErupt,injdur,m_H2O_emiss_vol_daily,m_H2O_emiss_vol, &
     128         ponde_lonlat_vol,altemiss_vol,sigma_alt_vol,xlat_min_vol,xlat_max_vol,xlon_min_vol,xlon_max_vol
    129129    USE strataer_nuc_mod, ONLY: strataer_nuc_init
    130130    USE strataer_emiss_mod, ONLY: strataer_emiss_init
    131 
     131   
    132132    USE lmdz_xios, ONLY: xios_update_calendar, xios_context_finalize
    133133    USE lmdz_xios, ONLY: xios_get_field_attr, xios_field_is_active, xios_context
Note: See TracChangeset for help on using the changeset viewer.