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