Changeset 4581 for palm/trunk/SOURCE/flow_statistics.f90

Timestamp:

Jun 29, 2020 8:49:58 AM (4 years ago)

Author:

suehring

Message:

mesoscale nesting: omit explicit pressure forcing via geostrophic wind components; chemistry: enable profile output of vertical fluxes; urban-surface: bugfix in initialization in case of cyclic_fill

File:

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

palm/trunk/SOURCE/flow_statistics.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Formatting adjustment
+! 
+! 4551 2020-06-02 10:22:25Z suehring
 ! Bugfix in summation for statistical regions
 ! 
@@ -1883 +1886 @@
        IF ( air_chemistry  .AND.  max_pr_cs > 0 )  THEN
           IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
-             DO  i = 1, max_pr_cs
-                CALL MPI_ALLREDUCE( sums_l(nzb,pr_palm+max_pr_user+i,0),       &
-                                    sums(nzb,pr_palm+max_pr_user+i),           &
-                                    nzt+2-nzb, MPI_REAL, MPI_SUM, comm2d, ierr )
-             ENDDO
+          DO  i = 1, max_pr_cs
+             CALL MPI_ALLREDUCE( sums_l(nzb,pr_palm+max_pr_user+i,0),          &
+                                 sums(nzb,pr_palm+max_pr_user+i),              &
+                                 nzt+2-nzb, MPI_REAL, MPI_SUM, comm2d, ierr )
+          ENDDO
        ENDIF
 
        IF ( salsa  .AND.  max_pr_salsa > 0 )  THEN
           IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
-             DO  i = 1, max_pr_salsa
-                CALL MPI_ALLREDUCE( sums_l(nzb,pr_palm+max_pr_user+max_pr_cs+i,0),                 &
-                                    sums(nzb,pr_palm+max_pr_user+max_pr_user+i),                   &
-                                    nzt+2-nzb, MPI_REAL, MPI_SUM, comm2d, ierr )
-             ENDDO
+          DO  i = 1, max_pr_salsa
+             CALL MPI_ALLREDUCE( sums_l(nzb,pr_palm+max_pr_user+max_pr_cs+i,0),                    &
+                                 sums(nzb,pr_palm+max_pr_user+max_pr_user+i),                      &
+                                 nzt+2-nzb, MPI_REAL, MPI_SUM, comm2d, ierr )
+          ENDDO
        ENDIF