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

    r4110 r4168  
    2525! -----------------
    2626! $Id$
     27! Replace function get_topography_top_index by topo_top_ind
     28!
     29! 4110 2019-07-22 17:05:21Z suehring
    2730! pass integer flag array as well as boundary flags to WS scalar advection
    2831! routine
     
    254257        ONLY:  advc_flags_s,                                                  &
    255258               nbgp, nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt,    &
     259               topo_top_ind,                                                  &
    256260               wall_flags_0
    257261
     
    272276        ONLY:  bc_h,                                                          &
    273277               bc_v,                                                          &
    274                get_topography_top_index_ji,                                   &
    275278               surf_def_h,                                                    &
    276279               surf_def_v,                                                    &
     
    13651368          DO  i = nxlg, nxrg
    13661369             DO  j = nysg, nyng
    1367                 nz_s_shift = get_topography_top_index_ji( j, i, 's' )
     1370                nz_s_shift = topo_top_ind(j,i,0)
    13681371
    13691372                e(nz_s_shift:nzt+1,j,i)  =  e(0:nzt+1-nz_s_shift,j,i)
     
    13751378             DO  i = nxlg, nxrg
    13761379                DO  j = nysg, nyng
    1377                    nz_s_shift = get_topography_top_index_ji( j, i, 's' )
     1380                   nz_s_shift = topo_top_ind(j,i,0)
    13781381
    13791382                   diss(nz_s_shift:nzt+1,j,i) = diss(0:nzt+1-nz_s_shift,j,i)
     
    13941397             IF ( nxlg <= 0 .AND. nxrg >= 0 .AND.  &
    13951398                  nysg <= 0 .AND. nyng >= 0        )  THEN
    1396                 nz_s_shift_l = get_topography_top_index_ji( 0, 0, 's' )
     1399                nz_s_shift_l = topo_top_ind(0,0,0)
    13971400             ELSE
    13981401                nz_s_shift_l = 0
Note: See TracChangeset for help on using the changeset viewer.