- Timestamp:
- Mar 27, 2020 2:06:23 PM (5 years ago)
- Location:
- palm/trunk
- Files:
-
- 5 edited
-
SOURCE/salsa_mod.f90 (modified) (4 diffs)
-
TESTS/cases/urban_environment_salsa/INPUT/urban_environment_salsa_p3d (modified) (15 diffs)
-
TESTS/cases/urban_environment_salsa/MONITORING/urban_environment_salsa_rc.000 (modified) (3 diffs)
-
TESTS/cases/urban_environment_salsa/OUTPUT/urban_environment_salsa_pr.000.nc (modified) (previous)
-
TESTS/cases/urban_environment_salsa/OUTPUT/urban_environment_salsa_ts.000.nc (modified) (previous)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/salsa_mod.f90
r4457 r4478 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Bug fixes: 29 ! - call salsa_driver in salsa_init also for the ghost points 30 ! - decycle flags missing from advec_s_ws call in salsa_tendency 31 ! 32 ! 4457 2020-03-11 14:20:43Z raasch 28 33 ! use statement for exchange horiz added 29 34 ! … … 1810 1815 1811 1816 !-- Initalisation run of SALSA + calculate the vertical top index of the topography 1812 DO i = nxl , nxr1813 DO j = nys , nyn1817 DO i = nxlg, nxrg 1818 DO j = nysg, nyng 1814 1819 1815 1820 k_topo_top(j,i) = MAXLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,j,i), 12 ) ), & … … 8095 8100 IF ( timestep_scheme(1:5) == 'runge' ) THEN 8096 8101 IF ( ws_scheme_sca ) THEN 8097 CALL advec_s_ws( salsa_advc_flags_s, i, j, rs, id, flux_s, diss_s, flux_l, diss_l, & 8098 i_omp_start, tn, bc_dirichlet_l .OR. bc_radiation_l, & 8099 bc_dirichlet_n .OR. bc_radiation_n, & 8100 bc_dirichlet_r .OR. bc_radiation_r, & 8101 bc_dirichlet_s .OR. bc_radiation_s, monotonic_limiter_z ) 8102 CALL advec_s_ws( salsa_advc_flags_s, i, j, rs, id, & 8103 flux_s, diss_s, flux_l, diss_l, i_omp_start, tn, & 8104 bc_dirichlet_l .OR. bc_radiation_l .OR. decycle_salsa_lr, & 8105 bc_dirichlet_n .OR. bc_radiation_n .OR. decycle_salsa_ns, & 8106 bc_dirichlet_r .OR. bc_radiation_r .OR. decycle_salsa_lr, & 8107 bc_dirichlet_s .OR. bc_radiation_s .OR. decycle_salsa_ns ) 8102 8108 ELSE 8103 8109 CALL advec_s_pw( i, j, rs ) … … 8230 8236 IF ( timestep_scheme(1:5) == 'runge' ) THEN 8231 8237 IF ( ws_scheme_sca ) THEN 8232 CALL advec_s_ws( salsa_advc_flags_s, rs, id, bc_dirichlet_l .OR. bc_radiation_l, & 8233 bc_dirichlet_n .OR. bc_radiation_n, & 8234 bc_dirichlet_r .OR. bc_radiation_r, & 8235 bc_dirichlet_s .OR. bc_radiation_s ) 8238 CALL advec_s_ws( salsa_advc_flags_s, rs, id, & 8239 bc_dirichlet_l .OR. bc_radiation_l .OR. decycle_salsa_lr, & 8240 bc_dirichlet_n .OR. bc_radiation_n .OR. decycle_salsa_ns, & 8241 bc_dirichlet_r .OR. bc_radiation_r .OR. decycle_salsa_lr, & 8242 bc_dirichlet_s .OR. bc_radiation_s .OR. decycle_salsa_ns ) 8236 8243 ELSE 8237 8244 CALL advec_s_pw( rs ) -
palm/trunk/TESTS/cases/urban_environment_salsa/INPUT/urban_environment_salsa_p3d
r4342 r4478 20 20 !------------------------------------------------------------------------------ 21 21 fft_method = 'fftw', 22 monotonic_limiter_z = .T., 22 monotonic_limiter_z = .T., 23 23 ! 24 24 !-- Mode 25 !------------------------------------------------------------------------------ 25 !------------------------------------------------------------------------------ 26 26 humidity = .T., 27 27 q_surface = 0.001, … … 30 30 !-- Initialization 31 31 !------------------------------------------------------------------------------ 32 initializing_actions = 'inifor set_constant_profiles', !'set_constant_profiles', or 'inifor', 32 initializing_actions = 'inifor set_constant_profiles', !'set_constant_profiles', or 'inifor', 33 33 ! or both (met. from inifor, chem constant profiles) 34 34 … … 41 41 ! 42 42 !-- Wall/soil spinup 43 !------------------------------------------------------------------------------ 43 !------------------------------------------------------------------------------ 44 44 spinup_time = 1200.0, 45 45 spinup_pt_mean = 284.15, … … 50 50 ! 51 51 !-- Topography 52 !------------------------------------------------------------------------------ 52 !------------------------------------------------------------------------------ 53 53 topography = 'read_from_file', 54 54 … … 95 95 'e', 'e*', '#km', '#l', 96 96 '#t_soil', '#m_soil', 'rad_lw_in', 'rad_lw_out', 97 'salsa_Ntot', 'salsa_N_UFP', 'salsa_LDSA', 98 'salsa_PM0.1', 'salsa_PM2.5', 'salsa_PM10', 97 'salsa_Ntot', 'salsa_N_UFP', 'salsa_LDSA', 98 'salsa_PM0.1', 'salsa_PM2.5', 'salsa_PM10', 99 99 ! 100 100 !-- 2D/3D output settings 101 101 !------------------------------------------------------------------------------ 102 do3d_at_begin = .T., 102 do3d_at_begin = .T., 103 103 do2d_at_begin = .T., 104 104 105 dt_data_output = 60.0, 105 dt_data_output = 60.0, 106 106 dt_data_output_av = 60.0, 107 107 averaging_interval = 60.0, … … 113 113 ! 'Ntot', 'PM2.5', 'LDSA', 114 114 ! 's_OC','g_OCNV','g_OCSV', 115 ! 'N_bin3', 'm_bin4', 115 ! 'N_bin3', 'm_bin4', 116 116 117 117 … … 143 143 !-- LAND SURFACE MODEL PARAMETER NAMELIST 144 144 ! Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/lsmpar 145 !------------------------------------------------------------------------------- 146 &land_surface_parameters 145 !------------------------------------------------------------------------------- 146 &land_surface_parameters 147 147 148 148 constant_roughness = .T., … … 154 154 155 155 root_fraction = 0.49, 0.26, 0.25, 0.0, 0.0, 0.0, 0.0, 0.0, 156 soil_temperature = 277.15, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 156 soil_temperature = 277.15, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 277.0, 157 157 soil_moisture = 0.30, 0.30, 0.30, 0.30, 0.30, 0.30, 0.30, 0.30, 158 158 deep_soil_temperature = 277.0, 159 159 160 160 / ! end of lsm_par namelist 161 161 … … 171 171 172 172 / ! end of urban_surface_par namelist 173 173 174 174 175 175 !------------------------------------------------------------------------------- … … 178 178 !------------------------------------------------------------------------------- 179 179 &plant_canopy_parameters 180 180 181 181 canopy_mode = 'read_from_file', 182 182 canopy_drag_coeff = 0.3, 183 183 184 184 / ! end of canopy_par namelist 185 185 186 186 187 187 !------------------------------------------------------------------------------- 188 188 !-- SALSA PARAMETER NAMELIST 189 189 !------------------------------------------------------------------------------- 190 &salsa_parameters 190 &salsa_parameters 191 191 192 192 ! … … 196 196 197 197 ! 198 !-- If initializing_actions includes 'set_constant_profiles' , 198 !-- If initializing_actions includes 'set_constant_profiles' , 199 199 !-- define the initial aerosol size distribution by dpg, sigmag and n_lognorm 200 200 dpg = 13.5E-9, 54.0E-9, 864.1E-9, ! mean diameter per mode (in metres) 201 sigmag = 1.8, 2.16, 2.21, ! standard deviation per mode 201 sigmag = 1.8, 2.16, 2.21, ! standard deviation per mode 202 202 n_lognorm = 1.43E9, 4.45E8, 7.0E4, ! number concentration per mode (#/m3) 203 203 204 204 ! 205 !-- If initializing_actions includes 'set_constant_profiles', apply these initial gas concentrations 205 !-- If initializing_actions includes 'set_constant_profiles', apply these initial gas concentrations 206 206 H2SO4_init = 5.0E12, ! sulphuric acid (#/m3) 207 207 HNO3_init = 3.0E12, ! nitric acid (#/m3) 208 NH3_init = 6.0E12, ! ammonia (#/m3) 208 NH3_init = 6.0E12, ! ammonia (#/m3) 209 209 OCNV_init = 1.0E12, ! non-volatile organic gases (#/m3) 210 210 OCSV_init = 1.0E12, ! non-volatile organic gases (#/m3) 211 211 212 212 ! 213 !-- List of activated chemical components: 213 !-- List of activated chemical components: 214 214 !-- NOTE! Chemical species have to be activated here even if they are not initially present! 215 listspec = 'OC','NO','NH','','','','', 216 ! listspec = 'SO4','OC','BC','DU','SS','NO','NH', 215 listspec = 'OC','NO','NH','','','','', 216 ! listspec = 'SO4','OC','BC','DU','SS','NO','NH', 217 217 218 218 ! 219 219 !-- Sectional presentation of the particle size distribution 220 220 reglim = 3.0E-9, 1.0E-8, 2.5E-6, ! limits of the subranges (m) 221 nbin = 1, 7, ! number of bins per subrange 221 nbin = 1, 7, ! number of bins per subrange 222 222 !-- NOTE! Subrange 1 consists only of H2SO4 and/or OC 223 nf2a = 1.0, ! Number fraction allocated to subrange 2a223 nf2a = 1.0, ! Number fraction allocated to subrange 2a 224 224 225 225 ! 226 226 !-- Aerosol emissions: 227 227 salsa_emission_mode = 'read_from_file', ! 'no_emission','uniform' or 'read_from_file' 228 !-- NOTE! chemical components of the source have to be activated in 'listspec' 228 !-- NOTE! chemical components of the source have to be activated in 'listspec' 229 229 !-- and have to be in the same order 230 230 … … 237 237 ! 238 238 !-- Boundary conditions 239 !-- Decycle at the left & right boundaries: initial concentration 240 !-- is copied to ghost and first three layers at both boundaries at the left boundary 241 !-- and a zero gradient is set at the right boundary 242 !-- Decycle at the north & south boundaries: initial concentration 243 !-- is copied to ghost and first three layers at both boundaries 244 decycle_salsa_lr = .T., 239 !-- Decycle at the left & right boundaries: 240 decycle_salsa_lr = .T., 245 241 decycle_salsa_ns = .F., 246 decycle_method_salsa = 'dirichlet','dirichlet','dirichlet','dirichlet', 242 decycle_method_salsa = 'dirichlet','dirichlet','dirichlet','dirichlet', 247 243 bc_salsa_b = 'neumann', ! surface flux requires 'neumann' 248 244 bc_salsa_t = 'dirichlet', ! top 249 245 250 246 ! 251 247 !-- Switches for aerosol processes: … … 259 255 nlcoag = .F., ! Coagulation master switch 260 256 nsnucl = 0, ! Nucleation scheme (0 = off) 261 nj3 = 1, ! J3 parametrization for nucleation 257 nj3 = 1, ! J3 parametrization for nucleation 262 258 263 259 ! … … 265 261 depo_pcm_par = 'zhang2001', ! plant canopy 266 262 depo_surf_par = 'zhang2001', ! ground, roofs and walls 267 season_z01 = 5 ! Season for zhang2001, 263 season_z01 = 5 ! Season for zhang2001, 268 264 ! 1 = summer, 2 = autumn, 3 = late autumn 269 265 ! 4 = winter, 5 = transitional spring 270 266 ! 271 267 !-- Other switches: 272 advect_particle_water = .T., ! particle water: advect or calculate at each dt_salsa 268 advect_particle_water = .T., ! particle water: advect or calculate at each dt_salsa 273 269 feedback_to_palm = .F., ! feedback to flow due to condensation of water 274 270 nesting_salsa = .F., ! apply self-nesting for salsa variables … … 276 272 write_binary_salsa = .F., ! skip writing restart data even if it's done for the flow 277 273 nesting_offline_salsa = .F., ! apply offline nesting for salsa 278 274 279 275 / ! end of salsa_par namelist 280 281 -
palm/trunk/TESTS/cases/urban_environment_salsa/MONITORING/urban_environment_salsa_rc.000
r4446 r4478 1 1 2 2 ****************************** -------------------------------------------- 3 * PALM 6.0 Rev: 4 379M * atmosphere - run without 1D - prerun3 * PALM 6.0 Rev: 4477M * atmosphere - run without 1D - prerun 4 4 ****************************** -------------------------------------------- 5 5 6 Date: 2020-0 1-17 Run: urban_environment_salsa7 Time: 1 8:25:39Run-No.: 006 Date: 2020-03-27 Run: urban_environment_salsa 7 Time: 15:56:23 Run-No.: 00 8 8 Run on host: default 9 9 Number of PEs: 4 Processor grid (x,y): ( 2, 2) calculated … … 13 13 ----------------- 14 14 15 --> Use the Moeng_Wyngaardturbulence closure (LES mode).15 --> Use the 1.5-order turbulence closure (LES mode). 16 16 --> Use the boussinesq approximation for the model equations. 17 17 --> Solve perturbation pressure via FFT using fftw routines … … 190 190 --------------------------- 191 191 192 Mixing length is limited to 1.80 * z192 Mixing length is limited close to surfaces 193 193 194 194
Note: See TracChangeset
for help on using the changeset viewer.
