Changeset 229 for palm/trunk

Timestamp:

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

Author:

raasch

Message:

bugfixes concerning particle tails

Location:

palm/trunk

Files:

• DOC/tec/technical_documentation.html (modified) (2 diffs)
• SOURCE/CURRENT_MODIFICATIONS (modified) (1 diff)
• SOURCE/advec_particles.f90 (modified) (3 diffs)
• SOURCE/init_particles.f90 (modified) (2 diffs)

palm/trunk/DOC/tec/technical_documentation.html

@@ -12 +12 @@
 <br><table nosave="" cellpadding="0" cellspacing="0"> <caption>&nbsp; <br> </caption><tbody>
 </tbody><tbody> </tbody> <tbody> <tr nosave=""> <td nosave=""><b>Current
-model version:</b></td> <td><span style="font-weight: bold;">3.6</span></td> </tr>
+model version:</b></td> <td><span style="font-weight: bold;">3.6 (rev 228)</span></td> </tr>
 <tr nosave=""> <td nosave=""><b>Last
 change of this document</b>:&nbsp;</td> <td nosave=""><b>$Id$</b></td> </tr>
@@ -2434 +2434 @@
 init_3d_model, init_dvrp, local_stop, plant_canopy_model, poismg,
 prandtl_fluxes, pres, read_3d_binary, user_interface, wall_fluxes,
-write_3d_binary<br></td></tr><tr><td style="vertical-align: top;">02/02/09</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.6</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">This
+write_3d_binary<br></td></tr><tr><td style="vertical-align: top;">02/02/09</td><td style="vertical-align: top;">SR</td><td style="vertical-align: top;">3.6<br>(r228)</td><td style="vertical-align: top;">N</td><td style="vertical-align: top;">This
 release contains further additions, changes, and bugfixes for the
 SGI-ICE system, as well as important changes like the splitting of the

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, &