Changeset 4403 for palm/trunk/SOURCE/chemistry_model_mod.f90

Timestamp:

Feb 12, 2020 1:08:46 PM (4 years ago)

Author:

banzhafs

Message:

chemistry model: implemented on-demand emission read mode for LOD 2

File:

• palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -22 +22 @@
 ! Current revisions:
 ! -----------------
-!
+! 
 !
 ! Former revisions:
 ! -----------------
 ! $Id$
+! in subroutine chem_init (ECC)
+! - allows different init paths emission data for legacy
+!   mode emission and on-demand mode
+! in subroutine chem_init_internal (ECC)
+! - reads netcdf file only when legacy mode is activated
+!   (i.e., emiss_read_legacy_mode = .TRUE.)
+!   otherwise file is read once at the beginning to obtain
+!   header information, and emission data are extracted on
+!   an on-demand basis
+!
+! 4372 2020-01-14 10:20:35Z banzhafs
 ! chem_parin : added handler for new namelist item emiss_legacy_read_mode (ECC)
 ! added messages
@@ -1773 +1784 @@
  SUBROUTINE chem_init
 
+!
+!-- 20200203 (ECC)
+!-- introduced additional interfaces for on-demand emission update
+
+!    USE chem_emissions_mod,                                                    &
+!        ONLY:  chem_emissions_init
+        
     USE chem_emissions_mod,                                                    &
-        ONLY:  chem_emissions_init
+        ONLY:  chem_emissions_init, chem_emissions_header_init
         
     USE netcdf_data_input_mod,                                                 &
@@ -1791 +1809 @@
            ilu_permanent_crops + ilu_savanna + ilu_semi_natural_veg + ilu_tropical_forest +    &
            ilu_urban ) == 0 )  CONTINUE
-         
-    IF ( emissions_anthropogenic )  CALL chem_emissions_init
+
+!
+!-- 20200203 (ECC)
+!-- calls specific emisisons initialization subroutines
+!-- for legacy mode and on-demand mode         
+
+!    IF ( emissions_anthropogenic )  CALL chem_emissions_init
+
+    IF  ( emissions_anthropogenic )  THEN
+
+       IF  ( emiss_read_legacy_mode )  THEN
+          CALL chem_emissions_init
+       ELSE
+          CALL chem_emissions_header_init
+       ENDIF
+
+    ENDIF
+
+
 !
 !-- Chemistry variables will be initialized if availabe from dynamic 
@@ -1835 +1870 @@
     INTEGER(iwp) ::  lpr_lev           !< running index for chem spcs profile level
 
+!
+!-- 20200203 ECC
+!-- reads netcdf data only under legacy mode
+
+!    IF ( emissions_anthropogenic )  THEN
+!       CALL netcdf_data_input_chemistry_data( chem_emis_att, chem_emis )
+!    ENDIF
+
     IF ( emissions_anthropogenic )  THEN
-       CALL netcdf_data_input_chemistry_data( chem_emis_att, chem_emis )
+       IF ( emiss_read_legacy_mode )  THEN
+          CALL netcdf_data_input_chemistry_data( chem_emis_att, chem_emis )
+       ENDIF
     ENDIF
+
 !
 !-- Allocate memory for chemical species
@@ -2476 +2522 @@
                               'Reverting to legacy emission read mode'
              CALL message ( 'parin_chem', 'CM0466', 0, 0, 0, 6, 0 )
+
              emiss_read_legacy_mode = .TRUE.
 
@@ -2482 +2529 @@
              message_string = 'New emission read mode activated'           // &
                               CHAR(10)  //  '                    '         // &
-                              'LOD 2 emissions will be read on an '        // &           
-                              'hourly basis according to'                  // &             
-                              CHAR(10)  //  '                    '         // &
-                              'indicated timestamps'
+                              'LOD 2 emissions will be updated on-demand ' // &           
+                              'according to indicated timestamps'
              CALL message ( 'parin_chem', 'CM0467', 0, 0, 0, 6, 0 )