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