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

Timestamp:

Apr 8, 2019 7:35:54 PM (5 years ago)

Author:

knoop

Message:

Added chem_non_transport_physics module interface to chemistry_model_mod and moved respective calls into it

File:

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

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -314 +314 @@
         ONLY:  photolysis_control
 
+    USE cpulog,                                                                                    &
+        ONLY:  cpu_log, log_point_s
+
     USE statistics
 
@@ -486 +489 @@
     END INTERFACE chem_actions
 
+    INTERFACE chem_non_transport_physics
+       MODULE PROCEDURE chem_non_transport_physics
+       MODULE PROCEDURE chem_non_transport_physics_ij
+    END INTERFACE chem_non_transport_physics
+
     INTERFACE chem_prognostic_equations
        MODULE PROCEDURE chem_prognostic_equations
@@ -583 +591 @@
          chem_init_profiles, chem_integrate, chem_parin,                      &
          chem_actions, chem_prognostic_equations, chem_rrd_local,             &
-         chem_statistics, chem_swap_timelevel, chem_wrd_local, chem_depo 
+         chem_statistics, chem_swap_timelevel, chem_wrd_local, chem_depo,     &
+         chem_non_transport_physics
 
  CONTAINS
@@ -2349 +2358 @@
 
     END SUBROUTINE chem_actions_ij
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Call for all grid points
+!------------------------------------------------------------------------------!
+ SUBROUTINE chem_non_transport_physics()
+
+
+    INTEGER(iwp) ::  i  !<
+    INTEGER(iwp) ::  j  !<
+
+!
+!-- Calculation of chemical reactions and deposition.
+    IF ( chem_gasphase_on ) THEN
+
+       IF ( intermediate_timestep_count == 1 .OR. call_chem_at_all_substeps )  THEN
+
+          CALL cpu_log( log_point_s(19), 'chem.reactions', 'start' )
+          !$OMP PARALLEL PRIVATE (i,j)
+          !$OMP DO schedule(static,1)
+          DO  i = nxl, nxr
+             DO  j = nys, nyn
+                CALL chem_integrate( i, j )
+             ENDDO
+          ENDDO
+          !$OMP END PARALLEL
+          CALL cpu_log( log_point_s(19), 'chem.reactions', 'stop' )
+
+          IF ( deposition_dry )  THEN
+             CALL cpu_log( log_point_s(24), 'chem.deposition', 'start' )
+             DO  i = nxl, nxr
+                DO  j = nys, nyn
+                   CALL chem_depo( i, j )
+                ENDDO
+             ENDDO
+             CALL cpu_log( log_point_s(24), 'chem.deposition', 'stop' )
+          ENDIF
+
+       ENDIF
+
+    ENDIF
+
+ END SUBROUTINE chem_non_transport_physics
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Call for grid points i,j
+!------------------------------------------------------------------------------!
+
+ SUBROUTINE chem_non_transport_physics_ij( i, j )
+
+
+    INTEGER(iwp), INTENT(IN) ::  i  !< grid index in x-direction
+    INTEGER(iwp), INTENT(IN) ::  j  !< grid index in y-direction
+
+!
+!-- Calculation of chemical reactions and deposition.
+    IF ( chem_gasphase_on ) THEN
+
+       IF ( intermediate_timestep_count == 1 .OR. call_chem_at_all_substeps )  THEN
+
+          CALL cpu_log( log_point_s(19), 'chem.reactions', 'start' )
+          CALL chem_integrate( i, j )
+          CALL cpu_log( log_point_s(19), 'chem.reactions', 'stop' )
+
+          IF ( deposition_dry )  THEN
+             CALL cpu_log( log_point_s(24), 'chem.deposition', 'start' )
+             CALL chem_depo( i, j )
+             CALL cpu_log( log_point_s(24), 'chem.deposition', 'stop' )
+          ENDIF
+
+       ENDIF
+
+    ENDIF
+
+ END SUBROUTINE chem_non_transport_physics_ij