Changeset 4114 for palm/trunk/SOURCE


Ignore:
Timestamp:
Jul 23, 2019 2:09:27 PM (5 years ago)
Author:
schwenkel
Message:

Bugfix: Added working precision

File:
1 edited

Legend:

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

    r4054 r4114  
    2525! -----------------
    2626! $Id$
     27! Bugfix: Added working precision for if statement
     28!
     29! 4054 2019-06-27 07:42:18Z raasch
    2730! bugfix for calculating the minimum particle time step
    2831!
     
    13381341!
    13391342!--                            In case of stretching the actual k index is found iteratively
    1340                                IF ( dz_stretch_level .NE. -9999999.9  .OR.           &
    1341                                     dz_stretch_level_start(1) .NE. -9999999.9 ) THEN
     1343                               IF ( dz_stretch_level .NE. -9999999.9_wp  .OR.           &
     1344                                    dz_stretch_level_start(1) .NE. -9999999.9_wp ) THEN
    13421345                                  kp = MINLOC( ABS( tmp_particle%z - zu ), DIM = 1 ) - 1
    13431346                               ELSE
     
    73527355!
    73537356!--    In case of stretching the actual k index must be found
    7354        IF ( dz_stretch_level .NE. -9999999.9  .OR.         &
    7355             dz_stretch_level_start(1) .NE. -9999999.9 ) THEN
     7357       IF ( dz_stretch_level .NE. -9999999.9_wp  .OR.         &
     7358            dz_stretch_level_start(1) .NE. -9999999.9_wp ) THEN
    73567359          kp = MINLOC( ABS( particle_array(n)%z - zu ), DIM = 1 ) - 1
    73577360       ELSE
     
    75237526                IF( zw(k)   < particles_before_move(n)%z ) k = k + 1
    75247527                IF( zw(k-1) > particles_before_move(n)%z ) k = k - 1
    7525                
     7528
    75267529!--             For lpm_exchange_horiz to work properly particles need to be moved to the outermost gridboxes
    75277530!--             of the respective processor. If the particle index is inside the processor the following lines
     
    75317534                j = MIN ( j , nyn )
    75327535                j = MAX ( j , nys )
    7533                
     7536
    75347537                k = MIN ( k , nzt )
    75357538                k = MAX ( k , nzb+1 )
    7536                
     7539
    75377540!
    75387541!--             Check, if particle has moved to another grid cell.
Note: See TracChangeset for help on using the changeset viewer.