Changeset 4168 for palm/trunk/SOURCE/multi_agent_system_mod.f90

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:

• palm/trunk/SOURCE/multi_agent_system_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/multi_agent_system_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Replace function get_topography_top_index by topo_top_ind
+! 
+! 3987 2019-05-22 09:52:13Z kanani
 ! Introduce alternative switch for debug output during timestepping
 ! 
@@ -131 +134 @@
     USE indices,                                                               &
         ONLY:  nx, nxl, nxlg, nxr, nxrg, ny, nyn, nyng, nys, nysg, nzb,        &
+               topo_top_ind,                                                   &
                wall_flags_0
 
@@ -374 +378 @@
  SUBROUTINE multi_agent_system
 
-    USE biometeorology_mod,                                                   &
-        ONLY:  bio_calc_ipt, bio_calculate_mrt_grid, bio_get_thermal_index_input_ij
+    USE biometeorology_mod,                                                    &
+        ONLY:  bio_calc_ipt,                                                   &
+               bio_calculate_mrt_grid,                                         &
+               bio_get_thermal_index_input_ij
 
 
@@ -3018 +3024 @@
            ONLY:  coupling_char, initializing_actions, io_blocks, io_group
 
-       USE surface_mod,                                                        &
-           ONLY:  get_topography_top_index, get_topography_top_index_ji
-
        USE arrays_3d,                                                          &
            ONLY:  zu, zw
@@ -3077 +3080 @@
        DO il = nxlg, nxrg
           DO jl = nysg, nyng
-             top_top_s(jl,il) = get_topography_top_index_ji(jl,il,'s') + 1
-             top_top_w(jl,il) = get_topography_top_index_ji(jl,il,'w')
+             top_top_s(jl,il) = topo_top_ind(jl,il,0) + 1
+             top_top_w(jl,il) = topo_top_ind(jl,il,3)
           ENDDO
        ENDDO