Skip to content

Commit 9c87dcb

Browse files
mcgrattamcgratta
authored
Merge pull request #761 from mcgratta/master
Add LNG_MEP scripts and output
2 parents f7d5fcb + 6895d3f commit 9c87dcb

File tree

10 files changed

+508
-110
lines changed

10 files changed

+508
-110
lines changed

‎LNG_MEP/MEPFires_database_spreadsheet_v2.xlsx‎

5.86 MB
Binary file not shown.

‎LNG_MEP/README.md‎

Lines changed: 24 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,24 @@
1+
## Model Evaluation Protocol (MEP) for Fire Models Involving Fuels at Liquefied Natural Gas (LNG) Facilities
2+
3+
Under the sponsorship of the Pipeline and Hazardous Materials Safety Administration (PHMSA) of the U.S. Department of Transportation, Sandia National Laboratories has compiled a collection of experimental datasets for use in evaluating models that predict the consequences of accidental fires at LNG facilities. Details can be found in the report written by Anay Luketa, _Model Evaluation Protocol for Fire Models Involving Fuels at Liquefied Natural Gas Facilities (Version 2)_, Sandia Report SAND2026-19309, April 2026.
4+
5+
The experiments included in this database are also included in the FDS validation suite and the results are reported in the FDS Validation Guide. In addition, Sandia has created an Excel spreadsheet with the basic parameters of each experiment, which is included in this directory under the name `MEPFires_database_spreadsheet_v2.xlsx`. Modelers are expected to report their results in the various sheets that are included in the file.
6+
7+
This directory contains a Python script (`csv2xlsx.py`) that reads the comma-delimited (`.csv`) FDS output files that reside in the `firemodels/out` repository, and transfers snippets of data into the Excel file. A file called `csv2xlsx_parameters.csv` contains a list of the data snippets to be transferred from the FDS output files to the Excel file. Note that the original Excel file is not changed after running the script; but rather a copy is created called `MEPFires_database_spreadsheet_v2_filled.xlsx` which contains the model results.
8+
9+
The columns of `csv2xlsx_parameters.csv` are as follows:
10+
11+
A. The FDS output file
12+
13+
B. The starting row of the FDS data
14+
15+
C. The starting column of the FDS data
16+
17+
D. The number of values, with a negative sign added if the data is to be read across rather than down
18+
19+
E. The name of the tab in the Excel file where the data is to be copied
20+
21+
F. The starting row of the data
22+
23+
G. The starting column of the data
24+

‎LNG_MEP/csv2xlsx.py‎

Lines changed: 49 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,49 @@
1+
"""csv2xlsx cut and paste
2+
3+
This script opens a series of comma-delimited .csv files and copies portions of columns into specified
4+
worksheets of a specified Excel (.xlsx) file. It saves the Excel file with a new name. The names of the
5+
new and old Excel files are currently hardwired in the script below.
6+
7+
The names of the csv files are listed in the file called 'csv2xlxs_parameters.csv'. This file contains the full
8+
path to the csv files, followed by the row/column pair where the data starts, the length of the column
9+
snippet, the name of the worksheet in the Excel file, and the row/column pair where the snippet is to be
10+
pasted.
11+
12+
Note that the module used to save the new Excel file does not support wmf (Windows MetaFile) images. Any
13+
wmf images in the original Excel file will be dropped from the new Excel file.
14+
15+
This script requires the Excel parsing module 'openpyxl' and the csv reader/writer 'csv'.
16+
"""
17+
18+
import csv
19+
from openpyxl import load_workbook
20+
21+
wb = load_workbook('MEPFires_database_spreadsheet_v2.xlsx')
22+
23+
outdir = '../'
24+
param_file = open('csv2xlsx_parameters.csv','r')
25+
param_reader = csv.reader(param_file,delimiter=',')
26+
params = list(param_reader)
27+
28+
for i in range(0,len(params)):
29+
30+
orig_row = int(params[i][1]) - 1
31+
orig_col = int(params[i][2]) - 1
32+
csv_file = open(outdir+params[i][0],'r')
33+
data = list(csv.reader(csv_file,delimiter=','))
34+
new_row = int(params[i][5])
35+
new_col = int(params[i][6])
36+
length = int(params[i][3])
37+
print(i,new_row,new_col)
38+
39+
ws = wb[params[i][4]]
40+
41+
for counter in range(0,abs(length)):
42+
if length > 0:
43+
_ = ws.cell(row=new_row+counter, column=new_col, value=round(float(data[orig_row+counter][orig_col]),2))
44+
else:
45+
_ = ws.cell(row=new_row+counter, column=new_col, value=round(float(data[orig_row][orig_col+counter]),2))
46+
47+
csv_file.close()
48+
49+
wb.save('MEPFires_database_spreadsheet_v2_filled.xlsx')

