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

palm/trunk/SOURCE/surface_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Remove functions get_topography_top_index. These are now replaced by 
+! precalculated arrays because of too much CPU-time consumption
+! 
+! 4159 2019-08-15 13:31:35Z suehring
 ! Surface classification revised and adjusted to changes in init_grid
 ! 
@@ -646 +650 @@
 
     PRIVATE
-
+    
     INTERFACE init_bc
        MODULE PROCEDURE init_bc
@@ -698 +702 @@
 !
 !-- Public subroutines and functions
-    PUBLIC get_topography_top_index,                                                               &
-           get_topography_top_index_ji,                                                            &
-           init_bc,                                                                                &
+    PUBLIC init_bc,                                                                                &
            init_single_surface_properties,                                                         &
            init_surfaces,                                                                          &
@@ -710 +712 @@
 
 #if defined( _OPENACC )
-    PUBLIC enter_surface_arrays, exit_surface_arrays
+    PUBLIC enter_surface_arrays,                                                                   &
+           exit_surface_arrays
 #endif
 
@@ -3141 +3144 @@
        
     END SUBROUTINE init_single_surface_properties
-
-!------------------------------------------------------------------------------!
-! Description:
-! ------------
-!> Determines topography-top index at given (j,i)-position.  
-!------------------------------------------------------------------------------!
-    FUNCTION get_topography_top_index_ji( j, i, grid )
-
-       IMPLICIT NONE 
-
-       CHARACTER(LEN=*) ::  grid                         !< flag to distinquish between staggered grids
-       INTEGER(iwp)     ::  i                            !< grid index in x-dimension
-       INTEGER(iwp)     ::  ibit                         !< bit position where topography information is stored on respective grid
-       INTEGER(iwp)     ::  j                            !< grid index in y-dimension
-       INTEGER(iwp)     ::  get_topography_top_index_ji  !< topography top index
-
-       SELECT CASE ( TRIM( grid ) )
-
-          CASE ( 's'     )
-             ibit = 12
-          CASE ( 'u'     )
-             ibit = 14
-          CASE ( 'v'     )
-             ibit = 16
-          CASE ( 'w'     )
-             ibit = 18
-          CASE ( 's_out' )
-             ibit = 24
-          CASE DEFAULT
-!
-!--          Set default to scalar grid
-             ibit = 12 
-
-       END SELECT
-
-       get_topography_top_index_ji = MAXLOC(                                   &
-                                     MERGE( 1, 0,                              &
-                                            BTEST( wall_flags_0(:,j,i), ibit ) &
-                                          ), DIM = 1                           &
-                                           ) - 1
-
-       RETURN
-
-    END FUNCTION get_topography_top_index_ji
-
-!------------------------------------------------------------------------------!
-! Description:
-! ------------
-!> Determines topography-top index at each (j,i)-position.  
-!------------------------------------------------------------------------------!
-    FUNCTION get_topography_top_index( grid )
-
-       IMPLICIT NONE 
-
-       CHARACTER(LEN=*) ::  grid                      !< flag to distinquish between staggered grids
-       INTEGER(iwp)     ::  ibit                      !< bit position where topography information is stored on respective grid
-       INTEGER(iwp), DIMENSION(nys:nyn,nxl:nxr) ::  get_topography_top_index  !< topography top index
-
-       SELECT CASE ( TRIM( grid ) )
-
-          CASE ( 's'     )
-             ibit = 12
-          CASE ( 'u'     )
-             ibit = 14
-          CASE ( 'v'     )
-             ibit = 16
-          CASE ( 'w'     )
-             ibit = 18
-          CASE ( 's_out' )
-             ibit = 24
-          CASE DEFAULT
-!
-!--          Set default to scalar grid
-             ibit = 12 
-
-       END SELECT
-
-       get_topography_top_index(nys:nyn,nxl:nxr) = MAXLOC(                     &
-                         MERGE( 1, 0,                                          &
-                                 BTEST( wall_flags_0(:,nys:nyn,nxl:nxr), ibit )&
-                              ), DIM = 1                                       &
-                                                         ) - 1
-
-       RETURN
-
-    END FUNCTION get_topography_top_index
 
 !------------------------------------------------------------------------------!