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

Timestamp:

Oct 7, 2020 1:25:11 PM (4 years ago)

Author:

schwenkel

Message:

Move exchange_horiz from time_integration to modules

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Move exchange_horiz from time_integration to modules
+!
+! 4671 2020-09-09 20:27:58Z pavelkrc
 ! Implementation of downward facing USM and LSM surfaces
 !
@@ -2934 +2937 @@
 !> routine for exchange horiz of chemical quantities
 !--------------------------------------------------------------------------------------------------!
- SUBROUTINE chem_exchange_horiz_bounds
+ SUBROUTINE chem_exchange_horiz_bounds( location )
 
     USE exchange_horiz_mod,                                                                        &
@@ -2940 +2943 @@
 
    INTEGER(iwp) ::  lsp       !<
-
-!
-!--    Loop over chemical species
-       CALL cpu_log( log_point_s(84), 'chem.exch-horiz', 'start' )
-       DO  lsp = 1, nvar
-          CALL exchange_horiz( chem_species(lsp)%conc, nbgp,                                       &
-                               alternative_communicator = communicator_chem )
-       ENDDO
-
-       CALL chem_boundary_conditions( horizontal_conditions_only = .TRUE. )
-
-       CALL cpu_log( log_point_s(84), 'chem.exch-horiz', 'stop' )
-
+   INTEGER(iwp) ::  n
+
+   CHARACTER (LEN=*), INTENT(IN) ::  location !< call location string
+
+   SELECT CASE ( location )
+
+       CASE ( 'before_prognostic_equation' )
+!
+!--       Loop over chemical species
+          CALL cpu_log( log_point_s(84), 'chem.exch-horiz', 'start' )
+          DO  lsp = 1, nvar
+             CALL exchange_horiz( chem_species(lsp)%conc, nbgp,                                    &
+                                  alternative_communicator = communicator_chem )
+          ENDDO
+
+          CALL chem_boundary_conditions( horizontal_conditions_only = .TRUE. )
+
+          CALL cpu_log( log_point_s(84), 'chem.exch-horiz', 'stop' )
+
+       CASE ( 'after_prognostic_equation' )
+
+          IF ( air_chemistry )  THEN
+             DO  n = 1, nvar
+                CALL exchange_horiz( chem_species(n)%conc_p, nbgp,                                 &
+                                     alternative_communicator = communicator_chem )
+             ENDDO
+          ENDIF
+
+       CASE ( 'after_anterpolation' )
+
+          IF ( air_chemistry )  THEN
+             DO  n = 1, nvar
+                CALL exchange_horiz( chem_species(n)%conc, nbgp,                                   &
+                                     alternative_communicator = communicator_chem )
+             ENDDO
+          ENDIF
+
+    END SELECT
 
  END SUBROUTINE chem_exchange_horiz_bounds