Description of variables and model setup used in 'Impact of Arctic stratospheric polar vortex on the Mediterranean precipitation' by Zhang et al., submitted to JC, 2023. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Description of variables: cyclonelabel.nc: monthly cyclone frequency (%) in February for each member ensemble in the SPV (cyclonelabel_SPV) and WPV (cyclonelabel_WPV) run. frontogenesis.nc: monthly frontogenesis function (K/100 km/day) in February for each member ensemble in the SPV (frontogenesis_SPV) and WPV (frontogenesis_WPV) run. TNflux.nc: TN flux for quasi-stationary waves associated with the WPV in x (fx) and y (fy) direction (m**2/s**2). hadv_vort.nc: monthly cyclone horizontal geostrophic absolute vorticity advection (1/s/day) in February for each member ensemble in the SPV (hadv_vort_SPV) and WPV (hadv_vort_WPV) run. t_hadv.nc: monthly vertical motion contributed by horizontal temperature advection (K/day) in February for each member ensemble in the SPV (t_hadv_SPV) and WPV (t_hadv_WPV) run. omega_tadv.nc: monthly vertical motion contributed by horizontal temperature advection (Pa/s) in February for each member ensemble in the SPV (omega_tadv_SPV) and WPV (omega_tadv_SPV) run. omega_vradv.nc: monthly vertical motion contributed by horizontal absolute vorticity advection (Pa/s) in February for each member ensemble in the SPV (omega_vradv_SPV) and WPV (omega_vradv_SPV) run. omega_Q.nc: monthly vertical motion contributed by diabatic processes (Pa/s) in February for each member ensemble in the SPV (omega_Q_SPV) and WPV (omega_Q_SPV) run. U.nc: monthly zonal wind (m/s) in February for each member ensemble in the SPV (U_SPV) and WPV (U_WPV) run. V.nc: monthly meridional wind (m/s) in February for each member ensemble in the SPV (V_SPV) and WPV (V_WPV) run. T.nc: monthly temperature (K) in February for each member ensemble in the SPV (T_SPV) and WPV (T_WPV) run. Z3.nc: monthly geopotential height (gpm) in February for each member ensemble in the SPV (Z3_SPV) and WPV (Z3_WPV) run. OMEGA.nc: monthly vertical motion (Pa/s) in February for each member ensemble in the SPV (OMEGA_SPV) and WPV (OMEGA_WPV) run. PRECT.nc: monthly total precipitation rate (m/s) in February for each member ensemble in the SPV (PRECT_SPV) and WPV (PRECT_WPV) run. EPdays.nc: monthly extreme precipitation days (count of days) in February for each member ensemble in the SPV (EPdays_SPV) and WPV (EPdays_WPV) run. uqint.nc: monthly integrated water flux in x direction (kg/m/s) in February for each member ensemble in the SPV (uqint_SPV) and WPV (uqint_WPV) run. vqint.nc: monthly integrated water flux in y direction (kg/m/s) in February for each member ensemble in the SPV (vqint_SPV) and WPV (vqint_WPV) run. uvqintdiv.nc: monthly integrated water flux divergence (kg/m**2/s) in February for each member ensemble in the SPV (uvqintdiv_SPV) and WPV (uvqintdiv_WPV) run. uvadv.nc: monthly integrated horizontal moisture advection (kg/m**2/s) in February for each member ensemble in the SPV (uvadv_SPV) and WPV (uvadv_WPV) run. wadv.nc: monthly integrated vertical moisture advection (kg/m**2/s) in February for each member ensemble in the SPV (wadv_SPV) and WPV (wadv_WPV) run. waadvqc.nc: monthly dynamic component of the vertical moisture advection (kg/m**2/s) in February for each member ensemble in the SPV (waadvqc_SPV) and WPV (waadvqc_WPV) run. wcadvqa.nc: monthly thermodynamic component of the vertical moisture advection (kg/m**2/s) in February for each member ensemble in the SPV (wcadvqa_SPV) and WPV (wcadvqa_WPV) run. waadvqc_omegaterms.nc: monthly dynamic component of the vertical moisture advection (kg/m**2/s) induced by (waadvqc_omega_vradvWPVSPVdiff) horizontal vorticity advection, (waadvqc_omega_tadvWPVSPVdiff) temperature advection , and (waadvqc_omega_QWPVSPVdiff) diabatic heating between WPV and SPV runs. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! WACCM4 model in cesm version 1.2.2 is used in our study. To impose thermal forcing in January, we mainly added the following code in line 1901 of the physpkg.F90 file: !!!for SPV run: allocate(ptend%localheating(state%psetcols,pver), stat=ierr) ptend%localheating = 0._r8 if ( ierr /= 0 ) call endrun('physics_ptend_alloc error: allocation error for ptend%localheating') calday = get_curr_calday(); if (calday <= 20 ) then do i = 1,ncol stat_lat=state%lat(i)*180.0/3.1415926 !stat_lon=state%lon(i)*180.0/3.1415926 do k = 1,pver if (hypm(k).gt.15000.or. hypm(k).lt.100. )then else if (hypm(k).gt.6000.)then ptend%localheating(i,k)=-2./2.*(1+tanh((stat_lat-60.)/5.))*(hypm(k)-15000.)/(6000.-15000.) else if (hypm(k).lt.1000.)then ptend%localheating(i,k)=-2./2.*(1+tanh((stat_lat-60.)/5.))*(100.-hypm(k))/(100.-1000.) else ! if (hypm(k).le.6000.and.hypm(k).gt.1000.)then ptend%localheating(i,k)=-2./2.*(1+tanh((stat_lat-60.)/5.)) end if end do end do end if call outfld('localheating', ptend%localheating, pcols, lchnk) ptend%localheating=ptend%localheating*cpair/86400. !K/day ->1 J/kg/s =1004/(60*60*24) K/day = 1004 *K/s ptend%s(:ncol,:)=ptend%s(:ncol,:)+ptend%localheating !!!for WPV run: allocate(ptend%localheating(state%psetcols,pver), stat=ierr) ptend%localheating = 0._r8 if ( ierr /= 0 ) call endrun('physics_ptend_alloc error: allocation error for ptend%localheating') calday = get_curr_calday(); if (calday <= 20 ) then do i = 1,ncol stat_lat=state%lat(i)*180.0/3.1415926 !stat_lon=state%lon(i)*180.0/3.1415926 do k = 1,pver if (hypm(k).gt.15000.or. hypm(k).lt.100. )then else if (hypm(k).gt.6000.)then ptend%localheating(i,k)=2./2.*(1+tanh((stat_lat-60.)/5.))*(hypm(k)-15000.)/(6000.-15000.) else if (hypm(k).lt.1000.)then ptend%localheating(i,k)=2./2.*(1+tanh((stat_lat-60.)/5.))*(100.-hypm(k))/(100.-1000.) else ! if (hypm(k).le.6000.and.hypm(k).gt.1000.)then ptend%localheating(i,k)=2./2.*(1+tanh((stat_lat-60.)/5.)) end if end do end do end if call outfld('localheating', ptend%localheating, pcols, lchnk) ptend%localheating=ptend%localheating*cpair/86400. !K/day ->1 J/kg/s =1004/(60*60*24) K/day = 1004 *K/s ptend%s(:ncol,:)=ptend%s(:ncol,:)+ptend%localheating !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Please contact the authors for more details. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! The namelist settings in our study shown below are identical for all 80 ensemble members, with the exception of the initial value file (specified by ncdata): &aerodep_flx_nl aerodep_flx_cycle_yr = 2000 aerodep_flx_datapath = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/chem/trop_mozart_aero/aero' aerodep_flx_file = 'aerosoldep_monthly_1849-2006_1.9x2.5_c090803.nc' aerodep_flx_specifier = 'BCDEPWET', 'BCPHODRY', 'BCPHIDRY', 'OCDEPWET', 'OCPHODRY', 'OCPHIDRY', 'DSTX01DD', 'DSTX02DD', 'DSTX03DD', 'DSTX04DD', 'DSTX01WD', 'DSTX02WD', 'DSTX03WD', 'DSTX04WD' aerodep_flx_type = 'CYCLICAL' / &cam_inparm absems_data = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/rad/abs_ems_factors_fastvx.c030508.nc' avgflag_pertape = 'A','A' bnd_topo = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/topo/USGS-gtopo30_1.9x2.5_remap_c050602.nc' dtime = 1800 empty_htapes = .TRUE. fincl1 = 'SST', 'CLDTOT', 'ICEFRAC', 'FSNS', 'FLNS', 'QREFHT', 'TREFHT', 'PHIS', 'Z3', 'U', 'V', 'T', 'OMEGA','QRL','QRS', 'Q', 'PS', 'PSL', 'TS', 'PRECC', 'PRECL', 'PRECT', 'PRECSL', 'LHFLX', 'SHFLX', 'localheating' fincl2 = "PSL","Z3" fincl3 = '' fincl4 = '' fincl5 = '' inithist = 'NONE' iradsw = -1 mfilt = 365,999 ncdata= '/home/zhangzhy2019/soft/cesm1_2_22/i/amip_ctl.cam.i.0019-01-01-00000.nc' nhtfrq = -24,-6 nlte_use_mo = .false. phys_loadbalance = 2 print_energy_errors = .false. qbo_cyclic = .false. qbo_forcing_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/waccm/qbo/qbocyclic28months.nc' qbo_use_forcing = .false. / &carma_nl carma_model = 'none' / &chem_inparm aer_wetdep_list = '' bndtvg = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/ggas/noaamisc.r8.nc' clim_soilw_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/chem/trop_mozart/dvel/clim_soilw.nc' depvel_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/chem/trop_mozart/dvel/depvel_monthly.nc' depvel_lnd_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/chem/trop_mozart/dvel/regrid_vegetation.nc' ghg_chem = .true. h2orates = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/waccm/phot/xh2o_c080826.nc' season_wes_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/chem/trop_mozart/dvel/season_wes.nc' sol_facti_cloud_borne = 1.0D0 solar_parms_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/waccm/phot/wa_smax_c100517.nc' / &chem_surfvals_nl co2vmr = 368.9e-6 flbc_cycle_yr = 2000 flbc_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/waccm/lb/LBC_1765-2005_1.9x2.5_CMIP5_za_c091204.nc' flbc_type = 'CYCLICAL' / &cldfrc_nl cldfrc_dp1 = 0.10D0 cldfrc_dp2 = 500.0D0 cldfrc_freeze_dry = .true. cldfrc_ice = .false. cldfrc_icecrit = 0.95D0 cldfrc_iceopt = 1 cldfrc_premib = 750.0D2 cldfrc_premit = 75000.0D0 cldfrc_rhminh = 0.800D0 cldfrc_rhminl = 0.90D0 cldfrc_rhminl_adj_land = 0.100D0 cldfrc_sh1 = 0.04D0 cldfrc_sh2 = 500.0D0 / &cldsed_nl cldsed_ice_stokes_fac = 0.5D0 / &cldwat_nl cldwat_conke = 5.0e-6 cldwat_icritc = 9.5e-6 cldwat_icritw = 2.0e-4 cldwat_r3lcrit = 10.0e-6 / &clubb_his_nl clubb_history = .false. clubb_rad_history = .false. / &dyn_fv_inparm del2coef = 3.e+5 div24del2flag = 2 fft_flt = 1 filtcw = 0 nsplit = 8 nspltrac = 2 nspltvrm = 2 / &gw_drag_nl gw_drag_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/waccm/gw/newmfspectra40_dc25.nc' / &hkconv_nl hkconv_c0 = 1.0e-4 hkconv_cmftau = 1800.0D0 / &phys_ctl_nl cam_chempkg = 'waccm_ghg' cam_physpkg = 'cam4' conv_water_in_rad = 0 deep_scheme = 'ZM' do_clubb_sgs = .false. do_tms = .true. eddy_scheme = 'HB' history_aero_optics = .false. history_aerosol = .false. history_amwg = .true. history_budget = .false. history_eddy = .false. macrop_scheme = 'rk' microp_scheme = 'RK' radiation_scheme = 'camrt' shallow_scheme = 'Hack' srf_flux_avg = 0 use_subcol_microp = .false. waccmx_opt = 'off' / &physconst_nl tms_orocnst = 1.0D0 tms_z0fac = 0.100D0 / &prescribed_aero_nl prescribed_aero_cycle_yr = 2000 prescribed_aero_datapath = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/chem/trop_mozart_aero/aero' prescribed_aero_file = 'aero_1.9x2.5_L26_1850-2005_c091112.nc' prescribed_aero_specifier = 'sulf:SO4', 'bcar1:CB1', 'bcar2:CB2', 'ocar1:OC1', 'ocar2:OC2', 'sslt1:SSLT01', 'sslt2:SSLT02', 'sslt3:SSLT03', 'sslt4:SSLT04', 'dust1:DST01', 'dust2:DST02', 'dust3:DST03', 'dust4:DST04' prescribed_aero_type = 'CYCLICAL' / &prescribed_ozone_nl prescribed_ozone_cycle_yr = 0 prescribed_ozone_datapath = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/ozone' prescribed_ozone_file = 'waccm_ozone_c121126.nc' prescribed_ozone_name = 'O3' prescribed_ozone_type = 'CYCLICAL' / &rad_cnst_nl rad_climate = 'A:Q:H2O', 'N:O2:O2', 'N:CO2:CO2', 'N:ozone:O3', 'A:N2O:N2O', 'A:CH4:CH4', 'A:CFC11:CFC11', 'A:CFC12:CFC12', 'N:sulf:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/sulfate_camrt_c080918.nc', 'N:dust1:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/dust1_camrt_c080918.nc', 'N:dust2:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/dust2_camrt_c080918.nc', 'N:dust3:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/dust3_camrt_c080918.nc', 'N:dust4:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/dust4_camrt_c080918.nc', 'N:bcar1:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/bcpho_camrt_c080918.nc', 'N:bcar2:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/bcphi_camrt_c080918.nc', 'N:ocar1:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/ocpho_camrt_c080918.nc', 'N:ocar2:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/ocphi_camrt_c080918.nc', 'N:SSLTA:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/ssam_camrt_c080918.nc', 'N:SSLTC:/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/physprops/sscm_camrt_c080918.nc' / &ref_pres_nl trop_cloud_top_press = 1.D2 / &solar_inparm solar_data_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/solar/spectral_irradiance_Lean_1610-2009_ann_c100405.nc' solar_data_type = 'FIXED' solar_data_ymd = 20000101 solar_htng_spctrl_scl = .true. / &spmd_fv_inparm npr_yz = 32,3,3,32 / &tropopause_nl tropopause_climo_file = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/cam/chem/trop_mozart/ub/clim_p_trop.nc' / &vert_diff_nl diff_cnsrv_mass_check = .false. do_iss = .false. eddy_lbulk_max = 40.D3 eddy_leng_max = 40.D3 eddy_max_bot_pressure = 100.D3 kv_freetrop_scale = 1.D0 kv_top_pressure = 0.D0 kv_top_scale = 1.D0 / &waccm_forcing_nl waccm_forcing_datapath = '/home/zhangzhy2019/soft/cesm1_2_22/inputdata/atm/waccm/ub' waccm_forcing_file = 'ghg_forcing_2000_c110321.nc' / &wetdep_inparm gas_wetdep_list = '' / &zmconv_nl zmconv_c0_lnd = 0.0035D0 zmconv_c0_ocn = 0.0035D0 zmconv_ke = 1.0E-6 /