Ignore:
Timestamp:
Aug 16, 2019 1:50:17 PM (5 years ago)
Author:
suehring
Message:

Replace get_topography_top_index functions by pre-calculated arrays in order to save computational resources

File:
1 edited

Legend:

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

    r4145 r4168  
    2525! -----------------
    2626! $Id$
     27! Replace function get_topography_top_index by topo_top_ind
     28!
     29! 4145 2019-08-06 09:55:22Z schwenkel
    2730! Some reformatting
    2831!
     
    227230    USE indices,                                                               &
    228231        ONLY:  nx, nxl, nxlg, nxrg, nxr, ny, nyn, nys, nyng, nysg, nz, nzb,    &
    229                nzb_max, nzt, wall_flags_0,nbgp, ngp_2dh_outer
     232               nzb_max, nzt,nbgp, ngp_2dh_outer,                               &
     233               topo_top_ind,                                                   &
     234               wall_flags_0
    230235
    231236    USE kinds
     
    259264
    260265    USE surface_mod,                                                           &
    261         ONLY:  get_topography_top_index_ji, surf_def_h, surf_lsm_h, surf_usm_h,&
    262                bc_h
     266        ONLY:  bc_h,                                                           &
     267               surf_def_h,                                                     &
     268               surf_lsm_h,                                                     &
     269               surf_usm_h
    263270
    264271#if defined( __parallel )  &&  !defined( __mpifh )
     
    14021409!--                            Determine surface level. Therefore, check for
    14031410!--                            upward-facing wall on w-grid.
    1404                                k_surf = get_topography_top_index_ji( jp, ip, 'w' )
     1411                               k_surf = topo_top_ind(jp,ip,3)
    14051412                               IF ( seed_follows_topography )  THEN
    14061413!
     
    34693476!--          above topography (Prandtl-layer height)
    34703477!--          Determine vertical index of topography top
    3471              k_wall = get_topography_top_index_ji( jp, ip, 's' )
     3478             k_wall = topo_top_ind(jp,ip,0)
    34723479
    34733480             IF ( constant_flux_layer  .AND.  zv(n) - zw(k_wall) < z_p )  THEN
     
    35583565!
    35593566!--          Determine vertical index of topography top
    3560              k_wall = get_topography_top_index_ji( jp,ip, 's' )
     3567             k_wall = topo_top_ind(jp,ip,0)
    35613568
    35623569             IF ( constant_flux_layer  .AND.  zv(n) - zw(k_wall) < z_p )  THEN
Note: See TracChangeset for help on using the changeset viewer.