| 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 = 1023, ! Number of gridboxes in x-direction (nx+1) |
|---|
| 10 | ny = 1023, ! Number of gridboxes in y-direction (ny+1) |
|---|
| 11 | nz = 768, ! Number of gridboxes in z-direction (nz) |
|---|
| 12 | |
|---|
| 13 | dx = 0.86E-3, ! Size of single gridbox in x-direction |
|---|
| 14 | dy = 0.86E-3, ! Size of single gridbox in y-direction |
|---|
| 15 | dz = 0.86E-3, ! Size of single gridbox in z-direction |
|---|
| 16 | |
|---|
| 17 | ! |
|---|
| 18 | !-- initialization |
|---|
| 19 | !------------------------------------------------------------------------------- |
|---|
| 20 | initializing_actions = 'set_constant_profiles', ! initial conditions |
|---|
| 21 | |
|---|
| 22 | ug_surface = 0.0, ! u-comp of geostrophic wind at surface |
|---|
| 23 | vg_surface = 0.0, ! v-comp of geostrophic wind at surface |
|---|
| 24 | |
|---|
| 25 | pt_surface = 300.0, ! Initial surface potential temperature |
|---|
| 26 | pt_reference = 288.0, ! Reference temperature to be used in the buoyancy term (in K) |
|---|
| 27 | |
|---|
| 28 | reference_state = 'single_value', ! A constant single potential temperature value is used as reference state |
|---|
| 29 | |
|---|
| 30 | ! |
|---|
| 31 | !-- boundary conditions |
|---|
| 32 | !------------------------------------------------------------------------------- |
|---|
| 33 | bc_lr = 'cyclic', |
|---|
| 34 | bc_ns = 'cyclic', |
|---|
| 35 | |
|---|
| 36 | bc_uv_b = 'dirichlet', ! (no-slip condition) |
|---|
| 37 | bc_uv_t = 'neumann', ! (free-slip condition) |
|---|
| 38 | |
|---|
| 39 | bc_pt_b = 'dirichlet', |
|---|
| 40 | bc_pt_t = 'neumann', |
|---|
| 41 | |
|---|
| 42 | constant_flux_layer = .FALSE., ! Parameter to switch on a constant flux layer at the bottom boundary. |
|---|
| 43 | ! Constant "eddy diffusivities" are only allowed with the constant flux layer switched off |
|---|
| 44 | |
|---|
| 45 | |
|---|
| 46 | ! |
|---|
| 47 | !-- numerics |
|---|
| 48 | !------------------------------------------------------------------------------- |
|---|
| 49 | fft_method = 'fftw', ! http://www.fftw.org/ |
|---|
| 50 | momentum_advec = 'ws-scheme', ! 5th order upwind scheme |
|---|
| 51 | scalar_advec = 'ws-scheme', ! 5th order upwind scheme |
|---|
| 52 | |
|---|
| 53 | ! |
|---|
| 54 | !-- physics |
|---|
| 55 | !------------------------------------------------------------------------------- |
|---|
| 56 | prandtl_number = 1.0, ! Kinematic viscosity/molecular diffusivity for momentum and molecular diffusivity for heat are the same |
|---|
| 57 | km_constant = 1.76E-5, ! Constant "eddy diffusivities" are used. In DNS it's the kinematic viscosity |
|---|
| 58 | omega = 0.0, ! Angular velocity of the rotating system (Earth, in rad/s) - No coriolis force |
|---|
| 59 | latitude = 52.37, ! Geographical latitude (in degrees north) of Hannover |
|---|
| 60 | |
|---|
| 61 | / ! end of initialization parameter namelist |
|---|
| 62 | |
|---|
| 63 | !------------------------------------------------------------------------------- |
|---|
| 64 | !-- RUNTIME PARAMETER NAMELIST |
|---|
| 65 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par |
|---|
| 66 | !------------------------------------------------------------------------------- |
|---|
| 67 | &runtime_parameters |
|---|
| 68 | ! |
|---|
| 69 | !-- run steering |
|---|
| 70 | !------------------------------------------------------------------------------- |
|---|
| 71 | end_time = 10.0, ! simulation time of the 3D model |
|---|
| 72 | |
|---|
| 73 | dt_max = 0.1, ! Maximum allowed value of the time step (in s) |
|---|
| 74 | |
|---|
| 75 | termination_time_needed = 600.0, ! CPU time needed for terminal actions |
|---|
| 76 | |
|---|
| 77 | npex = 32, npey = 16, ! Virtual processor topology. |
|---|
| 78 | |
|---|
| 79 | create_disturbances = .TRUE., ! randomly perturbate horiz. velocity |
|---|
| 80 | dt_disturb = 10, ! interval for random perturbations in s |
|---|
| 81 | disturbance_energy_limit = 0.1E-6, ! upper limit for perturbation energy E* |
|---|
| 82 | disturbance_amplitude = 0.2, ! Maximum perturbation amplitude of the random perturbations in m/s |
|---|
| 83 | |
|---|
| 84 | data_output_2d_on_each_pe = .FALSE., ! don't do 2D output on each MPI rank |
|---|
| 85 | |
|---|
| 86 | ! |
|---|
| 87 | !-- data output |
|---|
| 88 | !------------------------------------------------------------------------------- |
|---|
| 89 | netcdf_data_format = 2, ! use NetCDF3 |
|---|
| 90 | |
|---|
| 91 | dt_run_control = 0.0, ! output interval for run control |
|---|
| 92 | dt_data_output = 1.0, ! output interval for general data |
|---|
| 93 | dt_data_output_av = 1.0, ! output interval for averaged data |
|---|
| 94 | dt_dopr = 1.0, ! output interval for profile data |
|---|
| 95 | |
|---|
| 96 | do3d_at_begin = .TRUE., |
|---|
| 97 | do2d_at_begin = .TRUE., |
|---|
| 98 | |
|---|
| 99 | data_output = 'u', 'u_xy', 'u_xz', |
|---|
| 100 | 'v', 'v_xy', 'v_xz', |
|---|
| 101 | 'w', 'w_xy', 'w_xz', |
|---|
| 102 | 'theta', 'theta_xy', 'theta_xz', |
|---|
| 103 | 'p', 'p_xy', 'p_xz', 'tsurf*_xy', |
|---|
| 104 | |
|---|
| 105 | data_output_pr = '#u', '#v', 'w', |
|---|
| 106 | 'wtheta', 'w"theta"', 'w*theta*', |
|---|
| 107 | 'e', 'e*', |
|---|
| 108 | 'u*2', 'v*2', 'w*2', 'theta*2', |
|---|
| 109 | '#theta', '#km', '#kh', |
|---|
| 110 | 'hyp', 'p', 'rho', |
|---|
| 111 | |
|---|
| 112 | section_xy = 1, 10, 20, 30, 40, |
|---|
| 113 | 50, 60, 70, 80, 90, 100, |
|---|
| 114 | 200, 300, 400, 500, 600, |
|---|
| 115 | 700, ! grid index for 2D XY cross sections |
|---|
| 116 | section_xz = 512, -1, ! grid index for 2D XZ cross sections |
|---|
| 117 | section_yz = 512, -1, ! grid index for 2D YZ cross sections |
|---|
| 118 | |
|---|
| 119 | |
|---|
| 120 | / ! end of runtime parameter namelist |
|---|
| 121 | |
|---|
| 122 | !------------------------------------------------------------------------------- |
|---|
| 123 | !-- RUNTIME PARAMETER NAMELIST |
|---|
| 124 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/userpar |
|---|
| 125 | !------------------------------------------------------------------------------- |
|---|
| 126 | &user_parameters |
|---|
| 127 | |
|---|
| 128 | / ! end of user parameter namelist |
|---|
