| 1 | !------------------------------------------------------------------------------- |
|---|
| 2 | !-- INITIALIZATION PARAMETER NAMELIST |
|---|
| 3 | ! https://palm.muk.uni-hannover.de/trac/wiki/doc/app/initialization_parameters |
|---|
| 4 | !------------------------------------------------------------------------------- |
|---|
| 5 | &initialization_parameters |
|---|
| 6 | ! |
|---|
| 7 | !-- grid parameters |
|---|
| 8 | !------------------------------------------------------------------------------- |
|---|
| 9 | nx = 39, ! Number of gridboxes in x-direction (nx+1) |
|---|
| 10 | ny = 39, ! Number of gridboxes in y-direction (ny+1) |
|---|
| 11 | nz = 40, ! Number of gridboxes in z-direction (nz) |
|---|
| 12 | |
|---|
| 13 | dx = 50.0, ! Size of single gridbox in x-direction |
|---|
| 14 | dy = 50.0, ! Size of single gridbox in y-direction |
|---|
| 15 | dz = 50.0, ! Size of single gridbox in z-direction |
|---|
| 16 | |
|---|
| 17 | dz_stretch_level = 1200.0, ! Height (in m) where stretching starts |
|---|
| 18 | dz_stretch_factor = 1.08, ! dz(k+1) = dz(k) * dz_stretch_factor |
|---|
| 19 | ! |
|---|
| 20 | !-- initialization |
|---|
| 21 | !------------------------------------------------------------------------------- |
|---|
| 22 | initializing_actions = 'set_constant_profiles', ! initial conditions |
|---|
| 23 | |
|---|
| 24 | ug_surface = 0.0, ! u-comp of geostrophic wind at surface |
|---|
| 25 | vg_surface = 0.0, ! v-comp of geostrophic wind at surface |
|---|
| 26 | |
|---|
| 27 | pt_surface = 300.0, ! initial surface potential temp |
|---|
| 28 | pt_vertical_gradient = 0.0, 1.0, ! piecewise temp gradients |
|---|
| 29 | pt_vertical_gradient_level = 0.0, 800.0, ! height level of temp gradients |
|---|
| 30 | |
|---|
| 31 | ! |
|---|
| 32 | !-- boundary conditions |
|---|
| 33 | !------------------------------------------------------------------------------- |
|---|
| 34 | surface_heatflux = 0.1, ! sensible heat flux at the bottom surface |
|---|
| 35 | bc_pt_b = 'neumann', ! required with surface_heatflux |
|---|
| 36 | |
|---|
| 37 | ! |
|---|
| 38 | !-- numerics |
|---|
| 39 | !------------------------------------------------------------------------------- |
|---|
| 40 | fft_method = 'temperton-algorithm', ! build-in fft method |
|---|
| 41 | |
|---|
| 42 | / ! end of initialization parameter namelist |
|---|
| 43 | |
|---|
| 44 | !------------------------------------------------------------------------------- |
|---|
| 45 | !-- RUNTIME PARAMETER NAMELIST |
|---|
| 46 | ! https://palm.muk.uni-hannover.de/trac/wiki/doc/app/runtime_parameters |
|---|
| 47 | !------------------------------------------------------------------------------- |
|---|
| 48 | &runtime_parameters |
|---|
| 49 | ! |
|---|
| 50 | !-- run steering |
|---|
| 51 | !------------------------------------------------------------------------------- |
|---|
| 52 | end_time = 10800.0, ! simulation time of the 3D model |
|---|
| 53 | |
|---|
| 54 | create_disturbances = .TRUE., ! randomly perturbate horiz. velocity |
|---|
| 55 | dt_disturb = 150.0, ! interval for random perturbations |
|---|
| 56 | disturbance_energy_limit = 0.015, ! upper limit for perturbation energy |
|---|
| 57 | |
|---|
| 58 | data_output_2d_on_each_pe = .FALSE., ! don't do 2D output on each MPI rank |
|---|
| 59 | |
|---|
| 60 | ! |
|---|
| 61 | !-- data output |
|---|
| 62 | !------------------------------------------------------------------------------- |
|---|
| 63 | netcdf_data_format = 2, ! use NetCDF3 |
|---|
| 64 | |
|---|
| 65 | dt_run_control = 0.0, ! output interval for run control |
|---|
| 66 | dt_data_output = 900.0, ! output interval for general data |
|---|
| 67 | dt_data_output_av = 1800.0, ! output interval for averaged data |
|---|
| 68 | dt_do3d = 900.0, ! output interval for 3D data |
|---|
| 69 | dt_do2d_xy = 900.0, ! output interval for 2D XY data |
|---|
| 70 | dt_do2d_xz = 900.0, ! output interval for 2D XZ data |
|---|
| 71 | dt_do2d_yz = 900.0, ! output interval for 2D YZ data |
|---|
| 72 | dt_dopr = 900.0, ! output interval for profile data |
|---|
| 73 | |
|---|
| 74 | data_output = 'u', 'u_av','v', 'v_av','w', 'w_av','theta', 'theta_av','u_xy', 'u_xz', 'u_yz', 'u_xz_av', 'v_xy', 'v_xz', 'v_yz', 'v_xz_av', 'w_xy', 'w_xz', 'w_yz', 'w_xz_av', 'theta_xy', 'theta_xz', 'theta_yz', |
|---|
| 75 | |
|---|
| 76 | data_output_pr = 'u', 'u*2', 'v', 'v*2', 'w', 'w*2', 'w*3', '#theta', 'theta*2', 'wtheta', 'w"theta"', 'w*theta*', |
|---|
| 77 | |
|---|
| 78 | section_xy = 2, |
|---|
| 79 | 10, ! grid index for 2D XY cross sections |
|---|
| 80 | section_xz = 5, |
|---|
| 81 | 20, ! grid index for 2D XZ cross sections |
|---|
| 82 | section_yz = 20, |
|---|
| 83 | 25, ! grid index for 2D YZ cross sections |
|---|
| 84 | |
|---|
| 85 | averaging_interval = 900.0, ! averaging interval general data |
|---|
| 86 | dt_averaging_input = 10.0, ! averaging general data sampling rate |
|---|
| 87 | |
|---|
| 88 | averaging_interval_pr = 900.0, ! averaging interval profile data |
|---|
| 89 | dt_averaging_input_pr = 0.0, ! averaging profile data sampling rate |
|---|
| 90 | |
|---|
| 91 | / ! end of runtime parameter namelist |
|---|
| 92 | |
|---|
