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

Timestamp:

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

Author:

banzhafs

Message:

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

File:

• palm/trunk/SOURCE/time_integration.f90 (modified) (4 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Allowing both existing and on-demand emission read modes
+!
+! 4360 2020-01-07 11:25:50Z suehring
 ! Bugfix, hour_call_emis uninitialized at first call of time_integration
 ! 
@@ -208 +211 @@
 
     USE chem_emissions_mod,                                                                        &
-        ONLY:  chem_emissions_setup
+        ONLY:  chem_emissions_setup, chem_emissions_update_on_demand
 
     USE chem_gasphase_mod,                                                                         &
@@ -214 +217 @@
 
     USE chem_modules,                                                                              &
-        ONLY:  bc_cs_t_val, chem_species, cs_name, emissions_anthropogenic, n_matched_vars
+        ONLY:  bc_cs_t_val, chem_species, cs_name,                                                 &
+               emissions_anthropogenic, emiss_read_legacy_mode,                                    &
+               n_matched_vars
 
     USE chemistry_model_mod,                                                                       &
@@ -1162 +1167 @@
        ENDIF
 
-!
-!--    If required, consider chemical emissions
-       IF ( air_chemistry  .AND.  emissions_anthropogenic )  THEN
-!
-!--       Call emission routine only once an hour
-          !> @todo Move check of hour_call_emis into chem_emissions_setup
-          CALL get_date_time( time_since_reference_point, hour=hour )
-          IF ( hour_call_emis /= hour )  THEN
-          !> @question Which of these two if-statements is correct?
-          ! IF ( hour_of_year  >  hour_call_emis )  THEN
-             CALL chem_emissions_setup( chem_emis_att, chem_emis, n_matched_vars )
-             hour_call_emis = hour
-          ENDIF
-       ENDIF
+       
+!
+!-- 20200203 (ECC)
+!-- allows for emission update mode in legacy mode as well as on-demand mode
+!-- note that under on-demand mode emission update is no longer restricted to
+!-- an hourly frequency, but whenever the simulation time corresponds to an
+!-- inrement in emission timestamp value
+
+!
+!-- If required, consider chemical emissions
+
+       IF  ( air_chemistry .AND. emissions_anthropogenic )  THEN
+
+          IF  ( emiss_read_legacy_mode )  THEN
+!
+!-- get hourly index and updates emission data when the hour is passed
+
+             CALL get_date_time( time_since_reference_point, hour=hour )
+
+             IF  ( hour_call_emis /= hour )   THEN
+
+                CALL chem_emissions_setup( chem_emis_att, chem_emis, n_matched_vars )
+                hour_call_emis = hour
+
+             ENDIF
+
+          ELSE
+
+             CALL chem_emissions_update_on_demand
+
+          ENDIF
+
+       ENDIF
+
+
 !
 !--    If required, consider aerosol emissions for the salsa model