Changeset 2317 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Jul 20, 2017 5:27:19 PM (7 years ago)

Author:

suehring

Message:

get topograpyh top index via function call

File:

• palm/trunk/SOURCE/init_grid.f90 (modified) (3 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Get topography top index via Function call 
 ! 
 ! Former revisions:
@@ -278 +278 @@
 
     USE surface_mod,                                                           &
-        ONLY:  init_bc
+        ONLY:  get_topography_top_index, init_bc
 
     IMPLICIT NONE
@@ -1713 +1713 @@
 !--    nxl to nxr and nys to nyn on south and right model boundary, hence, 
 !--    use intrinsic lbound and ubound functions to infer array bounds.
-       DO  i = lbound(zu_s_inner, 1), ubound(zu_s_inner, 1)
-          DO  j = lbound(zu_s_inner, 2), ubound(zu_s_inner, 2)
+       DO  i = LBOUND(zu_s_inner, 1), UBOUND(zu_s_inner, 1)
+          DO  j = LBOUND(zu_s_inner, 2), UBOUND(zu_s_inner, 2)
 !
 !--          Topography height on scalar grid. Therefore, determine index of 
-!--          upward-facing surface element on scalar grid (bit 12). 
-             zu_s_inner(i,j) = zu( MAXLOC( MERGE(                              &
-                                         1, 0, BTEST( wall_flags_0(:,j,i), 12 )&
-                                                ), DIM = 1                     &
-                                         ) - 1                                 &
-                                 )
+!--          upward-facing surface element on scalar grid. 
+             zu_s_inner(i,j) = zu( get_topography_top_index( j, i, 's' ) )
+
+             write(9,*)  get_topography_top_index( j, i, 's' ), MAXLOC(                                                    &
+                          MERGE( 1, 0,                                         &
+                                 BTEST( wall_flags_0(:,j,i), 12 )  &
+                               ), DIM = 1                                      &
+                          ) - 1
 !
 !--          Topography height on w grid. Therefore, determine index of 
-!--          upward-facing surface element on w grid (bit 18).
-             zw_w_inner(i,j) = zw( MAXLOC( MERGE(                              &
-                                         1, 0, BTEST( wall_flags_0(:,j,i), 18 )&
-                                                ), DIM = 1                     &
-                                         ) - 1                                 &
-                                 )
+!--          upward-facing surface element on w grid.
+             zw_w_inner(i,j) = zw( get_topography_top_index( j, i, 's' ) )
           ENDDO
        ENDDO
-
+       flush(9)
 
     ENDIF