!-------------------------------------------------------------------------------
!-- INITIALIZATION PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar
!-------------------------------------------------------------------------------
 &initialization_parameters
!-- Grid
!------------------------------------------------------------------------------
          nx = 39,                                  ! (nx+1)*dx for cyclic bc / (nx+1) for non-cyclic
          ny = 39,                                  ! (ny+1)*dy for cyclic bc / (ny+1) for non-cyclic
          nz = 40,                                  ! (nz+2)
          dx = 2.0, dy = 2.0, dz = 2.0,
          
          topography     = 'single_street_canyon',
          canyon_height  = 40.0,
          canyon_width_x = 40.0,
         
!-- Numerics
!------------------------------------------------------------------------------      
          fft_method     = 'temperton-algorithm',   ! The fast Fourier transformation is used for solving the perturbation pressure equation with a direct method and for calculating power spectra
          momentum_advec =           'pw-scheme',        
          scalar_advec   =           'pw-scheme',          
         
!
!-- Initialization
!------------------------------------------------------------------------------         
          initializing_actions = 'set_constant_profiles',

          ug_surface       = 1.0, 
          vg_surface       = 0.0,
          roughness_length = 0.1,
         
!
!-- Physics
!------------------------------------------------------------------------------
          omega = 0.0,
         
!
!-- Boundary conditions
!------------------------------------------------------------------------------
          bc_pt_b              = 'neumann',
          bc_uv_t              = 'neumann', 
          
!
!-- Mode
!------------------------------------------------------------------------------
          neutral = .T.,                    ! strictly neutral flow
          conserve_volume_flow =       .T., ! maintain constant bulk velocity
         
/ ! end of initialization namelist

!-------------------------------------------------------------------------------
!-- RUNTIME PARAMETER NAMELIST
!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
!-------------------------------------------------------------------------------
 &runtime_parameters
 
!
!-- Run steering
!------------------------------------------------------------------------------ 
          end_time                 = 1200.0,  
          create_disturbances      =    .T.,          ! Switch to impose random perturbations to the horizontal velocity field

! 
!-- General output settings
!------------------------------------------------------------------------------         
          netcdf_data_format = 2,

!
!-- Run-control/timeseries output settings
!------------------------------------------------------------------------------         
          dt_run_control = 6.0,                     ! Temporal interval at which run control output is to be made
          dt_dots = 6.0,                            ! Temporal interval at which time series data shall be output

!
!-- Profile output settings
!------------------------------------------------------------------------------
          skip_time_dopr            =      0.0,     ! No output of vertical profile data before this interval has passed 
          dt_dopr                   =   1200.0,     ! Temporal interval at which data of vertical profiles shall be output
          averaging_interval_pr     =   1200.0,
          dt_averaging_input_pr     =      6.0,     ! Temporal interval of data which are subject to temporal averaging
          data_output_pr            = '#u', 'w"u"', 'w*u*', 'wu', 'u*2', 
                                      '#v', 'w"v"', 'w*v*', 'wv', 'v*2',
                                      'w*2', 'e', 'e*', '#km', '#l',

!
!-- 2D/3D output settings
!------------------------------------------------------------------------------         
          do3d_at_begin          =    .F.,
         
          skip_time_data_output  =    0.0,
          dt_data_output         = 1200.0,
          averaging_interval     = 1200.0,
          dt_averaging_input     =    6.0,

          section_xz             =    1,                        !j*dy
         
          data_output = 'u',       'v',       'w',
                       'u_xz_av', 'v_xz_av', 'w_xz_av',
                       'u_av',    'v_av',    'w_av',

/ ! end of runtime namelist

