!-------------------------------------------------------------------------------
!-- INITIALIZATION PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar
!-------------------------------------------------------------------------------
&initialization_parameters
!
!-- grid parameters
!-------------------------------------------------------------------------------
    nx                         = 23, ! Number of gridboxes in x-direction (nx+1)
    ny                         = 23, ! Number of gridboxes in y-direction (ny+1)
    nz                         = 40, ! Number of gridboxes in z-direction (nz)

    dx                         = 50.0, ! Size of single gridbox in x-direction
    dy                         = 50.0, ! Size of single gridbox in y-direction
    dz                         = 50.0, ! Size of single gridbox in z-direction

    dz_stretch_level           = 1000.0, ! Height (in m) where stretching starts
    dz_stretch_factor          = 1.08,   ! dz(k+1) = dz(k) * dz_stretch_factor

!
!-- initialization and vertical profiles
!-------------------------------------------------------------------------------
    initializing_actions       = 'set_constant_profiles initialize_bubble',
                                         ! initial conditions (with warm air bubble)

    ug_surface                 = 0.0,    ! u-comp of geostrophic wind at surface
    surface_pressure           = 1015.4, ! surface pressure

    pt_surface                 = 297.9,            ! temperature at surface
    pt_vertical_gradient       = 0.0, 0.58588957,  ! vertical gradient of temperature
    pt_vertical_gradient_level = 0.0, 740.0,       ! height level of temp gradients

    q_surface                  = 0.016,                                    ! mixing ratio at surface
    q_vertical_gradient        = -2.97297E-4, -4.5238095E-4, -8.108108E-5, ! gradient for mix. ratio
    q_vertical_gradient_level  = 0.0, 740.0, 3260.0,                       ! height lev. for gradients

!
!--  passive scalar
!-------------------------------------------------------------------------------
     passive_scalar            = .TRUE.,      ! enable passive scalar
     s_surface                 = 0.0008,      ! equivalent to 420ppm background
     surface_scalarflux        = 0.002,       ! kg emission per square meter and second
     bc_s_b                    = 'neumann',   ! Neumann boundary conditions
     scalar_advec              = 'bc-scheme', ! Bott chlond scheme for scalar advection
     loop_optimization         = 'vector',    ! loops over k,j,i

!
!--  bulk microphysics
!-------------------------------------------------------------------------------
     humidity                  = .TRUE.,     ! prog. equation for total water mixing ratio

!
!-- numerics
!-------------------------------------------------------------------------------
    fft_method                 = 'temperton-algorithm',   ! build-in fft method

/ ! end of initialization parameter namelist


!-------------------------------------------------------------------------------
!-- RUNTIME PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
!-------------------------------------------------------------------------------
&runtime_parameters
!
!-- run steering
!-------------------------------------------------------------------------------
    end_time                   = 1200.0, ! simulation time of the 3D model

    create_disturbances        = .TRUE., ! randomly perturbate horiz. velocity
    dt_disturb                 = 150.0,  ! interval for random perturbations
    disturbance_energy_limit   = 0.01,   ! upper limit for perturbation energy

    data_output_2d_on_each_pe  = .FALSE., ! don't do 2D output on each MPI rank

!
!-- data output
!-------------------------------------------------------------------------------
    netcdf_data_format         = 2, ! use NetCDF3
    dt_run_control             = 0.0,    ! output interval for run control
    dt_data_output             = 400.0,  ! output interval for general data
    dt_dopr                    = 400.0,  ! output interval for profile data

    data_output                = 'w_xy',
                                 'w_xz',
                                 'pt_xy',
                                 'pt_xz',
                                 'ql_yz',
                                 'w_yz',
                                 'u_yz',
                                 'v_yz',
                                 's_yz',
                                 'nc_yz',

    data_output_pr             = 'u',
                                 'v',
                                 'w',
                                 'ql',
                                 '#pt',
                                 'w"pt"',
                                 'w*pt*',
                                 'wpt',
                                 'w*2',
                                 'pt*2',

    section_xy                 = 2, 10,  ! grid index for 2D XY cross sections
    section_xz                 = 20,     ! grid index for 2D XZ cross sections
    section_yz                 = -1, 10, ! grid index for 2D YZ cross sections

    averaging_interval_pr      = 100.0, ! averaging interval profile data
    dt_averaging_input_pr      = 10.0,  ! averaging profile data sampling rate

/ ! end of runtime parameter namelist


!-------------------------------------------------------------------------------
!-- BULK CLOUD PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/bcmpar
!-------------------------------------------------------------------------------
&bulk_cloud_parameters

     bulk_cloud_model          = .TRUE.,     ! enable bulk cloud model
     cloud_scheme              = 'morrison', ! two-moment liquid water microphysics
                                             ! including activation, diffusion, autoconversion,
                                             ! accretion, sedimentation, precipitation

     aerosol_bulk              = 'nacl',     ! chemical composition of aerosol

     na_init                   = 1.5E8,   ! background aerosol concetration
     dry_aerosol_radius        = 0.15E-6, ! mean dry aerosol radius
     sigma_bulk                = 1.61,    ! width of dry aerosol spectrum

     curvature_solution_effects_bulk   = .TRUE., ! considering aerosol properties for activation
     collision_turbulence              = .TRUE., ! parametize turbulence effects for collision
     cloud_water_sedimentation         = .TRUE., ! enable cloud water sedimentation

     call_microphysics_at_all_substeps = .TRUE., ! microphysics every timestep

/ ! end of bulk cloud parameter namelist