!-------------------------------------------------------------------------------
!-- INITIALIZATION PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar
!-------------------------------------------------------------------------------
 &initialization_parameters

!
!-- Grid
!------------------------------------------------------------------------------
    nx                   = 68,
    ny                   = 68,
    nz                   = 30,

    dx                   = 2.0,
    dy                   = 2.0,
    dz                   = 2.0,

    dz_stretch_level     = 300.0,
    dz_stretch_factor    = 1.08,
    dz_max               = 250.0,

!
!-- Numerics
!------------------------------------------------------------------------------
    fft_method           = 'fftw',
    psolver              = 'multigrid'
    cycle_mg             = 'w',
    mg_cycles            = 2,
    mg_switch_to_pe0_level = -1,
    bc_p_t               = 'neumann',

!
!-- Mode
!------------------------------------------------------------------------------
    humidity                   = .T.,
    q_surface                  = 0.010,
    q_vertical_gradient        = 0.0, -1.0e-3,             ! gradient for mix. ratio
    q_vertical_gradient_level  = 0.0, 30, 825              ! height lev. for gradients

!
!-- Initialization
!------------------------------------------------------------------------------
    initializing_actions = 'set_constant_profiles',   !'set_constant_profiles', or 'inifor',
                                                      ! or both (met. from inifor, chem constant profiles)

    pt_surface                 = 292.5,
    pt_vertical_gradient       =   0.0,  12.0,   0.5,
    pt_vertical_gradient_level =   0.0, 800.0, 825.0,

    ug_surface           = 1.0,
    vg_surface           = 0.0,

    origin_date_time     = '2019-06-01 09:00:00 +00'

!
!-- Topography
!------------------------------------------------------------------------------
    topography           = 'read_from_file',

!
!-- Physics
!------------------------------------------------------------------------------
    longitude            = 13.4,
    latitude             = 52.5,

!
!-- Wall/soil spinup
!------------------------------------------------------------------------------
    spinup_time         = .000001, !86400.0, !120.0,
    spinup_pt_mean      = 292.5
    spinup_pt_amplitude = 0.0,             ! give a dirunal cycle
    dt_spinup           = .000001, !60.0,
    data_output_during_spinup = .T.,

 / ! 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              = 0.0,
    create_disturbances   = .F.,
    npey                  = 1,
    npex                  = 1,
!
!-- Run-control/timeseries output settings
!------------------------------------------------------------------------------
    dt_run_control        =  0.0,
    dt_dots               =  60.0,

!
!-- Profile output settings
!------------------------------------------------------------------------------
    skip_time_dopr        =  0.0,
    dt_dopr               = 60.0,
    averaging_interval_pr = 60.0,
    dt_averaging_input_pr =  0.0,

    data_output_pr        = 'wtheta', 'w"theta"', 'w*theta*',
                            'w*2', '#theta', 'theta*2',
                            'u',
                            'v',
                            'q',

                            !'rad_lw_hr', 'rad_sw_hr',

!
!-- 2D/3D output settings
!------------------------------------------------------------------------------
    do3d_at_begin      = .T.,
    do2d_at_begin      = .T.,

    dt_data_output     = 180.0,
    dt_data_output_av  = 180.0,
    averaging_interval = 180.0,
    dt_averaging_input =  0.0,

    section_xy         = 0,2,10,
    section_xz         = 6,

    data_output = 'rad_lw_in', 'rad_lw_out',
                  'rad_sw_in', 'rad_sw_out',
                  'rad_sw_in*_xy', 'ta_2m*_xy', 'theta_2m*_xy', 'wspeed_10m*_xy',
                  'hr', 'rh','wspeed',
                  'q', 'theta',
                  'u', 'v', 'w',
                  'tsurf*_xy',
                  !'rad_lw_hr','rad_sw_hr',

!
!-- Debug
!---------------------------------------------------------------------------------
!    debug_output = .T.,
!    debug_output_timestep = .T.,

 / ! end of d3par namelist

