Skip to content

Commit

Permalink
FDS Source: use Z_LS to build terrain files for both OBST and GEOM cases
Browse files Browse the repository at this point in the history
  • Loading branch information
ericvmueller committed Jan 6, 2025
1 parent 6cee736 commit eb3f4de
Showing 1 changed file with 1 addition and 5 deletions.
6 changes: 1 addition & 5 deletions Source/dump.f90
Original file line number Diff line number Diff line change
Expand Up @@ -988,11 +988,7 @@ SUBROUTINE INITIALIZE_MESH_DUMPS(NM)
ALLOCATE(OUT_OF_MESH(0:M%IBAR,0:M%JBAR)) ; OUT_OF_MESH = .FALSE.
DO J=0,M%JBAR
DO I=0,M%IBAR
IF (CC_IBM) THEN
Z_TERRAIN(I,J) = REAL(M%GEOM_ZMAX(I,J),FB)
ELSE
Z_TERRAIN(I,J) = REAL(0.25_EB*(M%Z_LS(I,J)+M%Z_LS(I+1,J)+M%Z_LS(I,J+1)+M%Z_LS(I+1,J+1)),FB)
ENDIF
Z_TERRAIN(I,J) = REAL(0.25_EB*(M%Z_LS(I,J)+M%Z_LS(I+1,J)+M%Z_LS(I,J+1)+M%Z_LS(I+1,J+1)),FB)
IF (Z_TERRAIN(I,J)<M%ZS .OR. Z_TERRAIN(I,J)>M%ZF) OUT_OF_MESH(I,J) = .TRUE.
ENDDO
ENDDO
Expand Down

0 comments on commit eb3f4de

Please sign in to comment.