Changeset 262 for palm/trunk/SOURCE/advec_particles.f90

Timestamp:

Mar 18, 2009 8:32:37 AM (15 years ago)

Author:

raasch

Message:

further updates for dvr output, bugfix in advec_particles concerning particle boundary condition

File:

• palm/trunk/SOURCE/advec_particles.f90 (modified) (2 diffs)

palm/trunk/SOURCE/advec_particles.f90

@@ -5 +5 @@
 ! -----------------
 ! Output of messages replaced by message handling routine
+! Bugfix: error in check, if particles moved further than one subdomain length.
+!         This check must not be applied for newly released particles
 ! Bugfix: several tail counters are initialized, particle_tail_coordinates is
 ! only written to file if its third index is > 0
@@ -3066 +3068 @@
 !
 !--       Stop if particles have moved further than the length of one 
-!--       PE subdomain
-          IF ( ABS(particles(n)%speed_x) >                                  & 
-               ((nxr-nxl+2)*dx)/(particles(n)%age-particles(n)%age_m)  .OR. &
-               ABS(particles(n)%speed_y) >                                  &
-               ((nyn-nys+2)*dy)/(particles(n)%age-particles(n)%age_m) )  THEN
-
-               WRITE( message_string, * )  'particle too fast.  n = ',  n 
-               CALL message( 'advec_particles', 'PA0148', 2, 2, 0, 6, 0 )
+!--       PE subdomain (newly released particles have age = age_m!)
+          IF ( particles(n)%age /= particles(n)%age_m )  THEN
+             IF ( ABS(particles(n)%speed_x) >                                  &
+                  ((nxr-nxl+2)*dx)/(particles(n)%age-particles(n)%age_m)  .OR. &
+                  ABS(particles(n)%speed_y) >                                  &
+                  ((nyn-nys+2)*dy)/(particles(n)%age-particles(n)%age_m) )  THEN
+
+                  WRITE( message_string, * )  'particle too fast.  n = ',  n 
+                  CALL message( 'advec_particles', 'PA0148', 2, 2, 0, 6, 0 )
+             ENDIF
           ENDIF