Changeset 229 for palm/trunk/SOURCE

Timestamp:

Feb 3, 2009 9:55:14 AM (16 years ago)

Author:

raasch

Message:

bugfixes concerning particle tails

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (1 diff)
• advec_particles.f90 (modified) (3 diffs)
• init_particles.f90 (modified) (2 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -10 +10 @@
 ------
 
+Bugfix: several tail counters are initialized, particle_tail_coordinates is
+only written to file if its third index is > 0, arrays for tails are allocated
+with a minimum size of 10 tails if there is no tail initially (init_particles,
+advec_particles)
+ 

palm/trunk/SOURCE/advec_particles.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Bugfix: several tail counters are initialized, particle_tail_coordinates is
+! only written to file if its third index is > 0
 ! TEST: PRINT statements on unit 9 (commented out)
 !
@@ -159 +161 @@
        WRITE ( 85 )  maximum_number_of_tailpoints, maximum_number_of_tails, &
                      number_of_tails
-       WRITE ( 85 )  particle_tail_coordinates
+       IF ( maximum_number_of_tails > 0 )  THEN
+          WRITE ( 85 )  particle_tail_coordinates
+       ENDIF
        CALL close_file( 85 )
 
@@ -823 +827 @@
        trrp_count_recv   = 0
        trsp_count_recv   = 0
+       trlpt_count_recv  = 0
+       trnpt_count_recv  = 0
+       trrpt_count_recv  = 0
+       trspt_count_recv  = 0
        IF ( use_particle_tails )  THEN
           tail_mask     = .TRUE.
-          deleted_tails = 0
        ENDIF
+       deleted_tails = 0
 
 

palm/trunk/SOURCE/init_particles.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Bugfix: arrays for tails are allocated with a minimum size of 10 tails if
+! there is no tail initially
 !
 ! Former revisions:
@@ -457 +458 @@
           ENDDO
           maximum_number_of_tails = factor * INT( value )
+!
+!--       Create a minimum number of tails in case that there is no tail
+!--       initially (otherwise, index errors will occur when adressing the
+!--       arrays below)
+          IF ( maximum_number_of_tails == 0 )  maximum_number_of_tails = 10
 
           ALLOCATE( particle_tail_coordinates(maximum_number_of_tailpoints,5, &