!-------------------------------------------------------------------------------
!-- INITIALIZATION PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar
!-------------------------------------------------------------------------------
&initialization_parameters
!
!-- grid parameters
!-------------------------------------------------------------------------------
    nx                         = 39, ! Number of gridboxes in x-direction (nx+1)
    ny                         = 39, ! 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.585,       ! 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.97E-4, -4.52E-4, -8.1E-5,   ! gradient for mix. ratio
    q_vertical_gradient_level  = 0.0, 740.0, 3260.0,            ! height lev. for gradients

    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',
                                 'theta_xy',
                                 'theta_xz',
                                 'ql_yz',
                                 'qr_yz',
                                 'nr_yz',
                                 'w_yz',
                                 'u_yz',
                                 'v_yz',

    data_output_pr             = 'u',
                                 'v',
                                 'ql',
                                 '#theta',
                                 'w"theta"',
                                 'w*theta*',
                                 'wtheta',
                                 'w*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

    cloud_scheme              = 'seifert_beheng', ! two-moment liquid water microphysics
                                                  ! including autoconversion,
                                                  ! accretion, sedimentation, precipitation
    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
