Ignore:
Timestamp:
Oct 8, 2015 7:32:13 AM (6 years ago)
Author:
raasch
Message:

bugfix concerning vertical index calculation for particles in case of ocean runs

File:
1 edited

Legend:

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

    r1683 r1685  
    1919! Current revisions:
    2020! ------------------
    21 !
     21! bugfix concerning vertical index calculation in case of ocean
    2222!
    2323! Former revisions:
     
    5858        ONLY:  deleted_tails, grid_particles, new_tail_id,                     &
    5959               number_of_particles, number_of_tails, offset_ocean_nzt,         &
    60                offset_ocean_nzt_m1, particles, particle_type, prt_count,       &
     60               particles, particle_type, prt_count,                            &
    6161               particle_tail_coordinates, tail_mask, use_particle_tails
    6262
     
    202202
    203203      USE control_parameters,                                                  &
    204           ONLY: dz,  atmos_ocean_sign
     204          ONLY: dz
    205205
    206206      USE indices,                                                             &
     
    214214      IMPLICIT NONE
    215215
    216       INTEGER(iwp), INTENT(IN)    :: ip
    217       INTEGER(iwp), INTENT(IN)    :: jp
    218       INTEGER(iwp), INTENT(IN)    :: kp
    219 
    220       INTEGER(iwp)                :: i
    221       INTEGER(iwp)                :: j
    222       INTEGER(iwp)                :: k
    223       INTEGER(iwp)                :: n
    224       INTEGER(iwp)                :: nn
    225       INTEGER(iwp)                :: m
    226       INTEGER(iwp)                :: sort_index
    227       INTEGER(iwp)                :: is
    228       INTEGER(iwp)                :: kk
    229 
    230       INTEGER(iwp),DIMENSION(0:7) :: sort_count
     216      INTEGER(iwp), INTENT(IN) :: ip
     217      INTEGER(iwp), INTENT(IN) :: jp
     218      INTEGER(iwp), INTENT(IN) :: kp
     219
     220      INTEGER(iwp)             :: i
     221      INTEGER(iwp)             :: is
     222      INTEGER(iwp)             :: j
     223      INTEGER(iwp)             :: k
     224      INTEGER(iwp)             :: n
     225      INTEGER(iwp)             :: nn
     226      INTEGER(iwp)             :: m
     227      INTEGER(iwp)             :: sort_index
     228
     229      INTEGER(iwp), DIMENSION(0:7) :: sort_count
    231230
    232231      TYPE(particle_type), DIMENSION(number_of_particles,0:7) :: sort_particles
     
    242241             i = particles(n)%x * ddx
    243242             j = particles(n)%y * ddy
    244              k = ( particles(n)%z + 0.5_wp * dz * atmos_ocean_sign ) / dz +    &
    245                  offset_ocean_nzt
    246              kk= particles(n)%z / dz + 1 + offset_ocean_nzt_m1
     243             k = ( particles(n)%z + 0.5_wp * dz ) / dz + offset_ocean_nzt
    247244             IF ( i == ip )  sort_index = sort_index+4
    248245             IF ( j == jp )  sort_index = sort_index+2
Note: See TracChangeset for help on using the changeset viewer.