Changeset 2801 for palm/trunk/SOURCE/lpm_write_exchange_statistics.f90

Timestamp:

Feb 14, 2018 4:01:55 PM (6 years ago)

Author:

thiele

Message:

Introduce particle transfer in nested models

File:

• palm/trunk/SOURCE/lpm_write_exchange_statistics.f90 (modified) (5 diffs)

palm/trunk/SOURCE/lpm_write_exchange_statistics.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Introduce particle transfer in nested models.
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -69 +72 @@
  SUBROUTINE lpm_write_exchange_statistics
  
+    USE MPI
 
     USE control_parameters,                                                    &
@@ -82 +86 @@
                trsp_count_sum, trsp_count_recv_sum
 
+    USE pmc_particle_interface,                                                &
+        ONLY:  pmcp_g_print_number_of_particles
+
     USE pegrid
 
@@ -89 +96 @@
     INTEGER(iwp) :: jp         !<
     INTEGER(iwp) :: kp         !<
+    INTEGER(iwp) :: tot_number_of_particles
+
+
 
 !
@@ -110 +120 @@
                         trsp_count_recv_sum, pnorth, trnp_count_sum,     &
                         trnp_count_recv_sum
-    CALL close_file( 80 )
 #else
     WRITE ( 80, 8000 )  current_timestep_number+1, simulated_time+dt_3d, &
                         number_of_particles
 #endif
+    CALL close_file( 80 )
+
+    IF ( number_of_particles > 0 ) THEN
+        WRITE(9,*) 'number_of_particles ', number_of_particles, current_timestep_number + 1, simulated_time + dt_3d
+    ENDIF
+
+#if defined( __parallel )
+    CALL MPI_ALLREDUCE( number_of_particles, tot_number_of_particles, 1, MPI_INTEGER,      &
+                                    MPI_SUM, comm2d, ierr)
+#else
+    tot_number_of_particles = number_of_particles
+#endif
+
+    CALL pmcp_g_print_number_of_particles (simulated_time+dt_3d, tot_number_of_particles)
 
 !