To Extract Slice file data to csv

[2000-lang]

Hello. I am unable to extract instantaneous values for slice files. The model suggested in manual is averaging by default.

For example: For below fds input, if I want soot values of slice Y=18 as per user manual fds2ascii, its giving averaged value for all seconds in that plane, but I need values for 0,5,10,15,20 seconds separately. Any leads would be helpful.

&HEAD CHID='StreetCanyonFire', TITLE='Longitudinal Fire Pollutant Dispersion' /

!--- DOMAIN AND MESH SETUP ---

&MESH IJK=45,9,40, XB=0.0,22.5, 0.0, 4.5, 0.0,40.0 /
&MESH IJK=45,9,40, XB=0.0,22.5, 4.5, 9.0, 0.0,40.0 /
&MESH IJK=45,9,40, XB=0.0,22.5, 9.0, 13.5, 0.0,40.0 /
&MESH IJK=45,9,40, XB=0.0,22.5, 13.5, 18.0, 0.0,40.0 /
&MESH IJK=45,9,40, XB=0.0,22.5, 18.0, 22.5, 0.0,40.0 /
&MESH IJK=45,9,40, XB=0.0,22.5, 22.5, 27.0, 0.0,40.0 /
&MESH IJK=45,9,40, XB=0.0,22.5, 27.0, 31.5, 0.0,40.0 /
&MESH IJK=45,9,40, XB=0.0,22.5, 31.5, 36.0, 0.0,40.0 /

!--- SPECIES DEFINITIONS ---
&SPEC ID='PROPANE' /
&SPEC ID='CARBON MONOXIDE' /
&SPEC ID='CARBON DIOXIDE' /
&SPEC ID='SOOT' /

!--- MATERIAL/SURFACE DEFINITIONS ---
&MATL ID='CONCRETE', SPECIFIC_HEAT=0.88, DENSITY=2100, CONDUCTIVITY=1.4 /
&SURF ID='CONCRETE', MATL_ID='CONCRETE', THICKNESS=0.2 /
&SURF ID='OPEN', COLOR='INVISIBLE' /
&SURF ID='FIRE_POOL', HRRPUA=3500., TMP_FRONT=1200., RAMP_Q='ULTRAFAST', COLOR='RED' /

!--- OBSTACLES SIDE WALLS (adjust to canyon) ---
&OBST XB=0.0,2.25, 0.0,36.0, 0.0,18.0, SURF_ID='CONCRETE' /
&OBST XB=20.25,22.5, 0.0,36.0, 0.0,18.0, SURF_ID='CONCRETE' /

!--- REACTION DEFINITION ---

!--- REAC (Gasoline/Octane) ---
&REAC FUEL='PROPANE', SOOT_YIELD=0.04, CO_YIELD=0.06 /

!--- FIRE SOURCE ---
&OBST XB=10.0,12.5, 15.0,21.0, 0.0,2.8, SURF_ID='FIRE_POOL' /

!------ FIRE GROWTH (t-squared, ramp to 600s) ------
&RAMP ID='ULTRAFAST', T=0.0, F=0.0 /
&RAMP ID='ULTRAFAST', T=20.0, F=0.17 /
&RAMP ID='ULTRAFAST', T=40.0, F=0.33 /
&RAMP ID='ULTRAFAST', T=60.0, F=0.5 /
&RAMP ID='ULTRAFAST', T=80.0, F=0.67 /
&RAMP ID='ULTRAFAST', T=100.0, F=0.83 /
&RAMP ID='ULTRAFAST', T=120.0, F=1.0 /

!------ VENTILATION/WIND ------
&SURF ID='WIND_IN', VEL=6.0, TMP_FRONT=303.0 /
&VENT XB=0.0,22.5, 0.0,0.0, 0.0,40.0, SURF_ID='WIND_IN' /
&VENT XB=0.0,22.5, 36.0,36.0, 0.0,40.0, SURF_ID='OPEN' /
&VENT XB=0.0,22.5, 0.0,36.0, 40.0,40.0, SURF_ID='OPEN' /

&TIME T_END=20. /
&DUMP DT_HRR=5.0, DT_SMOKE3D=5.0 /

!--- SLCF: FIELD (SLICE) OUTPUTS ---
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE', PBY=18.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE', PBY=18.0 /
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT', PBY=18.0 /

!--- DEVCs: LONGITUDINAL PROFILES (x positions, y=18, z=5,10,15m) ---
! --- At z=5m ---
&DEVC ID='CO_x0_z5', XYZ=11.25,6.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE' /
&DEVC ID='CO_x2_z5', XYZ=11.25,12.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE' /
&DEVC ID='CO_x4_z5', XYZ=11.25,17.5,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE' /
&DEVC ID='CO_x6_z5', XYZ=11.25,25.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE' /
&DEVC ID='CO_x8_z5', XYZ=11.25,29.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON MONOXIDE' /

&DEVC ID='CO2_x0_z5', XYZ=11.25,6.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE' /
&DEVC ID='CO2_x2_z5', XYZ=11.25,12.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE' /
&DEVC ID='CO2_x4_z5', XYZ=11.25,17.5,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE' /
&DEVC ID='CO2_x6_z5', XYZ=11.25,25.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE' /
&DEVC ID='CO2_x8_z5', XYZ=11.25,29.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='CARBON DIOXIDE' /

&DEVC ID='SOOT_x0_z5', XYZ=11.25,6.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT' /
&DEVC ID='SOOT_x2_z5', XYZ=11.25,12.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT' /
&DEVC ID='SOOT_x4_z5', XYZ=11.25,17.5,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT' /
&DEVC ID='SOOT_x6_z5', XYZ=11.25,25.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT' /
&DEVC ID='SOOT_x8_z5', XYZ=11.25,29.0,10.0, QUANTITY='VOLUME FRACTION', SPEC_ID='SOOT' /

&TAIL /

[rmcdermo]

I've converted this to a Discussion. By now there are many ways of handling this, but the trend is toward Python-based helper scripts. I use this one that was based on Simo's Matlab version: slread.py

There are also more complete packages in

https://github.com/firemodels/fds/wiki/Third-Party-Tools

For example, pyFDStools by @johodges .

If you do this kind of thing regularly, I suggest you get handy with one of these tools.