Changeset 3611
- Timestamp:
- Dec 7, 2018 2:14:11 PM (6 years ago)
- Location:
- palm/trunk/SOURCE
- Files:
-
- 3 edited
-
chem_emissions_mod.f90 (modified) (31 diffs)
-
chem_modules.f90 (modified) (3 diffs)
-
chemistry_model_mod.f90 (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/chem_emissions_mod.f90
r3589 r3611 22 22 ! Current revisions: 23 23 ! ------------------ 24 ! 24 ! 25 25 ! 26 26 ! Former revisions: 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Code update to comply PALM coding rules 30 ! removed unnecessary informative messsages/location 31 ! messages and corrected comments on PM units from to kg 32 ! bug fix: spcs_name replaced by nvar in DO loops 33 ! 34 ! 3591 2018-11-30 17:37:32Z suehring 29 35 ! - Removed salsa dependency. 30 36 ! - Enabled PARAMETRIZED mode for default surfaces when LSM is not applied but 31 37 ! salsa is (M. Kurppa) 32 ! 38 ! 33 39 ! 3582 2018-11-29 19:16:36Z suehring 34 40 ! resler: … … 67 73 ! Authors: 68 74 ! -------- 69 ! @author Emmanuele Russo (F u-Berlin)75 ! @author Emmanuele Russo (FU-Berlin) 70 76 ! @author Sabine Banzhaf (FU-Berlin) 71 77 ! @author Martijn Schaap (FU-Berlin, TNO Utrecht) … … 79 85 !> @todo revise indices of files read from the netcdf: preproc_emission_data and expert_emission_data 80 86 !> @todo for now emission data may be passed on a singular vertical level: need to be more flexible 87 !> @todo fill/activate restart option in chem_emissions_init 88 !> @todo discuss dt_emis 81 89 !> @note <Enter notes on the module> 82 90 !> @bug <Enter known bugs here> … … 86 94 87 95 USE arrays_3d, & 88 ONLY:rho_air96 ONLY: rho_air 89 97 90 98 USE control_parameters, & 91 ONLY: initializing_actions, end_time, message_string,&92 intermediate_timestep_count, dt_3d99 ONLY: end_time, message_string, initializing_actions, & 100 intermediate_timestep_count, dt_3d 93 101 94 102 USE indices … … 96 104 USE kinds 97 105 98 #if defined 99 USE NETCDF106 #if defined( __netcdf ) 107 USE netcdf 100 108 #endif 101 109 102 110 USE netcdf_data_input_mod, & 103 ONLY: chem_emis_att_type, chem_emis_val_type111 ONLY: chem_emis_att_type, chem_emis_val_type 104 112 105 113 USE date_and_time_mod, & 106 ONLY: time_default_indices, time_preprocessed_indices, & 107 month_of_year, day_of_month, hour_of_day, & 108 index_mm, index_dd, index_hh 109 114 ONLY: day_of_month, hour_of_day, & 115 index_mm, index_dd, index_hh, & 116 month_of_year, hour_of_day, & 117 time_default_indices, time_preprocessed_indices 118 110 119 USE chem_gasphase_mod, & 111 ONLY: spc_names120 ONLY: spc_names 112 121 113 122 USE chem_modules 114 123 115 124 USE statistics, & 116 ONLY: weight_pres 117 125 ONLY: weight_pres 126 127 118 128 IMPLICIT NONE 119 129 130 ! 120 131 !-- Declare all global variables within the module 121 132 122 CHARACTER (LEN=80) :: filename_emis !> Variable for the name of the netcdf input file 123 124 INTEGER(iwp) :: i !> index 1st selected dimension (some dims are not spatial) 125 INTEGER(iwp) :: j !> index 2nd selected dimension 126 INTEGER(iwp) :: i_start !> Index to start read variable from netcdf along one dims 127 INTEGER(iwp) :: i_end !> Index to end read variable from netcdf in one dims 128 INTEGER(iwp) :: j_start !> Index to start read variable from netcdf in additional dims 129 INTEGER(iwp) :: j_end !> Index to end read variable from netcdf in additional dims 130 INTEGER(iwp) :: z_start !> Index to start read variable from netcdf in additional dims 131 INTEGER(iwp) :: z_end !> Index to end read variable from netcdf in additional dims 132 INTEGER(iwp) :: dt_emis !> Time Step Emissions 133 INTEGER(iwp) :: len_index !> length of index (used for several indices) 134 INTEGER(iwp) :: len_index_voc !> length of voc index 135 INTEGER(iwp) :: len_index_pm !> length of PMs index 136 REAL(wp) :: con_factor !> Units Conversion Factor 137 138 139 ! --------------------------------------------------------------- 140 ! Set Values of molecules, mols, etc 141 ! --------------------------------------------------------------- 142 143 !> Avogadro number 144 REAL, PARAMETER :: Avog = 6.02205e23 ! mlc/mol 145 146 !> Dobson units: 147 REAL, PARAMETER :: Dobs = 2.68668e16 ! (mlc/cm2) / DU 148 149 !> sesalt composition: 150 ! (Seinfeld and Pandis, "Atmospheric Chemistry and Physics", 151 ! table 7.8 "Composition of Sea-Salt", p. 444) 152 REAL, PARAMETER :: massfrac_Cl_in_seasalt = 0.5504 ! (kg Cl )/(kg seasalt) 153 REAL, PARAMETER :: massfrac_Na_in_seasalt = 0.3061 ! (kg Na )/(kg seasalt) 154 REAL, PARAMETER :: massfrac_SO4_in_seasalt = 0.0768 ! (kg SO4)/(kg seasalt) 155 156 !> other numbers 157 REAL, PARAMETER :: xm_seasalt = 74.947e-3 ! kg/mol : NaCl, SO4, .. 158 REAL, PARAMETER :: rho_seasalt = 2.2e3 ! kg/m3 159 160 !> * amonium sulphate 161 162 REAL, PARAMETER :: xm_NH4HSO4 = xm_NH4 + xm_H + xm_SO4 ! kg/mol 163 REAL, PARAMETER :: rho_NH4HSO4a = 1.8e3 ! kg/m3 164 165 ! --------------------------------------------------------------- 166 ! gas 167 ! --------------------------------------------------------------- 168 169 !> gas constant 170 REAL, PARAMETER :: Rgas = 8.3144 ! J/mol/K 171 172 !> gas constant for dry air 173 REAL, PARAMETER :: Rgas_air = Rgas / xm_air ! J/kg/K 174 175 !> water vapour 176 REAL, PARAMETER :: Rgas_h2o = Rgas / xm_h2o ! J/kg/K 177 178 !> standard pressure 179 REAL, PARAMETER :: p0 = 1.0e5 ! Pa 180 181 !> sea-level pressure 182 REAL, PARAMETER :: p_sealevel = 1.01325e5 ! Pa <-- suggestion Bram Bregman 183 184 !> global mean pressure 185 REAL, PARAMETER :: p_global = 98500.0 ! Pa 186 187 !> specific heat of dry air at constant pressure 188 REAL, PARAMETER :: cp_air = 1004.0 ! J/kg/K 189 190 !> Latent heat of evaporation 191 REAL, PARAMETER :: lvap = 2.5e6 ! [J kg-1] 192 193 !> Latent heat of condensation at 0 deg Celcius 194 ! (heat (J) necesarry to evaporate 1 kg of water) 195 REAL, PARAMETER :: Lc = 22.6e5 ! J/kg 196 197 !> kappa = R/cp = 2/7 198 REAL, PARAMETER :: kappa = 2.0/7.0 199 200 !> Von Karman constant (dry_dep) 201 REAL, PARAMETER :: vkarman = 0.4 202 203 !> Boltzmann constant: 204 REAL(wp), PARAMETER :: kbolz = 1.38066e-23_wp ! J/K 205 206 !> Inverse Reference Pressure (1/Pa) 207 REAL(wp), PARAMETER :: pref_i = 1.0_wp / 100000.0_wp 208 209 !> 210 REAL(wp), PARAMETER :: r_cp = 0.286_wp !< R / cp (exponent for potential temperature) 133 CHARACTER (LEN=80) :: filename_emis !< Variable for the name of the netcdf input file 134 135 INTEGER(iwp) :: i !< index 1st selected dimension (some dims are not spatial) 136 INTEGER(iwp) :: j !< index 2nd selected dimension 137 INTEGER(iwp) :: i_start !< Index to start read variable from netcdf along one dims 138 INTEGER(iwp) :: i_end !< Index to end read variable from netcdf in one dims 139 INTEGER(iwp) :: j_start !< Index to start read variable from netcdf in additional dims 140 INTEGER(iwp) :: j_end !< Index to end read variable from netcdf in additional dims 141 INTEGER(iwp) :: z_start !< Index to start read variable from netcdf in additional dims 142 INTEGER(iwp) :: z_end !< Index to end read variable from netcdf in additional dims 143 INTEGER(iwp) :: dt_emis !< Time Step Emissions 144 INTEGER(iwp) :: len_index !< length of index (used for several indices) 145 INTEGER(iwp) :: len_index_voc !< length of voc index 146 INTEGER(iwp) :: len_index_pm !< length of PMs index 147 148 REAL(wp) :: con_factor !< Units Conversion Factor 149 150 REAL(wp), PARAMETER :: Rgas = 8.3144 !< gas constant in J/mol/K 151 REAL(wp), PARAMETER :: pref_i = 1.0_wp / 100000.0_wp !< Inverse Reference Pressure (1/Pa) 152 REAL(wp), PARAMETER :: r_cp = 0.286_wp !< R / cp (exponent for potential temperature) 211 153 212 154 213 155 SAVE 214 156 215 !-- Interfaces Part 216 157 158 ! 217 159 !-- Checks Input parameters 218 160 INTERFACE chem_emissions_check_parameters 219 161 MODULE PROCEDURE chem_emissions_check_parameters 220 162 END INTERFACE chem_emissions_check_parameters 221 222 !-- Reading of NAMELIST parameters 223 ! INTERFACE chem_emissions_parin 224 ! MODULE PROCEDURE chem_emissions_parin 225 ! END INTERFACE chem_emissions_parin 226 227 !-- Output of information to the header file 228 ! INTERFACE chem_emissions_header 229 ! MODULE PROCEDURE chem_emissions_header 230 ! END INTERFACE chem_emissions_header 231 163 ! 232 164 !-- Matching Emissions actions 233 165 INTERFACE chem_emissions_match 234 166 MODULE PROCEDURE chem_emissions_match 235 167 END INTERFACE chem_emissions_match 236 168 ! 237 169 !-- Initialization actions 238 170 INTERFACE chem_emissions_init 239 171 MODULE PROCEDURE chem_emissions_init 240 172 END INTERFACE chem_emissions_init 241 173 ! 242 174 !-- Setup of Emissions 243 175 INTERFACE chem_emissions_setup … … 246 178 247 179 248 !-- Public Interfaces 249 PUBLIC chem_emissions_init, chem_emissions_match,chem_emissions_setup250 180 181 PUBLIC chem_emissions_init, chem_emissions_match, chem_emissions_setup 182 ! 251 183 !-- Public Variables 252 253 PUBLIC con_factor, len_index,len_index_voc,len_index_pm 184 PUBLIC con_factor, len_index, len_index_pm, len_index_voc 254 185 255 186 CONTAINS … … 263 194 264 195 265 !TBD: Where Should we put the call to chem_emissions_check_parameters? In chem_init or in check_parameters?266 267 196 IMPLICIT NONE 268 197 269 INTEGER(iwp) :: tmp 270 271 TYPE(chem_emis_att_type) :: emt 272 ! 273 !-- Call routine for controlling chem_emissions namelist 274 ! CALL chem_emissions_parin 275 276 !TBD: In case the given emission values do not coincide with the passed names we should think of a solution: 277 ! I would personally do that if the name of the species differ from the number of emission values, I would 278 ! print a warning that says that in correspondance of that particular species the value is zero. 279 ! An alternative would be to initialize the cs_surface_flux values to a negative number. 280 281 !-- Check Emission Species Number equal to number of passed names for the chemistry species: 282 IF ( SIZE(emt%species_name) /= emt%nspec ) THEN 283 284 message_string = 'Numbers of input emission species names and number of species' // & 198 INTEGER(iwp) :: tmp 199 200 TYPE(chem_emis_att_type) :: emt 201 202 ! 203 !-- Check Emission Species Number equal to number of passed names for the chemistry species: 204 IF ( SIZE(emt%species_name) /= emt%nspec ) THEN 205 206 message_string = 'Numbers of input emission species names and number of species' // & 285 207 'for which emission values are given do not match' 286 208 CALL message( 'chem_emissions_check_parameters', 'CM0437', 2, 2, 0, 6, 0 ) 287 288 209 289 210 ENDIF 290 291 !-- Check Emission Species Number equals to number of passed names for the species 292 !IF ( SIZE(emt%species_name) /= SIZE(emt%species_index) ) THEN 293 ! message_string = 'Number of input emission species names and ' // & 294 ! ' number of passed species indices do not match' 295 ! CALL message( 'chem_emissions_check_parameters', 'CM0101', 2, 2, 0, 6, 0 ) 296 297 !ENDIF 298 299 !-- Check Emission Categories 300 ! IF ( SIZE(chem_emis%cat_name) /= SIZE(chem_emis%cat_index) ) THEN 301 ! WRITE( message_string, * ) & 302 ! 'Number of input emission categories name and categories index does not match\\' 303 ! CALL message( 'chem_emissions_check_parameters', 'CM0101', 1, 2, 0, 6, 0 ) 304 ! ENDIF 305 306 307 308 !TBD: Check which other consistency controls should be included 309 310 !TBD: Include check for spatial dimension of netcdf file. If they do not match with the ones 311 ! of the simulation, the model should print an error. 312 313 END SUBROUTINE chem_emissions_check_parameters 211 212 213 END SUBROUTINE chem_emissions_check_parameters 314 214 315 215 !------------------------------------------------------------------------------! … … 317 217 ! ------------ 318 218 !> Matching the chemical species indices. The routine checks what are the indices of the emission input species 319 ! and the corresponding ones of the model species. The routine gives as output a vector containing the number 320 ! of common species: it is important to note that while the model species are distinct, their values could be 321 ! given to a single species in input: for example, in the case of NO2 and NO, values may be passed in input as NOx values. 219 !> and the corresponding ones of the model species. The routine gives as output a vector containing the number 220 !> of common species: it is important to note that while the model species are distinct, their values could be 221 !> given to a single species in input: for example, in the case of NO2 and NO, values may be passed in input as 222 !> NOx values. 322 223 !------------------------------------------------------------------------------! 323 SUBROUTINE chem_emissions_match( emt_att,len_index)324 325 326 INTEGER(iwp), INTENT(INOUT) :: len_index!< Variable where to store the number of common species between the input dataset and the model species327 328 TYPE(chem_emis_att_type), INTENT(INOUT) :: emt_att!< Chemistry Emission Array containing information for all the input chemical emission species224 SUBROUTINE chem_emissions_match( emt_att,len_index ) 225 226 227 INTEGER(iwp), INTENT(INOUT) :: len_index !< Variable where to store the number of common species between the input dataset and the model species 228 229 TYPE(chem_emis_att_type), INTENT(INOUT) :: emt_att !< Chemistry Emission Array containing information for all the input chemical emission species 329 230 330 INTEGER(iwp) :: ind_mod, ind_inp !< Parameters for cycling through chemical model and input species 331 INTEGER(iwp) :: nspec_emis_inp !< Variable where to store the number of the emission species in input 332 333 INTEGER(iwp) :: ind_voc !< Indices to check whether a split for voc should be done 334 335 INTEGER(iwp) :: ispec !> index for cycle over effective number of emission species 336 337 338 !> Tell the user where we are!! 231 INTEGER(iwp) :: ind_mod, ind_inp !< Parameters for cycling through chemical model and input species 232 INTEGER(iwp) :: nspec_emis_inp !< Variable where to store the number of the emission species in input 233 INTEGER(iwp) :: ind_voc !< Indices to check whether a split for voc should be done 234 INTEGER(iwp) :: ispec !< index for cycle over effective number of emission species 235 236 339 237 CALL location_message( 'Matching input emissions and model chemistry species', .FALSE. ) 340 238 341 ! > Number of input emission species.342 239 ! 240 !-- Number of input emission species. 343 241 nspec_emis_inp=emt_att%nspec 344 242 345 ! > Check the emission mode: DEFAULT, PRE-PROCESSED or PARAMETERIZED !TBD: Add option for capital or small letters346 SELECT CASE(TRIM(mode_emis))347 348 !> PRE-PROCESSED case: in this case the input species have to coincide with the ones of the model (except VOCs for which we have the VOC split): NO and NO2 in input and not NOx. 349 CASE ("PRE-PROCESSED")350 351 CALL location_message( 'chem_emissions PRE-PROCESSED_mode_matching', .FALSE.)352 353 len_index =0354 len_index_voc =0355 356 !> The first condition is that both the number of model and input emissions species are not null357 IF ( (SIZE(spc_names) .GT. 0) .AND. (nspec_emis_inp .GT. 0)) THEN358 359 !> Cycle over model species360 DO ind_mod = 1, SIZE(spc_names)361 ! > Cycle over Input species243 ! 244 !-- Check the emission mode: DEFAULT, PRE-PROCESSED or PARAMETERIZED 245 SELECT CASE( TRIM( mode_emis ) ) 246 247 ! 248 !-- PRE-PROCESSED mode 249 CASE ( "PRE-PROCESSED" ) 250 251 len_index = 0 252 len_index_voc = 0 253 254 IF ( nvar > 0 .AND. (nspec_emis_inp > 0) ) THEN 255 ! 256 !-- Cycle over model species 257 DO ind_mod = 1, nvar 258 ! 259 !-- Cycle over input species 362 260 DO ind_inp = 1, nspec_emis_inp 363 261 364 !> In the PRE-PROCESSED mode each emission species is given input values, made exception for the VOCs, having the total number of VOCs in input: the model then executes a split through the different VOC species 365 ! > Check for VOC Species: In input in this case we only have one species (VOC) 366 IF (TRIM(emt_att%species_name(ind_inp))=="VOC") THEN 367 !> Cycle over the input voc species: they have to be one of the VOCs of the mechanism used. A list of VOC species for different mechanisms is provided in the module documentation 368 DO ind_voc= 1,emt_att%nvoc 369 !> Determine the indices of the coinciding species in the two cases and assign them to matching arrays 370 IF (TRIM(emt_att%voc_name(ind_voc))==TRIM(spc_names(ind_mod))) THEN 371 len_index=len_index+1 372 len_index_voc=len_index_voc+1 262 ! 263 !-- Check for VOC Species 264 IF ( TRIM( emt_att%species_name(ind_inp) ) == "VOC" ) THEN 265 DO ind_voc = 1, emt_att%nvoc 266 267 IF ( TRIM( emt_att%voc_name(ind_voc) ) == TRIM( spc_names(ind_mod) ) ) THEN 268 len_index = len_index + 1 269 len_index_voc = len_index_voc + 1 373 270 ENDIF 374 271 END DO 375 272 ENDIF 376 !> Other Species 377 IF (TRIM(emt_att%species_name(ind_inp))==TRIM(spc_names(ind_mod))) THEN 378 len_index=len_index+1 273 ! 274 !-- Other Species 275 IF ( TRIM( emt_att%species_name(ind_inp) ) == TRIM( spc_names(ind_mod) ) ) THEN 276 len_index = len_index + 1 379 277 ENDIF 380 278 ENDDO 381 279 ENDDO 382 280 383 !> Allocate array for storing the indices of the matched species 384 IF (len_index>0) THEN 385 386 ALLOCATE(match_spec_input(len_index)) 387 388 ALLOCATE(match_spec_model(len_index)) 389 390 IF (len_index_voc>0) THEN 391 392 ALLOCATE(match_spec_voc_model(len_index_voc)) !> Contains indices of the VOC model species 393 394 ALLOCATE(match_spec_voc_input(len_index_voc)) !> In input there is only one value for VOCs in the DEFAULT mode. This array contains the indices of the different values of VOC compositions of the input variable VOC_composition 281 ! 282 !-- Allocate array for storing the indices of the matched species 283 IF ( len_index > 0 ) THEN 284 285 ALLOCATE( match_spec_input(len_index) ) 286 287 ALLOCATE( match_spec_model(len_index) ) 288 289 IF ( len_index_voc > 0 ) THEN 290 ! 291 !-- contains indices of the VOC model species 292 ALLOCATE( match_spec_voc_model(len_index_voc) ) 293 ! 294 !-- contains the indices of different values of VOC composition of input variable VOC_composition 295 ALLOCATE( match_spec_voc_input(len_index_voc) ) 395 296 396 297 ENDIF 397 298 398 !> Repeat the same cycle of above but passing the species indices to the newly declared arrays 399 len_index=0 400 401 !> Cycle over model species 402 DO ind_mod = 1, SIZE(spc_names) 403 !> Cycle over Input species 299 ! 300 !-- pass the species indices to declared arrays 301 len_index = 0 302 303 ! 304 !-- Cycle over model species 305 DO ind_mod = 1, nvar 306 ! 307 !-- Cycle over Input species 404 308 DO ind_inp = 1, nspec_emis_inp 405 406 !> Determine the indices of the coinciding species in the two cases and assign them to matching arrays407 408 ! > VOCs409 IF ( TRIM(emt_att%species_name(ind_inp))=="VOC" .AND. ALLOCATED(match_spec_voc_input) ) THEN410 DO ind_voc= 1, emt_att%nvoc411 IF ( TRIM(emt_att%voc_name(ind_voc))==TRIM(spc_names(ind_mod)))THEN412 len_index =len_index+1413 len_index_voc =len_index_voc+1309 ! 310 !-- VOCs 311 IF ( TRIM(emt_att%species_name(ind_inp) ) == "VOC" .AND. & 312 ALLOCATED(match_spec_voc_input) ) THEN 313 314 DO ind_voc= 1, emt_att%nvoc 315 IF ( TRIM( emt_att%voc_name(ind_voc) ) == TRIM( spc_names(ind_mod) ) ) THEN 316 len_index = len_index + 1 317 len_index_voc = len_index_voc + 1 414 318 415 match_spec_input(len_index) =ind_inp416 match_spec_model(len_index) =ind_mod417 418 match_spec_voc_input(len_index_voc) =ind_voc419 match_spec_voc_model(len_index_voc) =ind_mod319 match_spec_input(len_index) = ind_inp 320 match_spec_model(len_index) = ind_mod 321 322 match_spec_voc_input(len_index_voc) = ind_voc 323 match_spec_voc_model(len_index_voc) = ind_mod 420 324 ENDIF 421 325 END DO 422 326 ENDIF 423 327 424 !>Other Species 425 IF (TRIM(emt_att%species_name(ind_inp))==TRIM(spc_names(ind_mod))) THEN 426 len_index=len_index+1 427 match_spec_input(len_index)=ind_inp 428 match_spec_model(len_index)=ind_mod 328 ! 329 !-- Other Species 330 IF ( TRIM( emt_att%species_name(ind_inp) ) == TRIM( spc_names(ind_mod) ) ) THEN 331 len_index = len_index + 1 332 match_spec_input(len_index) = ind_inp 333 match_spec_model(len_index) = ind_mod 429 334 ENDIF 430 335 END DO 431 336 END DO 432 337 433 !> In the case there are no species matching, the emission module should not be called434 338 ELSE 435 436 message_string = 'Given Chemistry Emission Species' // & 437 ' DO NOT MATCH' // & 438 ' model chemical species:' // & 439 ' Chemistry Emissions routine is not called' 339 ! 340 !-- in case there are no species matching 341 message_string = 'Non of given emission species' // & 342 ' matches' // & 343 ' model chemical species:' // & 344 ' Emission routine is not called' 440 345 CALL message( 'chem_emissions_matching', 'CM0438', 0, 0, 0, 6, 0 ) 441 ENDIF !> IF (len_index>0)346 ENDIF 442 347 443 348 ELSE 444 349 445 ! In this case either spc_names is null or nspec_emis_inp is not allocated446 350 ! 351 !-- either spc_names is zero or nspec_emis_inp is not allocated 447 352 message_string = 'Array of Emission species not allocated:' // & 448 ' Either no emission species are provided as input or' // 449 ' no chemical species are used by PALM:' // 450 ' Chemistry Emissionsroutine is not called'353 ' Either no emission species are provided as input or' // & 354 ' no chemical species are used by PALM:' // & 355 ' Emission routine is not called' 451 356 CALL message( 'chem_emissions_matching', 'CM0439', 0, 2, 0, 6, 0 ) 452 357 453 ENDIF !> IF ( (SIZE(spc_names) .GT. 0) .AND. ALLOCATED(nspec_emis_inp)) 454 455 !> ------------------------------------------------------------------ 456 !> DEFAULT case 457 358 ENDIF 359 360 ! 361 !-- DEFAULT mode 458 362 CASE ("DEFAULT") 459 363 460 CALL location_message( 'chem_emissions DEFAULT_mode_matching', .FALSE. ) 461 462 len_index=0 !> index for TOTAL number of species 463 len_index_voc=0 !> index for TOTAL number of VOCs 464 len_index_pm=3 !> index for TOTAL number of PM Types: PM1, PM2.5, PM10. It is needed because the number of emission inputs and the one of PMs in the model may not be the same. 465 466 !> The first condition is that both the number of model and input emissions species are not null 467 IF ( (SIZE(spc_names) .GT. 0) .AND. (nspec_emis_inp .GT. 0) ) THEN 468 469 !> Cycle over model species 470 DO ind_mod = 1, SIZE(spc_names) 471 !> Cycle over input species 364 len_index = 0 !< index for TOTAL number of species 365 len_index_voc = 0 !< index for TOTAL number of VOCs 366 len_index_pm = 3 !< index for TOTAL number of PMs: PM1, PM2.5, PM10. 367 368 IF ( nvar > 0 .AND. nspec_emis_inp > 0 ) THEN 369 370 ! 371 !-- Cycle over model species 372 DO ind_mod = 1, nvar 373 ! 374 !-- Cycle over input species 472 375 DO ind_inp = 1, nspec_emis_inp 473 376 474 ! > Check for VOC Species: In input in this case we only have one species (VOC)475 IF (TRIM(emt_att%species_name(ind_inp))=="VOC") THEN476 !> Cycle over the voc species given in input: they have to be one of the VOCs of the mechanism used. A list of VOC species for different mechanisms is provided in the module documentation477 DO ind_voc= 1, emt_att%nvoc478 !> Determine the indices of the coinciding species in the two cases and assign them to matching arrays479 IF ( TRIM(emt_att%voc_name(ind_voc))==TRIM(spc_names(ind_mod)))THEN480 len_index =len_index+1481 len_index_voc =len_index_voc+1377 ! 378 !-- Check for VOC Species 379 IF ( TRIM( emt_att%species_name(ind_inp) ) == "VOC" ) THEN 380 DO ind_voc= 1, emt_att%nvoc 381 382 IF ( TRIM( emt_att%voc_name(ind_voc) ) == TRIM( spc_names(ind_mod) ) ) THEN 383 len_index = len_index + 1 384 len_index_voc = len_index_voc + 1 482 385 ENDIF 386 483 387 END DO 484 388 ENDIF 485 389 486 !> PMs: For PMs there is only one input species name for all the PM types. This variable has 3 dimensions, one for each PM type. 487 IF (TRIM(emt_att%species_name(ind_inp))=="PM") THEN 488 !>pm1 489 IF (TRIM(spc_names(ind_mod))=="PM1") THEN 490 len_index=len_index+1 491 !>pm2.5 492 ELSE IF (TRIM(spc_names(ind_mod))=="PM25") THEN 493 len_index=len_index+1 494 !>pm10 495 ELSE IF (TRIM(spc_names(ind_mod))=="PM10") THEN 496 len_index=len_index+1 390 ! 391 !-- PMs: There is one input species name for all PM 392 !-- This variable has 3 dimensions, one for PM1, PM2.5 and PM10 393 IF ( TRIM( emt_att%species_name(ind_inp) ) == "PM" ) THEN 394 ! 395 !-- PM1 396 IF ( TRIM( spc_names(ind_mod) ) == "PM1" ) THEN 397 len_index = len_index + 1 398 ! 399 !-- PM2.5 400 ELSEIF ( TRIM( spc_names(ind_mod) ) == "PM25" ) THEN 401 len_index = len_index + 1 402 ! 403 !-- PM10 404 ELSEIF ( TRIM( spc_names(ind_mod) ) == "PM10" ) THEN 405 len_index = len_index + 1 497 406 ENDIF 498 407 ENDIF 499 408 500 !> NOx: for NOx by definition we have 2 splits. The Emission Input Name in this case is only one: NOx, while in the model we can have NO2 and NO 501 IF (TRIM(emt_att%species_name(ind_inp))=="NOx") THEN 502 !>no 503 IF (TRIM(spc_names(ind_mod))=="NO") THEN 504 len_index=len_index+1 505 !>no2 506 ELSE IF (TRIM(spc_names(ind_mod))=="NO2") THEN 507 len_index=len_index+1 409 ! 410 !-- NOx: NO2 and NO 411 IF ( TRIM( emt_att%species_name(ind_inp) ) == "NOx" ) THEN 412 ! 413 !-- NO 414 IF ( TRIM( spc_names(ind_mod) ) == "NO" ) THEN 415 len_index = len_index + 1 416 ! 417 !-- NO2 418 ELSEIF ( TRIM( spc_names(ind_mod) ) == "NO2" ) THEN 419 len_index = len_index + 1 508 420 ENDIF 509 421 ENDIF 510 422 511 !>SOX: same as for NOx, but with SO2 and SO4 512 IF (TRIM(emt_att%species_name(ind_inp))=="SOx") THEN 513 !>no 514 IF (TRIM(spc_names(ind_mod))=="SO2") THEN 515 len_index=len_index+1 516 !>no2 517 ELSE IF (TRIM(spc_names(ind_mod))=="SO4") THEN 518 len_index=len_index+1 423 ! 424 !-- SOx: SO2 and SO4 425 IF ( TRIM( emt_att%species_name(ind_inp) ) == "SOx" ) THEN 426 ! 427 !-- SO2 428 IF ( TRIM( spc_names(ind_mod) ) == "SO2" ) THEN 429 len_index = len_index + 1 430 ! 431 !-- SO4 432 ELSEIF ( TRIM( spc_names(ind_mod) ) == "SO4" ) THEN 433 len_index = len_index + 1 519 434 ENDIF 520 435 ENDIF 521 436 522 !> Other Species 523 IF (TRIM(emt_att%species_name(ind_inp))==TRIM(spc_names(ind_mod))) THEN 524 len_index=len_index+1 437 ! 438 !-- Other Species 439 IF ( TRIM( emt_att%species_name(ind_inp) ) == TRIM( spc_names(ind_mod) ) ) THEN 440 len_index = len_index + 1 525 441 ENDIF 526 END DO !> number of emission input species 527 END DO !> number of model species 528 529 530 !> Allocate Arrays 531 IF (len_index>0) THEN 532 533 ALLOCATE(match_spec_input(len_index)) 534 ALLOCATE(match_spec_model(len_index)) 535 536 IF (len_index_voc>0) THEN 537 ALLOCATE(match_spec_voc_model(len_index_voc)) !> contains indices of the VOC model species 538 ALLOCATE(match_spec_voc_input(len_index_voc)) !> In input there is only one value for VOCs in the DEFAULT mode. 539 ! This array contains the indices of the different values of VOC compositions of the input variable VOC_composition 442 END DO 443 END DO 444 445 446 ! 447 !-- Allocate arrays 448 IF ( len_index > 0 ) THEN 449 450 ALLOCATE( match_spec_input(len_index) ) 451 ALLOCATE( match_spec_model(len_index) ) 452 453 IF ( len_index_voc > 0 ) THEN 454 ! 455 !-- Contains indices of the VOC model species 456 ALLOCATE( match_spec_voc_model(len_index_voc) ) 457 ! 458 !-- Contains the indices of different values of VOC composition 459 !-- of input variable VOC_composition 460 ALLOCATE( match_spec_voc_input(len_index_voc) ) 540 461 ENDIF 541 462 542 ! > ASSIGN INDICES543 ! > Repeat the same cycles of above, but passing the species indices to the newlydeclared arrays544 len_index =0545 len_index_voc =0463 ! 464 !-- Pass the species indices to declared arrays 465 len_index = 0 466 len_index_voc = 0 546 467 547 DO ind_mod = 1, SIZE(spc_names)468 DO ind_mod = 1, nvar 548 469 DO ind_inp = 1, nspec_emis_inp 549 470 550 ! > VOCs 551 IF ( TRIM(emt_att%species_name(ind_inp))=="VOC" .AND. ALLOCATED(match_spec_voc_input) ) THEN 552 DO ind_voc= 1,emt_att%nvoc 553 IF (TRIM(emt_att%voc_name(ind_voc))==TRIM(spc_names(ind_mod))) THEN 554 len_index=len_index+1 555 len_index_voc=len_index_voc+1 471 ! 472 !-- VOCs 473 IF ( TRIM( emt_att%species_name(ind_inp) ) == "VOC" .AND. & 474 ALLOCATED(match_spec_voc_input) ) THEN 475 DO ind_voc= 1, emt_att%nvoc 476 IF ( TRIM( emt_att%voc_name(ind_voc) ) == TRIM( spc_names(ind_mod) ) ) THEN 477 len_index = len_index + 1 478 len_index_voc = len_index_voc + 1 556 479 557 match_spec_input(len_index) =ind_inp558 match_spec_model(len_index) =ind_mod559 560 match_spec_voc_input(len_index_voc) =ind_voc561 match_spec_voc_model(len_index_voc) =ind_mod480 match_spec_input(len_index) = ind_inp 481 match_spec_model(len_index) = ind_mod 482 483 match_spec_voc_input(len_index_voc) = ind_voc 484 match_spec_voc_model(len_index_voc) = ind_mod 562 485 ENDIF 563 486 END DO 564 487 ENDIF 565 488 566 !> PMs:In this case the Inputs have only PM while the model has three different possible types: PM1, PM2.5, PM10. We need an additional index for matching each PM index in the model. 567 IF (TRIM(emt_att%species_name(ind_inp))=="PM") THEN 568 !>pm1 569 IF (TRIM(spc_names(ind_mod))=="PM1") THEN 570 len_index=len_index+1 571 572 match_spec_input(len_index)=ind_inp 573 match_spec_model(len_index)=ind_mod 489 ! 490 !-- PMs 491 IF ( TRIM( emt_att%species_name(ind_inp) ) == "PM" ) THEN 492 ! 493 !-- PM1 494 IF ( TRIM( spc_names(ind_mod) ) == "PM1" ) THEN 495 len_index = len_index + 1 496 497 match_spec_input(len_index) = ind_inp 498 match_spec_model(len_index) = ind_mod 499 ! 500 !-- PM2.5 501 ELSEIF ( TRIM( spc_names(ind_mod) ) == "PM25" ) THEN 502 len_index = len_index + 1 503 504 match_spec_input(len_index) = ind_inp 505 match_spec_model(len_index) = ind_mod 506 ! 507 !-- PM10 508 ELSEIF ( TRIM( spc_names(ind_mod) ) == "PM10" ) THEN 509 len_index = len_index + 1 510 511 match_spec_input(len_index) = ind_inp 512 match_spec_model(len_index) = ind_mod 574 513 575 !match_spec_pm(1)=ind_mod576 !>pm2.5577 ELSE IF (TRIM(spc_names(ind_mod))=="PM25") THEN578 len_index=len_index+1579 580 match_spec_input(len_index)=ind_inp581 match_spec_model(len_index)=ind_mod582 583 !match_spec_pm(2)=ind_mod584 !>pm10585 ELSE IF (TRIM(spc_names(ind_mod))=="PM10") THEN586 len_index=len_index+1587 588 match_spec_input(len_index)=ind_inp589 match_spec_model(len_index)=ind_mod590 591 !match_spec_pm(3)=ind_mod592 514 ENDIF 593 515 ENDIF 594 516 595 !> NOx - The same as for PMs but here the model species are only 2: NO and NO2 596 IF (TRIM(emt_att%species_name(ind_inp))=="NOx") THEN 597 !>no 598 IF (TRIM(spc_names(ind_mod))=="NO") THEN 599 len_index=len_index+1 600 601 match_spec_input(len_index)=ind_inp 602 match_spec_model(len_index)=ind_mod 517 ! 518 !-- NOx 519 IF ( TRIM( emt_att%species_name(ind_inp) ) == "NOx" ) THEN 520 ! 521 !-- NO 522 IF ( TRIM( spc_names(ind_mod) ) == "NO" ) THEN 523 len_index = len_index + 1 524 525 match_spec_input(len_index) = ind_inp 526 match_spec_model(len_index) = ind_mod 527 ! 528 !-- NO2 529 ELSEIF ( TRIM( spc_names(ind_mod) ) == "NO2" ) THEN 530 len_index = len_index + 1 531 532 match_spec_input(len_index) = ind_inp 533 match_spec_model(len_index) = ind_mod 603 534 604 !match_spec_nox(1)=ind_mod605 !>no2606 ELSE IF (TRIM(spc_names(ind_mod))=="NO2") THEN607 len_index=len_index+1608 609 match_spec_input(len_index)=ind_inp610 match_spec_model(len_index)=ind_mod611 612 ! match_spec_nox(2)=ind_mod613 535 ENDIF 614 536 ENDIF 615 537 616 !> SOx 617 IF (TRIM(emt_att%species_name(ind_inp))=="SOx") THEN 618 !>so2 619 IF (TRIM(spc_names(ind_mod))=="SO2") THEN 620 len_index=len_index+1 621 622 match_spec_input(len_index)=ind_inp 623 match_spec_model(len_index)=ind_mod 538 ! 539 !-- SOx 540 IF ( TRIM( emt_att%species_name(ind_inp) ) == "SOx" ) THEN 541 ! 542 !-- SO2 543 IF ( TRIM( spc_names(ind_mod) ) == "SO2" ) THEN 544 len_index = len_index + 1 545 546 match_spec_input(len_index) = ind_inp 547 match_spec_model(len_index) = ind_mod 624 548 625 ! match_spec_sox(1)=ind_mod626 ! >so4627 ELSE IF (TRIM(spc_names(ind_mod))=="SO4")THEN628 len_index =len_index+1629 630 match_spec_input(len_index) =ind_inp631 match_spec_model(len_index) =ind_mod549 ! 550 !-- SO4 551 ELSEIF ( TRIM( spc_names(ind_mod) ) == "SO4" ) THEN 552 len_index = len_index + 1 553 554 match_spec_input(len_index) = ind_inp 555 match_spec_model(len_index) = ind_mod 632 556 633 ! match_spec_sox(2)=ind_mod634 557 ENDIF 635 558 ENDIF 636 559 637 !> Other Species 638 IF (TRIM(emt_att%species_name(ind_inp))==TRIM(spc_names(ind_mod))) THEN 639 len_index=len_index+1 560 ! 561 !-- Other Species 562 IF ( TRIM( emt_att%species_name(ind_inp) ) == TRIM( spc_names(ind_mod) ) ) THEN 563 len_index = len_index + 1 640 564 641 match_spec_input(len_index) =ind_inp642 match_spec_model(len_index) =ind_mod565 match_spec_input(len_index) = ind_inp 566 match_spec_model(len_index) = ind_mod 643 567 ENDIF 644 568 END DO … … 647 571 ELSE 648 572 649 message_string = ' Given Chemistry Emission Species' //&650 ' DO NOT MATCH' // &651 ' model chemical species' 652 ' Chemistry Emissionsroutine is not called'573 message_string = 'Non of given Emission Species' // & 574 ' matches' // & 575 ' model chemical species' // & 576 ' Emission routine is not called' 653 577 CALL message( 'chem_emissions_matching', 'CM0440', 0, 0, 0, 6, 0 ) 654 578 … … 658 582 659 583 message_string = 'Array of Emission species not allocated: ' // & 660 ' Either no emission species are provided as input or' // 661 ' no chemical species are used by PALM:' // 662 ' Chemistry Emissionsroutine is not called'584 ' Either no emission species are provided as input or' // & 585 ' no chemical species are used by PALM:' // & 586 ' Emission routine is not called' 663 587 CALL message( 'chem_emissions_matching', 'CM0441', 0, 2, 0, 6, 0 ) 664 588 665 589 ENDIF 666 667 !-- CASE parameterized: In the parameterized case the user gives in input the exact species names of the chemical mechanism: no split is required for NOx, SOx, PMs and VOCs, but units of emissions inputs must be in (mole/m**)/s, made exception for PMs. 668 590 591 ! 592 !-- PARAMETERIZED mode 669 593 CASE ("PARAMETERIZED") 670 594 671 len_index=0 672 673 !spc_names have to be greater than zero for proceeding 674 IF ( (SIZE(spc_names) .GT. 0) .AND. (nspec_emis_inp .GT. 0) ) THEN 675 676 677 !cycle over model species 678 DO ind_mod = 1, SIZE(spc_names) 679 ind_inp=1 680 DO WHILE ( TRIM( surface_csflux_name( ind_inp ) ) /= 'novalue' ) !<'novalue' is the default 681 IF (TRIM(surface_csflux_name( ind_inp ))==TRIM(spc_names(ind_mod))) THEN 682 len_index=len_index+1 595 len_index = 0 596 597 IF ( nvar > 0 .AND. nspec_emis_inp > 0 ) THEN 598 599 ! 600 !-- Cycle over model species 601 DO ind_mod = 1, nvar 602 ind_inp = 1 603 DO WHILE ( TRIM( surface_csflux_name(ind_inp) ) /= 'novalue' ) !< 'novalue' is the default 604 IF ( TRIM( surface_csflux_name(ind_inp) ) == TRIM( spc_names(ind_mod) ) ) THEN 605 len_index = len_index + 1 683 606 ENDIF 684 ind_inp=ind_inp+1 685 607 ind_inp = ind_inp + 1 686 608 ENDDO 687 609 ENDDO 688 610 689 !Proceed only if there are matched species 690 IF ( len_index .GT. 0 ) THEN 691 692 !Allocation of Arrays of the matched species 693 ALLOCATE(match_spec_input(len_index)) 694 695 ALLOCATE(match_spec_model(len_index)) 696 697 !Assign corresponding indices of input and model matched species 698 len_index=0 699 700 DO ind_mod = 1, SIZE(spc_names) 611 IF ( len_index > 0 ) THEN 612 613 ! 614 !-- Allocation of Arrays of the matched species 615 ALLOCATE( match_spec_input(len_index) ) 616 617 ALLOCATE( match_spec_model(len_index) ) 618 619 ! 620 !-- Pass the species indices to declared arrays 621 len_index = 0 622 623 DO ind_mod = 1, nvar 701 624 ind_inp = 1 702 DO WHILE ( TRIM( surface_csflux_name( ind_inp ) ) /= 'novalue' ) !<'novalue' is the default703 IF ( TRIM( surface_csflux_name( ind_inp ) ) == TRIM(spc_names(ind_mod)))THEN704 len_index =len_index+1705 match_spec_input(len_index) =ind_inp706 match_spec_model(len_index) =ind_mod625 DO WHILE ( TRIM( surface_csflux_name(ind_inp) ) /= 'novalue' ) 626 IF ( TRIM( surface_csflux_name(ind_inp) ) == TRIM( spc_names(ind_mod) ) ) THEN 627 len_index = len_index + 1 628 match_spec_input(len_index) = ind_inp 629 match_spec_model(len_index) = ind_mod 707 630 ENDIF 708 ind_inp =ind_inp+1631 ind_inp = ind_inp + 1 709 632 END DO 710 633 END DO 711 634 712 ! also, Add AN if to check that when the surface_csflux are passed to the namelist, also the street factor values are provided 713 DO ispec = 1 , len_index 714 715 IF ( emiss_factor_main(match_spec_input(ispec)) .LT. 0 .AND. & 716 emiss_factor_side(match_spec_input(ispec)) .LT. 0 ) THEN 717 718 message_string = 'PARAMETERIZED emissions mode selected:' // & 635 ! 636 !-- Check 637 DO ispec = 1, len_index 638 639 IF ( emiss_factor_main(match_spec_input(ispec) ) < 0 .AND. & 640 emiss_factor_side(match_spec_input(ispec) ) < 0 ) THEN 641 642 message_string = 'PARAMETERIZED emissions mode selected:' // & 719 643 ' EMISSIONS POSSIBLE ONLY ON STREET SURFACES' // & 720 644 ' but values of scaling factors for street types' // & … … 730 654 ELSE 731 655 732 message_string = ' Given ChemistryEmission Species' // &733 ' DO NOT MATCH' //&734 ' model chemical species' //&735 ' Chemistry Emissionsroutine is not called'656 message_string = 'Non of given Emission Species' // & 657 ' matches' // & 658 ' model chemical species' // & 659 ' Emission routine is not called' 736 660 CALL message( 'chem_emissions_matching', 'CM0443', 0, 0, 0, 6, 0 ) 737 661 ENDIF … … 739 663 ELSE 740 664 741 message_string = 'Array of Emission species not allocated: ' // &665 message_string = 'Array of Emission species not allocated: ' // & 742 666 ' Either no emission species are provided as input or' // & 743 667 ' no chemical species are used by PALM.' // & 744 ' Chemistry Emissionsroutine is not called'668 ' Emission routine is not called' 745 669 CALL message( 'chem_emissions_matching', 'CM0444', 0, 2, 0, 6, 0 ) 746 670 … … 748 672 749 673 750 !-- CASE when emission module is switched on but mode_emis is not specified or it is given the wrong name 674 ! 675 !-- If emission module is switched on but mode_emis is not specified or it is given the wrong name 751 676 CASE DEFAULT 752 677 753 message_string = ' Chemistry Emissions Module switched ON, but' //&678 message_string = 'Emission Module switched ON, but' // & 754 679 ' either no emission mode specified or incorrectly given :' // & 755 680 ' please, pass the correct value to the namelist parameter "mode_emis"' … … 758 683 END SELECT 759 684 760 END SUBROUTINE chem_emissions_match 761 685 END SUBROUTINE chem_emissions_match 686 687 762 688 !------------------------------------------------------------------------------! 763 689 ! Description: 764 690 ! ------------ 765 691 !> Initialization: 766 !> Netcdf reading, arrays allocation and first calculation of cssws fluxes at timestep 0767 !> 692 !> Netcdf reading, arrays allocation and first calculation of cssws 693 !> fluxes at timestep 0 768 694 !------------------------------------------------------------------------------! 769 SUBROUTINE chem_emissions_init( emt_att,emt,nspec_out)695 SUBROUTINE chem_emissions_init( emt_att, emt, nspec_out ) 770 696 771 697 USE surface_mod, & 772 ONLY: surf_lsm_h,surf_def_h,surf_usm_h698 ONLY: surf_def_h, surf_lsm_h, surf_usm_h 773 699 774 700 IMPLICIT NONE 775 701 776 TYPE(chem_emis_att_type),INTENT(INOUT) :: emt_att !> variable where to store all the information of 777 ! emission inputs whose values do not depend 778 ! on the considered species 779 780 TYPE(chem_emis_val_type),INTENT(INOUT), ALLOCATABLE, DIMENSION(:) :: emt !> variable where to store emission inputs values, 781 ! depending on the considered species 702 CHARACTER (LEN=80) :: units !< units of inputs 782 703 783 INTEGER(iwp),INTENT(INOUT) :: nspec_out !> number of outputs of chemistry emission routines 784 785 CHARACTER (LEN=80) :: units !> units of chemistry inputs 786 787 INTEGER(iwp) :: ispec !> Index to go through the emission chemical species 788 789 790 !-- Actions for initial runs : TBD: needs to be updated 791 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 704 INTEGER(iwp), INTENT(INOUT) :: nspec_out !< number of outputs 705 706 INTEGER(iwp) :: ispec !< running index 707 708 TYPE(chem_emis_att_type), INTENT(INOUT) :: emt_att !< variable where to store all the information 709 !< of emission inputs whose values do not 710 !< depend on the considered species 711 712 TYPE(chem_emis_val_type), INTENT(INOUT), ALLOCATABLE, DIMENSION(:) :: emt !< variable where to store emission input 713 !< values depending on the considered species 714 715 ! 716 !-- Actions for initial runs 717 ! IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 792 718 !-- ... 793 719 ! 794 720 ! 795 721 !-- Actions for restart runs 722 ! ELSE 723 !-- ... 724 ! 725 ! ENDIF 726 727 728 CALL location_message( 'Starting initialization of emission arrays', .FALSE. ) 729 730 ! 731 !-- Matching 732 CALL chem_emissions_match( emt_att, nspec_out ) 733 734 IF ( nspec_out == 0 ) THEN 735 736 emission_output_required = .FALSE. 737 796 738 ELSE 797 !-- ... 798 799 ENDIF 800 801 802 CALL location_message( 'Starting initialization of chemistry emissions arrays', .FALSE. ) 803 804 805 !-- Match Input and Model emissions 806 CALL chem_emissions_match(emt_att,nspec_out) 807 808 !> If nspec_out==0, then do not use emission module: The corresponding message is already printed in the matching routine 809 IF ( nspec_out == 0 ) THEN 810 811 emission_output_required=.FALSE. 812 813 ELSE 814 815 emission_output_required=.TRUE. 816 817 818 !----------------------------------------------------------------- 819 !> Set molecule masses' 820 ALLOCATE(emt_att%xm(nspec_out)) 821 ! loop over emisisons: 739 740 emission_output_required = .TRUE. 741 742 743 ! 744 !-- Set molecule masses' 745 ALLOCATE( emt_att%xm(nspec_out) ) 822 746 823 747 DO ispec = 1, nspec_out 824 ! switch: 825 SELECT CASE ( TRIM(spc_names(match_spec_model(ispec))) ) 826 CASE ( 'SO2','SO4' ) ; emt_att%xm(ispec) = xm_S + xm_O * 2 ! kg/mole 827 CASE ( 'NO','NO2' ) ; emt_att%xm(ispec) = xm_N + xm_O * 2 ! kg/mole NO2 equivalent 828 CASE ( 'NH3' ) ; emt_att%xm(ispec) = xm_N + xm_H * 3 ! kg/mole 829 CASE ( 'CO' ) ; emt_att%xm(ispec) = xm_C + xm_O ! kg/mole 830 CASE ( 'CO2' ) ; emt_att%xm(ispec) = xm_C + xm_O * 2 ! kg/mole 831 CASE ( 'CH4' ) ; emt_att%xm(ispec) = xm_C + xm_H * 4 ! kg/mole 832 CASE ( 'HNO3' ); emt_att%xm(ispec) = xm_H + xm_N + xm_O*3 !kg/mole 748 SELECT CASE ( TRIM( spc_names(match_spec_model(ispec)) ) ) 749 CASE ( 'SO2' ); emt_att%xm(ispec) = xm_S + xm_O * 2 !< kg/mole 750 CASE ( 'SO4' ); emt_att%xm(ispec) = xm_S + xm_O * 4 !< kg/mole 751 CASE ( 'NO' ); emt_att%xm(ispec) = xm_N + xm_O !< kg/mole 752 CASE ( 'NO2' ); emt_att%xm(ispec) = xm_N + xm_O * 2 !< kg/mole 753 CASE ( 'NH3' ); emt_att%xm(ispec) = xm_N + xm_H * 3 !< kg/mole 754 CASE ( 'CO' ); emt_att%xm(ispec) = xm_C + xm_O !< kg/mole 755 CASE ( 'CO2' ); emt_att%xm(ispec) = xm_C + xm_O * 2 !< kg/mole 756 CASE ( 'CH4' ); emt_att%xm(ispec) = xm_C + xm_H * 4 !< kg/mole 757 CASE ( 'HNO3' ); emt_att%xm(ispec) = xm_H + xm_N + xm_O*3 !< kg/mole 833 758 CASE DEFAULT 834 !-- TBD: check if this hase to be removed 835 !emt_att%xm(ispec) = 1.0_wp 759 emt_att%xm(ispec) = 1.0_wp 836 760 END SELECT 837 761 ENDDO 838 762 839 763 840 ! > ASSIGN EMISSION VALUES for the different emission modes841 SELECT CASE(TRIM(mode_emis)) !TBD: Add the option for CApital or small letters842 843 844 !> PRE-PROCESSED case 845 CASE ("PRE-PROCESSED")846 847 IF ( .NOT. ALLOCATED( emis_distribution) ) ALLOCATE(emis_distribution(nzb:nzt+1,0:ny,0:nx,nspec_out))848 849 CALL location_message( 'emis_distribution array allocated in PRE-PROCESSED mode', .FALSE. )850 851 ! > Calculate the values of the emissions at the first time step852 CALL chem_emissions_setup(emt_att,emt,nspec_out)853 854 !> Default case 855 CASE ("DEFAULT")856 857 IF ( .NOT. ALLOCATED( emis_distribution) ) ALLOCATE( emis_distribution( 1, 0:ny, 0:nx, nspec_out ) )858 859 CALL location_message( 'emis_distribution array allocated in DEFAULT mode', .FALSE. )860 861 ! > Calculate the values of the emissions at the first time step862 CALL chem_emissions_setup(emt_att,emt,nspec_out)863 864 !> PARAMETERIZED case 865 CASE ("PARAMETERIZED")866 867 CALL location_message( 'emis_distribution array allocated in PARAMETERIZED mode', .FALSE.)868 869 ! For now for PAR and DEF values only, first vertical level of emis_distribution is allocated, while for PRE-PROCESSED all.870 IF ( .NOT. ALLOCATED( emis_distribution) ) ALLOCATE(emis_distribution(1,0:ny,0:nx,nspec_out)) 871 872 ! > Calculate the values of theemissions at the first time step873 CALL chem_emissions_setup( emt_att,emt,nspec_out)764 ! 765 !-- assign emission values 766 SELECT CASE ( TRIM( mode_emis ) ) 767 768 769 ! 770 !-- PRE-PROCESSED case 771 CASE ( "PRE-PROCESSED" ) 772 773 IF ( .NOT. ALLOCATED( emis_distribution) ) ALLOCATE( emis_distribution(nzb:nzt+1,0:ny,0:nx,nspec_out) ) 774 775 ! 776 !-- Get emissions at the first time step 777 CALL chem_emissions_setup( emt_att, emt, nspec_out ) 778 779 ! 780 !-- Default case 781 CASE ( "DEFAULT" ) 782 783 IF ( .NOT. ALLOCATED( emis_distribution) ) ALLOCATE( emis_distribution(1,0:ny,0:nx,nspec_out) ) 784 785 ! 786 !-- Get emissions at the first time step 787 CALL chem_emissions_setup( emt_att, emt, nspec_out ) 788 789 ! 790 !-- PARAMETERIZED case 791 CASE ( "PARAMETERIZED" ) 792 793 IF ( .NOT. ALLOCATED( emis_distribution) ) ALLOCATE( emis_distribution(1,0:ny,0:nx,nspec_out) ) 794 795 ! 796 !-- Get emissions at the first time step 797 CALL chem_emissions_setup( emt_att, emt, nspec_out) 874 798 875 799 END SELECT … … 887 811 !-------------------------------------------------------------------------------! 888 812 889 SUBROUTINE chem_emissions_setup( emt_att,emt,nspec_out)813 SUBROUTINE chem_emissions_setup( emt_att, emt, nspec_out ) 890 814 891 815 USE arrays_3d, & 892 816 ONLY: dzw 893 USE grid_variables, 817 USE grid_variables, & 894 818 ONLY: dx, dy 895 USE indices, 896 ONLY: nnx, nny,nnz897 USE surface_mod, 898 ONLY: surf_ lsm_h,surf_def_h,surf_usm_h899 USE netcdf_data_input_mod, 819 USE indices, & 820 ONLY: nnx, nny, nnz 821 USE surface_mod, & 822 ONLY: surf_def_h, surf_lsm_h, surf_usm_h 823 USE netcdf_data_input_mod, & 900 824 ONLY: street_type_f 901 USE arrays_3d, 825 USE arrays_3d, & 902 826 ONLY: hyp, pt 903 827 828 904 829 IMPLICIT NONE 905 830 906 !--- IN/OUT 907 908 TYPE(chem_emis_att_type),INTENT(INOUT) :: emt_att !> variable where to store all the information of 909 ! emission inputs whose values do not depend 910 ! on the considered species 911 912 TYPE(chem_emis_val_type),INTENT(INOUT), ALLOCATABLE, DIMENSION(:) :: emt !> variable where to store emission inputs values, 913 ! depending on the considered species 914 915 INTEGER,INTENT(IN) :: nspec_out !> Output of routine chem_emis_matching with number 916 ! of matched species 917 !--- 918 919 REAL(wp),ALLOCATABLE,DIMENSION(:,:) :: delta_emis !> Term to add to the emission_distribution array 920 ! for each of the categories at each time step 921 ! in the default mode 922 923 CHARACTER(LEN=80) :: units !> Units of the emissions 924 925 INTEGER(iwp) :: icat !> Index for number of categories 926 INTEGER(iwp) :: ispec !> index for number of species 927 INTEGER(iwp) :: i_pm_comp !> index for number of PM components 928 INTEGER(iwp) :: ivoc !> Index for number of components of VOCs 929 INTEGER(iwp) :: time_index !> Index for time 831 832 TYPE(chem_emis_att_type), INTENT(INOUT) :: emt_att !< variable to store emission information 833 834 TYPE(chem_emis_val_type), INTENT(INOUT), ALLOCATABLE, DIMENSION(:) :: emt !< variable to store emission input values, 835 !< depending on the considered species 836 837 INTEGER,INTENT(IN) :: nspec_out !< Output of matching routine with number 838 !< of matched species 839 840 CHARACTER(LEN=80) :: units !< Units of the emission data 841 842 INTEGER(iwp) :: i !< running index for grid in x-direction 843 INTEGER(iwp) :: j !< running index for grid in y-direction 844 INTEGER(iwp) :: k !< running index for grid in z-direction 845 INTEGER(iwp) :: m !< running index for horizontal surfaces 930 846 931 REAL(wp),ALLOCATABLE, DIMENSION(:) :: time_factor !> factor for time scaling of emissions 932 REAL(wp),ALLOCATABLE, DIMENSION(:,:) :: emis 933 934 REAL(wp), DIMENSION(24) :: par_emis_time_factor !< time factors 935 ! for the parameterized mode: these are fixed for each hour 936 ! of a single day. 937 REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: conv_to_ratio !> factor used for converting input 938 ! to adimensional concentration ratio 939 940 ! ------------------------------------------ 941 ! variables for the conversion of indices i and j according to surface_mod 942 943 INTEGER(iwp) :: i !> running index for grid in x-direction 944 INTEGER(iwp) :: j !> running index for grid in y-direction 945 INTEGER(iwp) :: k !> running index for grid in z-direction 946 INTEGER(iwp) :: m !> running index for horizontal surfaces 947 948 REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: tmp_temp 949 950 ! --- const ------------------------------- 951 !-CONVERSION FACTORS: TIME 952 ! number of sec per hour = 3600 953 REAL, PARAMETER :: s_per_hour = 3600.0 ! (s)/(hour) 954 ! number of sec per day = 86400 955 REAL, PARAMETER :: s_per_day = 86400.0 ! (s)/(day) 956 ! number of hours in a year of 365 days: 957 REAL, PARAMETER :: hour_per_year = 8760.0 !> TBD: What about leap years? 958 ! number of hours in a day: 959 REAL, PARAMETER :: hour_per_day = 24.0 960 961 ! conversion from hours to seconds (s/hour) = 1/3600.0 ~ 0.2777778 962 REAL, PARAMETER :: hour_to_s = 1/s_per_hour ! (hour)/(s) 963 ! conversion from day to seconds (s/day) = 1/86400.0 ~ 1.157407e-05 964 REAL, PARAMETER :: day_to_s = 1/s_per_day ! (day)/(s) 965 ! conversion from year to sec (s/year) = 1/31536000.0 ~ 3.170979e-08 966 REAL, PARAMETER :: year_to_s = 1/(s_per_hour*hour_per_year) ! (year)/(s) 967 968 !-CONVERSION FACTORS: WEIGHT 969 ! Conversion from tons to kg (kg/tons) = 100.0/1 ~ 100 970 REAL, PARAMETER :: tons_to_kg = 100 ! (tons)/(kg) 971 ! Conversion from g to kg (kg/g) = 1/1000.0 ~ 0.001 972 REAL, PARAMETER :: g_to_kg = 0.001 ! (g)/(kg) 973 ! Conversion from g to kg (kg/g) = 1/1000.0 ~ 0.001 974 REAL, PARAMETER :: miug_to_kg = 0.000000001 ! (microng)/(kg) 975 976 !< CONVERSION FACTORS: fraction to ppm 977 REAL(wp), PARAMETER :: ratio2ppm = 1.0e06_wp 847 INTEGER(iwp) :: icat !< Index for number of categories 848 INTEGER(iwp) :: ispec !< index for number of species 849 INTEGER(iwp) :: i_pm_comp !< index for number of PM components 850 INTEGER(iwp) :: ivoc !< Index for number of VOCs 851 INTEGER(iwp) :: time_index !< Index for time 852 853 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: delta_emis 854 REAL(wp), ALLOCATABLE, DIMENSION(:) :: time_factor !< factor for time scaling of emissions 855 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: emis 856 857 REAL(wp), DIMENSION(24) :: par_emis_time_factor !< time factors for the parameterized mode: 858 !< fixed houlry profile for example day 859 REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: conv_to_ratio !< factor used for converting input 860 !< to concentration ratio 861 REAL(wp), DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) :: tmp_temp 862 863 ! 864 !-- CONVERSION FACTORS: TIME 865 REAL(wp), PARAMETER :: s_per_hour = 3600.0 !< number of sec per hour (s)/(hour) 866 REAL(wp), PARAMETER :: s_per_day = 86400.0 !< number of sec per day (s)/(day) 867 REAL(wp), PARAMETER :: hour_per_year = 8760.0 !< number of hours in a year of 365 days 868 REAL(wp), PARAMETER :: hour_per_day = 24.0 !< number of hours in a day 869 870 REAL(wp), PARAMETER :: hour_to_s = 1/s_per_hour !< conversion from hours to seconds (s/hour) ~ 0.2777778 871 REAL(wp), PARAMETER :: day_to_s = 1/s_per_day !< conversion from day to seconds (s/day) ~ 1.157407e-05 872 REAL(wp), PARAMETER :: year_to_s = 1/(s_per_hour*hour_per_year) !< conversion from year to sec (s/year) ~ 3.170979e-08 873 ! 874 !-- CONVERSION FACTORS: WEIGHT 875 REAL(wp), PARAMETER :: tons_to_kg = 100 !< Conversion from tons to kg (kg/tons) 876 REAL(wp), PARAMETER :: g_to_kg = 0.001 !< Conversion from g to kg (kg/g) 877 REAL(wp), PARAMETER :: miug_to_kg = 0.000000001 !< Conversion from g to kg (kg/g) 878 ! 879 !-- CONVERSION FACTORS: fraction to ppm 880 REAL(wp), PARAMETER :: ratio2ppm = 1.0e06 978 881 !------------------------------------------------------ 979 882 980 IF ( emission_output_required ) THEN981 982 !> Set emis_dt !TBD: this is the same as dt_chem. We should consider the fact that dt_emis should be the timestep of input emissions or better defined, the timestep at which the emission routines are called: for now one hour. It should be made changeable.983 883 IF ( emission_output_required ) THEN 884 885 ! 886 !-- Set emis_dt 984 887 IF ( call_chem_at_all_substeps ) THEN 985 888 … … 993 896 994 897 995 ! --- CHECK UNITS 996 !>----------------------------------------------------- 997 !> Conversion of the units to the ones employed in PALM. 998 !> Possible temporal units of input emissions data are: year, hour and second; 999 !> the mass could be expressed in: tons, kilograms or grams. 1000 !> IN the PARAMETERIZED mode no conversion is possible: in this case INPUTS have to have fixed units. 1001 !>----------------------------------------------------- 1002 1003 IF (TRIM(mode_emis)=="DEFAULT" .OR. TRIM(mode_emis)=="PRE-PROCESSED" ) THEN 1004 1005 SELECT CASE(TRIM(emt_att%units)) 1006 !> kilograms 1007 CASE ( 'kg/m2/s','KG/M2/S') 1008 CALL location_message( 'Units of the emissions are consistent with' // & 1009 ' the ones employed in the PALM-4U model ', .FALSE. ) 898 ! 899 !-- Conversion of units to the ones employed in PALM 900 !-- In PARAMETERIZED mode no conversion is performed: in this case input units are fixed 901 902 IF ( TRIM( mode_emis ) == "DEFAULT" .OR. TRIM( mode_emis ) == "PRE-PROCESSED" ) THEN 903 904 SELECT CASE ( TRIM( emt_att%units ) ) 905 ! 906 !-- kilograms 907 CASE ( 'kg/m2/s', 'KG/M2/S' ) 908 1010 909 con_factor=1 1011 910 1012 CASE ('kg/m2/hour','KG/M2/HOUR') 1013 CALL location_message( 'Units of emission inputs need conversion', .FALSE. ) 911 CASE ('kg/m2/hour', 'KG/M2/HOUR' ) 1014 912 1015 913 con_factor=hour_to_s 1016 914 1017 CASE ( 'kg/m2/day','KG/M2/DAY' ) 1018 CALL location_message( 'Units of emission inputs need conversion', .FALSE. ) 1019 915 CASE ( 'kg/m2/day', 'KG/M2/DAY' ) 916 1020 917 con_factor=day_to_s 1021 918 1022 CASE ( 'kg/m2/year','KG/M2/YEAR' ) 1023 CALL location_message( 'Units of emission inputs need conversion', .FALSE. ) 1024 919 CASE ( 'kg/m2/year', 'KG/M2/YEAR' ) 920 1025 921 con_factor=year_to_s 1026 922 1027 !> Tons1028 CASE ( 'ton/m2/s','TON/M2/S' )1029 CALL location_message( 'Units of emission inputs need conversion', .FALSE. )1030 923 ! 924 !-- Tons 925 CASE ( 'ton/m2/s', 'TON/M2/S' ) 926 1031 927 con_factor=tons_to_kg 1032 928 1033 CASE ( 'ton/m2/hour','TON/M2/HOUR' ) 1034 CALL location_message( 'Units of emission inputs need conversion', .FALSE. ) 1035 929 CASE ( 'ton/m2/hour', 'TON/M2/HOUR' ) 930 1036 931 con_factor=tons_to_kg*hour_to_s 1037 932 1038 CASE ( 'ton/m2/year','TON/M2/YEAR' ) 1039 CALL location_message( 'Units of emission inputs need conversion', .FALSE. ) 1040 933 CASE ( 'ton/m2/year', 'TON/M2/YEAR' ) 934 1041 935 con_factor=tons_to_kg*year_to_s 1042 936 1043 !> Grams1044 CASE ( 'g/m2/s','G/M2/S' )1045 CALL location_message( 'Units of emission inputs need conversion', .FALSE.)1046 937 ! 938 !-- Grams 939 CASE ( 'g/m2/s', 'G/M2/S' ) 940 1047 941 con_factor=g_to_kg 1048 942 1049 CASE ( 'g/m2/hour','G/M2/HOUR' ) 1050 CALL location_message( 'Units of emission inputs need conversion', .FALSE. ) 943 CASE ( 'g/m2/hour', 'G/M2/HOUR' ) 1051 944 1052 945 con_factor=g_to_kg*hour_to_s 1053 946 1054 CASE ( 'g/m2/year','G/M2/YEAR' ) 1055 CALL location_message( 'Units of emission inputs need conversion', .FALSE. ) 1056 947 CASE ( 'g/m2/year', 'G/M2/YEAR' ) 948 1057 949 con_factor=g_to_kg*year_to_s 1058 950 1059 !> Micro Grams1060 CASE ( 'micrograms/m2/s','MICROGRAMS/M2/S' )1061 CALL location_message( 'Units of emission inputs need conversion', .FALSE.)1062 951 ! 952 !-- Micrograms 953 CASE ( 'micrograms/m2/s', 'MICROGRAMS/M2/S' ) 954 1063 955 con_factor=miug_to_kg 1064 956 1065 CASE ( 'micrograms/m2/hour','MICROGRAMS/M2/HOUR' ) 1066 CALL location_message( 'Units of emission inputs need conversion', .FALSE. ) 1067 957 CASE ( 'micrograms/m2/hour', 'MICROGRAMS/M2/HOUR' ) 958 1068 959 con_factor=miug_to_kg*hour_to_s 1069 960 1070 CASE ( 'micrograms/m2/year','MICROGRAMS/M2/YEAR' ) 1071 CALL location_message( 'Units of emission inputs need conversion', .FALSE. ) 1072 961 CASE ( 'micrograms/m2/year', 'MICROGRAMS/M2/YEAR' ) 962 1073 963 con_factor=miug_to_kg*year_to_s 1074 964 1075 965 CASE DEFAULT 1076 message_string = 'The Units of the provided input chemistry emission species' // &966 message_string = 'The Units of the provided emission input' // & 1077 967 ' are not the ones required by PALM-4U: please check ' // & 1078 ' chemistryemission module documentation.'968 ' emission module documentation.' 1079 969 CALL message( 'chem_emissions_setup', 'CM0446', 2, 2, 0, 6, 0 ) 1080 970 … … 1082 972 1083 973 1084 !> PRE-PROCESSED and parameterized mode1085 ELSE1086 1087 message_string = 'No Units conversion required for units of chemistry emissions' // &1088 ' of the PARAMETERIZED mode: units have to be provided in' // &1089 ' micromole/m**2/day for GASES and' // &1090 ' kg/m**2/day for PMs'1091 CALL message( 'chem_emissions_setup', 'CM0447', 0, 0, 0, 6, 0 )1092 1093 974 ENDIF 1094 975 1095 !> Conversion factors (if the units are kg/m**2/s) we have to convert them to ppm/s 1096 DO i=nxl,nxr 1097 DO j=nys,nyn 1098 !> Derive Temperature from Potential Temperature 976 ! 977 !-- Conversion factor to convert kg/m**2/s to ppm/s 978 DO i = nxl, nxr 979 DO j = nys, nyn 980 ! 981 !-- Derive Temperature from Potential Temperature 1099 982 tmp_temp(nzb:nzt+1,j,i) = pt(nzb:nzt+1,j,i) * ( hyp(nzb:nzt+1) * pref_i )**r_cp 1100 1101 !> We need to pass to cssws<- (ppm/s) * dz.1102 ! Input is Nmole/(m^2*s)1103 ! To go to ppm*dz basically we need to multiply the input by (m**2/N)*dz1104 ! (m**2/N)*dz == V/N1105 ! V/N=RT/P1106 1107 !> m**3/Nmole(J/mol)*K^-1 K Pa983 984 985 !> We need to pass to cssws <- (ppm/s) * dz 986 !> Input is Nmole/(m^2*s) 987 !> To go to ppm*dz multiply the input by (m**2/N)*dz 988 !> (m**2/N)*dz == V/N 989 !> V/N = RT/P 990 !> m**3/Nmole (J/mol)*K^-1 K Pa 1108 991 conv_to_ratio(nzb:nzt+1,j,i) = ( (Rgas * tmp_temp(nzb:nzt+1,j,i)) / ((hyp(nzb:nzt+1))) ) 1109 992 ENDDO 1110 993 ENDDO 1111 994 1112 !>------------------------------------------------ 1113 !> Start The Processing of the input data 1114 1115 emis_distribution(:,nys:nyn,nxl:nxr,:) = 0.0_wp 1116 1117 !>----------------------------------------------- 1118 !> Distinguish between DEFAULT, PRE-PROCESSED and PARAMETERIZED mode when calculating time_factors: only done for DEFAULT mode. For the PARAMETERIZED mode there is a time factor but it is fixed in the model 995 996 ! 997 !-- Initialize 998 emis_distribution(:,nys:nyn,nxl:nxr,:) = 0.0_wp 999 1119 1000 1120 !> PRE-PROCESSED MODE1121 IF (TRIM(mode_emis)=="PRE-PROCESSED") THEN1122 1123 !> Update time indices 1124 CALL time_preprocessed_indices(index_hh)1125 1126 CALL location_message( 'PRE-PROCESSED MODE: No time-factor specification required', .FALSE.)1127 1128 ELSEIF ( TRIM(mode_emis)=="DEFAULT")THEN1129 1130 CALL location_message( 'DEFAULT MODE: time-factor specification required', .FALSE. )1131 1132 !> Allocate array where to store temporary emission values1133 IF(.NOT. ALLOCATED(emis)) ALLOCATE(emis(nys:nyn,nxl:nxr))1134 1135 !> Allocate time factor per emitted component category1136 ALLOCATE(time_factor(emt_att%ncat))1137 1138 !> Read-in HOURLY emission time factor1139 IF (TRIM(time_fac_type)=="HOUR") THEN 1140 1141 !>Update time indices1142 CALL time_default_indices( month_of_year,day_of_month,hour_of_day,index_hh)1143 1144 !>Check if the index is less or equal to the temporal dimension of HOURLY emission files1145 IF ( index_hh .LE. SIZE(emt_att%hourly_emis_time_factor(1,:)))THEN1146 1147 !>Read-in the correspondant time factor1148 time_factor(:) = emt_att%hourly_emis_time_factor(:,index_hh)1001 ! 1002 !-- PRE-PROCESSED MODE 1003 IF ( TRIM( mode_emis ) == "PRE-PROCESSED" ) THEN 1004 1005 ! 1006 !-- Update time indices 1007 CALL time_preprocessed_indices( index_hh ) 1008 1009 ELSEIF ( TRIM( mode_emis ) == "DEFAULT" ) THEN 1010 1011 ! 1012 !-- Allocate array where to store temporary emission values 1013 IF ( .NOT. ALLOCATED(emis) ) ALLOCATE( emis(nys:nyn,nxl:nxr) ) 1014 ! 1015 !-- Allocate time factor per category 1016 ALLOCATE( time_factor(emt_att%ncat) ) 1017 ! 1018 !-- Read-in hourly emission time factor 1019 IF ( TRIM( time_fac_type ) == "HOUR" ) THEN 1020 1021 ! 1022 !-- Update time indices 1023 CALL time_default_indices( month_of_year, day_of_month, hour_of_day, index_hh ) 1024 ! 1025 !-- Check if the index is less or equal to the temporal dimension of HOURLY emission files 1026 IF ( index_hh <= SIZE( emt_att%hourly_emis_time_factor(1,:) ) ) THEN 1027 ! 1028 !-- Read-in the correspondant time factor 1029 time_factor(:) = emt_att%hourly_emis_time_factor(:,index_hh) 1149 1030 1150 1031 ELSE 1151 1032 1152 message_string = 'The "HOUR" time-factors (DEFAULT mode) of the chemistry emission species' // &1033 message_string = 'The "HOUR" time-factors in the DEFAULT mode ' // & 1153 1034 ' are not provided for each hour of the total simulation time' 1154 1035 CALL message( 'chem_emissions_setup', 'CM0448', 2, 2, 0, 6, 0 ) 1155 1036 1156 1037 ENDIF 1157 1158 !> Read-in MDH emissions time factors 1159 ELSEIF (TRIM(time_fac_type)=="MDH") THEN 1160 1161 !> Update time indices 1162 CALL time_default_indices(daytype_mdh,month_of_year,day_of_month,hour_of_day,index_mm,index_dd,index_hh) 1163 1164 1165 !> Check if the index is less or equal to the temporal dimension of MDH emission files 1166 IF ((index_hh+index_dd+index_mm) .LE. SIZE(emt_att%mdh_emis_time_factor(1,:))) THEN 1167 1168 !> Read-in the correspondant time factor 1169 time_factor(:)=emt_att%mdh_emis_time_factor(:,index_mm)*emt_att%mdh_emis_time_factor(:,index_dd)* & 1038 ! 1039 !-- Read-in MDH emissions time factors 1040 ELSEIF ( TRIM( time_fac_type ) == "MDH" ) THEN 1041 1042 ! 1043 !-- Update time indices 1044 CALL time_default_indices( daytype_mdh, month_of_year, day_of_month, & 1045 hour_of_day, index_mm, index_dd,index_hh ) 1046 1047 ! 1048 !-- Check if the index is less or equal to the temporal dimension of MDH emission files 1049 IF ( ( index_hh + index_dd + index_mm) <= SIZE( emt_att%mdh_emis_time_factor(1,:) ) ) THEN 1050 ! 1051 !-- Read-in the correspondant time factor 1052 time_factor(:) = emt_att%mdh_emis_time_factor(:,index_mm) * emt_att%mdh_emis_time_factor(:,index_dd) * & 1170 1053 emt_att%mdh_emis_time_factor(:,index_hh) 1171 1054 1172 1055 ELSE 1173 1056 1174 message_string = 'The "MDH" time-factors (DEFAULT mode) of the chemistry emission species' //&1057 message_string = 'The "MDH" time-factors in the DEFAULT mode ' // & 1175 1058 ' are not provided for each hour/day/month of the total simulation time' 1176 1059 CALL message( 'chem_emissions_setup', 'CM0449', 2, 2, 0, 6, 0 ) … … 1179 1062 1180 1063 ELSE 1181 !> condition when someone used the DEFAULT mode but forgets to indicate the time-factor in the namelist 1182 1183 message_string = 'The time-factor (DEFAULT mode) of the chemistry emission species' // & 1184 ' is not provided in the NAMELIST' 1064 1065 message_string = 'In the DEFAULT mode the time factor' // & 1066 ' has to be defined in the NAMELIST' 1185 1067 CALL message( 'chem_emissions_setup', 'CM0450', 2, 2, 0, 6, 0 ) 1186 1068 1187 1069 ENDIF 1188 1070 1189 !> PARAMETERIZED MODE 1190 ELSEIF (TRIM(mode_emis)=="PARAMETERIZED") THEN 1191 CALL location_message( 'PARAMETERIZED MODE: time-factor specification is fixed: ' // & 1192 ' 24 values for every day of the year ', .FALSE. ) 1071 ! 1072 !-- PARAMETERIZED MODE 1073 ELSEIF ( TRIM( mode_emis ) == "PARAMETERIZED" ) THEN 1193 1074 1194 !Assign Constant Values of time factors, diurnal time profile for traffic sector: 1195 par_emis_time_factor( : ) = & 1196 (/ 0.009, 0.004, 0.004, 0.009, 0.029, 0.039, 0.056, 0.053, 0.051, 0.051, 0.052, 0.055, & 1075 ! 1076 !-- assign constant values of time factors, diurnal time profile for traffic sector 1077 par_emis_time_factor( : ) = & 1078 (/ 0.009, 0.004, 0.004, 0.009, 0.029, 0.039, 0.056, 0.053, 0.051, 0.051, 0.052, 0.055, & 1197 1079 0.059, 0.061, 0.064, 0.067, 0.069, 0.069, 0.049, 0.039, 0.039, 0.029, 0.024, 0.019 /) 1198 1080 1199 !> in this case allocate time factor each hour in a day1200 IF (.NOT. ALLOCATED(time_factor)) ALLOCATE(time_factor(1)) 1201 1202 ! >Pass the values of time factors in the parameterized mode to the time_factor variable. in this case index_hh==hour_of_day1203 index_hh =hour_of_day1081 IF ( .NOT. ALLOCATED( time_factor ) ) ALLOCATE( time_factor(1) ) 1082 1083 ! 1084 !-- Get time-factor for specific hour 1085 index_hh = hour_of_day 1204 1086 1205 1087 time_factor(1) = par_emis_time_factor(index_hh) … … 1207 1089 ENDIF 1208 1090 1209 !-------------------------------- 1210 !-- EMISSION DISTRIBUTION Calculation 1211 1212 !> PARAMETERIZED CASE 1213 IF ( mode_emis == "PARAMETERIZED" ) THEN 1214 1215 DO ispec=1,nspec_out 1216 1217 !> Values are still micromoles/(m**2*s). Units in this case are always micromoles/m**2*day (or kilograms/m**2*day for PMs) 1218 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=surface_csflux(match_spec_input(ispec))*time_factor(1)*hour_to_s 1091 1092 ! 1093 !-- Emission distribution calculation 1094 1095 ! 1096 !-- PARAMETERIZED case 1097 IF ( TRIM( mode_emis ) == "PARAMETERIZED" ) THEN 1098 1099 DO ispec = 1, nspec_out 1100 1101 ! 1102 !-- Units are micromoles/m**2*day (or kilograms/m**2*day for PMs) 1103 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = surface_csflux(match_spec_input(ispec)) * & 1104 time_factor(1) * hour_to_s 1219 1105 1220 1106 ENDDO 1221 1107 1222 !> PRE-PROCESSED CASE 1223 ELSEIF (TRIM(mode_emis)=="PRE-PROCESSED") THEN 1224 1225 !> Start Cycle over Species 1226 DO ispec=1,nspec_out !> nspec_out represents the number of species in common between 1227 ! the emission input data and the chemistry mechanism used 1108 ! 1109 !-- PRE-PROCESSED case 1110 ELSEIF ( TRIM( mode_emis ) == "PRE-PROCESSED" ) THEN 1111 1112 ! 1113 !-- Cycle over species: 1114 !-- nspec_out represents the number of species in common between the emission input data 1115 !-- and the chemistry mechanism used 1116 DO ispec=1,nspec_out 1228 1117 1229 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emt(match_spec_input(ispec))% &1230 preproc_emission_data(index_hh,1,nys+1:nyn+1,nxl+1:nxr+1)* &1231 1118 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emt(match_spec_input(ispec))% & 1119 preproc_emission_data(index_hh,1,nys+1:nyn+1,nxl+1:nxr+1) * & 1120 con_factor 1232 1121 1233 1122 ENDDO 1234 1235 !TBD: At the moment the default mode considers a single vertical level (the surface). So we need to change it accordingly or eventually include the variable vertical_profile to be used in the default mode i we want to consider additional levels. 1236 1237 ELSE IF (TRIM(mode_emis)=="DEFAULT") THEN 1238 1239 !> Allocate array for the emission value corresponding to a specific category and time factor 1240 ALLOCATE(delta_emis(nys:nyn,nxl:nxr)) 1241 1242 1243 !> Start Cycle over categories 1123 1124 ! 1125 !-- DEFAULT case 1126 ELSEIF ( TRIM( mode_emis ) == "DEFAULT" ) THEN 1127 1128 ! 1129 !-- Allocate array for the emission value corresponding to a specific category and time factor 1130 ALLOCATE( delta_emis(nys:nyn,nxl:nxr) ) 1131 1132 ! 1133 !-- Cycle over categories 1244 1134 DO icat = 1, emt_att%ncat 1245 1135 1246 !> Start Cycle over Species 1247 DO ispec=1,nspec_out !> nspec_out represents the number of species in common between 1248 ! the emission input data and the chemistry mechanism used 1136 ! 1137 !-- Cycle over Species 1138 !-- nspec_out represents the number of species in common between the emission input data 1139 !-- and the chemistry mechanism used 1140 DO ispec = 1, nspec_out 1249 1141 1250 1142 emis(nys:nyn,nxl:nxr) = emt(match_spec_input(ispec))%default_emission_data(icat,nys+1:nyn+1,nxl+1:nxr+1) 1251 1143 1252 !TBD: The consideration of dt_emis of the input data is still missing. Basically the emissions could be provided every 10, 30 minutes and not always at one hour. This should be eventually solved in the date_and_time mode routine. 1253 1254 !> the time factors are 24 for each day. When multiplied by a daily value, they allow to have an hourly value. Then to convert it to seconds, we still have to divide this value by 3600. 1255 ! So given any units, we convert them to seconds and finally multiply them by 24 ((value/sec)*(24*3600)=value/day ---- (value/day)*time_factor=value/hour ---(value/hour)/(3600)=value/sec ) 1256 ! ((value/sec)*(24*3600)*time_factor)/3600=24*(value/sec)*time_factor 1257 1258 !> NOX Compositions 1259 IF (TRIM(spc_names(match_spec_model(ispec)))=="NO") THEN 1260 !> Kg/m2*s kg/m2*s 1261 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* & 1262 emt_att%nox_comp(icat,1)*con_factor*hour_per_day 1263 1264 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr) 1265 1266 ELSE IF (TRIM(spc_names(match_spec_model(ispec)))=="NO2") THEN 1267 1268 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* & 1269 emt_att%nox_comp(icat,2)*con_factor*hour_per_day 1270 1271 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr) 1272 1273 !> SOX Compositions 1274 1275 ELSE IF (TRIM(spc_names(match_spec_model(ispec)))=="SO2") THEN 1276 !> Kg/m2*s kg/m2*s 1277 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* & 1278 emt_att%sox_comp(icat,1)*con_factor*hour_per_day 1279 1280 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr) 1281 1282 ELSE IF (TRIM(spc_names(match_spec_model(ispec)))=="SO4") THEN 1283 !> Kg/m2*s kg/m2*s 1284 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* & 1285 emt_att%sox_comp(icat,2)*con_factor*hour_per_day 1286 1287 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr) 1288 1289 1290 !> PMs should be in kg/m**2/s, so conversion factor is here still required 1291 !> PM1 Compositions 1292 ELSE IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1") THEN 1293 1294 !> Cycle over the different pm components for PM1 type 1295 DO i_pm_comp= 1,SIZE(emt_att%pm_comp(1,:,1)) 1296 1297 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* & 1298 emt_att%pm_comp(icat,i_pm_comp,1)*con_factor*hour_per_day 1299 1300 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+ & 1144 1145 ! 1146 !-- NOx 1147 IF ( TRIM( spc_names(match_spec_model(ispec)) ) == "NO" ) THEN 1148 1149 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) * time_factor(icat) * & !< kg/m2*s 1150 emt_att%nox_comp(icat,1) * con_factor * hour_per_day 1151 1152 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emis_distribution(1,nys:nyn,nxl:nxr,ispec) + & 1301 1153 delta_emis(nys:nyn,nxl:nxr) 1302 1154 1155 ELSEIF ( TRIM( spc_names(match_spec_model(ispec)) ) == "NO2" ) THEN 1156 1157 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) * time_factor(icat) * & !< kg/m2*s 1158 emt_att%nox_comp(icat,2) * con_factor * hour_per_day 1159 1160 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emis_distribution(1,nys:nyn,nxl:nxr,ispec) + & 1161 delta_emis(nys:nyn,nxl:nxr) 1162 1163 ! 1164 !-- SOx 1165 ELSEIF ( TRIM( spc_names(match_spec_model(ispec)) ) == "SO2" ) THEN 1166 1167 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) * time_factor(icat) * & !< kg/m2*s 1168 emt_att%sox_comp(icat,1) * con_factor * hour_per_day 1169 1170 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emis_distribution(1,nys:nyn,nxl:nxr,ispec) + & 1171 delta_emis(nys:nyn,nxl:nxr) 1172 1173 ELSEIF ( TRIM( spc_names(match_spec_model(ispec)) ) == "SO4" ) THEN 1174 1175 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) * time_factor(icat) * & !< kg/m2*s 1176 emt_att%sox_comp(icat,2) * con_factor * hour_per_day 1177 1178 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emis_distribution(1,nys:nyn,nxl:nxr,ispec) + & 1179 delta_emis(nys:nyn,nxl:nxr) 1180 1181 1182 ! 1183 !-- PMs 1184 !-- PM1 1185 ELSEIF ( TRIM( spc_names(match_spec_model(ispec)) ) == "PM1" ) THEN 1186 1187 ! 1188 !-- Cycle over PM1 components 1189 DO i_pm_comp = 1, SIZE( emt_att%pm_comp(1,:,1) ) 1190 1191 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) * time_factor(icat) * & !< kg/m2*s 1192 emt_att%pm_comp(icat,i_pm_comp,1) * con_factor * hour_per_day 1193 1194 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emis_distribution(1,nys:nyn,nxl:nxr,ispec) + & 1195 delta_emis(nys:nyn,nxl:nxr) 1303 1196 ENDDO 1304 1197 1305 !> PM2.5 Compositions 1306 ELSE IF (TRIM(spc_names(match_spec_model(ispec)))=="PM25") THEN 1307 1308 !> Cycle over the different pm components for PM2.5 type 1309 DO i_pm_comp= 1,SIZE(emt_att%pm_comp(1,:,2)) 1310 1311 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* & 1312 emt_att%pm_comp(icat,i_pm_comp,2)*con_factor*hour_per_day 1313 1314 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+ & 1315 delta_emis(nys:nyn,nxl:nxr) 1198 ! 1199 !-- PM2.5 1200 ELSEIF ( TRIM( spc_names(match_spec_model(ispec)) ) == "PM25" ) THEN 1201 1202 ! 1203 !-- Cycle over PM2.5 components 1204 DO i_pm_comp = 1, SIZE( emt_att%pm_comp(1,:,2) ) 1205 1206 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) * time_factor(icat) * & !< kg/m2*s 1207 emt_att%pm_comp(icat,i_pm_comp,2) * con_factor * hour_per_day 1208 1209 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emis_distribution(1,nys:nyn,nxl:nxr,ispec) + & 1210 delta_emis(nys:nyn,nxl:nxr) 1316 1211 1317 1212 ENDDO 1318 1213 1319 !> PM10 Compositions 1320 ELSE IF (TRIM(spc_names(match_spec_model(ispec)))=="PM10") THEN 1321 1322 !> Cycle over the different pm components for PM10 type 1323 DO i_pm_comp= 1,SIZE(emt_att%pm_comp(1,:,3)) 1324 1325 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* & 1326 emt_att%pm_comp(icat,i_pm_comp,3)*con_factor*hour_per_day 1214 ! 1215 !-- PM10 1216 ELSEIF ( TRIM( spc_names(match_spec_model(ispec)) ) == "PM10" ) THEN 1217 1218 ! 1219 !-- Cycle over PM10 components 1220 DO i_pm_comp = 1, SIZE( emt_att%pm_comp(1,:,3) ) 1221 1222 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) * time_factor(icat) * & !< kg/m2*s 1223 emt_att%pm_comp(icat,i_pm_comp,3) * con_factor * hour_per_day 1327 1224 1328 1225 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+ & … … 1331 1228 ENDDO 1332 1229 1333 !> VOCs Compositions: for VOCs, the input value is provided in kg/m**2*s but the composition is provided in mole/kg: a conversion factor for the input that could be eventually provided in, for example, tons/(m**2*s) is still required 1334 ELSE IF (SIZE(match_spec_voc_input) .GT. 0) THEN 1335 1336 !TBD: maybe we can avoid the cycle 1337 DO ivoc= 1,SIZE(match_spec_voc_input) 1338 1339 IF (TRIM(spc_names(match_spec_model(ispec)))==TRIM(emt_att%voc_name(ivoc))) THEN 1340 1341 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr)*time_factor(icat)* & 1342 emt_att%voc_comp(icat,match_spec_voc_input(ivoc))*con_factor*hour_per_day 1343 1344 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+ & 1345 delta_emis(nys:nyn,nxl:nxr) 1230 ! 1231 !-- VOCs 1232 ELSEIF ( SIZE( match_spec_voc_input ) > 0 ) THEN 1233 1234 DO ivoc = 1, SIZE( match_spec_voc_input ) 1235 1236 IF ( TRIM( spc_names(match_spec_model(ispec)) ) == TRIM( emt_att%voc_name(ivoc) ) ) THEN 1237 1238 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) * time_factor(icat) * & 1239 emt_att%voc_comp(icat,match_spec_voc_input(ivoc)) * & 1240 con_factor * hour_per_day 1241 1242 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emis_distribution(1,nys:nyn,nxl:nxr,ispec) + & 1243 delta_emis(nys:nyn,nxl:nxr) 1346 1244 1347 1245 ENDIF … … 1349 1247 ENDDO 1350 1248 1351 !> General case (other species) 1249 ! 1250 !-- any other species 1352 1251 ELSE 1353 1252 1354 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) *time_factor(icat)*con_factor*hour_per_day1355 1356 emis_distribution(1,nys:nyn,nxl:nxr,ispec)=emis_distribution(1,nys:nyn,nxl:nxr,ispec)+delta_emis(nys:nyn,nxl:nxr)1357 1358 ENDIF ! IF (spc_names==)1359 1360 !> for every species and category set emis to 0 so to avoid overwriting. TBD: discuss whether necessary1361 1253 delta_emis(nys:nyn,nxl:nxr) = emis(nys:nyn,nxl:nxr) * time_factor(icat) * & 1254 con_factor * hour_per_day 1255 1256 emis_distribution(1,nys:nyn,nxl:nxr,ispec) = emis_distribution(1,nys:nyn,nxl:nxr,ispec) + & 1257 delta_emis(nys:nyn,nxl:nxr) 1258 1259 ENDIF 1260 1362 1261 emis(:,:)= 0 1363 1262 1364 1263 ENDDO 1365 1366 !> for every category set delta_emis to 0 so to avoid overwriting. TBD: discuss whether necessary 1367 1368 delta_emis(:,:)=0 1369 1264 1265 delta_emis(:,:)=0 1266 1370 1267 ENDDO 1371 1268 1372 ENDIF !> mode_emis1373 1374 !------------------------------------------------------------------------------------------------------- 1375 ! --- Cycle to transform x,y coordinates to the one of surface_mod and to assign emission values to cssws1376 !-- -----------------------------------------------------------------------------------------------------1377 1269 ENDIF 1270 1271 1272 ! 1273 !-- Cycle to transform x,y coordinates to the one of surface_mod and to assign emission values to cssws 1274 ! 1378 1275 !-- PARAMETERIZED mode 1379 !> For the PARAMETERIZED mode units of inputs are always micromoles/(m**2*s). In this case we do not need the molar mass for conversion into ppms 1380 IF (TRIM(mode_emis)=="PARAMETERIZED") THEN 1276 ! 1277 !-- Units of inputs are micromoles/(m**2*s) 1278 IF ( TRIM( mode_emis ) == "PARAMETERIZED" ) THEN 1381 1279 1382 1280 IF ( street_type_f%from_file ) THEN 1383 1281 1384 !> Streets are lsm surfaces, hence, no usm surface treatment required 1385 IF (surf_lsm_h%ns .GT. 0 ) THEN 1282 ! 1283 !-- Streets are lsm surfaces, hence, no usm surface treatment required 1284 IF ( surf_lsm_h%ns > 0 ) THEN 1386 1285 DO m = 1, surf_lsm_h%ns 1387 1286 i = surf_lsm_h%i(m) … … 1389 1288 k = surf_lsm_h%k(m) 1390 1289 1391 1392 1290 IF ( street_type_f%var(j,i) >= main_street_id .AND. street_type_f%var(j,i) < max_street_id ) THEN 1393 1291 1394 !> Cycle over already matched species 1395 DO ispec=1,nspec_out 1396 1397 !> PMs are already in mass units: kilograms 1398 IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1" & 1399 .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25" & 1400 .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM10") THEN 1401 1402 ! kg/(m^2*s) *kg/m^3 1292 ! 1293 !-- Cycle over matched species 1294 DO ispec = 1, nspec_out 1295 1296 ! 1297 !-- PMs are already in kilograms 1298 IF ( TRIM( spc_names(match_spec_model(ispec)) ) == "PM1 " & 1299 .OR. TRIM( spc_names(match_spec_model(ispec)) ) == "PM25" & 1300 .OR. TRIM( spc_names(match_spec_model(ispec)) )=="PM10") THEN 1301 1302 ! 1303 !-- kg/(m^2*s) * kg/m^3 1403 1304 surf_lsm_h%cssws(match_spec_model(ispec),m) = emiss_factor_main(match_spec_input(ispec)) * & 1404 ! kg/(m^2*s) 1405 emis_distribution(1,j,i,ispec)* & 1406 ! kg/m^3 1407 rho_air(k) 1408 1305 emis_distribution(1,j,i,ispec) * & !< in kg/(m^2*s) 1306 rho_air(k) !< in kg/m^3 1307 1308 ! 1309 !-- Other Species 1310 !-- Inputs are micromoles 1409 1311 ELSE 1410 1312 1411 !> Other Species: inputs are micromoles: they have to be converted in moles 1412 ! ppm/s *m *kg/m^3 1413 surf_lsm_h%cssws(match_spec_model(ispec),m) = emiss_factor_main(match_spec_input(ispec))* & 1414 ! micromoles/(m^2*s) 1415 emis_distribution(1,j,i,ispec) * & 1416 ! m^3/Nmole 1417 conv_to_ratio(k,j,i)* & 1418 ! kg/m^3 1419 rho_air(k) 1420 1313 ! 1314 !-- ppm/s *m *kg/m^3 1315 surf_lsm_h%cssws(match_spec_model(ispec),m) = emiss_factor_main(match_spec_input(ispec)) * & 1316 emis_distribution(1,j,i,ispec) * & !< in micromoles/(m^2*s) 1317 conv_to_ratio(k,j,i) * & !< in m^3/Nmole 1318 rho_air(k) !< in kg/m^3 1421 1319 1422 1320 ENDIF … … 1427 1325 ELSEIF ( street_type_f%var(j,i) >= side_street_id .AND. street_type_f%var(j,i) < main_street_id ) THEN 1428 1326 1429 !> Cycle over already matched species 1430 DO ispec=1,nspec_out 1431 1432 !> PMs are already in mass units: micrograms 1433 IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1" & 1434 .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25" & 1435 .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM10") THEN 1436 1437 ! kg/(m^2*s) *kg/m^3 1438 surf_lsm_h%cssws(match_spec_model(ispec),m)= emiss_factor_side(match_spec_input(ispec)) * & 1439 ! kg/(m^2*s) 1440 emis_distribution(1,j,i,ispec)* & 1441 ! kg/m^3 1442 rho_air(k) 1327 ! 1328 !-- Cycle over matched species 1329 DO ispec = 1, nspec_out 1330 1331 ! 1332 !-- PMs are already in kilograms 1333 IF ( TRIM( spc_names(match_spec_model(ispec)) ) == "PM1" & 1334 .OR. TRIM( spc_names(match_spec_model(ispec)) ) == "PM25" & 1335 .OR. TRIM( spc_names(match_spec_model(ispec)) ) == "PM10" ) THEN 1336 1337 ! 1338 !-- kg/(m^2*s) * kg/m^3 1339 surf_lsm_h%cssws(match_spec_model(ispec),m) = emiss_factor_side(match_spec_input(ispec)) * & 1340 emis_distribution(1,j,i,ispec) * & !< in kg/(m^2*s) 1341 rho_air(k) !< in kg/m^3 1342 ! 1343 !-- Other species 1344 !-- Inputs are micromoles 1443 1345 ELSE 1444 1346 1445 1446 !>Other Species: inputs are micromoles 1447 ! ppm/s *m *kg/m^3 1448 surf_lsm_h%cssws(match_spec_model(ispec),m) = emiss_factor_side(match_spec_input(ispec)) * & 1449 ! micromoles/(m^2*s) 1450 emis_distribution(1,j,i,ispec) * & 1451 ! m^3/Nmole 1452 conv_to_ratio(k,j,i)* & 1453 ! kg/m^3 1454 rho_air(k) 1347 ! 1348 !-- ppm/s *m *kg/m^3 1349 surf_lsm_h%cssws(match_spec_model(ispec),m) = emiss_factor_side(match_spec_input(ispec)) * & 1350 emis_distribution(1,j,i,ispec) * & !< in micromoles/(m^2*s) 1351 conv_to_ratio(k,j,i) * & !< in m^3/Nmole 1352 rho_air(k) !< in kg/m^3 1455 1353 ENDIF 1456 1354 … … 1459 1357 ELSE 1460 1358 1461 !> If no street type is defined, then assign null emissions to all the species 1359 ! 1360 !-- If no street type is defined, then assign zero emission to all the species 1462 1361 surf_lsm_h%cssws(:,m) = 0.0_wp 1463 1362 … … 1465 1364 1466 1365 ENDDO 1366 1467 1367 ENDIF 1468 1368 1469 1369 ENDIF 1470 1370 1471 !> For both DEFAULT and PRE-PROCESSED 1371 ! 1372 !-- For both DEFAULT and PRE-PROCESSED mode 1472 1373 ELSE 1473 1374 1474 1375 1475 DO ispec=1,nspec_out 1476 !TBD: for the PRE-PROCESSED mode in the future, values at different heights should be included! 1477 !> Default surface type 1478 IF (surf_def_h(0)%ns .GT. 0) THEN 1376 DO ispec = 1, nspec_out 1377 1378 ! 1379 !-- Default surfaces 1380 IF ( surf_def_h(0)%ns > 0 ) THEN 1479 1381 1480 1382 DO m = 1, surf_def_h(0)%ns … … 1485 1387 IF ( emis_distribution(1,j,i,ispec) > 0.0_wp ) THEN 1486 1388 1487 1488 !> Distinguish between PMs (no needing conversion in ppms), 1489 ! VOC (already converted to moles/(m**2*s) using conv_factors: they do not need molar masses for their conversion to PPMs ) and 1490 ! other Species (still expressed in Kg/(m**2*s) at this point) 1491 1492 !> PMs 1493 IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1" .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25" & 1494 .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM10") THEN 1389 ! 1390 !-- PMs 1391 IF ( TRIM( spc_names(match_spec_model(ispec)) ) == "PM1" & 1392 .OR. TRIM( spc_names(match_spec_model(ispec)) ) == "PM25" & 1393 .OR. TRIM( spc_names(match_spec_model(ispec)) ) == "PM10" ) THEN 1495 1394 1496 ! kg/(m^2*s) *kg/m^3 kg/(m^2*s)1497 surf_def_h(0)%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec)* &1498 ! kg/m^31499 rho_air(nzb)1395 ! 1396 !-- kg/(m^2*s) *kg/m^3 kg/(m^2*s) 1397 surf_def_h(0)%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec)* & 1398 rho_air(nzb) !< in kg/m^3 1500 1399 1501 1400 1502 1401 ELSE 1503 1402 1504 !> VOCs 1505 IF ( len_index_voc .GT. 0 .AND. emt_att%species_name(match_spec_input(ispec))=="VOC" ) THEN 1506 ! ( ppm/s) * m * kg/m^3 mole/(m^2/s) 1507 surf_def_h(0)%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1508 ! m^3/mole ppm 1509 conv_to_ratio(nzb,j,i)*ratio2ppm * & 1510 ! kg/m^3 1511 rho_air(nzb) 1512 1513 1514 !> OTHER SPECIES 1403 ! 1404 !-- VOCs 1405 IF ( len_index_voc > 0 .AND. emt_att%species_name(match_spec_input(ispec)) == "VOC" ) THEN 1406 ! 1407 !-- (ppm/s) * m * kg/m^3 mole/(m^2/s) 1408 surf_def_h(0)%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1409 conv_to_ratio(nzb,j,i) * & !< in m^3/mole 1410 ratio2ppm * & !< in ppm 1411 rho_air(nzb) !< in kg/m^3 1412 1413 1414 ! 1415 !-- Other species 1515 1416 ELSE 1516 1417 1517 ! ( ppm/s )*m * kg/m^3 kg/(m^2/s) 1518 surf_def_h(0)%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1519 ! mole/Kg 1520 (1/emt_att%xm(ispec))* & 1521 ! m^3/mole ppm 1522 conv_to_ratio(nzb,j,i)*ratio2ppm* & 1523 ! kg/m^3 1524 rho_air(nzb) 1418 ! 1419 !-- (ppm/s) * m * kg/m^3 kg/(m^2/s) 1420 surf_def_h(0)%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1421 ( 1.0_wp / emt_att%xm(ispec) ) * & !< in mole/kg 1422 conv_to_ratio(nzb,j,i) * & !< in m^3/mole 1423 ratio2ppm * & !< in ppm 1424 rho_air(nzb) !< in kg/m^3 1525 1425 1526 1426 … … 1533 1433 ENDDO 1534 1434 1535 END 1435 ENDIF 1536 1436 1537 !-- Land Surface Mode surface type 1538 IF (surf_lsm_h%ns .GT. 0) THEN 1437 ! 1438 !-- LSM surfaces 1439 IF ( surf_lsm_h%ns > 0 ) THEN 1539 1440 1540 1441 DO m = 1, surf_lsm_h%ns … … 1546 1447 IF ( emis_distribution(1,j,i,ispec) > 0.0_wp ) THEN 1547 1448 1548 !> Distinguish between PMs (no needing conversion in ppms), 1549 ! VOC (already converted to moles/(m**2*s) using conv_factors: they do not need molar masses for their conversion to PPMs ) and 1550 ! other Species (still expressed in Kg/(m**2*s) at this point) 1551 1552 !> PMs 1553 IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1" .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25" & 1554 .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM10") THEN 1449 ! 1450 !-- PMs 1451 IF ( TRIM( spc_names(match_spec_model(ispec)) ) == "PM1" & 1452 .OR. TRIM( spc_names(match_spec_model(ispec)) ) == "PM25" & 1453 .OR. TRIM( spc_names(match_spec_model(ispec)) ) == "PM10" ) THEN 1555 1454 1556 ! kg/(m^2*s) * kg/m^3 kg/(m^2*s)1557 surf_lsm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * &1558 ! kg/m^31559 rho_air(k)1455 ! 1456 !-- kg/(m^2*s) * kg/m^3 kg/(m^2*s) 1457 surf_lsm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1458 rho_air(k) !< in kg/m^3 1560 1459 1561 1460 ELSE 1562 1461 1563 !> VOCs 1564 IF ( len_index_voc .GT. 0 .AND. emt_att%species_name(match_spec_input(ispec))=="VOC" ) THEN 1565 ! ( ppm/s) * m * kg/m^3 mole/(m^2/s) 1566 surf_lsm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1567 ! m^3/mole ppm 1568 conv_to_ratio(k,j,i)*ratio2ppm* & 1569 ! kg/m^3 1570 rho_air(k) 1571 1572 !> OTHER SPECIES 1462 ! 1463 !-- VOCs 1464 IF ( len_index_voc > 0 .AND. emt_att%species_name(match_spec_input(ispec)) == "VOC" ) THEN 1465 ! 1466 !-- (ppm/s) * m * kg/m^3 mole/(m^2/s) 1467 surf_lsm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1468 conv_to_ratio(k,j,i) * & !< in m^3/mole 1469 ratio2ppm * & !< in ppm 1470 rho_air(k) !< in kg/m^3 1471 1472 ! 1473 !-- Other species 1573 1474 ELSE 1574 1575 ! ( ppm/s) * m * kg/m^3 kg/(m^2/s) 1576 surf_lsm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1577 ! mole/Kg 1578 (1/emt_att%xm(ispec))* & 1579 ! m^3/mole ppm 1580 conv_to_ratio(k,j,i)*ratio2ppm* & 1581 ! kg/m^3 1582 rho_air(k) 1475 ! 1476 !-- (ppm/s) * m * kg/m^3 kg/(m^2/s) 1477 surf_lsm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1478 ( 1.0_wp / emt_att%xm(ispec) ) * & !< in mole/kg 1479 conv_to_ratio(k,j,i) * & !< in m^3/mole 1480 ratio2ppm * & !< in ppm 1481 rho_air(k) !< in kg/m^3 1583 1482 1584 1483 ENDIF … … 1590 1489 ENDDO 1591 1490 1592 END IF 1593 1594 !-- Urban Surface Mode surface type 1595 IF (surf_usm_h%ns .GT. 0) THEN 1491 ENDIF 1492 1493 ! 1494 !-- USM surfaces 1495 IF ( surf_usm_h%ns > 0 ) THEN 1596 1496 1597 1497 … … 1604 1504 IF ( emis_distribution(1,j,i,ispec) > 0.0_wp ) THEN 1605 1505 1606 !> PMs 1607 IF (TRIM(spc_names(match_spec_model(ispec)))=="PM1" .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM25" & 1608 .OR. TRIM(spc_names(match_spec_model(ispec)))=="PM10") THEN 1506 ! 1507 !-- PMs 1508 IF ( TRIM( spc_names(match_spec_model(ispec)) ) == "PM1" & 1509 .OR. TRIM( spc_names(match_spec_model(ispec)) ) == "PM25" & 1510 .OR. TRIM( spc_names(match_spec_model(ispec)) ) == "PM10" ) THEN 1609 1511 1610 ! kg/(m^2*s) *kg/m^3 kg/(m^2*s)1611 surf_usm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec)* &1612 ! kg/m^31613 rho_air(k)1512 ! 1513 !-- kg/(m^2*s) *kg/m^3 kg/(m^2*s) 1514 surf_usm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec)* & 1515 rho_air(k) !< in kg/m^3 1614 1516 1615 1517 1616 1518 ELSE 1617 1618 !> VOCs 1619 IF ( len_index_voc .GT. 0 .AND. emt_att%species_name(match_spec_input(ispec))=="VOC" ) THEN 1620 ! ( ppm/s) * m * kg/m^3 mole/(m^2/s) 1621 surf_usm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1622 ! m^3/mole ppm 1623 conv_to_ratio(k,j,i)*ratio2ppm* & 1624 ! kg/m^3 1625 rho_air(k) 1626 1627 !> OTHER SPECIES 1519 1520 ! 1521 !-- VOCs 1522 IF ( len_index_voc > 0 .AND. emt_att%species_name(match_spec_input(ispec)) == "VOC" ) THEN 1523 ! 1524 !-- (ppm/s) * m * kg/m^3 mole/(m^2/s) 1525 surf_usm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1526 conv_to_ratio(k,j,i) * & !< in m^3/mole 1527 ratio2ppm * & !< in ppm 1528 rho_air(k) !< in kg/m^3 1529 1530 ! 1531 !-- Other species 1628 1532 ELSE 1629 1533 1630 1631 ! ( ppm/s) * m * kg/m^3 kg/(m^2/s) 1632 surf_usm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1633 ! mole/Kg 1634 (1/emt_att%xm(ispec))* & 1635 ! m^3/mole ppm 1636 conv_to_ratio(k,j,i)*ratio2ppm* & 1637 ! kg/m^3 1638 rho_air(k) 1534 ! 1535 !-- (ppm/s) * m * kg/m^3 kg/(m^2/s) 1536 surf_usm_h%cssws(match_spec_model(ispec),m) = emis_distribution(1,j,i,ispec) * & 1537 ( 1.0_wp / emt_att%xm(ispec) ) * & !< in mole/kg 1538 conv_to_ratio(k,j,i) * & !< in m^3/mole 1539 ratio2ppm* & !< in ppm 1540 rho_air(k) !< in kg/m^3 1639 1541 1640 1542 … … 1646 1548 1647 1549 ENDDO 1648 END IF 1550 1551 ENDIF 1552 1649 1553 ENDDO 1650 1554 1651 ENDIF 1652 1653 1654 !> At the end of each call to chem_emissions setup, the time_factor is deallocated (ALLOCATED ONLY in the DEFAULT mode) 1655 1555 ENDIF 1556 1557 ! 1558 !-- Deallocate time_factor in case of DEFAULT mode) 1656 1559 IF ( ALLOCATED ( time_factor ) ) DEALLOCATE( time_factor ) 1657 1560 1658 ENDIF !> emis_output_required 1659 1561 ENDIF 1660 1562 1661 1563 END SUBROUTINE chem_emissions_setup -
palm/trunk/SOURCE/chem_modules.f90
r3458 r3611 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Minor formatting 30 ! 31 ! 3458 2018-10-30 14:51:23Z kanani 29 32 ! from chemistry branch r3443: 30 33 ! ?? 31 ! 34 ! 32 35 ! 3298 2018-10-02 12:21:11Z kanani 33 36 ! - Minor formatting … … 77 80 PUBLIC spc_names 78 81 79 INTEGER(iwp), DIMENSION(99) 80 INTEGER(iwp) ::ibc_cs_b !< integer flag for bc_cs_b81 INTEGER(iwp) ::ibc_cs_t !< integer flag for bc_cs_t82 INTEGER(iwp) ::cs_pr_count = 083 INTEGER(iwp) ::max_pr_cs = 084 INTEGER(iwp) ::cs_vertical_gradient_level_ind(99,10) = -9999 !< grid index values of cs_vertical_gradient_level_ind(s)85 86 LOGICAL ::constant_top_csflux(99) = .TRUE. !< chem spcs at the top orig .TRUE.87 LOGICAL ::constant_csflux(99) = .TRUE. !< chem spcs at namelist parameter orig TRUE88 LOGICAL ::call_chem_at_all_substeps = .FALSE. !< namelist parameter89 LOGICAL ::chem_debug0 = .FALSE. !< namelist parameter flag for minimum print output90 LOGICAL ::chem_debug1 = .FALSE. !< namelist parameter flag for print output91 LOGICAL ::chem_debug2 = .FALSE. !< namelist parameter flag for further print output92 LOGICAL ::chem_gasphase_on = .TRUE. !< namelist parameter93 LOGICAL ::emission_output_required = .TRUE. !< Logical Variable for requiring Emission Outputs94 LOGICAL ::do_emis = .FALSE. !< Flag for turning on chemistry emissions95 LOGICAL ::cs_pr_namelist_found = .FALSE. !< Namelist parameter: Names of t96 LOGICAL ::do_depo = .FALSE. !< namelist parameter for activation of deposition calculation97 98 99 82 INTEGER(iwp), DIMENSION(99) :: cs_pr_index = 0 83 INTEGER(iwp) :: ibc_cs_b !< integer flag for bc_cs_b 84 INTEGER(iwp) :: ibc_cs_t !< integer flag for bc_cs_t 85 INTEGER(iwp) :: cs_pr_count = 0 86 INTEGER(iwp) :: max_pr_cs = 0 87 INTEGER(iwp) :: cs_vertical_gradient_level_ind(99,10) = -9999 !< grid index values of cs_vertical_gradient_level_ind(s) 88 89 LOGICAL :: constant_top_csflux(99) = .TRUE. !< chem spcs at the top orig .TRUE. 90 LOGICAL :: constant_csflux(99) = .TRUE. !< chem spcs at namelist parameter orig TRUE 91 LOGICAL :: call_chem_at_all_substeps = .FALSE. !< namelist parameter 92 LOGICAL :: chem_debug0 = .FALSE. !< namelist parameter flag for minimum print output 93 LOGICAL :: chem_debug1 = .FALSE. !< namelist parameter flag for print output 94 LOGICAL :: chem_debug2 = .FALSE. !< namelist parameter flag for further print output 95 LOGICAL :: chem_gasphase_on = .TRUE. !< namelist parameter 96 LOGICAL :: emission_output_required = .TRUE. !< Logical Variable for requiring Emission Outputs 97 LOGICAL :: do_emis = .FALSE. !< Flag for turning on chemistry emissions 98 LOGICAL :: cs_pr_namelist_found = .FALSE. !< Namelist parameter: Names of t 99 LOGICAL :: do_depo = .FALSE. !< namelist parameter for activation of deposition calculation 100 101 102 ! 100 103 !-- Namelist parameters for creating initial chemistry profiles 101 REAL(wp) :: wall_csflux (99,0:5) = 0.0_wp !< namelist parameter102 REAL(wp) :: cs_vertical_gradient (99,10) = 0.0_wp !< namelist parameter103 REAL(wp) :: cs_vertical_gradient_level (99,10) = -999999.9_wp !< namelist parameter104 REAL(wp) :: top_csflux ( 99 ) = 0.0_wp !< namelist parameter105 REAL(wp) :: cs_surface_initial_change(99) = 0.0_wp !< namelist parameter106 REAL(wp) :: surface_csflux(99 ) = 0.0_wp !< namelist parameter: fluxes where 'surface_csflux_name' is in the namelist107 108 REAL(wp), DIMENSION(:), ALLOCATABLE ::bc_cs_t_val109 REAL(wp), DIMENSION(:), ALLOCATABLE :: css!< scaling parameter for chem spcs110 REAL(wp), DIMENSION(99) :: cs_surface = 0.0_wp !< Namelist parameter: Surface conc of chem spcs'111 REAL(wp), DIMENSION(99,100) :: cs_heights = 9999999.9_wp !< Namelist parameter: Height lvls(m) for cs_profiles112 REAL(wp), DIMENSION(99,100) :: cs_profile = 9999999.9_wp !< Namelist parameter: Chem conc for each spcs defined104 REAL(wp) :: wall_csflux (99,0:5) = 0.0_wp !< namelist parameter 105 REAL(wp) :: cs_vertical_gradient (99,10) = 0.0_wp !< namelist parameter 106 REAL(wp) :: cs_vertical_gradient_level (99,10) = -999999.9_wp !< namelist parameter 107 REAL(wp) :: top_csflux ( 99 ) = 0.0_wp !< namelist parameter 108 REAL(wp) :: cs_surface_initial_change(99) = 0.0_wp !< namelist parameter 109 REAL(wp) :: surface_csflux(99 ) = 0.0_wp !< namelist parameter: fluxes where 'surface_csflux_name' is in the namelist 110 111 REAL(wp), DIMENSION(:), ALLOCATABLE :: bc_cs_t_val 112 REAL(wp), DIMENSION(:), ALLOCATABLE :: css !< scaling parameter for chem spcs 113 REAL(wp), DIMENSION(99) :: cs_surface = 0.0_wp !< Namelist parameter: Surface conc of chem spcs' 114 REAL(wp), DIMENSION(99,100) :: cs_heights = 9999999.9_wp !< Namelist parameter: Height lvls(m) for cs_profiles 115 REAL(wp), DIMENSION(99,100) :: cs_profile = 9999999.9_wp !< Namelist parameter: Chem conc for each spcs defined 113 116 114 117 115 118 #if defined( __nopointer ) 116 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: cs !< chem spcs117 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: cs_p !< prognostic value of chem spc118 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tcs_m !< weighted tendency of cs for previous sub-timestep (Runge-Kutta)119 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: cs !< chem spcs 120 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: cs_p !< prognostic value of chem spc 121 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: tcs_m !< weighted tendency of cs for previous sub-timestep (Runge-Kutta) 119 122 120 123 #else 121 ! use pointers cs, cs_p and tcs_m to point arrays cs_1, cs_2, and cs_3 122 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE, TARGET :: cs_1 !< pointer for swapping of timelevels for respective quantity 123 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE, TARGET :: cs_2 !< pointer for swapping of timelevels for respective quantity 124 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE, TARGET :: cs_3 !< pointer for swapping of timelevels for respective quantity 125 REAL(wp), DIMENSION(:,:,:), POINTER :: cs !< pointer: sgs chem spcs) 126 REAL(wp), DIMENSION(:,:,:), POINTER :: cs_p !< pointer: prognostic value of sgs chem spcs 127 REAL(wp), DIMENSION(:,:,:), POINTER :: tcs_m !< pointer: 124 ! 125 !-- Use pointers cs, cs_p and tcs_m to point arrays cs_1, cs_2, and cs_3 126 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE, TARGET :: cs_1 !< pointer for swapping of timelevels for respective quantity 127 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE, TARGET :: cs_2 !< pointer for swapping of timelevels for respective quantity 128 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE, TARGET :: cs_3 !< pointer for swapping of timelevels for respective quantity 129 REAL(wp), DIMENSION(:,:,:), POINTER :: cs !< pointer: sgs chem spcs) 130 REAL(wp), DIMENSION(:,:,:), POINTER :: cs_p !< pointer: prognostic value of sgs chem spcs 131 REAL(wp), DIMENSION(:,:,:), POINTER :: tcs_m !< pointer: 128 132 129 133 #endif 130 134 131 CHARACTER (LEN=20) :: bc_cs_b = 'dirichlet'!< namelist parameter132 CHARACTER (LEN=20) :: bc_cs_t = 'initial_gradient'!< namelist parameter133 CHARACTER (LEN=11), DIMENSION(99) :: cs_name = 'novalue'!< Namelist parameter: chem spcs names134 CHARACTER (LEN=11), DIMENSION(99) :: cs_profile_name = 'novalue'!< Namelist parameter: Names of the chem for profiles135 CHARACTER (LEN=11), DIMENSION(99) :: surface_csflux_name = 'novalue'!< Namelist parameter: chem species surface fluxes names136 !< active chem spcs, default is 'novalue') ????137 CHARACTER (LEN=80) :: mode_emis ='PARAMETERIZED'!< Mode of chemistry emissions: DEFAULT .OR. EXPERT .OR.138 !PARAMETERIZED139 CHARACTER (LEN=80) :: time_fac_type ='MDH'!< Type of time treatment in the emis DEFAULT mode: HOUR .OR. MDH140 CHARACTER (LEN=80) :: daytype_mdh ='workday'!< Type of day in the MDH case: workday, weekend, holiday141 CHARACTER (LEN=11), DIMENSION(99) :: data_output_pr_cs = 'novalue'!< Namelist parameter: Names of the che m for profile output142 !< by cs_name for each height lvls defined by cs_heights135 CHARACTER (LEN=20) :: bc_cs_b = 'dirichlet' !< namelist parameter 136 CHARACTER (LEN=20) :: bc_cs_t = 'initial_gradient' !< namelist parameter 137 CHARACTER (LEN=11), DIMENSION(99) :: cs_name = 'novalue' !< Namelist parameter: chem spcs names 138 CHARACTER (LEN=11), DIMENSION(99) :: cs_profile_name = 'novalue' !< Namelist parameter: Names of the chem for profiles 139 CHARACTER (LEN=11), DIMENSION(99) :: surface_csflux_name = 'novalue' !< Namelist parameter: chem species surface fluxes names 140 !< active chem spcs, default is 'novalue') ???? 141 CHARACTER (LEN=80) :: mode_emis ='PARAMETERIZED' !< Mode of chemistry emissions: DEFAULT .OR. EXPERT .OR. 142 !< PARAMETERIZED 143 CHARACTER (LEN=80) :: time_fac_type ='MDH' !< Type of time treatment in the emis DEFAULT mode: HOUR .OR. MDH 144 CHARACTER (LEN=80) :: daytype_mdh ='workday' !< Type of day in the MDH case: workday, weekend, holiday 145 CHARACTER (LEN=11), DIMENSION(99) :: data_output_pr_cs = 'novalue' !< Namelist parameter: Names of the che m for profile output 146 !< by cs_name for each height lvls defined by cs_heights 143 147 ! 144 148 !-- Namelist parameters for chem_emissions … … 150 154 REAL(wp) :: emiss_factor_main ( 99 ) = -9999.0_wp 151 155 REAL(wp) :: emiss_factor_side ( 99 ) = -9999.0_wp 152 156 ! 153 157 !-- Other Emissions Variables 154 158 INTEGER(iwp) :: nspec_out !< Output of routine chem_emis_matching with 155 159 !< number of matched species 156 REAL(wp),ALLOCATABLE, DIMENSION(:,:,:,:) :: emis_distribution !> Emissions Final Values (main module output) 157 158 INTEGER(iwp),ALLOCATABLE,DIMENSION(:) :: match_spec_input !< Index of Input chem species for matching routine 159 INTEGER(iwp),ALLOCATABLE,DIMENSION(:) :: match_spec_model !< Index of Model chem species for matching routine 160 INTEGER(iwp),ALLOCATABLE,DIMENSION(:) :: match_spec_voc_input !< index of VOC input components matching the model's VOCs 161 INTEGER(iwp),ALLOCATABLE,DIMENSION(:) :: match_spec_voc_model !< index of VOC model species matching the input VOCs comp. 162 INTEGER(iwp),DIMENSION(:) :: match_spec_pm(1:3) !< results of matching the input and model's PMs 163 INTEGER(iwp),DIMENSION(:) :: match_spec_nox(1:2) !< results of matching the input and model's NOx 164 INTEGER(iwp),DIMENSION(:) :: match_spec_sox(1:2) !< results of matching the input and model's SOx! 165 ! TBD: evaluate whether to make them allocatable 166 ! and allocate their dimension in the matching 167 ! routine according to the number of components 168 ! matching between the model and the input files 169 160 REAL(wp),ALLOCATABLE, DIMENSION(:,:,:,:) :: emis_distribution !> Emissions Final Values (main module output) 161 162 INTEGER(iwp),ALLOCATABLE,DIMENSION(:) :: match_spec_input !< Index of Input chem species for matching routine 163 INTEGER(iwp),ALLOCATABLE,DIMENSION(:) :: match_spec_model !< Index of Model chem species for matching routine 164 INTEGER(iwp),ALLOCATABLE,DIMENSION(:) :: match_spec_voc_input !< index of VOC input components matching the model's VOCs 165 INTEGER(iwp),ALLOCATABLE,DIMENSION(:) :: match_spec_voc_model !< index of VOC model species matching the input VOCs comp. 166 INTEGER(iwp),DIMENSION(:) :: match_spec_pm(1:3) !< results of matching the input and model's PMs 167 INTEGER(iwp),DIMENSION(:) :: match_spec_nox(1:2) !< results of matching the input and model's NOx 168 INTEGER(iwp),DIMENSION(:) :: match_spec_sox(1:2) !< results of matching the input and model's SOx! 169 170 171 ! 172 !-- Selected atomic/molecular weights: 173 REAL, PARAMETER :: xm_H = 1.00790e-3 !< kg/mol 174 REAL, PARAMETER :: xm_N = 14.00670e-3 !< kg/mol 175 REAL, PARAMETER :: xm_C = 12.01115e-3 !< kg/mol 176 REAL, PARAMETER :: xm_S = 32.06400e-3 !< kg/mol 177 REAL, PARAMETER :: xm_O = 15.99940e-3 !< kg/mol 178 REAL, PARAMETER :: xm_F = 18.99840e-3 !< kg/mol 179 REAL, PARAMETER :: xm_Na = 22.98977e-3 !< kg/mol 180 REAL, PARAMETER :: xm_Cl = 35.45300e-3 !< kg/mol 181 REAL, PARAMETER :: xm_Rn222 = 222.00000e-3 !< kg/mol 182 REAL, PARAMETER :: xm_Pb210 = 210.00000e-3 !< kg/mol 183 REAL, PARAMETER :: xm_Ca = 40.07800e-3 !< kg/mol 184 REAL, PARAMETER :: xm_K = 39.09800e-3 !< kg/mol 185 REAL, PARAMETER :: xm_Mg = 24.30500e-3 !< kg/mol 186 REAL, PARAMETER :: xm_Pb = 207.20000e-3 !< kg/mol 187 REAL, PARAMETER :: xm_Cd = 112.41000e-3 !< kg/mol 170 188 171 !-- Selected atomic/molecular weights: 189 REAL, PARAMETER :: xm_h2o = xm_H * 2 + xm_O !< kg/mol 190 REAL, PARAMETER :: xm_o3 = xm_O * 3 !< kg/mol 191 REAL, PARAMETER :: xm_N2O5 = xm_N * 2 + xm_O * 5 !< kg/mol 192 REAL, PARAMETER :: xm_HNO3 = xm_H + xm_N + xm_O * 3 !< kg/mol 193 REAL, PARAMETER :: xm_NH4 = xm_N + xm_H * 4 !< kg/mol 194 REAL, PARAMETER :: xm_SO4 = xm_S + xm_O * 4 !< kg/mol 195 REAL, PARAMETER :: xm_NO3 = xm_N + xm_O * 3 !< kg/mol 196 REAL, PARAMETER :: xm_CO2 = xm_C + xm_O * 2 !< kg/mol 172 197 173 REAL, PARAMETER :: xm_H = 1.00790e-3 ! kg/mol 174 REAL, PARAMETER :: xm_N = 14.00670e-3 ! kg/mol 175 REAL, PARAMETER :: xm_C = 12.01115e-3 ! kg/mol 176 REAL, PARAMETER :: xm_S = 32.06400e-3 ! kg/mol 177 REAL, PARAMETER :: xm_O = 15.99940e-3 ! kg/mol 178 REAL, PARAMETER :: xm_F = 18.99840e-3 ! kg/mol 179 REAL, PARAMETER :: xm_Na = 22.98977e-3 ! kg/mol 180 REAL, PARAMETER :: xm_Cl = 35.45300e-3 ! kg/mol 181 REAL, PARAMETER :: xm_Rn222 = 222.00000e-3 ! kg/mol 182 REAL, PARAMETER :: xm_Pb210 = 210.00000e-3 ! kg/mol 183 REAL, PARAMETER :: xm_Ca = 40.07800e-3 ! kg/mol 184 REAL, PARAMETER :: xm_K = 39.09800e-3 ! kg/mol 185 REAL, PARAMETER :: xm_Mg = 24.30500e-3 ! kg/mol 186 REAL, PARAMETER :: xm_Pb = 207.20000e-3 ! kg/mol 187 REAL, PARAMETER :: xm_Cd = 112.41000e-3 ! kg/mol 188 189 REAL, PARAMETER :: xm_h2o = xm_H * 2 + xm_O ! kg/mol 190 REAL, PARAMETER :: xm_o3 = xm_O * 3 ! kg/mol 191 REAL, PARAMETER :: xm_N2O5 = xm_N * 2 + xm_O * 5 ! kg/mol 192 REAL, PARAMETER :: xm_HNO3 = xm_H + xm_N + xm_O * 3 ! kg/mol 193 REAL, PARAMETER :: xm_NH4 = xm_N + xm_H * 4 ! kg/mol 194 REAL, PARAMETER :: xm_SO4 = xm_S + xm_O * 4 ! kg/mol 195 REAL, PARAMETER :: xm_NO3 = xm_N + xm_O * 3 ! kg/mol 196 REAL, PARAMETER :: xm_CO2 = xm_C + xm_O * 2 ! kg/mol 197 198 ! mass of air 199 REAL, PARAMETER :: xm_air = 28.964e-3 ! kg/mol 198 ! 199 !-- mass of air 200 REAL, PARAMETER :: xm_air = 28.964e-3 !< kg/mol 200 201 201 ! dummy weight, used for complex molecules: 202 ! 203 !-- dummy weight, used for complex molecules: 202 204 REAL, PARAMETER :: xm_dummy = 1000.0e-3 ! kg/mol 203 205 -
palm/trunk/SOURCE/chemistry_model_mod.f90
r3600 r3611 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Minor formatting 30 ! 31 ! 3600 2018-12-04 13:49:07Z banzhafs 29 32 ! Code update to comply PALM coding rules 30 33 ! Bug fix in par_dir_diff subroutine … … 825 828 INTEGER(iwp) :: ss !< 826 829 REAL(wp), DIMENSION(nzb:nzt+1) :: cs_pr_init 827 REAL(wp), DIMENSION(nzb:nzt+1, nysg:nyng,nxlg:nxrg) :: cs_3d830 REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) :: cs_3d 828 831 REAL(wp) :: flag !< flag to mask topography grid points 829 832 … … 1207 1210 LOGICAL :: two_d !< flag parameter that indicates 2D variables (horizontal cross sections) 1208 1211 REAL(wp) :: fill_value 1209 REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) :: local_pf !<1212 REAL(wp), DIMENSION(nxl:nxr,nys:nyn,nzb:nzt+1) :: local_pf 1210 1213 1211 1214 ! … … 1294 1297 1295 1298 REAL(wp) :: fill_value !< 1296 REAL(sp), DIMENSION(nxl:nxr, nys:nyn,nzb_do:nzt_do) :: local_pf1299 REAL(sp), DIMENSION(nxl:nxr,nys:nyn,nzb_do:nzt_do) :: local_pf 1297 1300 1298 1301 … … 2359 2362 2360 2363 INTEGER(iwp),INTENT(IN) :: i, j, i_omp_start, tn, ilsp 2361 REAL(wp), DIMENSION(nzb+1:nzt, 2362 REAL(wp), DIMENSION(nzb+1:nzt, 2363 REAL(wp), DIMENSION(nzb+1:nzt, nys:nyn, 0:threads_per_task-1):: flux_l_cs !<2364 REAL(wp), DIMENSION(nzb+1:nzt, nys:nyn, 0:threads_per_task-1):: diss_l_cs !<2365 REAL(wp), DIMENSION(0:nz+1) :: pr_init_cs !<2364 REAL(wp), DIMENSION(nzb+1:nzt,0:threads_per_task-1) :: flux_s_cs !< 2365 REAL(wp), DIMENSION(nzb+1:nzt,0:threads_per_task-1) :: diss_s_cs !< 2366 REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn,0:threads_per_task-1) :: flux_l_cs !< 2367 REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn,0:threads_per_task-1) :: diss_l_cs !< 2368 REAL(wp), DIMENSION(0:nz+1) :: pr_init_cs !< 2366 2369 2367 2370 ! … … 2473 2476 LOGICAL, INTENT(OUT) :: found 2474 2477 2475 REAL(wp), DIMENSION(nzb:nzt+1, nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_3d !< 3D array to temp store data2478 REAL(wp), DIMENSION(nzb:nzt+1,nys_on_file-nbgp:nyn_on_file+nbgp,nxl_on_file-nbgp:nxr_on_file+nbgp) :: tmp_3d !< 3D array to temp store data 2476 2479 2477 2480
Note: See TracChangeset
for help on using the changeset viewer.
