source: palm/trunk/TUTORIALS/cases/dispersion_eulerian_and_lpm_extended/INPUT/dispersion_eulerian_and_lpm_extended_p3d @ 4180

&initialization_parameters

!-- grid setup
    nx = 39,
    ny = 39,
    nz = 40,

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

!-- init-action
    initializing_actions = 'set_constant_profiles',

!-- advection schemes
    momentum_advec = 'ws-scheme',           ! default advection scheme
    scalar_advec   = 'ws-scheme',

!-- wind speed setup
    ug_surface = 1.0,
    vg_surface = 0.0,

!-- temperature setup
    neutral = .T.,                          ! strictly neutral flow

!-- tracer setup
    passive_scalar     = .T.,
    s_surface          = 0.0008,            !equivalent to 420ppm CO2 as background concentration
    surface_scalarflux = 0.02,              !kg emission per square meter and second
    bc_s_b             = 'neumann',

!-- boundary conditions
    bc_uv_t = 'neumann',                    ! free-slip boundary condition
    roughness_length = 0.1,                 ! applies to all surfaces locally

!-- special methods
    dp_external = .T.,                      ! use horizontal pressure gradient
    dpdxy = -0.0002, 0.0,                   ! for forcing

!-- special setups
    omega = 0.0,                            ! no Coriolis force

    topography = 'read_from_file',
    topography_grid_convention = 'cell_center',   ! default
/


&runtime_parameters

!-- simulation time
    end_time                = 10800.0,
    termination_time_needed =   900.0,

!-- others
    create_disturbances = .T.,

!-- global output settings
    skip_time_data_output =    0.0,
    dt_data_output        = 3600.0,

    averaging_interval    = 3600.0,
    dt_averaging_input    =    6.0,

!-- time series / run control
    dt_run_control        =  6.0,
    dt_dots               =  6.0,

!-- profiles
    skip_time_dopr        =    0.0,
    dt_dopr               = 3600.0,
    averaging_interval_pr = 3600.0,
    dt_averaging_input_pr =    6.0,

    do3d_at_begin         = .TRUE.,
    do2d_at_begin         = .TRUE.,

    dt_do2d_xy            = 3600.0,
    dt_do2d_xz            = 3600.0,

    data_output_pr = '#u', 'u*2', 'wu', 'w*u*', 'w"u"',
                     '#v', 'v*2', 'wv', 'w*v*', 'w"v"',
                     'w', 'w*2',
                     'e', 'e*',
                     '#km',
                     '#l',
                     '#s', 'ws', 'w*s*', 'w"s"', 

!-- output variables
    data_output = 'u', 'u_av',
                  'u_xy', 'u_xz', 'u_xz_av',
                  'v', 'v_av',
                  'w', 'w_av',
                  's', 's_av',
                  'ssws*_xy',
                  'pc', 'pc_av', 'pc_xy', 'pc_xz', 'pc_xz_av',
                  

    section_xy = -1,21,22,23,
    section_xz = -1,17,18,19,

/

&particle_parameters

!-- Boundary conditions for particles
    bc_par_lr = 'absorb',
    bc_par_b  = 'reflect',
    bc_par_t  = 'reflect', 

!-- Define point source
    psb = 45.0,  
    pst = 45.0, 
    pss = 40.0,
    psn = 40.0,
    psl = 25.0, 
    psr = 25.0, 

!-- Time after new particle set is released 
    dt_prel = 1.0,

!-- Time steps after array are deallocated
    step_dealloc = 50,

!-- Particles per point (default = 1)
    particles_per_point = 100,

!-- Maximum time for particles
    particle_maximum_age = 200.0,

!-- Add random fluctuations for initial particle position
    random_start_position = .FALSE.,

!-- SGS-velocities for particles 
    use_sgs_for_particles = .TRUE., 

!-- Particle time series interval
    dt_dopts = 10.0, 

/

&user_parameters
    emission_stripe_width_y = 6.0,  !in meters
    emission_stripe_orig_y  = 12.0, 62.0,  !in meters
    s_ts_pos_x = 17.0, 61.0, 75.0,
    s_ts_pos_y = 25.0, 15.0, 41.0,
           
/