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

Timestamp:

Aug 2, 2019 11:06:18 AM (5 years ago)

Author:

monakurppa

Message:

Several changes in the salsa aerosol module:

• Add "salsa_" before each salsa output variable
• Add a possibility to output the number (salsa_N_UFP) and mass concentration (salsa_PM0.1) of ultrafine particles, i.e. particles with a diameter smaller than 100 nm
• Implement aerosol emission mode "parameterized" which is based on the street type (similar to the chemistry module).
• Remove unnecessary nucleation subroutines.
• Add the z-dimension for gaseous emissions to correspond the implementation in the chemistry module

File:

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

palm/trunk/SOURCE/flow_statistics.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Allow profile output for salsa variables.
+! 
+! 4039 2019-06-18 10:32:41Z suehring
 ! Correct conversion to kinematic scalar fluxes in case of pw-scheme and 
 ! statistic regions
@@ -335 +338 @@
                 simulated_time, simulated_time_at_begin,                       &
                 use_subsidence_tendencies, use_surface_fluxes, use_top_fluxes, &
-                ws_scheme_mom, ws_scheme_sca
+                ws_scheme_mom, ws_scheme_sca, salsa, max_pr_salsa
 
     USE cpulog,                                                                &
@@ -1954 +1957 @@
                 ENDIF
              ENDIF
+             IF ( salsa )  THEN
+                IF ( max_pr_cs > 0 )  THEN
+                   sums_l(:,pr_palm+max_pr_user+max_pr_cs+1:pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0) =    &
+                      sums_l(:,pr_palm+max_pr_user+max_pr_cs+1:pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,0) + &
+                      sums_l(:,pr_palm+max_pr_user+max_pr_cs+1:pr_palm+max_pr_user+max_pr_cs+max_pr_salsa,i)
+
+                ENDIF
+             ENDIF
           ENDDO
        ENDIF
@@ -1974 +1985 @@
                 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
@@ -2051 +2071 @@
                                  sums(k, pr_palm+1:pr_palm+max_pr_user+max_pr_cs) / &
                                  ngp_2dh_s_inner(k,sr)
+             ENDDO
+          ENDIF  
+       ENDIF
+
+       IF ( salsa ) THEN
+          IF ( max_pr_salsa > 0 )  THEN
+             DO k = nzb, nzt+1
+                sums(k,pr_palm+max_pr_user+max_pr_cs+1:pr_palm+max_pr_user+max_pr_cs+max_pr_salsa) = &
+                  sums(k,pr_palm+max_pr_user+max_pr_cs+1:pr_palm+max_pr_user+max_pr_cs+max_pr_salsa) &
+                  / ngp_2dh_s_inner(k,sr)
              ENDDO
           ENDIF  
@@ -2202 +2232 @@
           ENDIF
        ENDIF
+
+       IF ( salsa )  THEN
+          IF ( max_pr_salsa > 0 )  THEN    ! salsa profiles
+             hom(:,1,pr_palm+max_pr_user+max_pr_cs+1:pr_palm+max_pr_user+max_pr_cs+max_pr_salsa, sr) = &
+                  sums(:,pr_palm+max_pr_user+max_pr_cs+1:pr_palm+max_pr_user+max_pr_cs+max_pr_salsa)
+          ENDIF
+       ENDIF
 !
 !--    Determine the boundary layer height using two different schemes.