Changeset 4628 for palm/trunk/SOURCE/mod_particle_attributes.f90

Timestamp:

Jul 29, 2020 7:23:03 AM (4 years ago)

Author:

raasch

Message:

extensions required for MPI-I/O of particle data to restart files

File:

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

palm/trunk/SOURCE/mod_particle_attributes.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-! 
-! 
+!
+!
 ! Former revisions:
 ! -----------------
 ! $Id$
+! extensions required for MPI-I/O of particle data to restart files
+! 
+! 4360 2020-01-07 11:25:50Z suehring
 ! Corrected "Former revisions" section
 ! 
@@ -47 +50 @@
     USE, INTRINSIC ::  ISO_C_BINDING
 
+    USE control_parameters,                                                                        &
+        ONLY: varnamelength
+
     USE kinds
+
+    CHARACTER(LEN=varnamelength), DIMENSION(50) ::  part_output = ' '  !< namelist parameter
 
     INTEGER(iwp) ::  dissipation_classes = 10                     !< namelist parameter (see documentation)
@@ -54 +62 @@
     INTEGER(iwp) ::  ibc_par_ns                                   !< particle north/south boundary condition dummy
     INTEGER(iwp) ::  ibc_par_t                                    !< particle top boundary condition dummy
+    INTEGER(iwp) ::  number_of_output_particles = 0               !< number of output particles
     INTEGER(iwp) ::  number_of_particles = 0                      !< number of particles for each grid box (3d array is saved on prt_count)
     INTEGER(iwp) ::  number_of_particle_groups = 1                !< namelist parameter (see documentation)
+    INTEGER(iwp) ::  part_inc = 1                                 !< increment of particles in output file
 
     INTEGER(iwp), PARAMETER ::  max_number_of_particle_groups = 10 !< maximum allowed number of particle groups
@@ -62 +72 @@
     
     LOGICAL ::  particle_advection = .FALSE.              !< parameter to steer the advection of particles
+    LOGICAL ::  unlimited_dimension = .TRUE.              !< umlimited dimension for particle output
     LOGICAL ::  use_sgs_for_particles = .FALSE.           !< namelist parameter (see documentation)    
     LOGICAL ::  wang_kernel = .FALSE.                     !< flag for collision kernel
 
-    REAL(wp) ::  alloc_factor = 20.0_wp                    !< namelist parameter (see documentation)
-    REAL(wp) ::  particle_advection_start = 0.0_wp         !< namelist parameter (see documentation)
+    REAL(wp) ::  alloc_factor = 20.0_wp                   !< namelist parameter (see documentation)
+    REAL(wp) ::  oversize = 100.0_wp                      !< reserve spare particles in output file (in % relative to initial number)
+    REAL(wp) ::  particle_advection_start = 0.0_wp        !< namelist parameter (see documentation)
+    REAL(wp) ::  part_percent = 100.0_wp                  !< percentage of particles in output file
 
-!
-!-- WARNING: For compatibility of derived types, the BIND attribute is required, and interoperable C
-!-- datatypes must be used. These type are hard wired here! So changes in working precision (wp, iwp)
-!-- will not affect the particle_type!
-!-- The main reason for introducing the interoperable datatypes was to avoid compiler warnings of
-!-- the gfortran compiler.
-!-- The BIND attribite is required because of C_F_POINTER usage in the pmc particle interface.
-    TYPE, BIND(C) ::  particle_type
-        REAL(C_DOUBLE) ::  aux1          !< auxiliary multi-purpose feature
-        REAL(C_DOUBLE) ::  aux2          !< auxiliary multi-purpose feature
-        REAL(C_DOUBLE) ::  radius        !< radius of particle
-        REAL(C_DOUBLE) ::  age           !< age of particle
-        REAL(C_DOUBLE) ::  age_m         !<
-        REAL(C_DOUBLE) ::  dt_sum        !<
-        REAL(C_DOUBLE) ::  e_m           !< interpolated sgs tke
-        REAL(C_DOUBLE) ::  origin_x      !< origin x-position of particle (changed cyclic bc)
-        REAL(C_DOUBLE) ::  origin_y      !< origin y-position of particle (changed cyclic bc)
-        REAL(C_DOUBLE) ::  origin_z      !< origin z-position of particle (changed cyclic bc)
-        REAL(C_DOUBLE) ::  rvar1         !<
-        REAL(C_DOUBLE) ::  rvar2         !<
-        REAL(C_DOUBLE) ::  rvar3         !<
-        REAL(C_DOUBLE) ::  speed_x       !< speed of particle in x
-        REAL(C_DOUBLE) ::  speed_y       !< speed of particle in y
-        REAL(C_DOUBLE) ::  speed_z       !< speed of particle in z
-        REAL(C_DOUBLE) ::  weight_factor !< weighting factor
-        REAL(C_DOUBLE) ::  x             !< x-position
-        REAL(C_DOUBLE) ::  y             !< y-position
-        REAL(C_DOUBLE) ::  z             !< z-position
-        INTEGER(C_INT) ::  class         !< radius class needed for collision
-        INTEGER(C_INT) ::  group         !< number of particle group
-        INTEGER(C_LONG_LONG) ::  id            !< particle ID (64 bit integer)
-        LOGICAL(C_BOOL) ::  particle_mask !< if this parameter is set to false the particle will be deleted
-        INTEGER(C_INT) ::  block_nr      !< number for sorting (removable?)
+    TYPE, PUBLIC ::  particle_type
+        REAL(wp)     ::  aux1          !< auxiliary multi-purpose feature
+        REAL(wp)     ::  aux2          !< auxiliary multi-purpose feature
+        REAL(wp)     ::  radius        !< radius of particle
+        REAL(wp)     ::  age           !< age of particle
+        REAL(wp)     ::  age_m         !<
+        REAL(wp)     ::  dt_sum        !<
+        REAL(wp)     ::  e_m           !< interpolated sgs tke
+        REAL(wp)     ::  origin_x      !< origin x-position of particle (changed cyclic bc)
+        REAL(wp)     ::  origin_y      !< origin y-position of particle (changed cyclic bc)
+        REAL(wp)     ::  origin_z      !< origin z-position of particle (changed cyclic bc)
+        REAL(wp)     ::  rvar1         !<
+        REAL(wp)     ::  rvar2         !<
+        REAL(wp)     ::  rvar3         !<
+        REAL(wp)     ::  speed_x       !< speed of particle in x
+        REAL(wp)     ::  speed_y       !< speed of particle in y
+        REAL(wp)     ::  speed_z       !< speed of particle in z
+        REAL(wp)     ::  weight_factor !< weighting factor
+        REAL(wp)     ::  x             !< x-position
+        REAL(wp)     ::  y             !< y-position
+        REAL(wp)     ::  z             !< z-position
+        INTEGER(iwp) ::  class         !< radius class needed for collision
+        INTEGER(iwp) ::  group         !< number of particle group
+        INTEGER(idp) ::  id            !< particle ID (64 bit integer)
+        LOGICAL      ::  particle_mask !< if this parameter is set to false the particle will be deleted
+        INTEGER(iwp) ::  block_nr      !< number for sorting (removable?)
+        INTEGER(iwp) ::  particle_nr=-1   !< particle number for particle IO (increment one
     END TYPE particle_type