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

    dx                         = 243.0, ! Size of single gridbox in x-direction
    dy                         = 243.0, ! Size of single gridbox in y-direction
    dz                         = 54.0, ! Size of single gridbox in z-direction
!
!-- initialization
!-------------------------------------------------------------------------------
    turbulence_closure         = '1.5-order',   
    initializing_actions       = 'set_constant_profiles', ! initial conditions


    omega = 0.0,  ! no coriolis force

    uv_heights = 0.0, 10.0, 38.0, 119.0, 338.0, 638.0, 760.0, 
                 938.0, 1238.0, 1449.0, 1538.0, 1838.0, 2438.0, 
                 2989.0, 3038.0, 4238.0, 5521.0, 6338.0, 7098.0, 7200.0,

    u_profile = 0.0, 0.98480773, 1.44888878, -0.76604444, -3.5999999, 
                4.58249569, 3.38289332, 2.56050014, 2.54558468, 2.37473798, 
                2.37473798, 0.17364843, 1.90324628, 3.55070424, 3.55070424, 
                9.76270866, 14.84330082, 14.00142002, 16.74173164, 16.74173164,

    v_profile = 0.0, -0.17364816, -0.38822848, -0.64278758, 1.6e-07, 0.40091625, 
                1.23127222, 0.45148513, -2.54558396, -1.99264121, -1.99264121, 
                -0.98480773, 0.88749814, 2.04999948, 2.04999948, 0.85412598, 
                1.29861999, 5.0960989, 2.95201826, 2.95201826,

    pt_surface = 273.649994,

    pt_vertical_gradient = 8.64461263, 4.92338958, 0.95464317, 
                           0.3187736, 0.3327503, 0.09285783,

    pt_vertical_gradient_level = 2.0, 38.0, 119.0, 760.0, 1449.0, 5492.0,

    q_surface = 0.004793,

    q_vertical_gradient = -0.00066591, -0.00017517, 5.641e-05, -0.0001911, 
                          -5.109e-05, -9.316e-05, -6.651e-05, -3.444e-05, -2.619e-05,

    q_vertical_gradient_level = 2.0, 60.56463242, 115.83196259, 390.80038452, 
                                925.44238281, 1845.13232422, 3072.58642578, 
                                4287.02050781, 5461.89697266,

    passive_scalar = .TRUE.,
!-- humidity
!-------------------------------------------------------------------------------
    humidity                   = .TRUE.,
    cloud_droplets             = .FALSE.,
    check_realistic_q          = .FALSE.,

!
!-- geophysical setup
!-------------------------------------------------------------------------------
    latitude                   = -43.5,
    longitude                  = 172.56,
    origin_date_time           = '2001-08-05 12:00:00 +00',

!
!-- boundary conditions
!-------------------------------------------------------------------------------
    bc_lr = 'nested',
    bc_ns = 'nested',
    bc_p_b = 'neumann',
    bc_p_t = 'nested',

    constant_flux_layer        = .TRUE.,
    bc_pt_b                    = 'dirichlet', ! required with surface_heatflux
    bc_pt_t                    = 'nested',
    bc_q_b                     = 'dirichlet',
    bc_q_t                     = 'nested',
    bc_s_b                     = 'dirichlet',
    bc_s_t                     = 'nested',
    bc_uv_b                    = 'dirichlet',
    bc_uv_t                    = 'nested',
    bc_e_b                     = 'neumann',

   ! rayleigh_damping_height    = 6000.0,
   ! rayleigh_damping_factor    = 0.01,

!
!-- numerics
!-------------------------------------------------------------------------------
    fft_method                 = 'fftw',   ! build-in fft method
    loop_optimization          = 'vector',
    approximation = 'boussinesq',
    random_heatflux = .F.,
    scalar_advec   = 'ws-scheme',
    momentum_advec = 'ws-scheme',
    roughness_length    =  0.1,
    z0h_factor           =  1.0,
    psolver = 'multigrid',
    reference_state            = 'horizontal_average', ! based on Cabauw
    
    topography = 'read_from_file',
