Changeset 1822 for palm/trunk/SOURCE/lpm_write_restart_file.f90

Timestamp:

Apr 7, 2016 7:49:42 AM (8 years ago)

Author:

hoffmann

Message:

changes in LPM and bulk cloud microphysics

File:

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

palm/trunk/SOURCE/lpm_write_restart_file.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! Tails removed. Unused variables removed.
 ! 
 ! Former revisions:
@@ -50 +50 @@
 !------------------------------------------------------------------------------!
  SUBROUTINE lpm_write_restart_file
- 
 
-    USE control_parameters,                                                    &
-        ONLY:  io_blocks, io_group
 
     USE indices,                                                               &
@@ -62 +59 @@
     USE particle_attributes,                                                   &
         ONLY:  bc_par_b, bc_par_lr, bc_par_ns, bc_par_t, grid_particles,       &
-               maximum_number_of_tails, maximum_number_of_tailpoints,          &
                number_of_particles, number_of_particle_groups,                 &
-               number_of_tails, particles, particle_groups,                    &
-               particle_tail_coordinates, prt_count, time_prel,                &
-               time_write_particle_data, uniform_particles,                    &
-               use_particle_tails, zero_particle
+               particles, particle_groups, prt_count, time_prel,               &
+               time_write_particle_data, uniform_particles
 
     USE pegrid
@@ -75 +69 @@
     CHARACTER (LEN=10) ::  particle_binary_version   !<
 
-    INTEGER(iwp) ::  i                               !<
     INTEGER(iwp) ::  ip                              !<
     INTEGER(iwp) ::  jp                              !<
@@ -97 +90 @@
     ENDIF
 
-!    DO  i = 0, io_blocks-1
-
-!       IF ( i == io_group )  THEN
+!
+!-- Write the version number of the binary format.
+!-- Attention: After changes to the following output commands the version
+!-- ---------  number of the variable particle_binary_version must be
+!--            changed! Also, the version number and the list of arrays
+!--            to be read in lpm_read_restart_file must be adjusted
+!--            accordingly.
+    particle_binary_version = '4.0'
+    WRITE ( 90 )  particle_binary_version
 
 !
-!--       Write the version number of the binary format.
-!--       Attention: After changes to the following output commands the version
-!--       ---------  number of the variable particle_binary_version must be
-!--                  changed! Also, the version number and the list of arrays
-!--                  to be read in lpm_read_restart_file must be adjusted
-!--                  accordingly.
-          particle_binary_version = '3.2'
-          WRITE ( 90 )  particle_binary_version
+!-- Write some particle parameters, the size of the particle arrays as
+!-- well as other dvrp-plot variables.
+    WRITE ( 90 )  bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,              &
+                  number_of_particle_groups, particle_groups, time_prel, &
+                  time_write_particle_data, uniform_particles
 
-!
-!--       Write some particle parameters, the size of the particle arrays as
-!--       well as other dvrp-plot variables.
-          WRITE ( 90 )  bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,              &
-                        maximum_number_of_tailpoints, maximum_number_of_tails, &
-                        number_of_particle_groups, number_of_tails,            &
-                        particle_groups, time_prel, time_write_particle_data,  &
-                        uniform_particles
+    WRITE ( 90 )  prt_count
+          
+    DO  ip = nxl, nxr
+       DO  jp = nys, nyn
+          DO  kp = nzb+1, nzt
+             number_of_particles = prt_count(kp,jp,ip)
+             particles => grid_particles(kp,jp,ip)%particles(1:number_of_particles)
+             IF ( number_of_particles <= 0 )  CYCLE
+             WRITE ( 90 )  particles
+          ENDDO
+       ENDDO
+    ENDDO
 
-          WRITE ( 90 )  prt_count
-          
-          DO  ip = nxl, nxr
-             DO  jp = nys, nyn
-                DO  kp = nzb+1, nzt
-                   number_of_particles = prt_count(kp,jp,ip)
-                   particles => grid_particles(kp,jp,ip)%particles(1:number_of_particles)
-                   IF ( number_of_particles <= 0 )  CYCLE
-                   WRITE ( 90 )  particles
-                ENDDO
-             ENDDO
-          ENDDO
-
-!
-!--       particle tails currently not available
-!          IF ( use_particle_tails )  THEN
-!             WRITE ( 90 )  particle_tail_coordinates
-!          ENDIF
-
-          CLOSE ( 90 )
-
-!       ENDIF
+    CLOSE ( 90 )
 
 #if defined( __parallel )
@@ -147 +126 @@
 #endif
 
- !   ENDDO
-
 
  END SUBROUTINE lpm_write_restart_file