Ignore:
Timestamp:
Oct 8, 2015 7:32:13 AM (9 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_exchange_horiz.f90

    r1683 r1685  
    1919! Current revisions:
    2020! ------------------
    21 !
     21! bugfix concerning vertical index offset in case of ocean
    2222!
    2323! Former revisions:
     
    8585               maximum_number_of_tailpoints, min_nr_particle,                  &
    8686               mpi_particle_type, number_of_tails, number_of_particles,        &
    87                offset_ocean_nzt, offset_ocean_nzt_m1, particles,               &
     87               offset_ocean_nzt, particles,                                    &
    8888               particle_tail_coordinates, particle_type, prt_count,            &
    8989               tail_mask, trlp_count_sum,                                      &
     
    10791079       ip = ( particle_array(n)%x + 0.5_wp * dx ) * ddx
    10801080       jp = ( particle_array(n)%y + 0.5_wp * dy ) * ddy
    1081        kp = particle_array(n)%z / dz + 1 + offset_ocean_nzt_m1
     1081       kp = particle_array(n)%z / dz + 1 + offset_ocean_nzt
    10821082
    10831083       IF ( ip >= nxl  .AND.  ip <= nxr  .AND.  jp >= nys  .AND.  jp <= nyn    &
Note: See TracChangeset for help on using the changeset viewer.