diff --git a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv
index 56ae60b0990..732ecc0bcc8 100644
--- a/Utilities/Matlab/FDS_verification_dataplot_inputs.csv
+++ b/Utilities/Matlab/FDS_verification_dataplot_inputs.csv
@@ -109,15 +109,18 @@ d,cyl_test,Flowfields/cyl_test_3_git.txt,Flowfields/cyl_test_3_devc.csv,2,3,Time
 d,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,Preconditioner,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_tight_no_precon_devc.csv,2,3,Time,iter,No Preconditioner,k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
 f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_devc.csv,2,3,Time,iter,UGLMAT,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
 f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_devc.csv,2,3,Time,iter,ULMAT,b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
+f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_hypre_devc.csv,2,3,Time,iter,UGLMAT HYPRE,m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
 f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,iter,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_hypre_devc.csv,2,3,Time,iter,ULMAT HYPRE,g-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (dancing\_eddies),Time (s),Pressure Iterations,0,2,1,0,100,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
 d,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,Preconditioner,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_tight_no_precon_devc.csv,2,3,Time,cpu,No Preconditioner,k--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
 f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_devc.csv,2,3,Time,cpu,UGLMAT,r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
 f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_devc.csv,2,3,Time,cpu,ULMAT,b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
+f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_hypre_devc.csv,2,3,Time,cpu,UGLMAT HYPRE,m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
 f,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,cpu,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_hypre_devc.csv,2,3,Time,cpu,ULMAT HYPRE,g-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,CPU Time (dancing\_eddies),Simulated Time (s),CPU Time (s),0,2,1,0,400,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight_precon,N/A,end,0,Pressure Solver,kd,k,TeX
 d,dancing_eddies,Pressure_Solver/dancing_eddies_tight_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,"FDS, 1 mesh",k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_tight_devc.csv,2,3,Time,pres,"FDS, 4 mesh, tol=1E-5 m/s",r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_tight),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_tight,Absolute Error,end,0.3,Pressure Solver,kd,k,TeX
 d,dancing_eddies,Pressure_Solver/dancing_eddies_default_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,"FDS, 1 mesh",k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_default_devc.csv,2,3,Time,pres,"FDS, 4 mesh, tol=5E-4 m/s",r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_default),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_default,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
 d,dancing_eddies,Pressure_Solver/dancing_eddies_uglmat_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,"FDS, 1 mesh",k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_devc.csv,2,3,Time,pres,"FDS, 4 mesh, uglmat",r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_uglmat),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_uglmat,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
 f,dancing_eddies,Pressure_Solver/dancing_eddies_uglmat_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_devc.csv,2,3,Time,pres,"FDS, 4 mesh, ulmat",b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_uglmat),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_uglmat,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
+f,dancing_eddies,Pressure_Solver/dancing_eddies_uglmat_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_uglmat_hypre_devc.csv,2,3,Time,pres,"FDS, 4 mesh, uglmat hypre",m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_uglmat),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_uglmat,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
 f,dancing_eddies,Pressure_Solver/dancing_eddies_uglmat_git.txt,Pressure_Solver/dancing_eddies_1mesh_devc.csv,2,3,Time,pres,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/dancing_eddies_ulmat_hypre_devc.csv,2,3,Time,pres,"FDS, 4 mesh, ulmat hypre",g-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure (dancing\_eddies\_uglmat),Time (s),Pressure (Pa),0,2,1,-0.3,0.2,1,no,0.05 0.90,SouthWest,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/dancing_eddies_uglmat,Absolute Error,end,0.01,Pressure Solver,kd,k,TeX
 d,device_restart,Restart/device_restart_base_case_git.txt,Restart/device_restart_base_case_devc.csv,2,3,Time,temp,Uninterrupted Case,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Restart/device_restart_a_devc.csv,2,3,Time,temp,Interrupted Case,ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Temperature (device\_restart),Time (s),Temperature (°C),0,60,1,0,100,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/device_restart,Absolute Error,end,5,Restart,kd,k,TeX
 d,devc_interpolation,Miscellaneous/devc_interpolation_git.txt,Miscellaneous/devc_interpolation.csv,1,2,z,u,{\it u} prescribed,k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Miscellaneous/devc_interpolation_velocity_line.csv,2,3,Z,U-VEL-1|U-VEL-2|U-VEL-3,U-VEL-1|U-VEL-2|U-VEL-3,r--|b--|g--,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Velocity (devc\_interpolation),{\it z} (m),Velocity (m/s),0,1,1,0,1,1,no,0.05 0.90,SouthEast,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/devc_interpolation_velocity,Absolute Error,end_1_3,0.01,Miscellaneous,kd,k,TeX
