Ignore:
Timestamp:
Aug 16, 2019 1:50:17 PM (2 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/multi_agent_system_mod.f90

    r3987 r4168  
    2525! -----------------
    2626! $Id$
     27! Replace function get_topography_top_index by topo_top_ind
     28!
     29! 3987 2019-05-22 09:52:13Z kanani
    2730! Introduce alternative switch for debug output during timestepping
    2831!
     
    131134    USE indices,                                                               &
    132135        ONLY:  nx, nxl, nxlg, nxr, nxrg, ny, nyn, nyng, nys, nysg, nzb,        &
     136               topo_top_ind,                                                   &
    133137               wall_flags_0
    134138
     
    374378 SUBROUTINE multi_agent_system
    375379
    376     USE biometeorology_mod,                                                   &
    377         ONLY:  bio_calc_ipt, bio_calculate_mrt_grid, bio_get_thermal_index_input_ij
     380    USE biometeorology_mod,                                                    &
     381        ONLY:  bio_calc_ipt,                                                   &
     382               bio_calculate_mrt_grid,                                         &
     383               bio_get_thermal_index_input_ij
    378384
    379385
     
    30183024           ONLY:  coupling_char, initializing_actions, io_blocks, io_group
    30193025
    3020        USE surface_mod,                                                        &
    3021            ONLY:  get_topography_top_index, get_topography_top_index_ji
    3022 
    30233026       USE arrays_3d,                                                          &
    30243027           ONLY:  zu, zw
     
    30773080       DO il = nxlg, nxrg
    30783081          DO jl = nysg, nyng
    3079              top_top_s(jl,il) = get_topography_top_index_ji(jl,il,'s') + 1
    3080              top_top_w(jl,il) = get_topography_top_index_ji(jl,il,'w')
     3082             top_top_s(jl,il) = topo_top_ind(jl,il,0) + 1
     3083             top_top_w(jl,il) = topo_top_ind(jl,il,3)
    30813084          ENDDO
    30823085       ENDDO
Note: See TracChangeset for help on using the changeset viewer.