| 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 = 7, ! number of gridboxes in x-direction (nx+1) |
|---|
| 10 | ny = 23, ! number of gridboxes in y-direction (ny+1) |
|---|
| 11 | nz = 60, ! 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 | ! |
|---|
| 18 | !-- initialization and vertical profiles |
|---|
| 19 | !------------------------------------------------------------------------------- |
|---|
| 20 | initializing_actions = 'set_constant_profiles initialize_bubble', |
|---|
| 21 | ! initial conditions (with warm air bubble) |
|---|
| 22 | |
|---|
| 23 | ug_surface = 0.0, ! u-comp of geostrophic wind at surface |
|---|
| 24 | surface_pressure = 1015.4, ! surface pressure |
|---|
| 25 | |
|---|
| 26 | pt_surface = 297.9, ! temperature at surface |
|---|
| 27 | pt_vertical_gradient = 0.0, 0.585, ! vertical gradient of temperature |
|---|
| 28 | pt_vertical_gradient_level = 0.0, 740.0, ! height level of temp gradients |
|---|
| 29 | |
|---|
| 30 | q_surface = 0.016, ! mixing ratio at surface |
|---|
| 31 | q_vertical_gradient = -2.97E-4, -4.52E-4, -8.1E-5, ! gradient for mix. ratio |
|---|
| 32 | q_vertical_gradient_level = 0.0, 740.0, 3260.0, ! height lev. for gradients |
|---|
| 33 | |
|---|
| 34 | ! |
|---|
| 35 | !-- humdity and microphysics |
|---|
| 36 | !------------------------------------------------------------------------------- |
|---|
| 37 | humidity = .TRUE., ! enables prog. equation for total water mixing ratio |
|---|
| 38 | cloud_droplets = .TRUE., ! instead the lcm is used for microphysics |
|---|
| 39 | |
|---|
| 40 | ! |
|---|
| 41 | !-- numerics |
|---|
| 42 | !------------------------------------------------------------------------------- |
|---|
| 43 | fft_method = 'temperton-algorithm', ! build-in fft method |
|---|
| 44 | |
|---|
| 45 | / ! end of initialization parameter namelist |
|---|
| 46 | |
|---|
| 47 | |
|---|
| 48 | !------------------------------------------------------------------------------- |
|---|
| 49 | !-- RUNTIME PARAMETER NAMELIST |
|---|
| 50 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par |
|---|
| 51 | !------------------------------------------------------------------------------- |
|---|
| 52 | &runtime_parameters |
|---|
| 53 | ! |
|---|
| 54 | !-- run steering |
|---|
| 55 | !------------------------------------------------------------------------------- |
|---|
| 56 | end_time = 950.0, ! simulation time of the 3D model |
|---|
| 57 | dt = 1.0, ! time step (for this setup unfortunately unphysically; |
|---|
| 58 | ! however, a smaller time step is unfeasible for this |
|---|
| 59 | ! test case ) |
|---|
| 60 | netcdf_data_format = 2, ! use NetCDF3 |
|---|
| 61 | |
|---|
| 62 | |
|---|
| 63 | create_disturbances = .TRUE., ! randomly perturbate horiz. velocity |
|---|
| 64 | dt_disturb = 150.0, ! interval for random perturbations |
|---|
| 65 | disturbance_energy_limit = 0.01, ! upper limit for perturbation energy |
|---|
| 66 | |
|---|
| 67 | data_output_2d_on_each_pe = .FALSE., ! don't do 2D output on each MPI rank |
|---|
| 68 | |
|---|
| 69 | ! |
|---|
| 70 | !-- data output |
|---|
| 71 | !------------------------------------------------------------------------------- |
|---|
| 72 | dt_run_control = 0.0, ! output interval for run control |
|---|
| 73 | dt_data_output = 400.0, ! output interval for general data |
|---|
| 74 | dt_dopr = 400.0, ! output interval for profile data |
|---|
| 75 | |
|---|
| 76 | data_output = 'w_xy', |
|---|
| 77 | 'w_xz', |
|---|
| 78 | 'theta_xy', |
|---|
| 79 | 'theta_xz', |
|---|
| 80 | 'ql_yz', |
|---|
| 81 | 'w_yz', |
|---|
| 82 | 'u_yz', |
|---|
| 83 | 'v_yz', |
|---|
| 84 | 'theta_yz', |
|---|
| 85 | |
|---|
| 86 | data_output_pr = 'u', |
|---|
| 87 | 'v', |
|---|
| 88 | 'w', |
|---|
| 89 | 'ql', |
|---|
| 90 | '#theta', |
|---|
| 91 | 'w"theta"', |
|---|
| 92 | 'w*theta*', |
|---|
| 93 | 'wtheta', |
|---|
| 94 | 'w*2', |
|---|
| 95 | 'theta*2', |
|---|
| 96 | |
|---|
| 97 | section_xy = 2, 10, ! grid index for 2D XY cross sections |
|---|
| 98 | section_xz = 20, ! grid index for 2D XZ cross sections |
|---|
| 99 | section_yz = -1, 5, ! grid index for 2D YZ cross sections |
|---|
| 100 | |
|---|
| 101 | averaging_interval_pr = 100.0, ! averaging interval profile data |
|---|
| 102 | dt_averaging_input_pr = 10.0, ! averaging profile data sampling rate |
|---|
| 103 | |
|---|
| 104 | / ! end of runtime parameter namelist |
|---|
| 105 | |
|---|
| 106 | !------------------------------------------------------------------------------- |
|---|
| 107 | !-- PARTICLE PARAMETER NAMELIST |
|---|
| 108 | ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/particle_parameters |
|---|
| 109 | !------------------------------------------------------------------------------- |
|---|
| 110 | &particle_parameters |
|---|
| 111 | ! |
|---|
| 112 | !-- initialize particles in model domain |
|---|
| 113 | !------------------------------------------------------------------------------- |
|---|
| 114 | psb = 20.0, ! bottom of particle source |
|---|
| 115 | pst = 2000.0, ! top of particle source |
|---|
| 116 | |
|---|
| 117 | pdx = 25.0, ! distance between particles along x |
|---|
| 118 | pdy = 25.0, ! distance between particles along y |
|---|
| 119 | pdz = 25.0, ! distance between particles along z |
|---|
| 120 | |
|---|
| 121 | random_start_position = .FALSE., ! add random start positions |
|---|
| 122 | number_of_particle_groups = 1, ! only one particle group is used |
|---|
| 123 | particle_advection_start = 700, ! particles are released after 700s |
|---|
| 124 | |
|---|
| 125 | ! |
|---|
| 126 | !-- boundary conditions |
|---|
| 127 | !------------------------------------------------------------------------------- |
|---|
| 128 | bc_par_b = 'absorb', ! bottom boundary condition for particles |
|---|
| 129 | bc_par_t = 'reflect', ! top boundary condition for particles |
|---|
| 130 | ! |
|---|
| 131 | !-- initialize particle model as lagrangian cloud model |
|---|
| 132 | !------------------------------------------------------------------------------- |
|---|
| 133 | density_ratio = 0.001, ! density ratio of particles (air=1, liquid water = 1000) |
|---|
| 134 | radius = 1.0E-6, ! initial radius of particles |
|---|
| 135 | |
|---|
| 136 | curvature_solution_effects = .FALSE., ! switch off koehler effects |
|---|
| 137 | collision_kernel = 'hall', ! enabales collision |
|---|
| 138 | number_concentration = 50.0E6, ! initialize particles with a weighting factor |
|---|
| 139 | ! such that a number number concentration of |
|---|
| 140 | ! 50 cm^-3 is obtained |
|---|
| 141 | use_sgs_for_particles = .FALSE., ! disable SGS velocities for particles |
|---|
| 142 | |
|---|
| 143 | ! |
|---|
| 144 | !-- particle output |
|---|
| 145 | !------------------------------------------------------------------------------- |
|---|
| 146 | dt_dopts = 60.0, ! time interval for particle timeseries |
|---|
| 147 | write_particle_statistics = .FALSE., ! write particle statistics in file |
|---|
| 148 | dt_write_particle_data = 10000.0, ! interval for particle output (Binary Files) |
|---|
| 149 | ! Attention: they might be very huge, so here |
|---|
| 150 | ! particles are only implicitly outputted |
|---|
| 151 | / ! end of particle parameter namelist |
|---|
| 152 | |
|---|
