[4397] | 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 = 19, ! Number of gridboxes in x-direction (nx+1) |
---|
| 10 | ny = 19, ! Number of gridboxes in y-direction (ny+1) |
---|
| 11 | nz = 50, ! Number of gridboxes in z-direction (nz) |
---|
[4222] | 12 | |
---|
[4397] | 13 | dx = 100.0, ! Size of single gridbox in x-direction |
---|
| 14 | dy = 100.0, ! Size of single gridbox in y-direction |
---|
| 15 | dz = 25.0, ! Size of single gridbox in z-direction |
---|
| 16 | |
---|
[4222] | 17 | dz_stretch_level = 100.0, ! Height (in m) where stretching starts |
---|
| 18 | dz_max = 50.0, ! maximum vertical grid size |
---|
[4397] | 19 | ! |
---|
| 20 | !-- initialization |
---|
| 21 | !------------------------------------------------------------------------------ |
---|
| 22 | initializing_actions = 'set_constant_profiles', ! initial conditions |
---|
| 23 | |
---|
| 24 | ug_surface = 2.5, ! u-comp of geostrophic wind at surface |
---|
| 25 | vg_surface = 0.0, ! v-comp of geostrophic wind at surface |
---|
| 26 | |
---|
| 27 | pt_surface = 293.0, ! initial surface potential temp |
---|
| 28 | |
---|
| 29 | pt_vertical_gradient = 0.0, |
---|
| 30 | 3.0, |
---|
| 31 | 0.0, ! piecewise temp gradients |
---|
| 32 | pt_vertical_gradient_level = 0.0, |
---|
| 33 | 800.0, |
---|
| 34 | 1600.0, ! height level of temp gradients |
---|
| 35 | |
---|
| 36 | q_surface = 0.01, ! mixing ratio --> rel. hum. = ? |
---|
| 37 | |
---|
| 38 | surface_pressure = 1000.0, ! in hPa |
---|
| 39 | |
---|
| 40 | origin_date_time = '2019-06-21 12:00:00 +00' |
---|
| 41 | ! |
---|
| 42 | !-- physics |
---|
| 43 | !------------------------------------------------------------------------------ |
---|
| 44 | latitude = 53.37, ! geographical latitude |
---|
| 45 | longitude = 9.73, ! geographical longitude |
---|
| 46 | ! |
---|
| 47 | !-- boundary conditions |
---|
| 48 | !------------------------------------------------------------------------------ |
---|
| 49 | constant_flux_layer = .TRUE., ! use a constant flux layer at |
---|
| 50 | ! the bottom boundary |
---|
| 51 | |
---|
| 52 | bc_pt_b = 'dirichlet', ! required when using LSM |
---|
| 53 | bc_q_b = 'dirichlet', ! required when using LSM |
---|
| 54 | ! |
---|
| 55 | !-- mode |
---|
| 56 | !------------------------------------------------------------------------------ |
---|
| 57 | humidity = .TRUE., ! use prognostic equation for humidity |
---|
| 58 | |
---|
| 59 | reference_state = 'horizontal_average', ! set reference state for |
---|
| 60 | ! buoyancy term |
---|
| 61 | ! |
---|
| 62 | !-- numerics |
---|
| 63 | !------------------------------------------------------------------------------ |
---|
| 64 | fft_method = 'temperton-algorithm', ! build-in fft method |
---|
| 65 | |
---|
| 66 | rayleigh_damping_height = 1600.0, ! starting height for Rayleigh damping |
---|
| 67 | rayleigh_damping_factor = 0.01, ! strength of Rayleigh damping |
---|
| 68 | |
---|
| 69 | / |
---|
| 70 | |
---|
| 71 | !------------------------------------------------------------------------------- |
---|
| 72 | !-- RUNTIME PARAMETER NAMELIST |
---|
| 73 | ! https://palm.muk.uni-hannover.de/trac/wiki/doc/app/runtime_parameters |
---|
| 74 | !------------------------------------------------------------------------------- |
---|
| 75 | &runtime_parameters |
---|
| 76 | ! |
---|
| 77 | !-- run steering |
---|
| 78 | !------------------------------------------------------------------------------- |
---|
| 79 | end_time = 216000.0, ! simulation time of the 3D model |
---|
| 80 | ! (2.5 days) |
---|
| 81 | |
---|
| 82 | create_disturbances = .T., ! randomly perturbate horiz. velocity |
---|
| 83 | dt_disturb = 120.0, ! interval for random perturbations |
---|
| 84 | disturbance_energy_limit = 0.0001, ! upper limit for perturbation energy |
---|
| 85 | ! |
---|
| 86 | !-- data output |
---|
| 87 | !------------------------------------------------------------------------------ |
---|
| 88 | netcdf_data_format = 2, ! use NetCDF3 |
---|
| 89 | |
---|
| 90 | dt_run_control = 900.0, ! output interval for run control |
---|
| 91 | dt_data_output = 21600.0, ! output interval for general data |
---|
| 92 | dt_data_output_av = 3600.0, ! output interval for averaged data |
---|
| 93 | |
---|
| 94 | averaging_interval = 1800.0, ! averaging interval general data |
---|
| 95 | |
---|
| 96 | data_output_pr = '#u', '#v', '#theta', '#q', '#thetav', |
---|
| 97 | '#km', '#kh', '#l', |
---|
| 98 | 'w', 'e', 'e*', 'p', |
---|
| 99 | 'w"u"', 'w*u*', 'w"v"', 'w*v*', |
---|
| 100 | 'w"theta"', 'w*theta*', 'w"thetav"', |
---|
| 101 | 'w*thetav*', |
---|
| 102 | 'w"q"', 'w*q*', 'w*e*', |
---|
| 103 | 'theta*2', 'q*2', 'u*2', 'v*2', 'w*2', |
---|
| 104 | '#t_soil', '#m_soil', ! Profile output |
---|
| 105 | |
---|
| 106 | section_xy = 0, ! surface variables only |
---|
| 107 | |
---|
| 108 | data_output = 'theta', 'theta_av', |
---|
| 109 | 'thetav', 'thetav_av', |
---|
| 110 | 'q', 'q_av', |
---|
| 111 | 'u', 'u_av', |
---|
| 112 | 'v', 'v_av', |
---|
| 113 | |
---|
| 114 | 't_soil', 't_soil_av', |
---|
| 115 | 'm_soil', 'm_soil_av', |
---|
| 116 | |
---|
| 117 | 'us*_xy', 'us*_xy_av', |
---|
| 118 | 't*_xy', 't*_xy_av', |
---|
| 119 | 'r_a*_xy', 'r_a*_xy_av', |
---|
| 120 | 'r_s*_xy', 'r_s*_xy_av', |
---|
| 121 | 'tsurf*_xy', 'tsurf*_xy_av' |
---|
| 122 | |
---|
| 123 | 'ghf*_xy', 'ghf*_xy_av', |
---|
| 124 | 'shf*_xy', 'shf*_xy_av', |
---|
| 125 | |
---|
| 126 | 'qsws_liq*_xy', 'qsws_liq*_xy_av', |
---|
| 127 | 'qsws_soil*_xy', 'qsws_soil*_xy_av', |
---|
| 128 | 'qsws_veg*_xy', 'qsws_veg*_xy_av', |
---|
| 129 | |
---|
| 130 | 'c_liq*_xy', 'c_liq*_xy_av', |
---|
| 131 | 'c_soil*_xy', 'c_soil*_xy_av', |
---|
| 132 | 'c_veg*_xy', 'c_veg*_xy_av', ! 2d and/or |
---|
| 133 | ! 3d output |
---|
| 134 | / |
---|
| 135 | |
---|
| 136 | !------------------------------------------------------------------------------- |
---|
| 137 | !-- LAND-SURFACE PARAMETER NAMELIST |
---|
| 138 | ! https://palm.muk.uni-hannover.de/trac/wiki/doc/app/land_surface_parameters |
---|
| 139 | !------------------------------------------------------------------------------- |
---|
| 140 | &land_surface_parameters |
---|
| 141 | soil_type = 3, ! type to be used in soil model |
---|
| 142 | soil_moisture = 0.35, 0.35, 0.35, 0.35, ! moisture at every |
---|
| 143 | 0.35, 0.35, 0.35, 0.35, ! soil level |
---|
| 144 | soil_temperature = 293.0, 293.0, 293.0, 293.0, ! temperature at |
---|
| 145 | 293.0, 293.0, 293.0, 293.0, ! every soil level |
---|
| 146 | deep_soil_temperature = 293.0, ! temperature of deep soil |
---|
| 147 | |
---|
| 148 | surface_type = 'vegetation', ! surface classification |
---|
| 149 | vegetation_type = 3, ! vegetation to be used in surface model |
---|
| 150 | |
---|
| 151 | vegetation_coverage = 0.95, ! part of surface covered by vegetation |
---|
| 152 | ! (0 - 1) |
---|
| 153 | c_surface = 0.0, ! surface heat capacity |
---|
| 154 | conserve_water_content = .TRUE., ! closed (T) or open (F) lower boundary |
---|
| 155 | |
---|
| 156 | / |
---|
| 157 | |
---|
| 158 | !------------------------------------------------------------------------------- |
---|
| 159 | !-- RADIATION PARAMETER NAMELIST |
---|
| 160 | ! https://palm.muk.uni-hannover.de/trac/wiki/doc/app/radiation_parameters |
---|
| 161 | !------------------------------------------------------------------------------- |
---|
| 162 | &radiation_parameters |
---|
| 163 | radiation_scheme = 'clear-sky', ! Radiation scheme to be used |
---|
| 164 | dt_radiation = 60.0, ! Time step of the radiation model |
---|
| 165 | |
---|
| 166 | / |
---|