‎LNG_MEP/csv2xlsx_parameters.csv‎

Lines changed: 321 additions & 0 deletions
Large diffs are not rendered by default.

‎Sandia_Crude_Oil/Tight_1_devc.csv‎

Lines changed: 9 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -4,8 +4,12 @@ Time,NW-5 Narrow,NE-5 Narrow,SW-5 Narrow,SE-5 Narrow,NW-5 Wide,SW-5 Wide,SE-5 Wi
44
5.0039436E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 7.3456631E-001, 4.7358862E+000, 0.0000000E+000, 0.0000000E+000
55
1.0005257E+001, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 1.9977767E+000, 7.6930487E+000, 0.0000000E+000, 0.0000000E+000
66
1.5002217E+001, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 0.0000000E+000, 2.4429347E+000, 7.2746015E+000, 0.0000000E+000, 0.0000000E+000
7-
2.0000533E+001, 3.8953414E-001, 7.1323025E-001, 4.7228516E-001, 3.6932019E-001, 1.0048591E-001, 5.3475540E-002, 6.1674106E-002, 9.0493226E-003, 5.3435040E-003, 1.1728330E-001, 3.5655257E-001, 1.6302175E-001, 1.0657898E-001, 1.1103471E-001, 3.8488820E-001, 2.1813011E-001, 8.1192108E-002, 1.5473391E-001, 3.5334984E-001, 3.4994583E-001, 9.4013192E-002, 1.9911982E+000, 9.2140454E+000, 1.3677130E+001, 1.4912292E+001
8-
2.5000071E+001, 6.5013430E+001, 1.0148553E+002, 8.3098976E+001, 6.2909361E+001, 1.6268733E+001, 8.5850998E+000, 1.0238768E+001, 1.2630733E+000, 8.2203610E-001, 2.0104525E+001, 5.3387090E+001, 3.1421804E+001, 1.3596531E+001, 1.6722492E+001, 5.1059371E+001, 2.6364234E+001, 1.1088429E+001, 2.4108082E+001, 4.4429973E+001, 4.4573576E+001, 8.1265006E+000, 2.9334383E+000, 8.6561746E+000, 5.5399045E+000, 1.0644126E+001
9-
3.0000505E+001, 6.4195341E+001, 1.0634325E+002, 8.5400516E+001, 6.2908626E+001, 1.6317453E+001, 8.7344407E+000, 1.0569504E+001, 1.2747990E+000, 8.1137452E-001, 2.3136838E+001, 5.6272215E+001, 3.2277006E+001, 1.4374526E+001, 1.9336048E+001, 5.2299439E+001, 2.6897936E+001, 1.1512172E+001, 2.5815203E+001, 4.4084935E+001, 4.5061377E+001, 7.9693073E+000, 3.5153155E+000, 8.2023085E+000, 5.6110304E+000, 1.0736744E+001
10-
3.5000985E+001, 6.5458595E+001, 1.0365769E+002, 8.3030547E+001, 6.2586339E+001, 1.5824869E+001, 8.5424785E+000, 9.9979658E+000, 1.3023844E+000, 8.3494577E-001, 2.0754290E+001, 5.4062378E+001, 3.3574435E+001, 1.4210844E+001, 1.8851654E+001, 5.2750103E+001, 2.9211036E+001, 1.2489011E+001, 2.5537057E+001, 4.6457353E+001, 4.6362158E+001, 9.6995364E+000, 2.2354883E+000, 8.9315234E+000, 5.6054248E+000, 1.0834983E+001
11-
4.0000200E+001, 6.4252883E+001, 9.4234473E+001, 7.9332097E+001, 6.1954870E+001, 1.5990172E+001, 8.5138740E+000, 9.8380201E+000, 1.4329775E+000, 8.4741002E-001, 2.1243452E+001, 5.6423667E+001, 3.2313132E+001, 1.4761171E+001, 1.7775856E+001, 5.3683234E+001, 2.6105249E+001, 1.1815669E+001, 2.3114971E+001, 4.2227893E+001, 4.3221956E+001, 7.9269464E+000, 1.1368004E+000, 9.2607447E+000, 5.6364733E+000, 1.0792427E+001
7+
2.0000533E+001, 3.8953414E-001, 7.1323025E-001, 4.7228516E-001, 3.6932019E-001, 1.0048591E-001, 5.3475540E-002, 6.1674106E-002, 9.0493226E-003, 5.3435040E-003, 1.1728330E-001, 3.5655257E-001, 1.6302175E-001, 1.0657898E-001, 1.1103471E-001, 3.8488820E-001, 2.1813011E-001, 8.1192108E-002, 1.5473391E-001, 3.5334984E-001, 3.4994583E-001, 9.4013192E-002, 1.9911982E+000, 9.2140454E+000, 2.3601312E+001, 5.5945106E+001
8+
2.5000071E+001, 6.5013430E+001, 1.0148553E+002, 8.3098976E+001, 6.2909361E+001, 1.6268733E+001, 8.5850998E+000, 1.0238768E+001, 1.2630733E+000, 8.2203610E-001, 2.0104525E+001, 5.3387090E+001, 3.1421804E+001, 1.3596531E+001, 1.6722492E+001, 5.1059371E+001, 2.6364234E+001, 1.1088429E+001, 2.4108082E+001, 4.4429973E+001, 4.4573576E+001, 8.1265006E+000, 2.9334383E+000, 8.6561746E+000, 7.9449223E+000, 4.5411656E+001
9+
3.0000505E+001, 6.4195341E+001, 1.0634325E+002, 8.5400516E+001, 6.2908626E+001, 1.6317453E+001, 8.7344407E+000, 1.0569504E+001, 1.2747990E+000, 8.1137452E-001, 2.3136838E+001, 5.6272215E+001, 3.2277006E+001, 1.4374526E+001, 1.9336048E+001, 5.2299439E+001, 2.6897936E+001, 1.1512172E+001, 2.5815203E+001, 4.4084935E+001, 4.5061377E+001, 7.9693073E+000, 3.5153155E+000, 8.2023085E+000, 8.0775574E+000, 4.5649551E+001
10+
3.5000985E+001, 6.5458595E+001, 1.0365769E+002, 8.3030547E+001, 6.2586339E+001, 1.5824869E+001, 8.5424785E+000, 9.9979658E+000, 1.3023844E+000, 8.3494577E-001, 2.0754290E+001, 5.4062378E+001, 3.3574435E+001, 1.4210844E+001, 1.8851654E+001, 5.2750103E+001, 2.9211036E+001, 1.2489011E+001, 2.5537057E+001, 4.6457353E+001, 4.6362158E+001, 9.6995364E+000, 2.2354883E+000, 8.9315234E+000, 8.1327250E+000, 4.5885036E+001
11+
4.0000200E+001, 6.4252883E+001, 9.4234473E+001, 7.9332097E+001, 6.1954870E+001, 1.5990172E+001, 8.5138740E+000, 9.8380201E+000, 1.4329775E+000, 8.4741002E-001, 2.1243452E+001, 5.6423667E+001, 3.2313132E+001, 1.4761171E+001, 1.7775856E+001, 5.3683234E+001, 2.6105249E+001, 1.1815669E+001, 2.3114971E+001, 4.2227893E+001, 4.3221956E+001, 7.9269464E+000, 1.1368004E+000, 9.2607447E+000, 8.1567565E+000, 4.5706261E+001
12+
4.5002569E+001, 6.3123084E+001, 1.0746820E+002, 8.3616899E+001, 6.2367000E+001, 1.6090065E+001, 8.4825455E+000, 1.0037913E+001, 1.2895400E+000, 8.2094707E-001, 2.2956375E+001, 5.9394823E+001, 3.0317639E+001, 1.4053115E+001, 1.8685205E+001, 5.5644173E+001, 2.4297841E+001, 1.1167887E+001, 2.4731814E+001, 4.1928126E+001, 4.1882458E+001, 8.1784513E+000, 2.5351524E+000, 8.3464519E+000, 8.1409563E+000, 4.5425821E+001
13+
5.0001883E+001, 5.9790499E+001, 9.5498819E+001, 8.0033206E+001, 6.0382784E+001, 1.5673474E+001, 8.2999189E+000, 9.9699441E+000, 1.3139709E+000, 8.1683835E-001, 2.0620116E+001, 5.3179047E+001, 2.9252900E+001, 1.3205804E+001, 1.6942120E+001, 4.9155577E+001, 2.3085146E+001, 9.9545664E+000, 2.2842243E+001, 3.8438729E+001, 4.0857749E+001, 6.0636843E+000, 3.0955225E+000, 8.5154404E+000, 8.1458621E+000, 4.5479319E+001
14+
5.5005374E+001, 6.4644975E+001, 1.0183265E+002, 8.3097970E+001, 6.4866260E+001, 1.5828748E+001, 8.5170093E+000, 1.0282018E+001, 1.2718971E+000, 8.2567698E-001, 2.0787689E+001, 5.3595885E+001, 3.0904758E+001, 1.1270161E+001, 1.8644559E+001, 5.2129787E+001, 2.5042823E+001, 9.4624996E+000, 2.5293646E+001, 4.4044493E+001, 4.2280833E+001, 7.0718507E+000, 3.4429698E+000, 8.8620516E+000, 8.1472432E+000, 4.5635414E+001
15+
6.0000000E+001, 6.7051615E+001, 9.1004632E+001, 8.0017463E+001, 6.3179737E+001, 1.6109471E+001, 8.5623515E+000, 9.9311129E+000, 1.3836203E+000, 8.3891413E-001, 1.9624455E+001, 5.4873564E+001, 3.4312193E+001, 1.3537322E+001, 1.7678391E+001, 5.2181980E+001, 2.8458217E+001, 1.1850623E+001, 2.2911631E+001, 4.3765014E+001, 4.6614699E+001, 9.1173102E+000, 1.7259989E+000, 9.6814377E+000, 8.1856904E+000, 4.5806658E+001

