source: palm/trunk/EXAMPLES/test_salsa/test_salsa_p3d @ 3899

!-------------------------------------------------------------------------------
!-- 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