JIANG/EPANET
1
0
Fork 0
Code Issues Pull Requests Actions 4 Packages Projects Releases Wiki Activity

feat: Add column headings for sections in INP file

Browse Source 
This commit includes the following changes:
- Column headings have been added for each section in the INP file.
- The alignment of the values has been updated to match these headings.

These enhancements improve readability and understanding of the INP file.
This commit is contained in:
Yu Chun Tsao
2024-06-15 23:48:14 +08:00
parent 109a93f7d5
commit 317990a8d1
1 changed files with 76 additions and 46 deletions
Download Patch File Download Diff File Expand all files Collapse all files

122
src/inpfile.c
View File

@@ -152,40 +152,46 @@ int saveinpfile(Project *pr, const char *fname)
// (Leave demands for [DEMANDS] section) // (Leave demands for [DEMANDS] section)
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_JUNCTIONS); fprintf(f, s_JUNCTIONS);
fprintf(f, "\n;%-31s\t%-12s\t%-12s\t%-31s",
"ID", "Elev", "Demand", "Pattern");
for (i = 1; i <= net->Njuncs; i++) for (i = 1; i <= net->Njuncs; i++)
{ {
node = &net->Node[i]; node = &net->Node[i];
fprintf(f, "\n %-31s %12.4f", node->ID, node->El * pr->Ucf[ELEV]); fprintf(f, "\n %-31s\t%-12.4f", node->ID, node->El * pr->Ucf[ELEV]);
if (node->Comment) fprintf(f, " ;%s", node->Comment); if (node->Comment) fprintf(f, "\t;%s", node->Comment);
} }
// Write [RESERVOIRS] section // Write [RESERVOIRS] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_RESERVOIRS); fprintf(f, s_RESERVOIRS);
fprintf(f, "\n;%-31s\t%-12s\t%-31s",
"ID", "Head", "Pattern");
for (i = 1; i <= net->Ntanks; i++) for (i = 1; i <= net->Ntanks; i++)
{ {
tank = &net->Tank[i]; tank = &net->Tank[i];
if (tank->A == 0.0) if (tank->A == 0.0)
{ {
node = &net->Node[tank->Node]; node = &net->Node[tank->Node];
sprintf(s, " %-31s %12.4f", node->ID, node->El * pr->Ucf[ELEV]); sprintf(s, " %-31s\t%-12.4f", node->ID, node->El * pr->Ucf[ELEV]);
if ((j = tank->Pat) > 0) sprintf(s1, " %s", net->Pattern[j].ID); if ((j = tank->Pat) > 0) sprintf(s1, "%s", net->Pattern[j].ID);
else strcpy(s1, " "); else strcpy(s1, " ");
fprintf(f, "\n%s %-31s", s, s1); fprintf(f, "\n%s\t%-31s", s, s1);
if (node->Comment) fprintf(f, " ;%s", node->Comment); if (node->Comment) fprintf(f, "\t;%s", node->Comment);
} }
} }
// Write [TANKS] section // Write [TANKS] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_TANKS); fprintf(f, s_TANKS);
fprintf(f, "\n;%-31s\t%-12s\t%-12s\t%-12s\t%-12s\t%-12s\t%-12s\t%-31s\t%-12s",
"ID", "Elevation", "InitLevel", "MinLevel", "MaxLevel", "Diameter", "MinVol", "VolCurve", "Overflow");
for (i = 1; i <= net->Ntanks; i++) for (i = 1; i <= net->Ntanks; i++)
{ {
tank = &net->Tank[i]; tank = &net->Tank[i];
if (tank->A > 0.0) if (tank->A > 0.0)
{ {
node = &net->Node[tank->Node]; node = &net->Node[tank->Node];
sprintf(s, " %-31s %12.4f %12.4f %12.4f %12.4f %12.4f %12.4f", sprintf(s, " %-31s\t%-12.4f\t%-12.4f\t%-12.4f\t%-12.4f\t%-12.4f\t%-12.4f",
node->ID, node->El * pr->Ucf[ELEV], node->ID, node->El * pr->Ucf[ELEV],
(tank->H0 - node->El) * pr->Ucf[ELEV], (tank->H0 - node->El) * pr->Ucf[ELEV],
(tank->Hmin - node->El) * pr->Ucf[ELEV], (tank->Hmin - node->El) * pr->Ucf[ELEV],
@@ -195,15 +201,17 @@ int saveinpfile(Project *pr, const char *fname)
if ((j = tank->Vcurve) > 0) sprintf(s1, "%s", net->Curve[j].ID); if ((j = tank->Vcurve) > 0) sprintf(s1, "%s", net->Curve[j].ID);
else if (tank->CanOverflow) strcpy(s1, "*"); else if (tank->CanOverflow) strcpy(s1, "*");
else strcpy(s1, " "); else strcpy(s1, " ");
fprintf(f, "\n%s %-31s", s, s1); fprintf(f, "\n%s\t%-31s", s, s1);
if (tank->CanOverflow) fprintf(f, " YES "); if (tank->CanOverflow) fprintf(f, "\t%-12s", "YES");
if (node->Comment) fprintf(f, " ;%s", node->Comment); if (node->Comment) fprintf(f, "\t;%s", node->Comment);
} }
} }
// Write [PIPES] section // Write [PIPES] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_PIPES); fprintf(f, s_PIPES);
fprintf(f, "\n;%-31s\t%-31s\t%-31s\t%-12s\t%-12s\t%-12s\t%-12s\t%-6s",
"ID", "Node1", "Node2", "Length", "Diameter", "Roughness", "MinorLoss", "Status");
for (i = 1; i <= net->Nlinks; i++) for (i = 1; i <= net->Nlinks; i++)
{ {
link = &net->Link[i]; link = &net->Link[i];
@@ -214,36 +222,38 @@ int saveinpfile(Project *pr, const char *fname)
if (hyd->Formflag == DW) kc = kc * pr->Ucf[ELEV] * 1000.0; if (hyd->Formflag == DW) kc = kc * pr->Ucf[ELEV] * 1000.0;
km = link->Km * SQR(d) * SQR(d) / 0.02517; km = link->Km * SQR(d) * SQR(d) / 0.02517;
sprintf(s, " %-31s %-31s %-31s %12.4f %12.4f %12.4f %12.4f", sprintf(s, " %-31s\t%-31s\t%-31s\t%-12.4f\t%-12.4f\t%-12.4f\t%-12.4f",
link->ID, net->Node[link->N1].ID, net->Node[link->N2].ID, link->ID, net->Node[link->N1].ID, net->Node[link->N2].ID,
link->Len * pr->Ucf[LENGTH], d * pr->Ucf[DIAM], kc, km); link->Len * pr->Ucf[LENGTH], d * pr->Ucf[DIAM], kc, km);
if (link->Type == CVPIPE) sprintf(s2, "CV"); if (link->Type == CVPIPE) sprintf(s2, "CV");
else if (link->Status == CLOSED) sprintf(s2, "CLOSED"); else if (link->Status == CLOSED) sprintf(s2, "CLOSED");
else strcpy(s2, " "); else strcpy(s2, " ");
fprintf(f, "\n%s %-6s", s, s2); fprintf(f, "\n%s\t%-6s", s, s2);
if (link->Comment) fprintf(f, " ;%s", link->Comment); if (link->Comment) fprintf(f, "\t;%s", link->Comment);
} }
} }
// Write [PUMPS] section // Write [PUMPS] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_PUMPS); fprintf(f, s_PUMPS);
fprintf(f, "\n;%-31s\t%-31s\t%-31s\t%-12s",
"ID", "Node1", "Node2", "Parameters");
for (i = 1; i <= net->Npumps; i++) for (i = 1; i <= net->Npumps; i++)
{ {
n = net->Pump[i].Link; n = net->Pump[i].Link;
link = &net->Link[n]; link = &net->Link[n];
pump = &net->Pump[i]; pump = &net->Pump[i];
sprintf(s, " %-31s %-31s %-31s", link->ID, net->Node[link->N1].ID, sprintf(s, " %-31s\t%-31s\t%-31s", link->ID, net->Node[link->N1].ID,
net->Node[link->N2].ID); net->Node[link->N2].ID);
// Pump has constant power // Pump has constant power
if (pump->Ptype == CONST_HP) sprintf(s1, " POWER %.4f", link->Km); if (pump->Ptype == CONST_HP) sprintf(s1, "\tPOWER %.4f", link->Km);
// Pump has a head curve // Pump has a head curve
else if ((j = pump->Hcurve) > 0) else if ((j = pump->Hcurve) > 0)
{ {
sprintf(s1, " HEAD %s", net->Curve[j].ID); sprintf(s1, "\tHEAD %s", net->Curve[j].ID);
} }
// Old format used for pump curve // Old format used for pump curve
@@ -260,25 +270,26 @@ int saveinpfile(Project *pr, const char *fname)
// Optional speed pattern // Optional speed pattern
if ((j = pump->Upat) > 0) if ((j = pump->Upat) > 0)
{ {
sprintf(s1, " PATTERN %s", net->Pattern[j].ID); sprintf(s1, "\tPATTERN %s", net->Pattern[j].ID);
strcat(s, s1); strcat(s, s1);
} }
// Optional speed setting // Optional speed setting
if (link->Kc != 1.0) if (link->Kc != 1.0)
{ {
sprintf(s1, " SPEED %.4f", link->Kc); sprintf(s1, "\tSPEED %.4f", link->Kc);
strcat(s, s1); strcat(s, s1);
} }
fprintf(f, "\n%s", s); fprintf(f, "\n%s", s);
if (link->Comment) fprintf(f, " ;%s", link->Comment); if (link->Comment) fprintf(f, "\t;%s", link->Comment);
} }
// Write [VALVES] section // Write [VALVES] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_VALVES); fprintf(f, s_VALVES);
fprintf(f, "\n;%-31s\t%-31s\t%-31s\t%-12s\t%-6s\t%-12s\t%-12s",
"ID", "Node1", "Node2", "Diameter", "Type", "Setting", "MinorLoss");
for (i = 1; i <= net->Nvalves; i++) for (i = 1; i <= net->Nvalves; i++)
{ {
n = net->Valve[i].Link; n = net->Valve[i].Link;
@@ -303,7 +314,7 @@ int saveinpfile(Project *pr, const char *fname)
} }
km = link->Km * SQR(d) * SQR(d) / 0.02517; km = link->Km * SQR(d) * SQR(d) / 0.02517;
sprintf(s, " %-31s %-31s %-31s %12.4f %5s", sprintf(s, " %-31s\t%-31s\t%-31s\t%-12.4f\t%-6s",
link->ID, net->Node[link->N1].ID, link->ID, net->Node[link->N1].ID,
net->Node[link->N2].ID, d * pr->Ucf[DIAM], net->Node[link->N2].ID, d * pr->Ucf[DIAM],
LinkTxt[link->Type]); LinkTxt[link->Type]);
@@ -311,22 +322,24 @@ int saveinpfile(Project *pr, const char *fname)
// For GPV, setting = head curve index // For GPV, setting = head curve index
if (link->Type == GPV && (j = ROUND(link->Kc)) > 0) if (link->Type == GPV && (j = ROUND(link->Kc)) > 0)
{ {
sprintf(s1, "%-31s %12.4f", net->Curve[j].ID, km); sprintf(s1, "%-31s\t%-12.4f", net->Curve[j].ID, km);
} }
// For PCV add loss curve if present // For PCV add loss curve if present
else if (link->Type == PCV && (j = net->Valve[i].Curve) > 0) else if (link->Type == PCV && (j = net->Valve[i].Curve) > 0)
{ {
sprintf(s1, "%12.4f %12.4f %-31s", kc, km, net->Curve[j].ID); sprintf(s1, "%-12.4f\t%-12.4f\t%-31s", kc, km, net->Curve[j].ID);
} }
else sprintf(s1, "%12.4f %12.4f", kc, km); else sprintf(s1, "%-12.4f\t%-12.4f", kc, km);
fprintf(f, "\n%s %s", s, s1); fprintf(f, "\n%s\t%s", s, s1);
if (link->Comment) fprintf(f, " ;%s", link->Comment); if (link->Comment) fprintf(f, "\t;%s", link->Comment);
} }
// Write [DEMANDS] section // Write [DEMANDS] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_DEMANDS); fprintf(f, s_DEMANDS);
fprintf(f, "\n;%-31s\t%-14s\t%-31s\t%-31s",
"Junction", "Demand", "Pattern", "Category");
ucf = pr->Ucf[DEMAND]; ucf = pr->Ucf[DEMAND];
for (i = 1; i <= net->Njuncs; i++) for (i = 1; i <= net->Njuncs; i++)
{ {
@@ -334,11 +347,11 @@ int saveinpfile(Project *pr, const char *fname)
for (demand = node->D; demand != NULL; demand = demand->next) for (demand = node->D; demand != NULL; demand = demand->next)
{ {
if (demand->Base == 0.0) continue; if (demand->Base == 0.0) continue;
sprintf(s, " %-31s %14.6f", node->ID, ucf * demand->Base); sprintf(s, " %-31s\t%-14.6f", node->ID, ucf * demand->Base);
if ((j = demand->Pat) > 0) sprintf(s1, " %-31s", net->Pattern[j].ID); if ((j = demand->Pat) > 0) sprintf(s1, "%-31s", net->Pattern[j].ID);
else strcpy(s1, " "); else strcpy(s1, " ");
fprintf(f, "\n%s %-31s", s, s1); fprintf(f, "\n%s\t%-31s", s, s1);
if (demand->Name) fprintf(f, " ;%s", demand->Name); if (demand->Name) fprintf(f, "\t;%s", demand->Name);
} }
} }
@@ -346,17 +359,21 @@ int saveinpfile(Project *pr, const char *fname)
// Write [EMITTERS] section // Write [EMITTERS] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_EMITTERS); fprintf(f, s_EMITTERS);
fprintf(f, "\n;%-31s\t%-14s",
"Junction", "Coefficient");
for (i = 1; i <= net->Njuncs; i++) for (i = 1; i <= net->Njuncs; i++)
{ {
node = &net->Node[i]; node = &net->Node[i];
if (node->Ke == 0.0) continue; if (node->Ke == 0.0) continue;
ke = pr->Ucf[FLOW] / pow(pr->Ucf[PRESSURE] * node->Ke, (1.0 / hyd->Qexp)); ke = pr->Ucf[FLOW] / pow(pr->Ucf[PRESSURE] * node->Ke, (1.0 / hyd->Qexp));
fprintf(f, "\n %-31s %14.6f", node->ID, ke); fprintf(f, "\n %-31s\t%-14.6f", node->ID, ke);
} }
// Write [STATUS] section // Write [STATUS] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_STATUS); fprintf(f, s_STATUS);
fprintf(f, "\n;%-31s\t%-12s",
"ID", "Status/Setting");
for (i = 1; i <= net->Nlinks; i++) for (i = 1; i <= net->Nlinks; i++)
{ {
link = &net->Link[i]; link = &net->Link[i];
@@ -364,7 +381,7 @@ int saveinpfile(Project *pr, const char *fname)
{ {
if (link->Status == CLOSED) if (link->Status == CLOSED)
{ {
fprintf(f, "\n %-31s %s", link->ID, StatTxt[CLOSED]); fprintf(f, "\n %-31s\t%s", link->ID, StatTxt[CLOSED]);
} }
// Write pump speed here for pumps with old-style pump curve input // Write pump speed here for pumps with old-style pump curve input
@@ -375,7 +392,7 @@ int saveinpfile(Project *pr, const char *fname)
if (pump->Hcurve == 0 && pump->Ptype != CONST_HP && if (pump->Hcurve == 0 && pump->Ptype != CONST_HP &&
link->Kc != 1.0) link->Kc != 1.0)
{ {
fprintf(f, "\n %-31s %-.4f", link->ID, link->Kc); fprintf(f, "\n %-31s\t%-.4f", link->ID, link->Kc);
} }
} }
} }
@@ -385,11 +402,11 @@ int saveinpfile(Project *pr, const char *fname)
{ {
if (link->Status == OPEN) if (link->Status == OPEN)
{ {
fprintf(f, "\n %-31s %s", link->ID, StatTxt[OPEN]); fprintf(f, "\n %-31s\t%s", link->ID, StatTxt[OPEN]);
} }
if (link->Status == CLOSED) if (link->Status == CLOSED)
{ {
fprintf(f, "\n%-31s %s", link->ID, StatTxt[CLOSED]); fprintf(f, "\n%-31s\t%s", link->ID, StatTxt[CLOSED]);
} }
} }
} }
@@ -398,26 +415,30 @@ int saveinpfile(Project *pr, const char *fname)
// (Use 6 pattern factors per line) // (Use 6 pattern factors per line)
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_PATTERNS); fprintf(f, s_PATTERNS);
fprintf(f, "\n;%-31s\t%-12s",
"ID", "Multipliers");
for (i = 1; i <= net->Npats; i++) for (i = 1; i <= net->Npats; i++)
{ {
if (net->Pattern[i].Comment) fprintf(f, "\n;%s", net->Pattern[i].Comment); if (net->Pattern[i].Comment) fprintf(f, "\n;%s", net->Pattern[i].Comment);
for (j = 0; j < net->Pattern[i].Length; j++) for (j = 0; j < net->Pattern[i].Length; j++)
{ {
if (j % 6 == 0) fprintf(f, "\n %-31s", net->Pattern[i].ID); if (j % 6 == 0) fprintf(f, "\n %-31s", net->Pattern[i].ID);
fprintf(f, " %12.4f", net->Pattern[i].F[j]); fprintf(f, "\t%-12.4f", net->Pattern[i].F[j]);
} }
} }
// Write [CURVES] section // Write [CURVES] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_CURVES); fprintf(f, s_CURVES);
fprintf(f, "\n;%-31s\t%-12s\t%-12s",
"ID", "X-Value", "Y-Value");
for (i = 1; i <= net->Ncurves; i++) for (i = 1; i <= net->Ncurves; i++)
{ {
if (net->Curve[i].Comment) fprintf(f, "\n;%s", net->Curve[i].Comment); if (net->Curve[i].Comment) fprintf(f, "\n;%s", net->Curve[i].Comment);
for (j = 0; j < net->Curve[i].Npts; j++) for (j = 0; j < net->Curve[i].Npts; j++)
{ {
curve = &net->Curve[i]; curve = &net->Curve[i];
fprintf(f, "\n %-31s %12.4f %12.4f", curve->ID, curve->X[j], curve->Y[j]); fprintf(f, "\n %-31s\t%-12.4f\t%-12.4f", curve->ID, curve->X[j], curve->Y[j]);
} }
} }
@@ -502,39 +523,42 @@ int saveinpfile(Project *pr, const char *fname)
// (Skip nodes with default quality of 0) // (Skip nodes with default quality of 0)
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_QUALITY); fprintf(f, s_QUALITY);
fprintf(f, "\n;%-31s\t%-14s", "ID", "InitQual");
for (i = 1; i <= net->Nnodes; i++) for (i = 1; i <= net->Nnodes; i++)
{ {
node = &net->Node[i]; node = &net->Node[i];
if (node->C0 == 0.0) continue; if (node->C0 == 0.0) continue;
fprintf(f, "\n %-31s %14.6f", node->ID, node->C0 * pr->Ucf[QUALITY]); fprintf(f, "\n %-31s\t%-14.6f", node->ID, node->C0 * pr->Ucf[QUALITY]);
} }
// Write [SOURCES] section // Write [SOURCES] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_SOURCES); fprintf(f, s_SOURCES);
fprintf(f, "\n;%-31s\t%-9s\t%-14s\t%-31s", "ID", "Type", "Quality", "Pattern");
for (i = 1; i <= net->Nnodes; i++) for (i = 1; i <= net->Nnodes; i++)
{ {
node = &net->Node[i]; node = &net->Node[i];
source = node->S; source = node->S;
if (source == NULL) continue; if (source == NULL) continue;
sprintf(s, " %-31s %-8s %14.6f", node->ID, SourceTxt[source->Type], sprintf(s, " %-31s\t%-9s\t%-14.6f", node->ID, SourceTxt[source->Type],
source->C0); source->C0);
if ((j = source->Pat) > 0) if ((j = source->Pat) > 0)
{ {
sprintf(s1, "%s", net->Pattern[j].ID); sprintf(s1, "%s", net->Pattern[j].ID);
} }
else strcpy(s1, ""); else strcpy(s1, "");
fprintf(f, "\n%s %s", s, s1); fprintf(f, "\n%s\t%s", s, s1);
} }
// Write [MIXING] section // Write [MIXING] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_MIXING); fprintf(f, s_MIXING);
fprintf(f, "\n;%-31s\t%-8s", "ID", "Model");
for (i = 1; i <= net->Ntanks; i++) for (i = 1; i <= net->Ntanks; i++)
{ {
tank = &net->Tank[i]; tank = &net->Tank[i];
if (tank->A == 0.0) continue; if (tank->A == 0.0) continue;
fprintf(f, "\n %-31s %-8s %12.4f", net->Node[tank->Node].ID, fprintf(f, "\n %-31s\t%-8s\t%12.4f", net->Node[tank->Node].ID,
MixTxt[tank->MixModel], tank->V1frac); MixTxt[tank->MixModel], tank->V1frac);
} }
@@ -558,6 +582,10 @@ int saveinpfile(Project *pr, const char *fname)
fprintf(f, "\n ROUGHNESS CORRELATION %-.6f", qual->Rfactor); fprintf(f, "\n ROUGHNESS CORRELATION %-.6f", qual->Rfactor);
} }
fprintf(f, "\n\n");
fprintf(f, s_REACTIONS);
fprintf(f, "\n;%-9s\t%-31s\t%-12s", "Type", "Pipe/Tank", "Coefficient");
// Pipe-specific parameters // Pipe-specific parameters
for (i = 1; i <= net->Nlinks; i++) for (i = 1; i <= net->Nlinks; i++)
{ {
@@ -565,11 +593,11 @@ int saveinpfile(Project *pr, const char *fname)
if (link->Type > PIPE) continue; if (link->Type > PIPE) continue;
if (link->Kb != qual->Kbulk) if (link->Kb != qual->Kbulk)
{ {
fprintf(f, "\n BULK %-31s %-.6f", link->ID, link->Kb * SECperDAY); fprintf(f, "\n %-9s\t%-31s\t%-.6f", "BULK", link->ID, link->Kb * SECperDAY);
} }
if (link->Kw != qual->Kwall) if (link->Kw != qual->Kwall)
{ {
fprintf(f, "\n WALL %-31s %-.6f", link->ID, link->Kw * SECperDAY); fprintf(f, "\n %-9s\t%-31s\t%-.6f", "WALL", link->ID, link->Kw * SECperDAY);
} }
} }
@@ -580,7 +608,7 @@ int saveinpfile(Project *pr, const char *fname)
if (tank->A == 0.0) continue; if (tank->A == 0.0) continue;
if (tank->Kb != qual->Kbulk) if (tank->Kb != qual->Kbulk)
{ {
fprintf(f, "\n TANK %-31s %-.6f", net->Node[tank->Node].ID, fprintf(f, "\n %-9s\t%-31s\t%-.6f", "TANK", net->Node[tank->Node].ID,
tank->Kb * SECperDAY); tank->Kb * SECperDAY);
} }
} }
@@ -793,23 +821,25 @@ int saveinpfile(Project *pr, const char *fname)
// Write [COORDINATES] section // Write [COORDINATES] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_COORDS); fprintf(f, s_COORDS);
fprintf(f, "\n;%-31s\t%-14s\t%-14s", "ID", "X-Coord", "Y-Coord");
for (i = 1; i <= net->Nnodes; i++) for (i = 1; i <= net->Nnodes; i++)
{ {
node = &net->Node[i]; node = &net->Node[i];
if (node->X == MISSING || node->Y == MISSING) continue; if (node->X == MISSING || node->Y == MISSING) continue;
fprintf(f, "\n %-31s %14.6f %14.6f", node->ID, node->X, node->Y); fprintf(f, "\n %-31s\t%-14.6f\t%-14.6f", node->ID, node->X, node->Y);
} }
// Write [VERTICES] section // Write [VERTICES] section
fprintf(f, "\n\n"); fprintf(f, "\n\n");
fprintf(f, s_VERTICES); fprintf(f, s_VERTICES);
fprintf(f, "\n;%-31s\t%-14s\t%-14s", "ID", "X-Coord", "Y-Coord");
for (i = 1; i <= net->Nlinks; i++) for (i = 1; i <= net->Nlinks; i++)
{ {
link = &net->Link[i]; link = &net->Link[i];
if (link->Vertices != NULL) if (link->Vertices != NULL)
{ {
for (j = 0; j < link->Vertices->Npts; j++) for (j = 0; j < link->Vertices->Npts; j++)
fprintf(f, "\n %-31s %14.6f %14.6f", fprintf(f, "\n %-31s\t%-14.6f\t%-14.6f",
link->ID, link->Vertices->X[j], link->Vertices->Y[j]); link->ID, link->Vertices->X[j], link->Vertices->Y[j]);
} }
} }