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

Timestamp:

Aug 27, 2020 8:54:43 AM (4 years ago)

Author:

pavelkrc

Message:

Radiation Transfer Model version 4.0

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -2411 +2411 @@
 !-- Pre-calculate topography top indices (former get_topography_top_index
 !-- function)
-    ALLOCATE( topo_top_ind(nysg:nyng,nxlg:nxrg,0:4) )
+    ALLOCATE( topo_top_ind(nysg:nyng,nxlg:nxrg,0:5) )
 !
 !-- Uppermost topography index on scalar grid
@@ -2437 +2437 @@
     topo_top_ind(:,:,4) = MAXLOC( MERGE( 1, 0, BTEST( wall_flags_total_0(:,:,:), ibit ) ), DIM=1 ) &
                           - 1
+!
+!-- Uppermost topography index including full-3D geometry
+    ibit = 12
+    DO  k = nzb, nzt+1
+       WHERE( BTEST( wall_flags_total_0(k,:,:), ibit ) )  topo_top_ind(:,:,5) = k
+    ENDDO
 
  END SUBROUTINE set_topo_flags