Ignore:
Timestamp:
Mar 18, 2009 8:32:37 AM (14 years ago)
Author:
raasch
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/advec_particles.f90

    r230 r262  
    55! -----------------
    66! Output of messages replaced by message handling routine
     7! Bugfix: error in check, if particles moved further than one subdomain length.
     8!         This check must not be applied for newly released particles
    79! Bugfix: several tail counters are initialized, particle_tail_coordinates is
    810! only written to file if its third index is > 0
     
    30663068!
    30673069!--       Stop if particles have moved further than the length of one
    3068 !--       PE subdomain
    3069           IF ( ABS(particles(n)%speed_x) >                                  &
    3070                ((nxr-nxl+2)*dx)/(particles(n)%age-particles(n)%age_m)  .OR. &
    3071                ABS(particles(n)%speed_y) >                                  &
    3072                ((nyn-nys+2)*dy)/(particles(n)%age-particles(n)%age_m) )  THEN
    3073 
    3074                WRITE( message_string, * )  'particle too fast.  n = ',  n
    3075                CALL message( 'advec_particles', 'PA0148', 2, 2, 0, 6, 0 )
     3070!--       PE subdomain (newly released particles have age = age_m!)
     3071          IF ( particles(n)%age /= particles(n)%age_m )  THEN
     3072             IF ( ABS(particles(n)%speed_x) >                                  &
     3073                  ((nxr-nxl+2)*dx)/(particles(n)%age-particles(n)%age_m)  .OR. &
     3074                  ABS(particles(n)%speed_y) >                                  &
     3075                  ((nyn-nys+2)*dy)/(particles(n)%age-particles(n)%age_m) )  THEN
     3076
     3077                  WRITE( message_string, * )  'particle too fast.  n = ',  n
     3078                  CALL message( 'advec_particles', 'PA0148', 2, 2, 0, 6, 0 )
     3079             ENDIF
    30763080          ENDIF
    30773081
Note: See TracChangeset for help on using the changeset viewer.