Ignore:
Timestamp:
Aug 16, 2007 2:30:26 PM (17 years ago)
Author:
raasch
Message:

preliminary update of bugfixes and extensions for non-cyclic BCs

File:
1 edited

Legend:

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

    r83 r106  
    44! Actual revisions:
    55! -----------------
    6 !
     6! variable iran replaced by iran_part
    77!
    88! Former revisions:
     
    354354       IF ( random_start_position )  THEN
    355355
    356           iran = iran + myid    ! Random positions should be different on
    357                                 ! different PEs
    358 
    359356          DO  n = 1, number_of_initial_particles
    360357             IF ( psl(particles(n)%group) /= psr(particles(n)%group) )  THEN
    361358                particles(n)%x = particles(n)%x + &
    362                                  ( random_function( iran ) - 0.5 ) * &
     359                                 ( random_function( iran_part ) - 0.5 ) * &
    363360                                 pdx(particles(n)%group)
    364361                IF ( particles(n)%x  <=  ( nxl - 0.5 ) * dx )  THEN
     
    370367             IF ( pss(particles(n)%group) /= psn(particles(n)%group) )  THEN
    371368                particles(n)%y = particles(n)%y + &
    372                                  ( random_function( iran ) - 0.5 ) * &
     369                                 ( random_function( iran_part ) - 0.5 ) * &
    373370                                 pdy(particles(n)%group)
    374371                IF ( particles(n)%y  <=  ( nys - 0.5 ) * dy )  THEN
     
    380377             IF ( psb(particles(n)%group) /= pst(particles(n)%group) )  THEN
    381378                particles(n)%z = particles(n)%z + &
    382                                  ( random_function( iran ) - 0.5 ) * &
     379                                 ( random_function( iran_part ) - 0.5 ) * &
    383380                                 pdz(particles(n)%group)
    384381             ENDIF
Note: See TracChangeset for help on using the changeset viewer.