/ ! 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                   = 172800.0, ! simulation time of the 3D model
    npex = 12,
    npey = 8,
    termination_time_needed    = 900.0,

    create_disturbances        = .TRUE., ! randomly perturbate horiz. velocity
    dt_disturb                 = 150.0,  ! interval for random perturbations
    disturbance_energy_limit   = 0.0001, ! 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         = 4, ! use NetCDF4 on Maui
    dt_max                     = 5.0,
    dt_run_control             = 60.0,    ! output interval for run control
    dt_data_output             = 60.0,  ! output interval for general data
    dt_data_output_av          = 60.0, ! output interval for averaged data
    dt_dopr                    = 60.0,  ! output interval for profile data
    dt_dots                    = 60.0,

    data_output    = 'e', 'p', 'q', 'rh', 'theta', 'u', 'v', 'w', 'wdir', 'wspeed', 'm_soil', 
                     't_soil', 'w_av', 'q_av',
                     'qsurf*_xy', 'qsws*_xy', 'shf*_xy', 'ssurf*_xy', 'theta_2m*_xy',
                     'tsurf*_xy', 'us*_xy', 'wspeed_10m*_xy', 'z0*_xy', 'z0h*_xy', 'ghf*_xy', 
                     'qsws_soil*_xy', 'qsws_liq*_xy', 'qsws_veg*_xy', 'rad_net*_xy',
                     'rad_lw_in*_xy','rad_lw_out*_xy', 'rad_sw_in*_xy', 'rad_sw_out*_xy',
                     'theta_xy', 'p_xy', 'u_xy', 'v_xy', 'w_xy', 'wdir_xy', 'wspeed_xy',
                     
      

    data_output_pr = 'e', 'e*', 'kh', 'km', '#q', 'q*2', 'rh', 'rho', '#theta', 'theta*2', 
                     '#u', 'u*2', 'ug', '#v', 'v*2', 'vg', 'w', 'w*2', 'w*3', 'w*e*', 'wq', 'w"q"',
                     'w*q*', 'wqv', 'w"qv"', 'w*qv*', 'wtheta', 'w"theta"', 'w*theta*', 'wu', 'w"u"',
                     'w*u*', 'wv', 'w"v"', 'w*v*', '#qv', 'p', '#t_soil', '#m_soil',

    section_xy                 = 0, 1, 2, 3, 4, 5, ! grid index for 2D XY cross sections
    dt_do2d_xy                 = 60.0,
    skip_time_do2d_xy          = 0.0,

  !  section_xz                 = -1,  ! grid index for 2D XZ cross sections
  !  dt_do2d_xz                 = 600.0,
  !  skip_time_do2d_xz          = 0.0,

    do2d_at_begin              = .FALSE.,

    do3d_at_begin              = .FALSE.,
    skip_time_do3d             = 0.0,

    averaging_interval         = 10.0, ! averaging interval general data
    dt_averaging_input         = 0.0,  ! averaging general data sampling rate

    averaging_interval_pr      = 0.0, ! averaging interval profile data
    dt_averaging_input_pr      = 0.0,  ! averaging profile data sampling rate

/ ! end of runtime parameter namelist

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

    vegetation_type            = 3,
    soil_type                  = 3,
    surface_type               = 'netcdf',
    constant_roughness         = .TRUE.,

!
!-- lsm initialization
!-------------------------------------------------------------------------------
   ! vegetation_coverage        = 1.0,

!    dz_soil                    = 0.01, 0.02, 0.04, 0.06, 0.14, 0.26, 0.54, 1.86,
    soil_moisture              = 0.31257257, 0.31257257, 0.31257257, 0.31257036, 
                                 0.31255273, 0.31187841, 0.29367608, 0.26859525, 

    soil_temperature           = 278.12188721, 278.12188721, 278.12188721, 278.21648711, 
                                 278.97328643, 280.00322266, 279.70466312, 280.58267212,

    deep_soil_temperature      = 285.859955, 

/ ! end of land surface parameter namelist

!-------------------------------------------------------------------------------
!-- RADIATION PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radpar
!-------------------------------------------------------------------------------
&radiation_parameters
    radiation_scheme           = 'rrtmg',
!
!-- rrtmg initialization
!-------------------------------------------------------------------------------
  !  skip_time_do_radiation     = 2100.0,
  !  albedo_type                = 5,
    dt_radiation               = 60.0,
    nrefsteps                  = 3,
    surface_reflections        = .FALSE.,
    radiation_interactions_on  = .FALSE.,

/ ! end of radiation parameter namelist

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

!
!-- switch on bulk cloud model and choose microphyics scheme
!-------------------------------------------------------------------------------
!    bulk_cloud_model          = .TRUE.,
!    cloud_scheme              = 'kessler',
!!    microphysics_ice_phase     = .TRUE.,

!
!-- cloud number concentratiion and aersol concentratiion
!-------------------------------------------------------------------------------
!    nc_const                   = 2.0E8,

!
!-- setting aerosol conditions
!-------------------------------------------------------------------------------
!    curvature_solution_effects_bulk = .TRUE.,
!    aerosol_bulk                    = 'nacl',
!    na_init                         = 8.42E8,
!    dry_aerosol_radius              = 0.021E-6,
!    sigma_bulk                      = 3.19,

!
!-- enabling turbulence effects for collision process, cloud water sedimentation
!-------------------------------------------------------------------------------
!    collision_turbulence       = .TRUE.,     ! parametize turbulence effects for collision
!    cloud_water_sedimentation  = .TRUE.,     ! enable cloud water sedimentation
!!    precipitation              = .TRUE.,

!    call_microphysics_at_all_substeps = .FALSE.,   ! microphysics every timestep


!/ ! end of bulk cloud parameter namelist

                                    
