| 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, |
|---|
| 29 | 1.0, ! piecewise temp gradients |
|---|
| 30 | pt_vertical_gradient_level = 0.0, |
|---|
| 31 | 800.0, ! height level of temp gradients |
|---|
| 32 | |
|---|
| 33 | ! |
|---|
| 34 | !-- boundary conditions |
|---|
| 35 | !------------------------------------------------------------------------------- |
|---|
| 36 | surface_heatflux = 0.1, ! sensible heat flux at the bottom surface |
|---|
| 37 | bc_pt_b = 'neumann', ! required with surface_heatflux |
|---|
| 38 | |
|---|
| 39 | ! |
|---|
| 40 | !-- numerics |
|---|
| 41 | !------------------------------------------------------------------------------- |
|---|
| 42 | fft_method = 'temperton-algorithm', ! build-in fft method |
|---|
| 43 | |
|---|
| 44 | / ! end of initialization parameter namelist |
|---|
| 45 | |
|---|
| 46 | !------------------------------------------------------------------------------- |
|---|
| 47 | !-- RUNTIME PARAMETER NAMELIST |
|---|
| 48 | ! https://palm.muk.uni-hannover.de/trac/wiki/doc/app/runtime_parameters |
|---|
| 49 | !------------------------------------------------------------------------------- |
|---|
| 50 | &runtime_parameters |
|---|
| 51 | ! |
|---|
| 52 | !-- run steering |
|---|
| 53 | !------------------------------------------------------------------------------- |
|---|
| 54 | end_time = 10800.0, ! simulation time of the 3D model |
|---|
| 55 | |
|---|
| 56 | create_disturbances = .TRUE., ! randomly perturbate horiz. velocity |
|---|
| 57 | dt_disturb = 150.0, ! interval for random perturbations |
|---|
| 58 | disturbance_energy_limit = 0.015, ! upper limit for perturbation energy |
|---|
| 59 | |
|---|
| 60 | data_output_2d_on_each_pe = .FALSE., ! don't do 2D output on each MPI rank |
|---|
| 61 | |
|---|
| 62 | ! |
|---|
| 63 | !-- data output |
|---|
| 64 | !------------------------------------------------------------------------------- |
|---|
| 65 | netcdf_data_format = 2, ! use NetCDF3 |
|---|
| 66 | |
|---|
| 67 | dt_run_control = 0.0, ! output interval for run control |
|---|
| 68 | dt_data_output = 900.0, ! output interval for general data |
|---|
| 69 | dt_data_output_av = 1800.0, ! output interval for averaged data |
|---|
| 70 | dt_do3d = 900.0, ! output interval for 3D data |
|---|
| 71 | dt_do2d_xy = 900.0, ! output interval for 2D XY data |
|---|
| 72 | dt_do2d_xz = 900.0, ! output interval for 2D XZ data |
|---|
| 73 | dt_do2d_yz = 900.0, ! output interval for 2D YZ data |
|---|
| 74 | dt_dopr = 900.0, ! output interval for profile data |
|---|
| 75 | |
|---|
| 76 | data_output = 'u', 'u_av', |
|---|
| 77 | 'v', 'v_av', |
|---|
| 78 | 'w', 'w_av', |
|---|
| 79 | 'theta', 'theta_av', |
|---|
| 80 | 'u_xy', 'u_xz', 'u_yz', 'u_xz_av', |
|---|
| 81 | 'v_xy', 'v_xz', 'v_yz', 'v_xz_av', |
|---|
| 82 | 'w_xy', 'w_xz', 'w_yz', 'w_xz_av', |
|---|
| 83 | 'theta_xy', 'theta_xz', 'theta_yz', |
|---|
| 84 | |
|---|
| 85 | data_output_pr = 'u', 'u*2', |
|---|
| 86 | 'v', 'v*2', |
|---|
| 87 | 'w', 'w*2', 'w*3', |
|---|
| 88 | '#theta', 'theta*2', |
|---|
| 89 | 'wtheta', 'w"theta"', 'w*theta*', |
|---|
| 90 | |
|---|
| 91 | section_xy = 2, |
|---|
| 92 | 10, ! grid index for 2D XY cross sections |
|---|
| 93 | section_xz = 5, |
|---|
| 94 | 20, ! grid index for 2D XZ cross sections |
|---|
| 95 | section_yz = 20, |
|---|
| 96 | 25, ! grid index for 2D YZ cross sections |
|---|
| 97 | |
|---|
| 98 | averaging_interval = 900.0, ! averaging interval general data |
|---|
| 99 | dt_averaging_input = 10.0, ! averaging general data sampling rate |
|---|
| 100 | |
|---|
| 101 | averaging_interval_pr = 900.0, ! averaging interval profile data |
|---|
| 102 | dt_averaging_input_pr = 0.0, ! averaging profile data sampling rate |
|---|
| 103 | |
|---|
| 104 | / ! end of runtime parameter namelist |
|---|
| 105 | |
|---|
