epyt-flow 0.14.2__py3-none-any.whl → 0.15.0b1__py3-none-any.whl

epyt_flow/EPANET/EPANET/SRC_engines/inpfile.c

@@ -1,810 +0,0 @@
-/*
-******************************************************************************
-Project:      OWA EPANET
-Version:      2.2
-Module:       inpfile.c
-Description:  saves network data to an EPANET formatted text file
-Authors:      see AUTHORS
-Copyright:    see AUTHORS
-License:      see LICENSE
-Last Updated: 10/29/2019
-******************************************************************************
-*/
-
-#include <stdlib.h>
-#include <stdio.h>
-#include <string.h>
-#include <math.h>
-
-#include "types.h"
-#include "funcs.h"
-#include "hash.h"
-#include "text.h"
-
-// Defined in enumstxt.h in EPANET.C
-extern char *LinkTxt[];
-extern char *FormTxt[];
-extern char *StatTxt[];
-extern char *FlowUnitsTxt[];
-extern char *PressUnitsTxt[];
-extern char *ControlTxt[];
-extern char *SourceTxt[];
-extern char *MixTxt[];
-extern char *TstatTxt[];
-extern char *RptFlagTxt[];
-extern char *SectTxt[];
-
-void saveauxdata(Project *pr, FILE *f)
-/*
-------------------------------------------------------------
-  Writes auxilary data from original input file to new file.
-------------------------------------------------------------
-*/
-{
-    int sect, newsect;
-    char *tok;
-    char write;
-    char line[MAXLINE + 1];
-    char s[MAXLINE + 1];
-    FILE *InFile = pr->parser.InFile;
-
-    // Re-open the input file
-    if (InFile == NULL)
-    {
-        InFile = fopen(pr->parser.InpFname, "rt");
-        if (InFile == NULL) return;
-    }
-    rewind(InFile);
-    sect = -1;
-
-    // Read each line of the input file
-    while (fgets(line, MAXLINE, InFile) != NULL)
-    {
-        strcpy(s, line);
-        tok = strtok(s, SEPSTR);
-        if (tok == NULL) continue;
-
-        // Check if line begins with a new section heading
-        if (*tok == '[')
-        {
-            newsect = findmatch(tok, SectTxt);
-            if (newsect >= 0)
-            {
-                sect = newsect;
-                if (sect == _END) break;
-
-                // Write section heading to file
-                switch (sect)
-                {
-                case _LABELS:
-                case _BACKDROP:
-                case _TAGS:
-                    fprintf(f, "\n%s", line);
-                }
-            }
-        }
-
-        // Write line of auxilary data to file
-        else
-        {
-            write = FALSE;
-            switch (sect)
-            {
-            case _TAGS:
-                if (*tok == ';' ||
-                    (match("NODE", tok) && findnode(&pr->network, strtok(NULL, SEPSTR))) ||
-                    (match("LINK", tok) && findlink(&pr->network, strtok(NULL, SEPSTR))))
-                    write = TRUE;
-                break;
-            case _LABELS:
-            case _BACKDROP:
-                write = TRUE; break;
-            default:
-                break;
-            }
-            if (write) fprintf(f, "%s", line);
-        }
-    }
-    fclose(InFile);
-    InFile = NULL;
-}
-
-int saveinpfile(Project *pr, const char *fname)
-/*
--------------------------------------------------
-  Writes network data to text file.
--------------------------------------------------
-*/
-{
-    Network *net = &pr->network;
-    Parser  *parser = &pr->parser;
-    Report  *rpt = &pr->report;
-    Outfile *out = &pr->outfile;
-    Hydraul *hyd = &pr->hydraul;
-    Quality *qual = &pr->quality;
-    Times   *time = &pr->times;
-
-    int i, j, n;
-    double d, kc, ke, km, ucf;
-    char s[MAXLINE + 1], s1[MAXLINE + 1], s2[MAXLINE + 1];
-    Pdemand demand;
-    Psource source;
-    FILE *f;
-    Slink *link;
-    Stank *tank;
-    Snode *node;
-    Spump *pump;
-    Scontrol *control;
-    Scurve *curve;
-
-    // Open the new text file
-    if ((f = fopen(fname, "wt")) == NULL) return 302;
-
-    // Write [TITLE] section
-    fprintf(f, s_TITLE);
-    for (i = 0; i < 3; i++)
-    {
-        if (strlen(pr->Title[i]) > 0) fprintf(f, "\n%s", pr->Title[i]);
-    }
-
-    // Write [JUNCTIONS] section
-    // (Leave demands for [DEMANDS] section)
-    fprintf(f, "\n\n");
-    fprintf(f, s_JUNCTIONS);
-    for (i = 1; i <= net->Njuncs; i++)
-    {
-        node = &net->Node[i];
-        fprintf(f, "\n %-31s %12.4f", node->ID, node->El * pr->Ucf[ELEV]);
-        if (node->Comment) fprintf(f, "  ;%s", node->Comment);
-    }
-
-    // Write [RESERVOIRS] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_RESERVOIRS);
-    for (i = 1; i <= net->Ntanks; i++)
-    {
-        tank = &net->Tank[i];
-        if (tank->A == 0.0)
-        {
-            node = &net->Node[tank->Node];
-            sprintf(s, " %-31s %12.4f", node->ID, node->El * pr->Ucf[ELEV]);
-            if ((j = tank->Pat) > 0) sprintf(s1, " %s", net->Pattern[j].ID);
-            else strcpy(s1, " ");
-            fprintf(f, "\n%s %-31s", s, s1);
-            if (node->Comment) fprintf(f, " ;%s", node->Comment);
-        }
-    }
-
-    // Write [TANKS] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_TANKS);
-    for (i = 1; i <= net->Ntanks; i++)
-    {
-        tank = &net->Tank[i];
-        if (tank->A > 0.0)
-        {
-            node = &net->Node[tank->Node];
-            sprintf(s, " %-31s %12.4f %12.4f %12.4f %12.4f %12.4f %12.4f",
-                    node->ID, node->El * pr->Ucf[ELEV],
-                    (tank->H0 - node->El) * pr->Ucf[ELEV],
-                    (tank->Hmin - node->El) * pr->Ucf[ELEV],
-                    (tank->Hmax - node->El) * pr->Ucf[ELEV],
-                    sqrt(4.0 * tank->A / PI) * pr->Ucf[ELEV],
-                    tank->Vmin * SQR(pr->Ucf[ELEV]) * pr->Ucf[ELEV]);
-            if ((j = tank->Vcurve) > 0) sprintf(s1, "%s", net->Curve[j].ID);
-            else if (tank->CanOverflow) strcpy(s1, "*");
-            else strcpy(s1, " ");
-            fprintf(f, "\n%s %-31s", s, s1);
-            if (tank->CanOverflow) fprintf(f, "  YES  ");
-            if (node->Comment) fprintf(f, " ;%s", node->Comment);
-        }
-    }
-
-    // Write [PIPES] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_PIPES);
-    for (i = 1; i <= net->Nlinks; i++)
-    {
-        link = &net->Link[i];
-        if (link->Type <= PIPE)
-        {
-            d = link->Diam;
-            kc = link->Kc;
-            if (hyd->Formflag == DW)  kc = kc * pr->Ucf[ELEV] * 1000.0;
-            km = link->Km * SQR(d) * SQR(d) / 0.02517;
-
-            sprintf(s, " %-31s %-31s %-31s %12.4f %12.4f %12.4f %12.4f",
-                    link->ID, net->Node[link->N1].ID, net->Node[link->N2].ID,
-                    link->Len * pr->Ucf[LENGTH], d * pr->Ucf[DIAM], kc, km);
-
-            if (link->Type == CVPIPE) sprintf(s2, "CV");
-            else if (link->Status == CLOSED) sprintf(s2, "CLOSED");
-            else strcpy(s2, " ");
-            fprintf(f, "\n%s %-6s", s, s2);
-            if (link->Comment) fprintf(f, " ;%s", link->Comment);
-        }
-    }
-
-    // Write [PUMPS] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_PUMPS);
-    for (i = 1; i <= net->Npumps; i++)
-    {
-        n = net->Pump[i].Link;
-        link = &net->Link[n];
-        pump = &net->Pump[i];
-        sprintf(s, " %-31s %-31s %-31s", link->ID, net->Node[link->N1].ID,
-                net->Node[link->N2].ID);
-
-        // Pump has constant power
-        if (pump->Ptype == CONST_HP) sprintf(s1, "  POWER %.4f", link->Km);
-
-        // Pump has a head curve
-        else if ((j = pump->Hcurve) > 0)
-        {
-            sprintf(s1, "  HEAD %s", net->Curve[j].ID);
-        }
-
-        // Old format used for pump curve
-        else
-        {
-            fprintf(f, "\n%s %12.4f %12.4f %12.4f          0.0 %12.4f", s,
-                   -pump->H0 * pr->Ucf[HEAD],
-                   (-pump->H0 - pump->R * pow(pump->Q0, pump->N)) * pr->Ucf[HEAD],
-                   pump->Q0 * pr->Ucf[FLOW], pump->Qmax * pr->Ucf[FLOW]);
-            continue;
-        }
-        strcat(s, s1);
-
-        // Optional speed pattern
-        if ((j = pump->Upat) > 0)
-        {
-            sprintf(s1, "  PATTERN  %s", net->Pattern[j].ID);
-            strcat(s, s1);
-        }
-
-        // Optional speed setting
-        if (link->Kc != 1.0)
-        {
-            sprintf(s1, "  SPEED %.4f", link->Kc);
-            strcat(s, s1);
-        }
-
-        fprintf(f, "\n%s", s);
-        if (link->Comment) fprintf(f, "  ;%s", link->Comment);
-
-    }
-
-    // Write [VALVES] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_VALVES);
-    for (i = 1; i <= net->Nvalves; i++)
-    {
-        n = net->Valve[i].Link;
-        link = &net->Link[n];
-        d = link->Diam;
-
-        // Valve setting
-        kc = link->Kc;
-        if (kc == MISSING) kc = 0.0;
-        switch (link->Type)
-        {
-          case FCV:
-            kc *= pr->Ucf[FLOW];
-            break;
-          case PRV:
-          case PSV:
-          case PBV:
-            kc *= pr->Ucf[PRESSURE];
-            break;
-          default:
-            break;
-        }
-        km = link->Km * SQR(d) * SQR(d) / 0.02517;
-
-        sprintf(s, " %-31s %-31s %-31s %12.4f %5s",
-                link->ID, net->Node[link->N1].ID,
-                net->Node[link->N2].ID, d * pr->Ucf[DIAM],
-                LinkTxt[link->Type]);
-
-        // For GPV, setting = head curve index
-        if (link->Type == GPV && (j = ROUND(link->Kc)) > 0)
-        {
-            sprintf(s1, "%-31s %12.4f", net->Curve[j].ID, km);
-        }
-        else sprintf(s1, "%12.4f %12.4f", kc, km);
-        fprintf(f, "\n%s %s", s, s1);
-        if (link->Comment) fprintf(f, " ;%s", link->Comment);
-    }
-
-
-    // Write [DEMANDS] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_DEMANDS);
-    ucf = pr->Ucf[DEMAND];
-    for (i = 1; i <= net->Njuncs; i++)
-    {
-        node = &net->Node[i];
-        for (demand = node->D; demand != NULL; demand = demand->next)
-        {
-            sprintf(s, " %-31s %14.6f", node->ID, ucf * demand->Base);
-            if ((j = demand->Pat) > 0) sprintf(s1, " %-31s", net->Pattern[j].ID);
-            else strcpy(s1, " ");
-            fprintf(f, "\n%s %-31s", s, s1);
-            if (demand->Name) fprintf(f, " ;%s", demand->Name);
-        }
-    }
-
-
-    // Write [EMITTERS] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_EMITTERS);
-    for (i = 1; i <= net->Njuncs; i++)
-    {
-        node = &net->Node[i];
-        if (node->Ke == 0.0) continue;
-        ke = pr->Ucf[FLOW] / pow(pr->Ucf[PRESSURE] * node->Ke, (1.0 / hyd->Qexp));
-        fprintf(f, "\n %-31s %14.6f", node->ID, ke);
-    }
-
-    // Write [STATUS] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_STATUS);
-    for (i = 1; i <= net->Nlinks; i++)
-    {
-        link = &net->Link[i];
-        if (link->Type <= PUMP)
-        {
-            if (link->Status == CLOSED)
-            {
-                fprintf(f, "\n %-31s %s", link->ID, StatTxt[CLOSED]);
-            }
-
-            // Write pump speed here for pumps with old-style pump curve input
-            else if (link->Type == PUMP)
-            {
-                n = findpump(net, i);
-                pump = &net->Pump[n];
-                if (pump->Hcurve == 0 && pump->Ptype != CONST_HP &&
-                    link->Kc != 1.0)
-                {
-                    fprintf(f, "\n %-31s %-.4f", link->ID, link->Kc);
-                }
-            }
-        }
-
-        // Write fixed-status PRVs & PSVs (setting = MISSING)
-        else if (link->Kc == MISSING)
-        {
-            if (link->Status == OPEN)
-            {
-                fprintf(f, "\n %-31s %s", link->ID, StatTxt[OPEN]);
-            }
-            if (link->Status == CLOSED)
-            {
-                fprintf(f, "\n%-31s %s", link->ID, StatTxt[CLOSED]);
-            }
-        }
-    }
-
-    // Write [PATTERNS] section
-    // (Use 6 pattern factors per line)
-    fprintf(f, "\n\n");
-    fprintf(f, s_PATTERNS);
-    for (i = 1; i <= net->Npats; i++)
-    {
-        if (net->Pattern[i].Comment) fprintf(f, "\n;%s", net->Pattern[i].Comment);
-        for (j = 0; j < net->Pattern[i].Length; j++)
-        {
-            if (j % 6 == 0) fprintf(f, "\n %-31s", net->Pattern[i].ID);
-            fprintf(f, " %12.4f", net->Pattern[i].F[j]);
-        }
-    }
-
-    // Write [CURVES] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_CURVES);
-    for (i = 1; i <= net->Ncurves; i++)
-    {
-        if (net->Curve[i].Comment) fprintf(f, "\n;%s", net->Curve[i].Comment);
-        for (j = 0; j < net->Curve[i].Npts; j++)
-        {
-            curve = &net->Curve[i];
-            fprintf(f, "\n %-31s %12.4f %12.4f", curve->ID, curve->X[j], curve->Y[j]);
-        }
-    }
-
-    // Write [CONTROLS] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_CONTROLS);
-    for (i = 1; i <= net->Ncontrols; i++)
-    {
-        // Check that controlled link exists
-        control = &net->Control[i];
-        if ((j = control->Link) <= 0) continue;
-        link = &net->Link[j];
-
-        // Get text of control's link status/setting
-        if (control->Setting == MISSING)
-        {
-            sprintf(s, " LINK %s %s ", link->ID, StatTxt[control->Status]);
-        }
-        else
-        {
-            kc = control->Setting;
-            switch (link->Type)
-            {
-              case PRV:
-              case PSV:
-              case PBV:
-                kc *= pr->Ucf[PRESSURE];
-                break;
-              case FCV:
-                kc *= pr->Ucf[FLOW];
-                break;
-              default:
-                break;
-            }
-            sprintf(s, " LINK %s %.4f", link->ID, kc);
-        }
-
-        switch (control->Type)
-        {
-          // Print level control
-          case LOWLEVEL:
-          case HILEVEL:
-            n = control->Node;
-            node = &net->Node[n];
-            kc = control->Grade - node->El;
-            if (n > net->Njuncs) kc *= pr->Ucf[HEAD];
-            else kc *= pr->Ucf[PRESSURE];
-            fprintf(f, "\n%s IF NODE %s %s %.4f", s, node->ID,
-                    ControlTxt[control->Type], kc);
-            break;
-
-          // Print timer control
-          case TIMER:
-            fprintf(f, "\n%s AT %s %.4f HOURS", s, ControlTxt[TIMER],
-                    control->Time / 3600.);
-            break;
-
-          // Print time-of-day control
-          case TIMEOFDAY:
-            fprintf(f, "\n%s AT %s %s", s, ControlTxt[TIMEOFDAY],
-                    clocktime(rpt->Atime, control->Time));
-            break;
-        }
-    }
-
-    // Write [RULES] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_RULES);
-    for (i = 1; i <= net->Nrules; i++)
-    {
-        fprintf(f, "\nRULE %s", pr->network.Rule[i].label);
-        writerule(pr, f, i);  // see RULES.C
-        fprintf(f, "\n");
-    }
-
-    // Write [QUALITY] section
-    // (Skip nodes with default quality of 0)
-    fprintf(f, "\n\n");
-    fprintf(f, s_QUALITY);
-    for (i = 1; i <= net->Nnodes; i++)
-    {
-        node = &net->Node[i];
-        if (node->C0 == 0.0) continue;
-        fprintf(f, "\n %-31s %14.6f", node->ID, node->C0 * pr->Ucf[QUALITY]);
-    }
-
-    // Write [SOURCES] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_SOURCES);
-    for (i = 1; i <= net->Nnodes; i++)
-    {
-        node = &net->Node[i];
-        source = node->S;
-        if (source == NULL) continue;
-        sprintf(s, " %-31s %-8s %14.6f", node->ID, SourceTxt[source->Type],
-                source->C0);
-        if ((j = source->Pat) > 0)
-        {
-            sprintf(s1, "%s", net->Pattern[j].ID);
-        }
-        else strcpy(s1, "");
-        fprintf(f, "\n%s %s", s, s1);
-    }
-
-    // Write [MIXING] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_MIXING);
-    for (i = 1; i <= net->Ntanks; i++)
-    {
-        tank = &net->Tank[i];
-        if (tank->A == 0.0) continue;
-        fprintf(f, "\n %-31s %-8s %12.4f", net->Node[tank->Node].ID,
-                MixTxt[tank->MixModel], (tank->V1max / tank->Vmax));
-    }
-
-    // Write [REACTIONS] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_REACTIONS);
-
-    // General parameters
-    fprintf(f, "\n ORDER  BULK            %-.2f", qual->BulkOrder);
-    fprintf(f, "\n ORDER  WALL            %-.0f", qual->WallOrder);
-    fprintf(f, "\n ORDER  TANK            %-.2f", qual->TankOrder);
-    fprintf(f, "\n GLOBAL BULK            %-.6f", qual->Kbulk * SECperDAY);
-    fprintf(f, "\n GLOBAL WALL            %-.6f", qual->Kwall * SECperDAY);
-
-    if (qual->Climit > 0.0)
-    {
-        fprintf(f, "\n LIMITING POTENTIAL     %-.6f", qual->Climit * pr->Ucf[QUALITY]);
-    }
-    if (qual->Rfactor != MISSING && qual->Rfactor != 0.0)
-    {
-        fprintf(f, "\n ROUGHNESS CORRELATION  %-.6f", qual->Rfactor);
-    }
-
-    // Pipe-specific parameters
-    for (i = 1; i <= net->Nlinks; i++)
-    {
-        link = &net->Link[i];
-        if (link->Type > PIPE) continue;
-        if (link->Kb != qual->Kbulk)
-        {
-            fprintf(f, "\n BULK   %-31s %-.6f", link->ID, link->Kb * SECperDAY);
-        }
-        if (link->Kw != qual->Kwall)
-        {
-            fprintf(f, "\n WALL   %-31s %-.6f", link->ID, link->Kw * SECperDAY);
-        }
-    }
-
-    // Tank parameters
-    for (i = 1; i <= net->Ntanks; i++)
-    {
-        tank = &net->Tank[i];
-        if (tank->A == 0.0) continue;
-        if (tank->Kb != qual->Kbulk)
-        {
-            fprintf(f, "\n TANK   %-31s %-.6f", net->Node[tank->Node].ID,
-                    tank->Kb * SECperDAY);
-        }
-    }
-
-    // Write [ENERGY] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_ENERGY);
-
-    // General parameters
-    if (hyd->Ecost != 0.0)
-    {
-        fprintf(f, "\n GLOBAL PRICE        %-.4f", hyd->Ecost);
-    }
-    if (hyd->Epat != 0)
-    {
-        fprintf(f, "\n GLOBAL PATTERN      %s", net->Pattern[hyd->Epat].ID);
-    }
-    fprintf(f, "\n GLOBAL EFFIC        %-.4f", hyd->Epump);
-    fprintf(f, "\n DEMAND CHARGE       %-.4f", hyd->Dcost);
-
-    // Pump-specific parameters
-    for (i = 1; i <= net->Npumps; i++)
-    {
-        pump = &net->Pump[i];
-        if (pump->Ecost > 0.0)
-        {
-            fprintf(f, "\n PUMP %-31s PRICE   %-.4f", net->Link[pump->Link].ID,
-                    pump->Ecost);
-        }
-        if (pump->Epat > 0.0)
-        {
-            fprintf(f, "\n PUMP %-31s PATTERN %s", net->Link[pump->Link].ID,
-                    net->Pattern[pump->Epat].ID);
-        }
-        if (pump->Ecurve > 0.0)
-        {
-            fprintf(f, "\n PUMP %-31s EFFIC   %s", net->Link[pump->Link].ID,
-                    net->Curve[pump->Ecurve].ID);
-        }
-    }
-
-    // Write [TIMES] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_TIMES);
-    fprintf(f, "\n DURATION            %s", clocktime(rpt->Atime, time->Dur));
-    fprintf(f, "\n HYDRAULIC TIMESTEP  %s", clocktime(rpt->Atime, time->Hstep));
-    fprintf(f, "\n QUALITY TIMESTEP    %s", clocktime(rpt->Atime, time->Qstep));
-    fprintf(f, "\n REPORT TIMESTEP     %s", clocktime(rpt->Atime, time->Rstep));
-    fprintf(f, "\n REPORT START        %s", clocktime(rpt->Atime, time->Rstart));
-    fprintf(f, "\n PATTERN TIMESTEP    %s", clocktime(rpt->Atime, time->Pstep));
-    fprintf(f, "\n PATTERN START       %s", clocktime(rpt->Atime, time->Pstart));
-    fprintf(f, "\n RULE TIMESTEP       %s", clocktime(rpt->Atime, time->Rulestep));
-    fprintf(f, "\n START CLOCKTIME     %s", clocktime(rpt->Atime, time->Tstart));
-    fprintf(f, "\n STATISTIC           %s", TstatTxt[rpt->Tstatflag]);
-
-    // Write [OPTIONS] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_OPTIONS);
-    fprintf(f, "\n UNITS               %s", FlowUnitsTxt[parser->Flowflag]);
-    fprintf(f, "\n PRESSURE            %s", PressUnitsTxt[parser->Pressflag]);
-    fprintf(f, "\n HEADLOSS            %s", FormTxt[hyd->Formflag]);
-    switch (out->Hydflag)
-    {
-        case USE:
-          fprintf(f, "\n HYDRAULICS USE      %s", out->HydFname);
-          break;
-        case SAVE:
-          fprintf(f, "\n HYDRAULICS SAVE     %s", out->HydFname);
-          break;
-    }
-    if (hyd->ExtraIter == -1)
-    {
-        fprintf(f, "\n UNBALANCED          STOP");
-    }
-    if (hyd->ExtraIter >= 0)
-    {
-        fprintf(f, "\n UNBALANCED          CONTINUE %d", hyd->ExtraIter);
-    }
-
-    switch (qual->Qualflag)
-    {
-        case CHEM:
-          fprintf(f, "\n QUALITY             %s %s",
-                  qual->ChemName, qual->ChemUnits);
-          break;
-        case TRACE:
-          fprintf(f, "\n QUALITY             TRACE %-31s",
-                  net->Node[qual->TraceNode].ID);
-          break;
-        case AGE:
-          fprintf(f, "\n QUALITY             AGE");
-          break;
-        case NONE:
-          fprintf(f, "\n QUALITY             NONE");
-          break;
-    }
-
-    fprintf(f, "\n DEMAND MULTIPLIER   %-.4f", hyd->Dmult);
-    fprintf(f, "\n EMITTER EXPONENT    %-.4f", 1.0 / hyd->Qexp);
-    fprintf(f, "\n VISCOSITY           %-.6f", hyd->Viscos / VISCOS);
-    fprintf(f, "\n DIFFUSIVITY         %-.6f", qual->Diffus / DIFFUS);
-    fprintf(f, "\n SPECIFIC GRAVITY    %-.6f", hyd->SpGrav);
-    fprintf(f, "\n TRIALS              %-d", hyd->MaxIter);
-    fprintf(f, "\n ACCURACY            %-.8f", hyd->Hacc);
-    fprintf(f, "\n TOLERANCE           %-.8f", qual->Ctol * pr->Ucf[QUALITY]);
-    fprintf(f, "\n CHECKFREQ           %-d", hyd->CheckFreq);
-    fprintf(f, "\n MAXCHECK            %-d", hyd->MaxCheck);
-    fprintf(f, "\n DAMPLIMIT           %-.8f", hyd->DampLimit);
-    if (hyd->HeadErrorLimit > 0.0)
-    {
-        fprintf(f, "\n HEADERROR           %-.8f",
-                hyd->HeadErrorLimit * pr->Ucf[HEAD]);
-    }
-    if (hyd->FlowChangeLimit > 0.0)
-    {
-        fprintf(f, "\n FLOWCHANGE          %-.8f",
-                hyd->FlowChangeLimit * pr->Ucf[FLOW]);
-    }
-    if (hyd->DemandModel == PDA)
-    {
-        fprintf(f, "\n DEMAND MODEL        PDA");
-        fprintf(f, "\n MINIMUM PRESSURE    %-.4f", hyd->Pmin * pr->Ucf[PRESSURE]);
-        fprintf(f, "\n REQUIRED PRESSURE   %-.4f", hyd->Preq * pr->Ucf[PRESSURE]);
-        fprintf(f, "\n PRESSURE EXPONENT   %-.4f", hyd->Pexp);
-    }
-
-    // Write [REPORT] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_REPORT);
-
-    // General options
-    fprintf(f, "\n PAGESIZE            %d", rpt->PageSize);
-    fprintf(f, "\n STATUS              %s", RptFlagTxt[rpt->Statflag]);
-    fprintf(f, "\n SUMMARY             %s", RptFlagTxt[rpt->Summaryflag]);
-    fprintf(f, "\n ENERGY              %s", RptFlagTxt[rpt->Energyflag]);
-    fprintf(f, "\n MESSAGES            %s", RptFlagTxt[rpt->Messageflag]);
-    if (strlen(rpt->Rpt2Fname) > 0)
-    {
-        fprintf(f, "\n FILE                %s", rpt->Rpt2Fname);
-    }
-
-    // Node reporting
-    switch (rpt->Nodeflag)
-    {
-        case 0:
-          fprintf(f, "\n NODES               NONE");
-          break;
-        case 1:
-          fprintf(f, "\n NODES               ALL");
-          break;
-        default:
-          j = 0;
-          for (i = 1; i <= net->Nnodes; i++)
-          {
-              node = &net->Node[i];
-              if (node->Rpt == 1)
-              {
-                  if (j % 5 == 0) fprintf(f, "\n NODES               ");
-                  fprintf(f, "%s ", node->ID);
-                  j++;
-             }
-          }
-    }
-
-    // Link reporting
-    switch (rpt->Linkflag)
-    {
-        case 0:
-          fprintf(f, "\n LINKS               NONE");
-          break;
-        case 1:
-          fprintf(f, "\n LINKS               ALL");
-          break;
-        default:
-          j = 0;
-          for (i = 1; i <= net->Nlinks; i++)
-          {
-              link = &net->Link[i];
-              if (link->Rpt == 1)
-              {
-                  if (j % 5 == 0) fprintf(f, "\n LINKS               ");
-                  fprintf(f, "%s ", link->ID);
-                  j++;
-              }
-          }
-    }
-
-    // Field formatting options
-    for (i = 0; i < FRICTION; i++)
-    {
-        SField *field = &rpt->Field[i];
-        if (field->Enabled == TRUE)
-        {
-            fprintf(f, "\n %-20sPRECISION %d", field->Name, field->Precision);
-            if (field->RptLim[LOW] < BIG)
-            {
-                fprintf(f, "\n %-20sBELOW %.6f", field->Name, field->RptLim[LOW]);
-            }
-            if (field->RptLim[HI] > -BIG)
-            {
-                fprintf(f, "\n %-20sABOVE %.6f", field->Name, field->RptLim[HI]);
-            }
-        }
-        else fprintf(f, "\n %-20sNO",field->Name);
-    }
-
-    // Write [COORDINATES] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_COORDS);
-    for (i = 1; i <= net->Nnodes; i++)
-    {
-        node = &net->Node[i];
-        if (node->X == MISSING || node->Y == MISSING) continue;
-        fprintf(f, "\n %-31s %14.6f %14.6f", node->ID, node->X, node->Y);
-    }
-
-    // Write [VERTICES] section
-    fprintf(f, "\n\n");
-    fprintf(f, s_VERTICES);
-    for (i = 1; i <= net->Nlinks; i++)
-    {
-        link = &net->Link[i];
-        if (link->Vertices != NULL)
-        {
-            for (j = 0; j < link->Vertices->Npts; j++)
-                fprintf(f, "\n %-31s %14.6f %14.6f",
-                    link->ID, link->Vertices->X[j], link->Vertices->Y[j]);
-        }
-    }
-
-    // Save auxilary data to new input file
-    fprintf(f, "\n");
-    saveauxdata(pr, f);
-
-    // Close the new input file
-    fprintf(f, "\n%s\n", s_END);
-    fclose(f);
-    return 0;
-}