!&bulk_cloud_parameters
!!
!!-- set microphyiscs scheme
!!-------------------------------------------------------------------------------
!    cloud_scheme              = 'seifert_beheng', ! two-moment liquid water microphysics
!                                                  ! including autoconversion,
!                                                  ! accretion, sedimentation, precipitation
!!
!!-- steering of liquid phase
!!-------------------------------------------------------------------------------
!    collision_turbulence              = .TRUE.,   ! parametize turbulence effects for collision
!    cloud_water_sedimentation         = .TRUE.,   ! enable cloud water sedimentation
!    nc_const                          = 1200.0E6,  ! cloud droplet concentration
!    call_microphysics_at_all_substeps = .FALSE.,  ! microphysics every sub-timestep disabled
!
!!
!!-- steering of ice phase
!!-------------------------------------------------------------------------------
!    !microphysics_ice_phase    = .TRUE.,           ! turn on ice microphysics
!    !in_init                   = 4000.0,           ! set ice nucleii concentration
!    !ice_crystal_sedimentation = .TRUE.,           ! turn on ice sedimentation
!    !snow                      = .TRUE.,           ! turn on prog. quantities for snow
!    !graupel                   = .TRUE.,           ! turn on prog. quantities for graupel
!    !start_ice_microphysics    = 0.0,              ! start with ice microphysics at model start
!
!/ ! end of bulk cloud parameter namelist

!-------------------------------------------------------------------------------
!-- RADIATION MODEL PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radpar
!-------------------------------------------------------------------------------
  &radiation_parameters

     radiation_scheme = 'tenstream',
     !radiation_scheme = 'rrtmg',
     albedo_type      = 5,
     constant_albedo  = .T.,
     use_broadband_albedo = .T.,
     dt_radiation     = 180.0,
     !l_ts_debug       = .T.,
     !radiation_interactions_on = .F.
  /


!-------------------------------------------------------------------------------
!-- LAND SURFACE MODEL PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/lsmpar
!-------------------------------------------------------------------------------
 &land_surface_parameters

!    surface_type = 'vegetation',

    constant_roughness     = .T.,
    aero_resist_kray       = .T.,

    vegetation_type        = 2,
    soil_type              = 3,
    conserve_water_content = .T.,

    root_fraction          = 0.49, 0.26, 0.25, 0.0, 0.0, 0.0, 0.0, 0.0,
    soil_temperature       = 277.15, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0,
    soil_moisture          = 0.30,   0.30,  0.30,  0.30,  0.30, 0.30,  0.30,  0.30,
    deep_soil_temperature  = 277.0,
 / ! end of lsm_par namelist


!-------------------------------------------------------------------------------
!-- URBAN SURFACE MODEL PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/usmpar
!-------------------------------------------------------------------------------
 &urban_surface_parameters

    usm_wall_mod = .T.,

 / ! end of urban_surface_par namelist

!-------------------------------------------------------------------------------
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/usmpar
!-------------------------------------------------------------------------------
 &surface_data_output_parameters

    data_output_surf = 'albedo',        'emissivity',
                       'rad_net',       'rad_net_av',
                       'shf',           'shf_av',
                       'qsws',          'qsws_av',
                       'q_surface',     'q_surface_av',
                       'theta_surface', 'theta_surface_av',
                       'ghf',           'ghf_av',
                       'rad_lw_dif',    'rad_lw_dif_av',
                       'rad_lw_in',     'rad_lw_in_av',
                       'rad_lw_out',    'rad_lw_out_av',
                       'rad_sw_dir',    'rad_sw_dir_av',
                       'rad_sw_dif',    'rad_sw_dif_av',
                       'rad_sw_in',     'rad_sw_in_av',
                       'rad_sw_out',    'rad_sw_out_av',
    dt_dosurf = .000001,
    averaging_interval_surf = 180,
    dt_dosurf_av = 180,
    to_netcdf = .TRUE., to_vtk = .False.,
    skip_time_dosurf = 0.0,
/
