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

palm/trunk/SOURCE/biometeorology_mod.f90

@@ -27 +27 @@
 ! -----------------
 ! $Id$
+! Replace function get_topography_top_index by topo_top_ind
+! 
+! 4144 2019-08-06 09:11:47Z raasch
 ! relational operators .EQ., .NE., etc. replaced by ==, /=, etc.
 ! 
@@ -178 +181 @@
     USE indices,                                                               &
         ONLY:  nxl, nxr, nys, nyn, nzb, nzt, nys, nyn, nxl, nxr, nxlg, nxrg,   &
-               nysg, nyng
+               nysg, nyng, topo_top_ind
 
     USE kinds  !< Set precision of INTEGER and REAL arrays according to PALM
@@ -192 +195 @@
                radiation_interactions, rad_sw_in,                              &
                rad_sw_out, rad_lw_in, rad_lw_out
-
-    USE surface_mod,                                                           &
-        ONLY:  get_topography_top_index_ji
 
     IMPLICIT NONE
@@ -1557 +1557 @@
              j = mrtbl(iy,l)
              k = mrtbl(iz,l)
-             IF ( k - get_topography_top_index_ji( j, i, 's' ) ==              &
-                   bio_cell_level + 1_iwp)  THEN
+             IF ( k - topo_top_ind(j,i,0) == bio_cell_level + 1_iwp)  THEN
 !
 !-- Averaging was done before, so we can just copy the result here
@@ -1583 +1582 @@
           j = mrtbl(iy,l)
           k = mrtbl(iz,l)
-          IF ( k - get_topography_top_index_ji( j, i, 's' ) ==                 &
-                bio_cell_level + 1_iwp)  THEN
+          IF ( k - topo_top_ind(j,i,0) == bio_cell_level + 1_iwp)  THEN
              IF ( mrt_include_sw )  THEN
                  tmrt_grid(j,i) = ( ( human_absorb * mrtinsw(l) +              &
@@ -1608 +1606 @@
 !------------------------------------------------------------------------------!
  SUBROUTINE bio_get_thermal_index_input_ij( average_input, i, j, ta, vp, ws,   &
-    pair, tmrt )
+                                            pair, tmrt )
 
     IMPLICIT NONE
@@ -1633 +1631 @@
 !-- Determine cell level closest to 1.1m above ground
 !   by making use of truncation due to int cast
-    k = INT( get_topography_top_index_ji(j, i, 's') + bio_cell_level )  !< Vertical cell center closest to 1.1m
+    k = INT( topo_top_ind(j,i,0) + bio_cell_level )  !< Vertical cell center closest to 1.1m
 
 !