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

Timestamp:

Apr 8, 2019 7:09:16 PM (5 years ago)

Author:

knoop

Message:

Added chem_actions module interface to chemistry_model_mod and moved photolysis_control call into it

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -300 +300 @@
          ONLY:  bc_lr_cyc, bc_ns_cyc, dt_3d, humidity, initializing_actions, message_string,       &
          omega, tsc, intermediate_timestep_count, intermediate_timestep_count_max,                 &
-         max_pr_user, timestep_scheme, use_prescribed_profile_data, ws_scheme_sca         
+         max_pr_user, timestep_scheme, use_prescribed_profile_data, ws_scheme_sca, air_chemistry
 
     USE arrays_3d,                                                                                 &
@@ -310 +310 @@
 
     USE chem_modules
+
+    USE chem_photolysis_mod,                                                                       &
+        ONLY:  photolysis_control
 
     USE statistics
@@ -478 +481 @@
     END INTERFACE chem_parin
 
+    INTERFACE chem_actions
+       MODULE PROCEDURE chem_actions
+       MODULE PROCEDURE chem_actions_ij
+    END INTERFACE chem_actions
+
     INTERFACE chem_prognostic_equations
        MODULE PROCEDURE chem_prognostic_equations
@@ -574 +582 @@
          chem_define_netcdf_grid, chem_header, chem_init, chem_init_arrays,   &
          chem_init_profiles, chem_integrate, chem_parin,                      &
-         chem_prognostic_equations, chem_rrd_local,                           &
+         chem_actions, chem_prognostic_equations, chem_rrd_local,             &
          chem_statistics, chem_swap_timelevel, chem_wrd_local, chem_depo 
 
@@ -1084 +1092 @@
     USE arrays_3d
     USE control_parameters,                                                    &
-        ONLY:  air_chemistry, data_output_pr, message_string
+        ONLY:  data_output_pr, message_string
     USE indices
     USE profil_parameter
@@ -2281 +2289 @@
 
  END SUBROUTINE chem_parin
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Call for all grid points
+!------------------------------------------------------------------------------!
+    SUBROUTINE chem_actions( location )
+
+
+    CHARACTER (LEN=*), INTENT(IN) ::  location !< call location string
+
+    SELECT CASE ( location )
+
+       CASE ( 'before_prognostic_equations' )
+!
+!--       Chemical reactions and deposition
+          IF ( chem_gasphase_on ) THEN
+!
+!--          If required, calculate photolysis frequencies -
+!--          UNFINISHED: Why not before the intermediate timestep loop?
+             IF ( intermediate_timestep_count ==  1 )  THEN
+                CALL photolysis_control
+             ENDIF
+
+          ENDIF
+
+       CASE DEFAULT
+          CONTINUE
+
+    END SELECT
+
+    END SUBROUTINE chem_actions
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Call for grid points i,j
+!------------------------------------------------------------------------------!
+
+    SUBROUTINE chem_actions_ij( i, j, location )
+
+
+    INTEGER(iwp),      INTENT(IN) ::  i         !< grid index in x-direction
+    INTEGER(iwp),      INTENT(IN) ::  j         !< grid index in y-direction
+    CHARACTER (LEN=*), INTENT(IN) ::  location  !< call location string
+    INTEGER(iwp)  ::  dummy  !< call location string
+
+    IF ( air_chemistry    )   dummy = i + j
+
+    SELECT CASE ( location )
+
+       CASE DEFAULT
+          CONTINUE
+
+    END SELECT
+
+
+    END SUBROUTINE chem_actions_ij