‎Sandia_Crude_Oil/Tight_1_line.csv‎

Lines changed: 45 additions & 45 deletions
Original file line numberDiff line numberDiff line change
@@ -1,50 +1,50 @@
11
m,kW
22
z_HRRPUL,HRRPUL
3-
1.0000000E-001, 6.0792330E+003
4-
3.0000000E-001, 6.8716413E+003
5-
5.0000000E-001, 4.4163806E+003
6-
7.0000000E-001, 3.3413156E+003
7-
9.0000000E-001, 2.7821482E+003
8-
1.1000000E+000, 2.5076621E+003
9-
1.3000000E+000, 2.3341797E+003
10-
1.5000000E+000, 2.2058345E+003
11-
1.7000000E+000, 2.0990427E+003
12-
1.9000000E+000, 2.0058213E+003
13-
2.1000000E+000, 1.9129193E+003
14-
2.3000000E+000, 1.8122552E+003
15-
2.5000000E+000, 1.7062831E+003
16-
2.7000000E+000, 1.5979429E+003
17-
2.9000000E+000, 1.4955761E+003
18-
3.1000000E+000, 1.3980059E+003
19-
3.3000000E+000, 1.2926495E+003
20-
3.5000000E+000, 1.1869170E+003
21-
3.7000000E+000, 1.0934537E+003
22-
3.9000000E+000, 1.0036475E+003
23-
4.1000000E+000, 9.2322914E+002
24-
4.3000000E+000, 8.4183337E+002
25-
4.5000000E+000, 7.5849267E+002
26-
4.7000000E+000, 6.7312815E+002
27-
4.9000000E+000, 5.8815172E+002
28-
5.1000000E+000, 5.0803770E+002
29-
5.3000000E+000, 4.3583474E+002
30-
5.5000000E+000, 3.7124663E+002
31-
5.7000000E+000, 3.1636054E+002
32-
5.9000000E+000, 2.6991044E+002
33-
6.1000000E+000, 2.2944237E+002
34-
6.3000000E+000, 1.9413905E+002
35-
6.5000000E+000, 1.6345578E+002
36-
6.7000000E+000, 1.3724127E+002
37-
6.9000000E+000, 1.1512073E+002
38-
7.1000000E+000, 9.7246402E+001
39-
7.3000000E+000, 8.1601492E+001
40-
7.5000000E+000, 6.7033686E+001
41-
7.7000000E+000, 5.3910410E+001
42-
7.9000000E+000, 4.3378828E+001
43-
8.1000000E+000, 3.4567037E+001
44-
8.3000000E+000, 2.7723882E+001
45-
8.5000000E+000, 2.1974080E+001
46-
8.7000000E+000, 1.7959926E+001
47-
8.9000000E+000, 1.0809763E+001
3+
1.0000000E-001, 6.0473758E+003
4+
3.0000000E-001, 6.8503660E+003
5+
5.0000000E-001, 4.4379217E+003
6+
7.0000000E-001, 3.3577145E+003
7+
9.0000000E-001, 2.7883215E+003
8+
1.1000000E+000, 2.5267203E+003
9+
1.3000000E+000, 2.3483102E+003
10+
1.5000000E+000, 2.2109834E+003
11+
1.7000000E+000, 2.1026407E+003
12+
1.9000000E+000, 2.0066074E+003
13+
2.1000000E+000, 1.9104795E+003
14+
2.3000000E+000, 1.8162946E+003
15+
2.5000000E+000, 1.7204373E+003
16+
2.7000000E+000, 1.6124789E+003
17+
2.9000000E+000, 1.5083478E+003
18+
3.1000000E+000, 1.4051038E+003
19+
3.3000000E+000, 1.2989200E+003
20+
3.5000000E+000, 1.1940242E+003
21+
3.7000000E+000, 1.0948292E+003
22+
3.9000000E+000, 9.9815829E+002
23+
4.1000000E+000, 9.0722693E+002
24+
4.3000000E+000, 8.2005873E+002
25+
4.5000000E+000, 7.3523596E+002
26+
4.7000000E+000, 6.5362394E+002
27+
4.9000000E+000, 5.7317063E+002
28+
5.1000000E+000, 4.9652933E+002
29+
5.3000000E+000, 4.2572208E+002
30+
5.5000000E+000, 3.6312654E+002
31+
5.7000000E+000, 3.0992443E+002
32+
5.9000000E+000, 2.6400661E+002
33+
6.1000000E+000, 2.2341493E+002
34+
6.3000000E+000, 1.8843805E+002
35+
6.5000000E+000, 1.5867672E+002
36+
6.7000000E+000, 1.3343598E+002
37+
6.9000000E+000, 1.1193885E+002
38+
7.1000000E+000, 9.4169571E+001
39+
7.3000000E+000, 7.8273608E+001
40+
7.5000000E+000, 6.4730681E+001
41+
7.7000000E+000, 5.2897279E+001
42+
7.9000000E+000, 4.2952847E+001
43+
8.1000000E+000, 3.4449618E+001
44+
8.3000000E+000, 2.7617331E+001
45+
8.5000000E+000, 2.1853069E+001
46+
8.7000000E+000, 1.7856227E+001
47+
8.9000000E+000, 1.0606370E+001
4848
NaN,NaN
4949
NaN,NaN
5050
NaN,NaN

0 commit comments

Comments
 (0)