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

Timestamp:

Aug 16, 2019 1:50:17 PM (4 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/nesting_offl_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/nesting_offl_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Replace function get_topography_top_index by topo_top_ind
+! 
+! 4125 2019-07-29 13:31:44Z suehring
 ! In order to enable netcdf parallel access, allocate dummy arrays for the 
 ! lateral boundary data on cores that actually do not belong to these 
@@ -119 +122 @@
     USE indices,                                                               &
         ONLY:  nbgp, nx, nxl, nxlg, nxlu, nxr, nxrg, ny, nys,                  &
-               nysv, nysg, nyn, nyng, nzb, nz, nzt, wall_flags_0
+               nysv, nysg, nyn, nyng, nzb, nz, nzt,                            &
+               topo_top_ind,                                                   &
+               wall_flags_0
 
     USE kinds
@@ -873 +878 @@
        
        USE kinds
-       
-       USE surface_mod,                                                        &
-           ONLY:  get_topography_top_index, get_topography_top_index_ji
 
        IMPLICIT NONE 
@@ -918 +920 @@
 !
 !--          Determine topography top index at current (j,i) index
-             k_surface = get_topography_top_index_ji( j, i, 's' )
+             k_surface = topo_top_ind(j,i,0)
 !
 !--          Pre-compute surface virtual temperature. Therefore, use 2nd 
@@ -978 +980 @@
        
           DO  i = nxl, nxr
-             k_surface = get_topography_top_index_ji( j, i, 's' )
+             k_surface = topo_top_ind(j,i,0)
  
              IF ( humidity )  THEN
@@ -1032 +1034 @@
 !--    turbulence generator accordingly. If Rayleigh damping would be applied
 !--    near buildings, etc., this would spoil the simulation results. 
-       topo_max_l = zw(MAXVAL( get_topography_top_index( 's' )))
+       topo_max_l = zw(MAXVAL( topo_top_ind(nys:nyn,nxl:nxr,0) ))
        
 #if defined( __parallel )