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

palm/trunk/SOURCE/bulk_cloud_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Replace function get_topography_top_index by topo_top_ind
+! 
+! 4110 2019-07-22 17:05:21Z suehring
 ! Pass integer flag array as well as boundary flags to WS scalar advection 
 ! routine
@@ -280 +283 @@
         ONLY:  advc_flags_s,                                                   &
                nbgp, nxl, nxlg, nxr, nxrg, nys, nysg, nyn, nyng, nzb, nzt,     &
+               topo_top_ind,                                                   &
                wall_flags_0
 
@@ -291 +295 @@
 
     USE surface_mod,                                                           &
-        ONLY :  bc_h, get_topography_top_index_ji, surf_bulk_cloud_model,      &
+        ONLY :  bc_h,                                                          &
+                surf_bulk_cloud_model,                                         &
                 surf_microphysics_morrison, surf_microphysics_seifert, &
                 surf_def_h, surf_def_v, surf_lsm_h, surf_lsm_v, surf_usm_h, surf_usm_v
@@ -3807 +3812 @@
 !
 !--          Determine vertical index of topography top
-             k_wall = get_topography_top_index_ji( j, i, 's' )
+             k_wall = topo_top_ind(j,i,0)
              DO  k = nzb+1, nzt
 !
@@ -3854 +3859 @@
 !
 !--    Determine vertical index of topography top
-       k_wall = get_topography_top_index_ji( j, i, 's' )
+       k_wall = topo_top_ind(j,i,0)
        DO  k = nzb+1, nzt
 !