UCLALES-SALSA simulations ========================= +Simulations used in Raatikainen, T., Prank, M., Ahola, J., Kokkola, H., Tonttila, J., and Romakkaniemi, S.: The effect of marine ice nucleating particles on mixed-phase clouds, Atmos. Chem. Phys. Discuss. [preprint], https://doi.org/10.5194/acp-2021-537, in review, 2021. The model --------- +The model is available from GitHub: https://github.com/UCLALES-SALSA/UCLALES-SALSA +These simulations are made with branch "isdac_poa", last updated Wed Mar 3 14:59:05 2021 +Direct link: https://github.com/UCLALES-SALSA/UCLALES-SALSA/tree/isdac_poa +3D outputs are merged using post-processing code ../script/combine.f90 +UCLALES-SALSA and the post-processing code were compiled using Fortran +Simulations were made with 64 PUs +Template NAMELISTs and soundings are provided Outputs ------- +UCLALES-SALSA produces NetCDF files: 1) Time series, which are averaged/integrated over domain (*.ts.nc) 2) Profiles, which are averaged over horizontal dimensions (*.ps.nc) 3) Analysis outputs, which are 3D snapshots of the domain (e.g., isdac_s0b0.nc and isdac_s0b2.nc) +Parameters are described (long name and unit) in each netCDF file. +Model outputs are compressed into a single zip file isdac.zip Default simulations ------------------- +Output file names depend on sea spray INP emissions and background INP concentration: isdac_sb.nc Sea spray INP emissions (NAMELIST variable: model/ssa_inp_frac) id title variable 0 off 0.0 2 on 0.0005 Background INP concentrations (NAMELIST variable: salsa/nf2a) id title variable 0 high 0.99990 2 med. 0.99995 3 low 0.99999 4 zero 1.0 +Eight simulations in total (title is the name used in the manuscript) Sensitivity experiments ----------------------- +The same model version as for the default simulations, but aerosol freezing is hard-coded 1) Aerosol freezing +Used NAMELIST_afr +Experiment name: isdac_s2b0_afr +New ts, ps and analysis outputs: 'nucl_na','nucl_ra' +Need a new executable where the hard-coded aerosol freezing is enabled 2) Aerosol and cloud sedimentation +Used NAMELIST_sed +Experiment name: isdac_s2b0_sed +Changed NAMELIST parameters: sed_aero = T and sed_cloud = T +New ts, ps and analysis outputs: 'sedi_3a','sedi_3c','sedi_na','sedi_ra','sedi_nc','sedi_rc' 3) Decoupled boundary layer +Used NAMELIST_dec +Experiment names: isdac_s0b0_dec and isdac_s2b0_dec +Replace sound_in with sound_in_decoupled Producing manucript figures =========================== +Using the standard netCDF outputs +Possible calculations are described for each figure Figure 1: --------- Note: not from netCDF, but computed using the specified input size distribution and SSA emission parameterizations Figures 2 and 9: ----------------------- Note: time series data used as is time = simulation time (s) zt = cloud top height (m) zb = cloud base height (m) iiNib = ice number concentration (#/kg) iwp_bar = ice water path (kg/m2) lwp_bar = liquid water path (kg/m2) Figure 3: --------- Note: time series data used as is time = simulation time (s) sedi_3i = precipitation (kg/m2/s) for_3t = subsidence (kg/m2/s) srfc_3a = surface (kg/m2/s) sedi_ni = precipitation (#/m2/s) forc_nx = subsidence (#/m2/s) srfc_nx = surface (#/m2/s) Figure 4: --------- Note: time series data used; calcuted the sum of three terms and the offset is from the first point time = simulation time (s) tcDUab+tcDUcb+tcDUib = total INP mass (kg/m2) Figure 5: --------- Note: profile data used; averaged over simulation time from 7*3600-2 s to 8*3600+2 s zt = altitude (m) advf_3t = advection (kg/m3/s) sedi_3i = precipitation (kg/m3/s) diff_3t = diffusion (kg/m3/s) forc_3t = subsidence (kg/m3/s) srfc_3a = surface (kg/m3/s) tcDUab = INP mass in aerosol (kg/kg) tcDUcb = INP mass in cloud (kg/kg) tcDUib = INP mass in ice (kg/kg) Figure 6: --------- Note: time series data used as is time = simulation time (s) forc_ri = subsidence (kg/m2/s) cond_ri = sublimation/deposition (kg/m2/s) nucl_ri = nucleation (kg/m2/s) sedi_ri = precipitation (kg/m2/s) forc_ni = subsidence (#/m2/s) diag_ni = sublimation (#/m2/s) nucl_ni = nucleation (#/m2/s) sedi_ni = precipitation (#/m2/s) Figure 7: --------- Note: profile data used; averaged over simulation time from 7*3600-2 s to 8*3600+2 s zt = altitude (m) advf_ri = advection (kg/m3/s) sedi_ri = precipitation (kg/m3/s) cond_ri = sublimation/deposition (kg/m3/s) nucl_ri = nucleation (kg/m3/s) diff_ri = diffusion (kg/m3/s) forc_ri = subsidence (kg/m3/s) advf_ni = advection (#/m3/s) sedi_ni = precipitation (#/m3/s) diag_ni = sublimation (#/m3/s) nucl_ni = nucleation (#/m3/s) diff_ni = diffusion (#/m3/s) forc_ni = subsidence (#/m3/s) l = cloud water mass (kg/kg) tcH2Oit = ice mass (kg/kg) Ni_ii = ice number (#/kg) Figure 8: --------- Note: 3D data from simulation time 8*3600 s used Note: Total dust is the sum of three terms Note: Vertical integrals use altitude (zt, m) and air density, which is calculated from potential temperature (theta, K) and pressure (p, Pa) as in the model Note: In-cloud vertical velocity is averaged over grid cells where liquid water mixing ratio (l, kg/kg) exceeds 1e-5 kg/kg. C_DUab+C_DUcb+C_DUib = total INP mass mixing ratio (kg/kg) nucl_ni = ice nucleation rate (#/kg/s) w = vertical velocity (m/s)