@@ -144,8 +147,10 @@ d,droplet_distributions,Sprinklers_and_Sprays/droplet_distributions_2_git.txt,Sp
 d,droplet_distributions,Sprinklers_and_Sprays/droplet_distributions_2_git.txt,Sprinklers_and_Sprays/droplet_distributions_2.csv,1,2,Time,Mass,Total Mass (Prescibed),ko,0,40000,,0,100000,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/droplet_distributions_2_devc.csv,2,3,Time,USER_M,PDPA Mass (Sampled),k-,0,4000,,0,100000,-1.00E+09,1.00E+09,0,User Defined Distribution Total Mass,Time (s),Mass (kg),0,11,1,0,0.1,1,no,0.05 0.90,West,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/droplet_distributions_2_4_mass,Relative Error,end,2.00E-02,Sprinklers and Sprays,ro,r,TeX
 d,duct_flow,Pressure_Solver/duct_flow_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Flow,Ideal (Flow),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_devc.csv,2,3,Time,flow_in|flow_out,FFT IBM (flow\_in)|FFT IBM (flow\_out),k-|r-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Volume Flow (duct\_flow),Time (s),Volume Flow (m³/s),0,60,1,0,1.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow,Relative Error,end_1_2,0.05,Pressure Solver,kd,k,TeX
 f,duct_flow,Pressure_Solver/duct_flow_uglmat_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Flow,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_uglmat_devc.csv,2,3,Time,flow_in|flow_out,UGLMAT (flow\_in)|UGLMAT (flow\_out),k-|b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Volume Flow (duct\_flow),Time (s),Volume Flow (m³/s),0,60,1,0,1.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow,Relative Error,end_1_2,0.05,Pressure Solver,kd,k,TeX
+f,duct_flow,Pressure_Solver/duct_flow_uglmat_hypre_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Flow,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_uglmat_hypre_devc.csv,2,3,Time,flow_in|flow_out,UGLMAT HYPRE (flow\_in)|UGLMAT HYPRE (flow\_out),k-|m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Volume Flow (duct\_flow),Time (s),Volume Flow (m³/s),0,60,1,0,1.5,1,no,0.05 0.90,SouthEast,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow,Relative Error,end_1_2,0.05,Pressure Solver,kd,k,TeX
 d,duct_flow,Pressure_Solver/duct_flow_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Iterations,Ideal (Iterations),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_devc.csv,2,3,Time,iter,FFT IBM (iter),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (duct\_flow),Time (s),Iterations,0,60,1,0,100,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
 f,duct_flow,Pressure_Solver/duct_flow_uglmat_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Iterations,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_uglmat_devc.csv,2,3,Time,iter,UGLMAT (iter),b-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (duct\_flow),Time (s),Iterations,0,60,1,0,100,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
+f,duct_flow,Pressure_Solver/duct_flow_uglmat_hypre_git.txt,Pressure_Solver/duct_flow.csv,1,2,Time,Iterations,blank,blank,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure_Solver/duct_flow_uglmat_hypre_devc.csv,2,3,Time,iter,UGLMAT HYPRE (iter),m-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Pressure Iterations (duct\_flow),Time (s),Iterations,0,60,1,0,100,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/duct_flow_iterations,N/A,end,0,Pressure Solver,kd,k,TeX
 d,e_coefficient,Sprinklers_and_Sprays/e_coefficient_git.txt,Sprinklers_and_Sprays/e_coefficient.csv,1,1,Time,HRR,Expected,ko,0,100000,,14.9,15.1,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/e_coefficient_devc.csv,2,3,Time,HRR_1,FDS,k-,0,100000,,14.9,15.1,-1.00E+09,1.00E+09,0,Heat Release Rate (OBST),Time (s),Heat Release Rate (kW),0,30,1,0,60,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/e_coefficient_obst,Absolute Error,max,3,Sprinklers and Sprays,bs,b,TeX
 d,e_coefficient,Sprinklers_and_Sprays/e_coefficient_git.txt,Sprinklers_and_Sprays/e_coefficient.csv,1,1,Time,HRR,Expected,ko,0,100000,,14.9,15.1,-1.00E+09,1.00E+09,0,Sprinklers_and_Sprays/e_coefficient_devc.csv,2,3,Time,HRR_2,FDS,k-,0,100000,,14.9,15.1,-1.00E+09,1.00E+09,0,Heat Release Rate (GEOM),Time (s),Heat Release Rate (kW),0,30,1,0,60,1,no,0.05 0.90,East,,1,linear,FDS_User_Guide/SCRIPT_FIGURES/e_coefficient_geom,Absolute Error,max,3,Sprinklers and Sprays,bs,b,TeX
 d,emissivity,Radiation/emissivity_git.txt,Radiation/emissivity.csv,1,2,Time,Emissivity,Analytical (Emissivity),ko,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Radiation/emissivity_devc.csv,2,3,Time,EMISSIVITY,FDS (EMISSIVITY),k-,0,100000,,0,100000,-1.00E+09,1.00E+09,0,Surface Emissivity (emissivity),Time (s),Emissivity,0,10,1,0,1,1,no,0.30 0.90,SouthWest,,1,linear,FDS_Verification_Guide/SCRIPT_FIGURES/emissivity,Absolute Error,end,0.1,Radiation,kd,k,TeX
