!------------------------------------------------------------------------------- !-- INITIALIZATION PARAMETER NAMELIST ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar !------------------------------------------------------------------------------- &initialization_parameters ! !-- Grid !------------------------------------------------------------------------------ nx = 19, ny = 19, nz = 60, dx = 2.0, dy = 2.0, dz = 2.0, ! !-- Numerics !------------------------------------------------------------------------------ psolver = 'multigrid', momentum_advec = 'ws-scheme', ! default advection scheme scalar_advec = 'ws-scheme', conserve_volume_flow = .T., ! !-- Mode !------------------------------------------------------------------------------ humidity = .T., q_surface = 0.008, ! !-- Initialization !------------------------------------------------------------------------------ initializing_actions = 'set_constant_profiles', pt_surface = 293.15, ug_surface = 1.0, vg_surface = 0.0, roughness_length = 0.05, neutral = .T., day_of_year_init = 65, ! March 06, to avoid extreme temperatures time_utc_init = 71700.0, !(19:55, UTC in seconds) ! !-- Boundary conditions !------------------------------------------------------------------------------ bc_lr = 'cyclic', bc_ns = 'cyclic', bc_uv_t = 'dirichlet', ! channel flow boundary condition ! !-- Topography !------------------------------------------------------------------------------ topography = 'read_from_file', ! !-- Physics !------------------------------------------------------------------------------ longitude = 13.4, ! Berlin longitude latitude = 52.5, ! Berlin latitude / ! end of inipar namelist !------------------------------------------------------------------------------- !-- RUNTIME PARAMETER NAMELIST ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par !------------------------------------------------------------------------------- &runtime_parameters ! !-- Run steering !------------------------------------------------------------------------------ end_time = 1000.0, create_disturbances = .T., ! !-- General output settings !------------------------------------------------------------------------------ netcdf_data_format = 4, ! use NetCDF4 ! !-- Run-control/timeseries output settings !------------------------------------------------------------------------------ dt_run_control = 0.0, dt_dots = 0.0, ! !-- Profile output settings !------------------------------------------------------------------------------ skip_time_dopr = 0.0, dt_dopr = 500.0, averaging_interval_pr = 500.0, dt_averaging_input_pr = 0.0, data_output_pr = '#u', 'u*2', 'wu', 'w*u*', 'w"u"', '#v', 'v*2', 'wv', 'w*v*', 'w"v"', 'w', 'w*2', ! !-- 2D/3D output settings !------------------------------------------------------------------------------ do3d_at_begin = .T., do2d_at_begin = .T., dt_data_output = 50.0, data_output = 'u', 'v', 'w', 'Ntot', 'LDSA', 'PM2.5', 's_SO4', 's_OC', 's_BC', 's_NH', 's_NO', 'g_H2SO4','g_HNO3','g_NH3','g_OCNV','g_OCSV', 'N_bin3', 'm_bin4', / ! end of d3par namelist !------------------------------------------------------------------------------- !-- PLANT CANOPY MODEL PARAMETER NAMELIST ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/canpar !------------------------------------------------------------------------------- &plant_canopy_parameters canopy_mode = 'read_from_file_3d', canopy_drag_coeff = 0.5, / ! end of canopy_par namelist !------------------------------------------------------------------------------- !-- SALSA PARAMETER NAMELIST !------------------------------------------------------------------------------- &salsa_parameters ! !-- Time stepping dt_salsa = 4.0, ! time step for calculating aerosol processes (s) skip_time_do_salsa = 100.0, ! starting time of SALSA (s) ! !-- Concentration initial types: 0 = based on the parameter file, 1 = read from a NetCDF file init_aerosol_type = 1, ! size distribution init_gases_type = 1, ! gases ! !-- If init_aerosol_type = 0, define the initial aerosol size distribution by dpg, sigmag and n_lognorm ! dpg = 13.5E-9, 54.0E-9, 864.1E-9, ! mean diameter per mode (in metres) ! sigmag = 1.8, 2.16, 2.21, ! standard deviation per mode ! n_lognorm = 1.43E9, 4.45E8, 7.0E4, ! number concentration per mode (#/m3) ! !-- If init_gases_type = 0, apply these initial gas concentrations ! H2SO4_init = 5.0E12, ! sulphuric acid (#/m3) ! HNO3_init = 3.0E15, ! nitric acid (#/m3) ! NH3_init = 6.0E15, ! ammonia (#/m3) ! OCNV_init = 1.0E14, ! non-volatile organic gases (#/m3) ! OCSV_init = 1.0E14, ! non-volatile organic gases (#/m3) ! !-- List of activated chemical components: !-- NOTE! Chemical species have to be activated here even if they are not initially present! listspec = 'SO4','OC','BC','NH','NO','','', ! List of actived aerosols ! listspec = 'SO4','OC','BC','DU','SS','NO','NH', ! !-- If isdtyp = 0, set the chemical composition of the initial particle size distribution ! mass_fracs_a = 0.1, 0.4, 0.3, 0.1, 0.1, 0.0, 0.0, ! mass fractions of chemical components: soluble ! mass_fracs_b = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, ! mass fractions of chemical components: insoluble !-- NOTE! Set everyhing to zero in mass_fracs_b if you do not want include insoluble species (default) ! !-- Sectional presentation of the particle size distribution reglim = 3.0E-9, 1.0E-8, 2.5E-6, ! limits of the subranges (m) nbin = 1, 7, ! number of bins per subrange !-- NOTE! Subrange 1 consists only of H2SO4 and/or OC nf2a = 1.0, ! Number fraction allocated to subrange 2a (b-bins will get 1-nf2a) ! !-- Aerosol emissions: salsa_emission_mode = 'read_from_file', ! 'no_emission','uniform' or 'read_from_file' !-- If salsa_source_mode = 'uniform', set emissions based on the variables below: ! aerosol_flux_dpg = 8.0E-9, 60.0E-9, 1000.0E-9, ! mean diameter per mode (metres) ! aerosol_flux_sigmag = 1.8, 2.16, 2.21, ! standard deviation per mode ! surface_aerosol_flux = 1.0E8, 1.0E9, 1.0E5, ! number concentration per mode (#/m3) ! aerosol_flux_mass_fracs_a = 0.3, 0.0, 0.6, 0.0, 0.0, 0.0, 0.0, ! order: listspec ! aerosol_flux_mass_fracs_b = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, ! order: listspec !-- NOTE! chemical components of the source have to be activated in 'listspec' !-- and have to be in the same order ! !-- Boundary conditions !-- Decycle at the left & right boundaries: initial concentration !-- is copied to ghost and first three layers at both boundaries at the left boundary !-- and a zero gradient is set at the right boundary !-- Decycle at the north & south boundaries: initial concentration !-- is copied to ghost and first three layers at both boundaries decycle_lr = .T., decycle_ns = .F., decycle_method = 'dirichlet','dirichlet','dirichlet','dirichlet', bc_salsa_b = 'neumann', !'dirichlet', ! surface flux requires 'neumann' bc_salsa_t = 'dirichlet', ! top ! !-- Switches for aerosol processes: nldistupdate = .T., ! update aerosol size distribution nldepo = .T., ! Deposition master switch nldepo_pcm = .T., ! Deposition on vegetation nldepo_surf = .T., ! Deposition on walls nlcnd = .F., ! Condensation master switch nlcndgas = .F., ! Condensation of precursor gases nlcndh2oae = .F., ! Condensation of H2O on aerosols nlcoag = .F., ! Coagulation master switch nsnucl = 0, ! Nucleation scheme (0 = off) nj3 = 1, ! J3 parametrization for nucleation ! !-- Deposition on vegetation (pcm = plant canopy model) depo_pcm_par = 'zhang2001', ! or 'petroff2010' ! !-- Deposition on on ground, walls and roofs depo_surf_par = 'zhang2001', ! or 'petroff2010' ! !-- Other switches: advect_particle_water = .T., ! particle water: advect or calculate at each dt_salsa feedback_to_palm = .F., ! feedback to flow due to condensation of water nest_salsa = .F., ! apply nesting for salsa read_restart_data_salsa = .F., ! skip reading restart data even if it's done for the flow write_binary_salsa = .F., ! skip writing restart data even if it's done for the flow / ! end of salsa_par namelist