 Timestamp:
 Apr 13, 2018 11:55:09 AM (4 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/lpm_advec.f90
r2886 r2969 25 25 !  26 26 ! $Id$ 27 ! Bugfix in Interpolation indices. 28 ! 29 ! 2886 20180314 11:51:53Z thiele 27 30 ! Bugfix in passive particle SGS Model: 28 31 ! Sometimes the added SGS velocities would lead to a violation of the CFL … … 388 391 ELSE 389 392 390 x = xv(n) + ( 0.5_wp  i )* dx391 y = yv(n)  j* dy393 x = xv(n)  i * dx 394 y = yv(n) + ( 0.5_wp  j ) * dy 392 395 aa = x**2 + y**2 393 396 bb = ( dx  x )**2 + y**2 … … 482 485 483 486 ELSE 484 x = xv(n)  i* dx485 y = yv(n) + ( 0.5_wp  j )* dy487 x = xv(n) + ( 0.5_wp  i ) * dx 488 y = yv(n)  j * dy 486 489 aa = x**2 + y**2 487 490 bb = ( dx  x )**2 + y**2 … … 517 520 IF ( vertical_particle_advection(particles(n)%group) ) THEN 518 521 519 x = xv(n)  i* dx520 y = yv(n)  j* dy522 x = xv(n) + ( 0.5_wp  i ) * dx 523 y = yv(n) + ( 0.5_wp  j ) * dy 521 524 aa = x**2 + y**2 522 525 bb = ( dx  x )**2 + y**2 … … 588 591 ! 589 592 ! Interpolate TKE 590 x = xv(n)  i* dx591 y = yv(n)  j* dy593 x = xv(n) + ( 0.5_wp  i ) * dx 594 y = yv(n) + ( 0.5_wp  j ) * dy 592 595 aa = x**2 + y**2 593 596 bb = ( dx  x )**2 + y**2
Note: See TracChangeset
for help on using the changeset viewer.