diff --git a/Verification/FDS_Cases.sh b/Verification/FDS_Cases.sh
index a2badee239c..d9c20a06828 100755
--- a/Verification/FDS_Cases.sh
+++ b/Verification/FDS_Cases.sh
@@ -402,6 +402,7 @@ $QFDS -p 8 -d Pressure_Solver duct_flow.fds
 $QFDS -p 8 -d Pressure_Solver duct_flow_ulmat.fds
 $QFDS -p 8 -d Pressure_Solver duct_flow_ulmat_hypre.fds
 $QFDS -p 8 -d Pressure_Solver duct_flow_uglmat.fds
+$QFDS -p 8 -d Pressure_Solver duct_flow_uglmat_hypre.fds
 $QFDS -p 5 -d Pressure_Solver hallways.fds
 $QFDS -p 8 -d Pressure_Solver tunnel_demo.fds
 $QFDS -p 8 -d Pressure_Solver tunnel_demo_glmat.fds
diff --git a/Verification/Pressure_Solver/duct_flow_uglmat_hypre.fds b/Verification/Pressure_Solver/duct_flow_uglmat_hypre.fds
new file mode 100644
index 00000000000..3e9dfa54fff
--- /dev/null
+++ b/Verification/Pressure_Solver/duct_flow_uglmat_hypre.fds
@@ -0,0 +1,86 @@
+Test pressure iteration routine for flow through a duct across meshes
+
+&HEAD CHID='duct_flow_uglmat_hypre', TITLE='Flow through a duct, test of UGLMAT with HYPRE' /
+
+&MESH IJK=16,16,16, XB=-.2,3.0,-.2,3.0,-.2,3.0, MULT_ID='mesh' /
+&MULT ID='mesh', DX=3.2, DY=3.2, DZ=3.2, I_UPPER=1, J_UPPER=1, K_UPPER=1 /
+
+MESH IJK=32,32,32, XB=-.2,6.2,-.2,6.2,-.2,6.2 /
+
+&RADI RADIATION=.FALSE. /
+&TIME T_END=60.0 /
+
+&MISC STRATIFICATION=.FALSE. /
+&PRES CHECK_POISSON=T, SOLVER='UGLMAT HYPRE' /
+
+&SURF ID='DUCT', COLOR='GRAY 30', TRANSPARENCY=0.5 /
+
+&OBST XB=1.0,1.1,1.0,2.0,1.0,2.0, SURF_ID='DUCT', COLOR='RED',TRANSPARENCY=0.5 /
+&VENT XB=1.0,1.0,1.0,2.0,1.0,2.0, SURF_ID='HVAC', ID='SUCK', COLOR='RED',TRANSPARENCY=0.5/
+&VENT XB=1.1,1.1,1.0,2.0,1.0,2.0, SURF_ID='HVAC', ID='BLOW', COLOR='RED',TRANSPARENCY=0.5/
+&HVAC TYPE_ID='NODE',ID='SUCK',VENT_ID='SUCK',DUCT_ID='FAN'/
+&HVAC TYPE_ID='NODE',ID='BLOW',VENT_ID='BLOW',DUCT_ID='FAN'/
+&HVAC TYPE_ID='DUCT',ID='FAN',NODE_ID='SUCK','BLOW',VOLUME_FLOW=1.,AREA=1./
+
+&OBST XB=1.1,5.0,1.0,1.0,1.0,2.0, SURF_ID='DUCT' /
+&OBST XB=1.1,4.0,2.0,2.0,1.0,2.0, SURF_ID='DUCT' /
+&OBST XB=1.1,5.0,1.0,2.0,1.0,1.0, SURF_ID='DUCT' /
+&OBST XB=1.1,5.0,1.0,2.0,2.0,2.0, SURF_ID='DUCT' /
+
+&OBST XB=5.0,5.0,1.0,5.0,1.0,2.0, SURF_ID='DUCT' /
+&OBST XB=4.0,4.0,2.0,4.0,1.0,2.0, SURF_ID='DUCT' /
+&OBST XB=4.0,5.0,2.0,5.0,1.0,1.0, SURF_ID='DUCT' /
+&OBST XB=4.0,5.0,2.0,5.0,2.0,2.0, SURF_ID='DUCT' /
+
+&OBST XB=1.0,4.0,4.0,4.0,1.0,2.0, SURF_ID='DUCT' /
+&OBST XB=1.0,5.0,5.0,5.0,1.0,2.0, SURF_ID='DUCT' /
+&OBST XB=1.0,4.0,4.0,5.0,1.0,1.0, SURF_ID='DUCT' /
+&OBST XB=2.0,4.0,4.0,5.0,2.0,2.0, SURF_ID='DUCT' /
+
+&OBST XB=1.0,1.0,4.0,5.0,1.0,5.0, SURF_ID='DUCT' /
+&OBST XB=2.0,2.0,4.0,5.0,2.0,4.0, SURF_ID='DUCT' /
+&OBST XB=1.0,2.0,4.0,4.0,2.0,4.0, SURF_ID='DUCT' /
+&OBST XB=1.0,2.0,5.0,5.0,2.0,4.0, SURF_ID='DUCT' /
+
+&OBST XB=1.0,5.0,1.0,1.0,4.0,5.0, SURF_ID='DUCT' /
+&OBST XB=1.0,4.0,2.0,2.0,4.0,5.0, SURF_ID='DUCT' /
+&OBST XB=1.0,5.0,1.0,2.0,4.0,4.0, SURF_ID='DUCT' /
+&OBST XB=1.0,5.0,1.0,2.0,5.0,5.0, SURF_ID='DUCT' /
+
+&OBST XB=5.0,5.0,1.0,5.0,4.0,5.0, SURF_ID='DUCT' /
+&OBST XB=4.0,4.0,2.0,4.0,4.0,5.0, SURF_ID='DUCT' /
+&OBST XB=4.0,5.0,2.0,5.0,4.0,4.0, SURF_ID='DUCT' /
+&OBST XB=4.0,5.0,2.0,5.0,5.0,5.0, SURF_ID='DUCT' /
+
+&OBST XB=1.0,4.0,4.0,4.0,4.0,5.0, SURF_ID='DUCT' /
+&OBST XB=1.0,5.0,5.0,5.0,4.0,5.0, SURF_ID='DUCT' /
+&OBST XB=2.0,4.0,4.0,5.0,4.0,4.0, SURF_ID='DUCT' /
+&OBST XB=1.0,4.0,4.0,5.0,5.0,5.0, SURF_ID='DUCT' /
+
+&VENT MB='XMIN', SURF_ID='OPEN' /
+&VENT MB='XMAX', SURF_ID='OPEN' /
+&VENT MB='YMIN', SURF_ID='OPEN' /
+&VENT MB='YMAX', SURF_ID='OPEN' /
+&VENT MB='ZMIN', SURF_ID='OPEN' /
+&VENT MB='ZMAX', SURF_ID='OPEN' /
+
+&SLCF PBY=1.5,QUANTITY='VELOCITY',VECTOR=.TRUE. /
+&SLCF PBY=1.5,QUANTITY='H' /
+
+&SLCF PBY=4.5,QUANTITY='VELOCITY',VECTOR=.TRUE. /
+&SLCF PBY=4.5,QUANTITY='H' /
+
+&SLCF PBZ=1.5,QUANTITY='VELOCITY',VECTOR=.TRUE. /
+&SLCF PBZ=1.5,QUANTITY='H' /
+
+&SLCF PBZ=4.5,QUANTITY='VELOCITY',VECTOR=.TRUE. /
+&SLCF PBZ=4.5,QUANTITY='H' /
+
+&BNDF QUANTITY='VELOCITY ERROR' /
+
+&DEVC XB=1.0,1.0,1.0,2.0,1.0,2.0, QUANTITY='U-VELOCITY', SPATIAL_STATISTIC='AREA INTEGRAL', QUANTITY_RANGE(1)=0, ID='flow_in' /
+&DEVC XB=1.0,1.0,1.0,2.0,4.0,5.0, QUANTITY='U-VELOCITY', SPATIAL_STATISTIC='AREA INTEGRAL', QUANTITY_RANGE(2)=0, ID='flow_out', CONVERSION_FACTOR=-1. /
+&DEVC XYZ=0.,0.,0., QUANTITY='MAXIMUM VELOCITY ERROR', ID='error' /
+&DEVC XYZ=0.,0.,0., QUANTITY='PRESSURE ITERATIONS', ID='iter' /
+
+&TAIL /