Changeset 4847 for palm/trunk/SOURCE/chemistry_model_mod.f90
- Timestamp:
- Jan 19, 2021 1:44:45 PM (4 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (31 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/chemistry_model_mod.f90
r4843 r4847 26 26 ! ----------------- 27 27 ! $Id$ 28 ! error message IDs changed from CM to the default PA 29 ! 30 ! 4843 2021-01-15 15:22:11Z raasch 28 31 ! local namelist parameter added to switch off the module although the respective module namelist 29 32 ! appears in the namelist file 30 33 ! 31 34 ! 4833 2021-01-07 08:57:54Z raasch 32 ! openmp bugfix for time measurements of non advective processes 35 ! openmp bugfix for time measurements of non advective processes 33 36 ! 34 37 ! 4828 2021-01-05 11:21:41Z Giersch … … 1119 1122 message_string = 'data_output_pr = ' // TRIM( data_output_pr(var_count) ) // & 1120 1123 ' is not implemented for air_chemistry = .FALSE.' 1121 CALL message( 'chem_check_ parameters', 'CM0433', 1, 2, 0, 6, 0 )1124 CALL message( 'chem_check_data_output_pr', 'PA0293', 1, 2, 0, 6, 0 ) 1122 1125 1123 1126 ENDIF … … 1129 1132 ' is currently not implemented. Please output resolved- and '// & 1130 1133 'subgrid-scale fluxes individually to obtain the total flux.' 1131 CALL message( 'chem_check_ parameters', 'CM0498', 1, 2, 0, 6, 0 )1134 CALL message( 'chem_check_data_output_pr', 'PA0487', 1, 2, 0, 6, 0 ) 1132 1135 ENDIF 1133 1136 ENDIF … … 1228 1231 IF ( chem_gasphase_on ) THEN 1229 1232 message_string = 'Chemical reactions: ON' 1230 CALL message( 'chem_check_parameters', ' CM0421', 0, 0, 0, 6, 0 )1231 ELSEIF ( .NOT. ( chem_gasphase_on) ) THEN1233 CALL message( 'chem_check_parameters', 'PA0517', 0, 0, 0, 6, 0 ) 1234 ELSEIF ( .NOT. ( chem_gasphase_on ) ) THEN 1232 1235 message_string = 'Chemical reactions: OFF' 1233 CALL message( 'chem_check_parameters', ' CM0422', 0, 0, 0, 6, 0 )1236 CALL message( 'chem_check_parameters', 'PA0518', 0, 0, 0, 6, 0 ) 1234 1237 ENDIF 1235 1238 ! … … 1237 1240 IF ( call_chem_at_all_substeps ) THEN 1238 1241 message_string = 'Chemistry is calculated at all meteorology time-step' 1239 CALL message( 'chem_check_parameters', ' CM0423', 0, 0, 0, 6, 0 )1240 ELSEIF ( . not. (call_chem_at_all_substeps) ) THEN1242 CALL message( 'chem_check_parameters', 'PA0522', 0, 0, 0, 6, 0 ) 1243 ELSEIF ( .NOT. ( call_chem_at_all_substeps ) ) THEN 1241 1244 message_string = 'Sub-time-steps are skipped for chemistry time-steps' 1242 CALL message( 'chem_check_parameters', ' CM0424', 0, 0, 0, 6, 0 )1245 CALL message( 'chem_check_parameters', 'PA0523', 0, 0, 0, 6, 0 ) 1243 1246 ENDIF 1244 1247 ! 1245 1248 !-- Check for photolysis scheme 1246 IF ( ( photolysis_scheme /= 'simple') .AND. (photolysis_scheme /= 'constant') ) THEN1249 IF ( ( photolysis_scheme /= 'simple' ) .AND. ( photolysis_scheme /= 'constant' ) ) THEN 1247 1250 message_string = 'Incorrect photolysis scheme selected, please check spelling' 1248 CALL message( 'chem_check_parameters', ' CM0425', 1, 2, 0, 6, 0 )1251 CALL message( 'chem_check_parameters', 'PA0593', 1, 2, 0, 6, 0 ) 1249 1252 ENDIF 1250 1253 ! … … 1254 1257 message_string = & 1255 1258 'Incorrect chemistry mechanism selected, check spelling in namelist and/or chem_gasphase_mod' 1256 CALL message( 'chem_check_parameters', ' CM0462', 1, 2, 0, 6, 0 )1259 CALL message( 'chem_check_parameters', 'PA0675', 1, 2, 0, 6, 0 ) 1257 1260 ENDIF 1258 1261 … … 1265 1268 ELSE 1266 1269 message_string = 'unknown boundary condition: bc_cs_b ="' // TRIM( bc_cs_b ) // '"' 1267 CALL message( 'chem_ boundary_conds', 'CM0429', 1, 2, 0, 6, 0 )1270 CALL message( 'chem_check_parameters', 'PA0730', 1, 2, 0, 6, 0 ) 1268 1271 ENDIF 1269 1272 ! … … 1279 1282 ELSE 1280 1283 message_string = 'unknown boundary condition: bc_c_t ="' // TRIM( bc_cs_t ) // '"' 1281 CALL message( 'che ck_parameters', 'CM0430', 1, 2, 0, 6, 0 )1284 CALL message( 'chem_check_parameters', 'PA0731', 1, 2, 0, 6, 0 ) 1282 1285 ENDIF 1283 1286 … … 1407 1410 message_string = 'Unused/incorrect input for initial surface value: ' // & 1408 1411 TRIM( cs_name(lsp_usr) ) 1409 CALL message( 'chem_check_parameters', ' CM0427', 1, 2, 0, 6, 0 )1412 CALL message( 'chem_check_parameters', 'PA0732', 1, 2, 0, 6, 0 ) 1410 1413 ENDIF 1411 1414 lsp_usr = lsp_usr + 1 … … 1426 1429 message_string = 'Unused/incorrect input of chemical species for surface emission fluxes: ' & 1427 1430 // TRIM( surface_csflux_name(lsp_usr) ) 1428 CALL message( 'chem_check_parameters', ' CM0428', 1, 2, 0, 6, 0 )1431 CALL message( 'chem_check_parameters', 'PA0733', 1, 2, 0, 6, 0 ) 1429 1432 ENDIF 1430 1433 lsp_usr = lsp_usr + 1 … … 2342 2345 IF ( cs_heights(1,1) /= 0.0_wp ) THEN 2343 2346 message_string = 'The surface value of cs_heights must be 0.0' 2344 CALL message( 'chem_check_parameters', ' CM0434', 1, 2, 0, 6, 0 )2347 CALL message( 'chem_check_parameters', 'PA0734', 1, 2, 0, 6, 0 ) 2345 2348 ENDIF 2346 2349 … … 2355 2358 IF ( npr_lev == 100 ) THEN 2356 2359 message_string = 'number of chem spcs exceeding the limit' 2357 CALL message( 'chem_check_parameters', ' CM0435', 1, 2, 0, 6, 0 )2360 CALL message( 'chem_check_parameters', 'PA0735', 1, 2, 0, 6, 0 ) 2358 2361 EXIT 2359 2362 ENDIF … … 2625 2628 ( mode_emis /= 'PRE-PROCESSED' ) ) THEN 2626 2629 message_string = 'Incorrect mode_emiss option select. Please check spelling' 2627 CALL message( 'chem_ check_parameters', 'CM0436', 1, 2, 0, 6, 0 )2630 CALL message( 'chem_parin', 'PA0736', 1, 2, 0, 6, 0 ) 2628 2631 ENDIF 2629 2632 ELSE … … 2632 2635 ( emiss_lod /= 2 ) ) THEN 2633 2636 message_string = 'Invalid value for emiss_lod (0, 1, or 2)' 2634 CALL message( 'chem_ check_parameters', 'CM0436', 1, 2, 0, 6, 0 )2637 CALL message( 'chem_parin', 'PA0737', 1, 2, 0, 6, 0 ) 2635 2638 ENDIF 2636 2639 ENDIF … … 2640 2643 ! IF ( (mode_emis /= 'PARAMETERIZED') .AND. ( mode_emis /= 'DEFAULT' ) .AND. ( mode_emis /= 'PRE-PROCESSED' ) ) THEN 2641 2644 ! message_string = 'Incorrect mode_emiss option select. Please check spelling' 2642 ! CALL message( 'chem_ check_parameters', 'CM0436', 1, 2, 0, 6, 0 )2645 ! CALL message( 'chem_parin', 'PA0738', 1, 2, 0, 6, 0 ) 2643 2646 ! ENDIF 2644 2647 … … 2661 2664 END SELECT 2662 2665 2663 message_string = 'Synchronizing mode_emis to defined emiss_lod' // & 2664 CHAR( 10 ) // ' ' // & 2665 'NOTE - mode_emis will be depreciated in future releases' // & 2666 CHAR( 10 ) // ' ' // & 2666 message_string = 'Synchronizing mode_emis to defined emiss_lod&' // & 2667 'NOTE - mode_emis will be depreciated in future releases&' // & 2667 2668 'please use emiss_lod to define emission mode' 2668 2669 2669 CALL message ( 'parin_chem', 'CM0463', 0, 0, 0, 6, 0 )2670 CALL message( 'chem_parin', 'PA0739', 0, 0, 0, 6, 0 ) 2670 2671 2671 2672 ELSE ! if emiss_lod is not set … … 2680 2681 END SELECT 2681 2682 2682 message_string = 'emiss_lod undefined. Using existing mod_emiss setting' // & 2683 CHAR( 10 ) // ' ' // & 2684 'NOTE - mode_emis will be depreciated in future releases' // & 2685 CHAR( 10 ) // ' ' // & 2683 message_string = 'emiss_lod undefined. Using existing mod_emiss setting&' // & 2684 'NOTE - mode_emis will be depreciated in future releases&' // & 2686 2685 'please use emiss_lod to define emission mode' 2687 2686 2688 CALL message ( 'parin_chem', 'CM0464', 0, 0, 0, 6, 0 )2687 CALL message( 'chem_parin', 'PA0740', 0, 0, 0, 6, 0 ) 2689 2688 ENDIF 2690 2689 … … 2702 2701 IF ( emiss_read_legacy_mode ) THEN !< notify legacy read mode 2703 2702 2704 message_string = 'Legacy emission read mode activated' // & 2705 CHAR( 10 ) // ' ' // & 2703 message_string = 'Legacy emission read mode activated&' // & 2706 2704 'All emissions data will be loaded ' // & 2707 2705 'prior to start of simulation' 2708 CALL message ( 'parin_chem', 'CM0465', 0, 0, 0, 6, 0 )2706 CALL message( 'chem_parin', 'PA0741', 0, 0, 0, 6, 0 ) 2709 2707 2710 2708 ELSE !< if new read mode selected … … 2713 2711 2714 2712 message_string = 'New emission read mode ' // & 2715 'currently unavailable for LODs 0 and 1.' // & 2716 CHAR( 10 ) // ' ' // & 2713 'currently unavailable for LODs 0 and 1.&' // & 2717 2714 'Reverting to legacy emission read mode' 2718 CALL message ( 'parin_chem', 'CM0466', 0, 0, 0, 6, 0 )2715 CALL message( 'chem_parin', 'PA0742', 0, 0, 0, 6, 0 ) 2719 2716 2720 2717 emiss_read_legacy_mode = .TRUE. … … 2722 2719 ELSE !< notify new read mode 2723 2720 2724 message_string = 'New emission read mode activated' // & 2725 CHAR( 10 ) // ' ' // & 2721 message_string = 'New emission read mode activated&' // & 2726 2722 'LOD 2 emissions will be updated on-demand ' // & 2727 2723 'according to indicated timestamps' 2728 CALL message ( 'parin_chem', 'CM0467', 0, 0, 0, 6, 0 )2724 CALL message( 'chem_parin', 'PA0743', 0, 0, 0, 6, 0 ) 2729 2725 2730 2726 ENDIF … … 3877 3873 luv_palm = surf_lsm_h(0)%vegetation_type(m) 3878 3874 IF ( luv_palm == ind_luv_user ) THEN 3879 message_string = 'No vegetation type defined. Please define vegetation type to enable deposition calculation' 3880 CALL message( 'chem_depo', 'CM0451', 1, 2, 0, 6, 0 ) 3875 message_string = 'No vegetation type defined. Please define vegetation type to' // & 3876 ' enable deposition calculation' 3877 CALL message( 'chem_depo', 'PA0744', 1, 2, 0, 6, 0 ) 3881 3878 ELSEIF ( luv_palm == ind_luv_b_soil ) THEN 3882 3879 luv_dep = 9 … … 3921 3918 lup_palm = surf_lsm_h(0)%pavement_type(m) 3922 3919 IF ( lup_palm == ind_lup_user ) THEN 3923 message_string = 'No pavement type defined. Please define pavement type to enable deposition calculation' 3924 CALL message( 'chem_depo', 'CM0452', 1, 2, 0, 6, 0 ) 3920 message_string = 'No pavement type defined. Please define pavement type to enable' // & 3921 ' deposition calculation' 3922 CALL message( 'chem_depo', 'PA0745', 1, 2, 0, 6, 0 ) 3925 3923 ELSEIF ( lup_palm == ind_lup_asph_conc ) THEN 3926 3924 lup_dep = 9 … … 3959 3957 luw_palm = surf_lsm_h(0)%water_type(m) 3960 3958 IF ( luw_palm == ind_luw_user ) THEN 3961 message_string = 'No water type defined. Please define water type to enable deposition calculation' 3962 CALL message( 'chem_depo', 'CM0453', 1, 2, 0, 6, 0 ) 3959 message_string = 'No water type defined. Please define water type to enable' // & 3960 ' deposition calculation' 3961 CALL message( 'chem_depo', 'PA0746', 1, 2, 0, 6, 0 ) 3963 3962 ELSEIF ( luw_palm == ind_luw_lake ) THEN 3964 3963 luw_dep = 13 … … 4457 4456 lug_palm = ind_luv_s_grass 4458 4457 IF ( lug_palm == ind_luv_user ) THEN 4459 message_string = 'No vegetation type defined. Please define vegetation type to enable deposition calculation' 4460 CALL message( 'chem_depo', 'CM0454', 1, 2, 0, 6, 0 ) 4458 message_string = 'No vegetation type defined. Please define vegetation type to' // & 4459 ' enable deposition calculation' 4460 CALL message( 'chem_depo', 'PA0747', 1, 2, 0, 6, 0 ) 4461 4461 ELSEIF ( lug_palm == ind_luv_b_soil ) THEN 4462 4462 lug_dep = 9 … … 4504 4504 luu_palm = ind_lup_conc 4505 4505 IF ( luu_palm == ind_lup_user ) THEN 4506 message_string = 'No pavement type defined. Please define pavement type to enable deposition calculation' 4507 CALL message( 'chem_depo', 'CM0455', 1, 2, 0, 6, 0 ) 4506 message_string = 'No pavement type defined. Please define pavement type to enable' // & 4507 ' deposition calculation' 4508 CALL message( 'chem_depo', 'PA0748', 1, 2, 0, 6, 0 ) 4508 4509 ELSEIF ( luu_palm == ind_lup_asph_conc ) THEN 4509 4510 luu_dep = 9 … … 4545 4546 lud_palm = ind_lup_metal 4546 4547 IF ( lud_palm == ind_lup_user ) THEN 4547 message_string = 'No pavement type defined. Please define pavement type to enable deposition calculation' 4548 CALL message( 'chem_depo', 'CM0456', 1, 2, 0, 6, 0 ) 4548 message_string = 'No pavement type defined. Please define pavement type to enable' // & 4549 ' deposition calculation' 4550 CALL message( 'chem_depo', 'PA0749', 1, 2, 0, 6, 0 ) 4549 4551 ELSEIF ( lud_palm == ind_lup_asph_conc ) THEN 4550 4552 lud_dep = 9 … … 5303 5305 CASE default 5304 5306 message_string = 'Component '// TRIM( compnam ) // ' not supported' 5305 CALL message( 'rc_special', ' CM0457', 1, 2, 0, 6, 0 )5307 CALL message( 'rc_special', 'PA0750', 1, 2, 0, 6, 0 ) 5306 5308 END SELECT 5307 5309 … … 5356 5358 CASE default 5357 5359 message_string = 'Component '// TRIM( compnam ) // ' not supported' 5358 CALL message( 'rc_gw', ' CM0458', 1, 2, 0, 6, 0 )5360 CALL message( 'rc_gw', 'PA0751', 1, 2, 0, 6, 0 ) 5359 5361 END SELECT 5360 5362 … … 5569 5571 CASE default 5570 5572 message_string = 'Component '// TRIM( compnam ) // ' not supported' 5571 CALL message( 'rc_gstom', ' CM0459', 1, 2, 0, 6, 0 )5573 CALL message( 'rc_gstom', 'PA0752', 1, 2, 0, 6, 0 ) 5572 5574 END SELECT 5573 5575 … … 5925 5927 ELSE 5926 5928 message_string = 'nwet can only be 0 or 1' 5927 CALL message( 'rc_gsoil_eff', ' CM0460', 1, 2, 0, 6, 0 )5929 CALL message( 'rc_gsoil_eff', 'PA0753', 1, 2, 0, 6, 0 ) 5928 5930 ENDIF 5929 5931 ENDIF … … 6101 6103 ! rc_eff = -999.0_wp 6102 6104 ! message_string = 'This should not be possible, check ccomp_tot' 6103 ! CALL message( 'rc_comp_point_rc_eff', ' CM0461', 1, 2, 0, 6, 0 )6105 ! CALL message( 'rc_comp_point_rc_eff', 'PA0754', 1, 2, 0, 6, 0 ) 6104 6106 ! ENDIF 6105 6107 !
Note: See TracChangeset
for help on using the changeset viewer.
