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

Timestamp:

Apr 11, 2014 5:15:14 PM (10 years ago)

Author:

hoffmann

Message:

new Lagrangian particle structure integrated

File:

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

palm/trunk/SOURCE/lpm_write_exchange_statistics.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! New particle structure integrated. 
 !
 ! Former revisions:
@@ -47 +47 @@
 !------------------------------------------------------------------------------!
 
-    USE control_parameters,                                                     &
+    USE control_parameters,                                                    &
         ONLY:  current_timestep_number, dt_3d, simulated_time
 
-    USE particle_attributes,                                                    &
-        ONLY:  maximum_number_of_particles, number_of_particles, trlp_count_sum,&
-               trlp_count_recv_sum, trnp_count_sum, trnp_count_recv_sum,        &
-               trrp_count_sum, trrp_count_recv_sum, trsp_count_sum,             &
-               trsp_count_recv_sum
+    USE indices,                                                               &
+        ONLY:  nxl, nxr, nys, nyn, nzb, nzt
+
+    USE particle_attributes,                                                   &
+        ONLY:  grid_particles, maximum_number_of_particles,                    &
+               number_of_particles, particles, prt_count,                      &
+               trlp_count_sum, trlp_count_recv_sum, trnp_count_sum,            &
+               trnp_count_recv_sum, trrp_count_sum, trrp_count_recv_sum,       &
+               trsp_count_sum, trsp_count_recv_sum
 
     USE pegrid
@@ -60 +64 @@
     IMPLICIT NONE
 
+    INTEGER(iwp) :: ip         !:
+    INTEGER(iwp) :: jp         !:
+    INTEGER(iwp) :: kp         !:
+
+!
+!-- Determine maximum number of particles (i.e., all possible particles that 
+!-- have been allocated) and the current number of particles
+    number_of_particles         = 0
+    maximum_number_of_particles = 0
+    DO  ip = nxl, nxr
+       DO  jp = nys, nyn
+          DO  kp = nzb+1, nzt
+             number_of_particles = number_of_particles                         &
+                                     + prt_count(kp,jp,ip)
+             maximum_number_of_particles = maximum_number_of_particles         &
+                                     + SIZE(grid_particles(kp,jp,ip)%particles)
+          ENDDO
+       ENDDO
+    ENDDO
 
     CALL check_open( 80 )
@@ -77 +100 @@
 !
 !-- Formats
-8000 FORMAT (I6,1X,F7.2,4X,I6,5X,4(I3,1X,I4,'/',I4,2X),6X,I6)
+8000 FORMAT (I6,1X,F7.2,4X,I10,5X,4(I3,1X,I4,'/',I4,2X),6X,I10)