Changeset 3797 for palm/trunk/SOURCE/prognostic_equations.f90

Timestamp:

Mar 15, 2019 11:15:38 AM (5 years ago)

Author:

forkel

Message:

Modifications for OpenMP version by Klaus Ketelsen

File:

• palm/trunk/SOURCE/prognostic_equations.f90 (modified) (2 diffs)

palm/trunk/SOURCE/prognostic_equations.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Call chem_integegrate in OpenMP loop   (ketelsen)
+! 
+! 
+! 3771 2019-02-28 12:19:33Z raasch
 ! preprocessor directivs fro rrtmg added
 ! 
@@ -527 +531 @@
              CALL photolysis_control
           ENDIF
-          DO  i = nxl, nxr
-             DO  j = nys, nyn
-
-                IF ( intermediate_timestep_count == 1 .OR.                        &
-                     call_chem_at_all_substeps )  THEN
-
-                   CALL cpu_log( log_point_s(19), 'chem.reactions', 'start' )  
+
+          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)
-                   CALL cpu_log( log_point_s(19), 'chem.reactions', 'stop' )
-
-                   IF ( do_depo )  THEN
-                      CALL cpu_log( log_point_s(24), 'chem.deposition', 'start' )
+                ENDDO
+             ENDDO
+             !$OMP END PARALLEL
+             CALL cpu_log( log_point_s(19), 'chem.reactions', 'stop' )
+
+             IF ( do_depo )  THEN
+                CALL cpu_log( log_point_s(24), 'chem.deposition', 'start' )
+                DO  i = nxl, nxr
+                   DO  j = nys, nyn
                       CALL chem_depo(i,j)
-                      CALL cpu_log( log_point_s(24), 'chem.deposition', 'stop' )
-                   ENDIF
-                ENDIF
-             ENDDO
-          ENDDO
+                   ENDDO
+                ENDDO
+                CALL cpu_log( log_point_s(24), 'chem.deposition', 'stop' )
+             ENDIF
+          ENDIF
        ENDIF
 !