Closed khalidmoinuddin closed 8 years ago
I'll take a look at it.
The problem is that the particles are not being properly distributed over the 4 meshes. I see in the 4 mesh case 130400 particles in each mesh. There should be a total of 130400 for the entire domain. I'll fix.
I fixed the obvious bug related to how FDS distributes particles within a cone. In short, we need to know how much of the cone volume is in each mesh and parcel out the particles accordingly. I pushed new code (#3863). Compile the new code and run your case again. I still see a difference in HRR. I also split the mesh into 4 meshes stacked one on top of the other and the agreement between single and 4 mesh case is very close.
Good to know the progress. Where can I get the new code?
Right here at GitHub. It was merged in #3863. Update you local repository and recompile.
Hi Kevin,
We compiled the code. The case file gives numerical instability after 28 time step. We noticed this with previous versions of the code as well (6.2.0, 6.4.0). This is both single and multimesh cases.
Regards
Khalid
From: Kevin McGrattan [mailto:notifications@github.com] Sent: Thursday, 19 May 2016 10:23 PM To: firemodels/fds-smv Cc: Khalid Moinuddin; Author Subject: Re: [firemodels/fds-smv] Unusual behavoiur in HRR using multiple meshes (#3856)
Right here at GitHub. It was merged in #3863https://github.com/firemodels/fds-smv/pull/3863. Update you local repository and recompile.
— You are receiving this because you authored the thread. Reply to this email directly or view it on GitHubhttps://github.com/firemodels/fds-smv/issues/3856#issuecomment-220308456
This email, including any attachment, is intended solely for the use of the intended recipient. It is confidential and may contain personal information or be subject to legal professional privilege. If you are not the intended recipient any use, disclosure, reproduction or storage of it is unauthorised. If you have received this email in error, please advise the sender via return email and delete it from your system immediately. Victoria University does not warrant that this email is free from viruses or defects and accepts no liability for any damage caused by such viruses or defects.
Copy the version number, Git hash, date, etc, of the version you compiled.
Hi Kevin,
The requested information is given below:
Current Date : May 25, 2016 14:18:17 Version : FDS 6.4.0 Revision : Git-r10-517-gc7aa429 Revision Date : Tue May 24 16:37:10 2016 -0400 Compilation Date : May 25, 2016 14:13:11
In the attachment I've copy-pasted all the output from running the cases.
Kind regards
Khalid
From: Kevin McGrattan [mailto:notifications@github.com] Sent: Tuesday, 24 May 2016 10:27 PM To: firemodels/fds-smv Cc: Khalid Moinuddin; Author Subject: Re: [firemodels/fds-smv] Unusual behavoiur in HRR using multiple meshes (#3856)
Copy the version number, Git hash, date, etc, of the version you compiled.
— You are receiving this because you authored the thread. Reply to this email directly or view it on GitHubhttps://github.com/firemodels/fds-smv/issues/3856#issuecomment-221253008
This email, including any attachment, is intended solely for the use of the intended recipient. It is confidential and may contain personal information or be subject to legal professional privilege. If you are not the intended recipient any use, disclosure, reproduction or storage of it is unauthorised. If you have received this email in error, please advise the sender via return email and delete it from your system immediately. Victoria University does not warrant that this email is free from viruses or defects and accepts no liability for any damage caused by such viruses or defects.
$ ./fds_linux_64 tree_2_m_14_pc.fds OpenMP thread 18 of 31 is running OpenMP thread 22 of 31 is running OpenMP thread 24 of 31 is running OpenMP thread 0 of 31 is running OpenMP thread 3 of 31 is running OpenMP thread 7 of 31 is running OpenMP thread 13 of 31 is running OpenMP thread 30 of 31 is running OpenMP thread 14 of 31 is running OpenMP thread 26 of 31 is running OpenMP thread 15 of 31 is running OpenMP thread 16 of 31 is running OpenMP thread 6 of 31 is running OpenMP thread 25 of 31 is running OpenMP thread 2 of 31 is running OpenMP thread 28 of 31 is running OpenMP thread 8 of 31 is running OpenMP thread 29 of 31 is running OpenMP thread 12 of 31 is running OpenMP thread 27 of 31 is running OpenMP thread 11 of 31 is running OpenMP thread 17 of 31 is running OpenMP thread 20 of 31 is running OpenMP thread 9 of 31 is running OpenMP thread 21 of 31 is running OpenMP thread 4 of 31 is running OpenMP thread 5 of 31 is running OpenMP thread 1 of 31 is running OpenMP thread 19 of 31 is running OpenMP thread 31 of 31 is running OpenMP thread 10 of 31 is running OpenMP thread 23 of 31 is running Completed Initialization Step 1 Completed Initialization Step 2 Completed Initialization Step 3 Completed Initialization Step 4
Fire Dynamics Simulator
Current Date : May 25, 2016 14:07:37 Version : FDS 6.4.0 Revision : Git-r10-517-gc7aa429 Revision Date : Tue May 24 16:37:10 2016 -0400 Compilation Date : May 25, 2016 14:02:05
MPI Disabled OpenMP Enabled; Number of OpenMP Threads: 32
Job TITLE : Test for burning vegetation mass conservation Job ID string : tree_2_m_14_pc
Time Step: 1, Simulation Time: 0.04 s Time Step: 2, Simulation Time: 0.08 s Time Step: 3, Simulation Time: 0.12 s Time Step: 4, Simulation Time: 0.16 s Time Step: 5, Simulation Time: 0.20 s Time Step: 6, Simulation Time: 0.24 s Time Step: 7, Simulation Time: 0.28 s Time Step: 8, Simulation Time: 0.31 s Time Step: 9, Simulation Time: 0.34 s Time Step: 10, Simulation Time: 0.37 s Time Step: 20, Simulation Time: 0.67 s Time Step: 30, Simulation Time: 0.93 s Time Step: 38, Simulation Time: 1.05 s
STOP: Numerical Instability (CHID: tree_2_m_14_pc)
$ ./fds_linux_64 tree_2_m_14_pc.fds Completed Initialization Step 1 Completed Initialization Step 2 Completed Initialization Step 3 Completed Initialization Step 4
Fire Dynamics Simulator
Current Date : May 25, 2016 14:10:43 Version : FDS 6.4.0 Revision : Git-r10-517-gc7aa429 Revision Date : Tue May 24 16:37:10 2016 -0400 Compilation Date : May 25, 2016 14:02:05
MPI Disabled OpenMP Enabled; Number of OpenMP Threads: 1
Job TITLE : Test for burning vegetation mass conservation Job ID string : tree_2_m_14_pc
Time Step: 1, Simulation Time: 0.04 s Time Step: 2, Simulation Time: 0.08 s Time Step: 3, Simulation Time: 0.12 s Time Step: 4, Simulation Time: 0.16 s Time Step: 5, Simulation Time: 0.20 s Time Step: 6, Simulation Time: 0.24 s Time Step: 7, Simulation Time: 0.28 s Time Step: 8, Simulation Time: 0.31 s Time Step: 9, Simulation Time: 0.34 s Time Step: 10, Simulation Time: 0.37 s Time Step: 20, Simulation Time: 0.67 s Time Step: 30, Simulation Time: 0.93 s Time Step: 38, Simulation Time: 1.05 s
STOP: Numerical Instability (CHID: tree_2_m_14_pc)
$mpirun -np 2 fds_mpi_gnu_linux_64 tree_5_m_26_pc.fds Mesh 1 is assigned to MPI Process 0 Mesh 2 is assigned to MPI Process 1 MPI process 1 of 1 MPI process 0 of 1 Completed Initialization Step 1 Completed Initialization Step 2 Completed Initialization Step 3 Completed Initialization Step 4
Fire Dynamics Simulator
Current Date : May 25, 2016 14:18:17 Version : FDS 6.4.0 Revision : Git-r10-517-gc7aa429 Revision Date : Tue May 24 16:37:10 2016 -0400 Compilation Date : May 25, 2016 14:13:11
MPI Enabled; Number of MPI Processes: 2 OpenMP Disabled
MPI version: 3.0 MPI library version: Open MPI v1.10.2, package: Open MPI duncans@EWBTEST-01 Distribution, ident: 1.10.2, repo rev: v1.10.1-145-g799148f, Jan 21, 2016
Job TITLE : NIST Douglas Fir Fires, 5 m tall, 26% moisture content Job ID string : tree_5_m_26_pc
Time Step: 1, Simulation Time: 0.03 s Time Step: 2, Simulation Time: 0.05 s Time Step: 3, Simulation Time: 0.08 s Time Step: 4, Simulation Time: 0.10 s Time Step: 5, Simulation Time: 0.13 s Time Step: 6, Simulation Time: 0.15 s Time Step: 7, Simulation Time: 0.18 s Time Step: 8, Simulation Time: 0.20 s Time Step: 9, Simulation Time: 0.23 s Time Step: 10, Simulation Time: 0.25 s Time Step: 16, Simulation Time: 0.37 s MPI process 1 has completed MPI process 0 has completed
STOP: Numerical Instability (CHID: tree_5_m_26_pc)
I assume you are running the NIST_Douglas_Firs case files in the repository. I am too, and I am using the same FDS version as you. I am not (yet) seeing numerical instabilities. Which Fortran compiler do you use? Do you use the makefile that we provide? What is the version of your compiler?
Hi Kevin
Following your comments, we tried compiling the debugging version: gnu_linux_64_db, following the standard process with the fds makefile without modifications. The only real difference we can see between these two versions is the level of compiler optimisation -O0 for the db version and -O2 for the standard version.
The simulations is running past the point which it previously failed. It reached 100 steps for 2.11s of simulated time before we stopped it. We'll let it run longer and if there is an issue we will let you know.
The information as requested is below.
I assume you are running the NIST_Douglas_Firs case files in the repository. Yes. tree_2_m_14_pc.fds and tree_5_m_26_pc.fds
Which Fortran compiler do you use? gfortran
Do you use the makefile that we provide? Yes, we used the script in the FDS_Compilation/gnu_linux_64 directory (or the corresponding script in the mpi directory script). This is simply a nice wrapper script to the standard makefile.
What is the version of your compiler?
GNU Fortran (Ubuntu 4.8.4-2ubuntu1~14.04.1) 4.8.4
The only packages that we can see linked in the makefile are openmp and MPI.
This is the openmp library
lrwxrwxrwx 1 root root 16 Jan 27 17:35 /usr/lib/x86_64-linux-gnu/libgomp.so.1 -> libgomp.so.1.0.0
MPI version: 3.0
MPI library version: Open MPI v1.10.2,
No other external libraries appear to be linked, which surprises us.
Kind regards
Khalid
From: Kevin McGrattan [mailto:notifications@github.com] Sent: Wednesday, 25 May 2016 11:47 PM To: firemodels/fds-smv Cc: Khalid Moinuddin; Author Subject: Re: [firemodels/fds-smv] Unusual behavoiur in HRR using multiple meshes (#3856)
I assume you are running the NIST_Douglas_Firs case files in the repository. I am too, and I am using the same FDS version as you. I am not (yet) seeing numerical instabilities. Which Fortran compiler do you use? Do you use the makefile that we provide? What is the version of your compiler?
— You are receiving this because you authored the thread. Reply to this email directly or view it on GitHubhttps://github.com/firemodels/fds-smv/issues/3856#issuecomment-221580844
This email, including any attachment, is intended solely for the use of the intended recipient. It is confidential and may contain personal information or be subject to legal professional privilege. If you are not the intended recipient any use, disclosure, reproduction or storage of it is unauthorised. If you have received this email in error, please advise the sender via return email and delete it from your system immediately. Victoria University does not warrant that this email is free from viruses or defects and accepts no liability for any damage caused by such viruses or defects.
If the debug version runs and the optimized version does not, that generally indicates a problem with the compiler. I suggest you upgrade to the latest version if you're not there already, or consider purchasing the Intel Fortran compiler. We have found that it is more reliable and faster than gfortran.
Thanks, Kevin. We will purchase the Intel Fortran compiler.
Regards
Khalid
From: Kevin McGrattan [mailto:notifications@github.com] Sent: Thursday, 26 May 2016 10:53 PM To: firemodels/fds-smv Cc: Khalid Moinuddin; Author Subject: Re: [firemodels/fds-smv] Unusual behavoiur in HRR using multiple meshes (#3856)
If the debug version runs and the optimized version does not, that generally indicates a problem with the compiler. I suggest you upgrade to the latest version if you're not there already, or consider purchasing the Intel Fortran compiler. We have found that it is more reliable and faster than gfortran.
— You are receiving this because you authored the thread. Reply to this email directly or view it on GitHubhttps://github.com/firemodels/fds-smv/issues/3856#issuecomment-221862095
This email, including any attachment, is intended solely for the use of the intended recipient. It is confidential and may contain personal information or be subject to legal professional privilege. If you are not the intended recipient any use, disclosure, reproduction or storage of it is unauthorised. If you have received this email in error, please advise the sender via return email and delete it from your system immediately. Victoria University does not warrant that this email is free from viruses or defects and accepts no liability for any damage caused by such viruses or defects.
I carried out simulation with a Fir Tree validation case (2m14PC). When I used a single mesh I got one result and when I used multi-mesh I got different result . For multi-mesh I symmetrically divided the domain into 4 domains. Both FDS 6.2.0 and 6.4.0 were used. The problem was worse with FDS 6.4.0
There is a known understanding that the FDS' MPI introduces discontinuity in the pressure solver eg see the FDS' roadmap. We understand` that the developers are working towards improvement of pressure solver. However, in the meantime is there any other way, the simulation can be improved?
There is a suggestion not to put flames at mesh/ sub-domain boundaries. It is difficult to adhere to such suggestion when big fire (several MW) is involved. It becomes even harder when we carry out grid convergence test. `` Is there anything else that could cause the problem? Because when in simulation, no flame is involved, we do get expected results eg an isothermal channel flow or a backward facing step case.
Case file with 4 mesh is given below
&HEAD CHID='tree_2_m_14_pc-4Proc', TITLE='Test for burning vegetation mass conservation' /
&MESH IJK=30,30,160, XB=-1.5,0.0,-1.5,0.0,0.0,8.0 / &MESH IJK=30,30,160, XB=-1.5,0.0,0.0,1.5,0.0,8.0 /
&MESH IJK=30,30,160, XB=0.0,1.5,-1.5,0.0,0.0,8.0 / &MESH IJK=30,30,160, XB=0.0,1.5,0.0,1.5,0.0,8.0 /
&TIME T_END=30. /
&REAC FUEL='CELLULOSE', C=6, H=10, O=5, SOOT_YIELD=0.015 /
&SPEC ID='WATER VAPOR' /
&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='ZMAX', SURF_ID='OPEN' / &VENT MB='ZMIN', SURF_ID='FOIL' /
&SURF ID='FOIL' COLOR = 'SILVER' MATL_ID='ALUMINUM','CONCRETE' THICKNESS(1:2) = 0.0002,0.01 /
&MATL ID='ALUMINUM' DENSITY = 2707. SPECIFIC_HEAT = 0.896 CONDUCTIVITY = 204. EMISSIVITY = 0.1 /
&MATL ID='CONCRETE' DENSITY = 1500. SPECIFIC_HEAT = 1. CONDUCTIVITY = 0.3 /
&SURF ID = 'foliage' MATL_ID(1,1:2) = 'vegetation','MOISTURE' MATL_MASS_FRACTION(1,1:2) = 0.877,0.123 THICKNESS = 0.0005 LENGTH = 0.05 GEOMETRY = 'CYLINDRICAL' /
&SURF ID = 'small roundwood' MATL_ID(1,1:2) = 'vegetation','MOISTURE' MATL_MASS_FRACTION(1,1:2) = 0.877,0.123 THICKNESS = 0.001 LENGTH = 0.10 GEOMETRY = 'CYLINDRICAL' /
&SURF ID = 'medium roundwood' MATL_ID(1,1:2) = 'vegetation','MOISTURE' MATL_MASS_FRACTION(1,1:2) = 0.877,0.123 THICKNESS = 0.002 LENGTH = 0.10 GEOMETRY = 'CYLINDRICAL' /
&SURF ID = 'large roundwood' MATL_ID(1,1:2) = 'vegetation','MOISTURE' MATL_MASS_FRACTION(1,1:2) = 0.877,0.123 THICKNESS = 0.003 LENGTH = 0.10 GEOMETRY = 'CYLINDRICAL' /
&MATL ID = 'vegetation' DENSITY = 514. CONDUCTIVITY = 2.0 SPECIFIC_HEAT= 1.2 N_REACTIONS = 1 REFERENCE_TEMPERATURE = 200. REFERENCE_RATE = 0.0005 HEATING_RATE = 1.6 NU_SPEC = 0.74 SPEC_ID = 'CELLULOSE' NU_MATL = 0.26 MATL_ID = 'CHAR' HEAT_OF_COMBUSTION = 17700. HEAT_OF_REACTION = 418. /
&MATL ID = 'CHAR' DENSITY = 300. CONDUCTIVITY = 2.0 SPECIFIC_HEAT= 1.2 N_REACTIONS = 1 REFERENCE_TEMPERATURE = 350. REFERENCE_RATE = 0.0002 HEATING_RATE = 1.6 NU_SPEC = 1.00 SPEC_ID = 'CELLULOSE' HEAT_OF_COMBUSTION = 17700. HEAT_OF_REACTION = 418. /
&MATL ID = 'MOISTURE' DENSITY = 1000. CONDUCTIVITY = 2.0 SPECIFIC_HEAT= 4.184 N_REACTIONS = 1 REFERENCE_TEMPERATURE = 100. REFERENCE_RATE = 0.002 HEATING_RATE = 1.6 NU_SPEC = 1.0 SPEC_ID = 'WATER VAPOR' HEAT_OF_REACTION= 2500. /
&PART ID='foliage', DRAG_LAW='CYLINDER', SAMPLING_FACTOR=10, SURF_ID='foliage', PROP_ID='needle image' QUANTITIES='PARTICLE TEMPERATURE','PARTICLE MASS','PARTICLE DIAMETER', STATIC=.TRUE., COLOR='FOREST GREEN' / &INIT PART_ID='foliage', XYZ=0.0,0.0,0.15, RADIUS=0.85, HEIGHT=1.9, SHAPE='CONE', N_PARTICLES=100000, MASS_PER_VOLUME=2.0 / &PROP ID='needle image', SMOKEVIEW_ID='TUBE', SMOKEVIEW_PARAMETERS='L=0.05','D=0.001','RANDXYZ=1' /
&PART ID='small roundwood', DRAG_LAW='CYLINDER', SAMPLING_FACTOR=10, SURF_ID='small roundwood', PROP_ID='small roundwood image' QUANTITIES='PARTICLE TEMPERATURE','PARTICLE MASS','PARTICLE DIAMETER', STATIC=.TRUE., COLOR='BROWN' / &INIT PART_ID='small roundwood', XYZ=0.0,0.0,0.15, RADIUS=0.85, HEIGHT=1.9, SHAPE='CONE', N_PARTICLES=10000, MASS_PER_VOLUME=0.4 / &PROP ID='small roundwood image', SMOKEVIEW_ID='TUBE', SMOKEVIEW_PARAMETERS='L=0.10','D=0.002','RANDXYZ=1' /
&PART ID='medium roundwood', DRAG_LAW='CYLINDER', SAMPLING_FACTOR=10, SURF_ID='medium roundwood', PROP_ID='medium roundwood image' QUANTITIES='PARTICLE TEMPERATURE','PARTICLE MASS','PARTICLE DIAMETER', STATIC=.TRUE., COLOR='BROWN' / &INIT PART_ID='medium roundwood', XYZ=0.0,0.0,0.15, RADIUS=0.85, HEIGHT=1.9, SHAPE='CONE', N_PARTICLES=10000, MASS_PER_VOLUME=0.3 / &PROP ID='medium roundwood image', SMOKEVIEW_ID='TUBE', SMOKEVIEW_PARAMETERS='L=0.10','D=0.004','RANDXYZ=1' /
&PART ID='large roundwood', DRAG_LAW='CYLINDER', SAMPLING_FACTOR=10, SURF_ID='large roundwood', PROP_ID='large roundwood image' QUANTITIES='PARTICLE TEMPERATURE','PARTICLE MASS','PARTICLE DIAMETER', STATIC=.TRUE., COLOR='BROWN' / &INIT PART_ID='large roundwood', XYZ=0.0,0.0,0.15, RADIUS=0.85, HEIGHT=1.9, SHAPE='CONE', N_PARTICLES=10000, MASS_PER_VOLUME=0.5 / &PROP ID='large roundwood image', SMOKEVIEW_ID='TUBE', SMOKEVIEW_PARAMETERS='L=0.10','D=0.006','RANDXYZ=1' /
&PART ID='ignitor', SAMPLING_FACTOR=1, SURF_ID='burner', STATIC=.TRUE. / &SURF ID='burner', TMP_FRONT=1300., TAU_T=0.05, RADIUS=0.01, GEOMETRY='SPHERICAL' / &INIT PART_ID='ignitor', XYZ=0.0,0.0,0.21, RADIUS=0.40, SHAPE='RING', N_PARTICLES=400 /
&DEVC XB=-1.5,0.0,-1.5,0.0,0.0,8.0, QUANTITY='MPUV', PART_ID='foliage', ID='foliage1', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=-1.5,0.0,0.0,1.5,0.0,8.0 , QUANTITY='MPUV', PART_ID='foliage', ID='foliage2', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=0.0,1.5,-1.5,0.0,0.0,8.0, QUANTITY='MPUV', PART_ID='foliage', ID='foliage3', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=0.0,1.5,0.0,1.5,0.0,8.0, QUANTITY='MPUV', PART_ID='foliage', ID='foliage4', STATISTICS='VOLUME INTEGRAL' /
&DEVC XB=-1.5,0.0,-1.5,0.0,0.0,8.0, QUANTITY='MPUV', PART_ID='small roundwood', ID='small roundwood1', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=-1.5,0.0,0.0,1.5,0.0,8.0 , QUANTITY='MPUV', PART_ID='small roundwood', ID='small roundwood2', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=0.0,1.5,-1.5,0.0,0.0,8.0, QUANTITY='MPUV', PART_ID='small roundwood', ID='small roundwood3', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=0.0,1.5,0.0,1.5,0.0,8.0, QUANTITY='MPUV', PART_ID='small roundwood', ID='small roundwood4', STATISTICS='VOLUME INTEGRAL' /
&DEVC XB=-1.5,0.0,-1.5,0.0,0.0,8.0, QUANTITY='MPUV', PART_ID='medium roundwood', ID='medium roundwood1', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=-1.5,0.0,0.0,1.5,0.0,8.0 , QUANTITY='MPUV', PART_ID='medium roundwood', ID='medium roundwood2', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=0.0,1.5,-1.5,0.0,0.0,8.0, QUANTITY='MPUV', PART_ID='medium roundwood', ID='medium roundwood3', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=0.0,1.5,0.0,1.5,0.0,8.0, QUANTITY='MPUV', PART_ID='medium roundwood', ID='medium roundwood4', STATISTICS='VOLUME INTEGRAL' /
&DEVC XB=-1.5,0.0,-1.5,0.0,0.0,8.0, QUANTITY='MPUV', PART_ID='large roundwood', ID='large roundwood1', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=-1.5,0.0,0.0,1.5,0.0,8.0 , QUANTITY='MPUV', PART_ID='large roundwood', ID='large roundwood2', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=0.0,1.5,-1.5,0.0,0.0,8.0, QUANTITY='MPUV', PART_ID='large roundwood', ID='large roundwood3', STATISTICS='VOLUME INTEGRAL' / &DEVC XB=0.0,1.5,0.0,1.5,0.0,8.0, QUANTITY='MPUV', PART_ID='large roundwood', ID='large roundwood4', STATISTICS='VOLUME INTEGRAL' /
&SLCF PBY=0., QUANTITY='RADIATION LOSS', CELL_CENTERED=.TRUE. / &SLCF PBY=0., QUANTITY='INTEGRATED INTENSITY', CELL_CENTERED=.TRUE. / &SLCF PBY=0., QUANTITY='TEMPERATURE', CELL_CENTERED=.TRUE. / &SLCF PBY=0., QUANTITY='TEMPERATURE', VECTOR=.TRUE. /
&BNDF QUANTITY='WALL TEMPERATURE' / &BNDF QUANTITY='GAUGE HEAT FLUX' /
&TAIL /