Changeset 4168 for palm/trunk/SOURCE/pmc_interface_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/pmc_interface_mod.f90 (modified) (7 diffs)

palm/trunk/SOURCE/pmc_interface_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Replace function get_topography_top_index by topo_top_ind
+! 
+! 4029 2019-06-14 14:04:35Z raasch
 ! nest_chemistry switch removed
 ! 
@@ -478 +481 @@
     USE indices,                                                               &
         ONLY:  nbgp, nx, nxl, nxlg, nxlu, nxr, nxrg, ny, nyn, nyng, nys, nysg, &
-               nysv, nz, nzb, nzt, wall_flags_0
+               nysv, nz, nzb, nzt, topo_top_ind, wall_flags_0
 
     USE bulk_cloud_model_mod,                                                  &
@@ -528 +531 @@
 
     USE surface_mod,                                                           &
-        ONLY:  get_topography_top_index_ji, surf_def_h, surf_lsm_h, surf_usm_h
+        ONLY:  surf_def_h, surf_lsm_h, surf_usm_h
 
     IMPLICIT NONE
@@ -2379 +2382 @@
           DO  j = nys, nyn
              sub_sum = 0.0_wp
-             k_wall = get_topography_top_index_ji( j, i, 'u' )
+             k_wall = topo_top_ind(j,i,1)
              DO   k = k_wall + 1, nzt
                 sub_sum = sub_sum + dzw(k)
@@ -2401 +2404 @@
           DO  j = nys, nyn
              sub_sum = 0.0_wp
-             k_wall = get_topography_top_index_ji( j, i, 'u' )
+             k_wall = topo_top_ind(j,i,1)
              DO   k = k_wall + 1, nzt
                 sub_sum = sub_sum + dzw(k)
@@ -2423 +2426 @@
           DO  i = nxl, nxr
              sub_sum = 0.0_wp
-             k_wall = get_topography_top_index_ji( j, i, 'v' )
+             k_wall = topo_top_ind(j,i,2)
              DO  k = k_wall + 1, nzt
                 sub_sum = sub_sum + dzw(k)
@@ -2445 +2448 @@
           DO  i = nxl, nxr
              sub_sum = 0.0_wp
-             k_wall = get_topography_top_index_ji( j, i, 'v' )
+             k_wall = topo_top_ind(j,i,2)
              DO  k = k_wall + 1, nzt
                 sub_sum = sub_sum + dzw(k)