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

palm/trunk/SOURCE/init_3d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Replace function get_topography_top_index by topo_top_ind
+! 
+! 4151 2019-08-09 08:24:30Z suehring
 ! Add netcdf directive around input calls (fix for last commit)
 ! 
@@ -719 +722 @@
                 surf_def_v,                                                    &
                 surf_lsm_h,                                                    &
-                surf_usm_h,                                                    &
-                get_topography_top_index_ji
+                surf_usm_h
 
 #if defined( _OPENACC )
@@ -1633 +1635 @@
           DO  i = nxlg, nxrg
              DO  j = nysg, nyng
-                nz_u_shift = get_topography_top_index_ji( j, i, 'u' )
-                nz_v_shift = get_topography_top_index_ji( j, i, 'v' )
-                nz_w_shift = get_topography_top_index_ji( j, i, 'w' )
-                nz_s_shift = get_topography_top_index_ji( j, i, 's' )
+                nz_u_shift = topo_top_ind(j,i,1)
+                nz_v_shift = topo_top_ind(j,i,2)
+                nz_w_shift = topo_top_ind(j,i,3)
+                nz_s_shift = topo_top_ind(j,i,0)
 
                 u(nz_u_shift:nzt+1,j,i)  = u(0:nzt+1-nz_u_shift,j,i)                
@@ -1677 +1679 @@
           IF ( complex_terrain )  THEN
              IF ( nxlg <= 0 .AND. nxrg >= 0 .AND. nysg <= 0 .AND. nyng >= 0 )  THEN
-                nz_u_shift_l = get_topography_top_index_ji( 0, 0, 'u' )
-                nz_v_shift_l = get_topography_top_index_ji( 0, 0, 'v' )
-                nz_w_shift_l = get_topography_top_index_ji( 0, 0, 'w' )
-                nz_s_shift_l = get_topography_top_index_ji( 0, 0, 's' )
+                nz_u_shift_l = topo_top_ind(j,i,1)
+                nz_v_shift_l = topo_top_ind(j,i,2)
+                nz_w_shift_l = topo_top_ind(j,i,3)
+                nz_s_shift_l = topo_top_ind(j,i,0)
              ELSE
                 nz_u_shift_l = 0