Ignore:
Timestamp:
Nov 20, 2017 12:40:38 PM (4 years ago)
Author:
schwenkel
Message:

enable particle advection with grid stretching and some formatation changes in lpm

File:
1 edited

Legend:

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

    r2609 r2628  
    2525! -----------------
    2626! $Id$
     27! Enabled particle advection with grid stretching.
     28!
     29! 2609 2017-11-14 14:14:44Z schwenkel
    2730! Integrated subroutine pack_and_sort into lpm_sort_in_subboxes
    2831!
     
    112115
    113116       USE cpulog,                                                              &
    114            ONLY:  cpu_log, log_point_s
     117          ONLY:  cpu_log, log_point_s
    115118
    116119       USE indices,                                                             &
    117            ONLY:  nxl, nxr, nys, nyn, nzb, nzt
     120          ONLY:  nxl, nxr, nys, nyn, nzb, nzt
    118121
    119122       USE kinds
    120123
    121124       USE control_parameters,                                                  &
    122            ONLY: dz
     125          ONLY: dz
    123126       
    124127       USE grid_variables,                                                      &
    125            ONLY: dx,dy,ddx, ddy
    126          
     128          ONLY: dx,dy,ddx, ddy
     129           
     130       USE arrays_3d,                                                           &
     131          ONLY:  zu
    127132       IMPLICIT NONE
    128133
     
    167172                      i = (particles(n)%x + 0.5_wp * dx) * ddx
    168173                      j = (particles(n)%y + 0.5_wp * dy) * ddy
    169                       k = ( particles(n)%z+ 0.5_wp *dz ) / dz+1 + offset_ocean_nzt
    170174                     
    171175                      IF ( i == ip )  sort_index = sort_index + 4
    172                       IF ( j == jp )  sort_index = sort_index + 2
    173                       IF ( k == kp ) sort_index = sort_index + 1
     176                      IF ( j == jp )  sort_index = sort_index + 2                     
     177                      IF ( zu(kp) > particles(n)%z ) sort_index = sort_index + 1
    174178                     
    175179                      sort_count(sort_index) = sort_count(sort_index) + 1
     
    343347
    344348
    345  END module lpm_pack_and_sort_mod
     349 END MODULE lpm_pack_and_sort_mod
Note: See TracChangeset for help on using the changeset viewer.