Merge pull request firemodels#13988 from mcgratta/master

FDS Source: ERROR for misuse of TUNNEL_PRECONDITIONER
gforney · Jan 6, 2025 · 3a1f30f · 3a1f30f
2 parents 1d89541 + 0d5bdf6
commit 3a1f30f
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 0 deletions.
diff --git a/Manuals/FDS_User_Guide/FDS_User_Guide.tex b/Manuals/FDS_User_Guide/FDS_User_Guide.tex
@@ -13976,6 +13976,7 @@ \chapter{Error Codes}
 373  \> \ct{Cannot have FISHPAK_BC>0.                                             } \> Section~\ref{sec:periodic}  \\
 374  \> \ct{Numerical Instability - FDS stopped.                                   } \> Section~\ref{info:Errors}  \\
 375  \> \ct{OBST ... is VARIABLE_THICKNESS or HT3D and needs a MATL_ID.      } \> Section~\ref{info:HT3D_Limitations}  \\
+376  \> \ct{Meshes must have the same y/z bounds for TUNNEL_PRECONDITIONER.  } \> Section~\ref{tunnel_solver}  \\
      \>                                                                               \>                                     \\
 381  \> \ct{Need more spectral band limits.                                        } \> Section~\ref{info:RADI_Wide_Band}  \\
 382  \> \ct{Spectral band limits should be given in ascending order.               } \> Section~\ref{info:RADI_Wide_Band}  \\

diff --git a/Source/read.f90 b/Source/read.f90
@@ -9807,6 +9807,12 @@ SUBROUTINE READ_PRES
    IF (MAX_PRESSURE_ITERATIONS<20) MAX_PRESSURE_ITERATIONS = 20
    TUNNEL_NXP = 0
    DO NM=1,NMESHES
+      IF (NM>1) THEN
+         IF (MESHES(NM)%YS/=MESHES(NM-1)%YS .OR. MESHES(NM)%YF/=MESHES(NM-1)%YF .OR. &
+             MESHES(NM)%ZS/=MESHES(NM-1)%ZS .OR. MESHES(NM)%ZF/=MESHES(NM-1)%ZF) THEN
+            CALL SHUTDOWN('ERROR(376): Meshes must have the same y/z bounds for TUNNEL_PRECONDITIONER.') ; RETURN
+         ENDIF
+      ENDIF
       TUNNEL_NXP = TUNNEL_NXP + MESHES(NM)%IBAR
    ENDDO
    ALLOCATE(TP_AA(TUNNEL_NXP))