void gage_validate(int j)
//
// Input: j = rain gage index
// Output: none
// Purpose: checks for valid rain gage parameters
//
// NOTE: assumes that any time series used by a rain gage has been
// previously validated.
//
{
int i, k;
int gageInterval;
// --- for gage with time series data:
if ( Gage[j].dataSource == RAIN_TSERIES )
{
// --- check gage's recording interval against that of time series
k = Gage[j].tSeries;
if ( Tseries[k].refersTo >= 0 )
{
report_writeErrorMsg(ERR_RAIN_GAGE_TSERIES, Gage[j].ID);
}
gageInterval = (int)(floor(Tseries[k].dxMin*SECperDAY + 0.5));
if ( gageInterval > 0 && Gage[j].rainInterval > gageInterval )
{
report_writeErrorMsg(ERR_RAIN_GAGE_INTERVAL, Gage[j].ID);
}
if ( Gage[j].rainInterval < gageInterval )
{
report_writeWarningMsg(WARN09, Gage[j].ID);
}
if ( Gage[j].rainInterval < WetStep )
{
report_writeWarningMsg(WARN01, Gage[j].ID);
WetStep = Gage[j].rainInterval;
}
// --- see if gage uses same time series as another gage
for (i=0; i<j; i++)
{
if ( Gage[i].dataSource == RAIN_TSERIES && Gage[i].tSeries == k )
{
Gage[j].coGage = i;
// --- check that both gages record same type of data
if ( Gage[j].rainType != Gage[i].rainType )
{
report_writeErrorMsg(ERR_RAIN_GAGE_FORMAT, Gage[j].ID);
}
return;
}
}
}
}
void node_validate(int j)
//
// Input: j = node index
// Output: none
// Purpose: validates a node's properties.
//
{
TDwfInflow* inflow;
// --- see if full depth was increased to accommodate conduit crown
if ( Node[j].fullDepth > Node[j].oldDepth && Node[j].oldDepth > 0.0 )
{
report_writeWarningMsg(WARN02, Node[j].ID);
}
// --- check that initial depth does not exceed max. depth
if ( Node[j].initDepth > Node[j].fullDepth + Node[j].surDepth )
report_writeErrorMsg(ERR_NODE_DEPTH, Node[j].ID);
if ( Node[j].type == DIVIDER ) divider_validate(j);
// --- initialize dry weather inflows
inflow = Node[j].dwfInflow;
while (inflow)
{
inflow_initDwfInflow(inflow);
inflow = inflow->next;
}
}
void node_validate(int j)
//
// Input: j = node index
// Output: none
// Purpose: validates a node's properties.
//
{
// --- see if full depth was increased to accommodate conduit crown //(5.0.014 - LR)
if ( Node[j].fullDepth > Node[j].oldDepth && Node[j].oldDepth > 0.0 ) //(5.0.014 - LR)
{ //(5.0.014 - LR)
report_writeWarningMsg(WARN02, Node[j].ID); //(5.0.015 - LR)
} //(5.0.014 - LR)
// --- check that initial depth does not exceed max. depth
if ( Node[j].initDepth > Node[j].fullDepth + Node[j].surDepth )
report_writeErrorMsg(ERR_NODE_DEPTH, Node[j].ID);
if ( Node[j].type == DIVIDER ) divider_validate(j);
}
void project_validate()
//
// Input: none
// Output: none
// Purpose: checks validity of project data.
//
{
int i;
int j;
int err;
// --- validate Curves and TimeSeries
for ( i=0; i<Nobjects[CURVE]; i++ )
{
err = table_validate(&Curve[i]);
if ( err ) report_writeErrorMsg(ERR_CURVE_SEQUENCE, Curve[i].ID);
}
for ( i=0; i<Nobjects[TSERIES]; i++ )
{
err = table_validate(&Tseries[i]);
if ( err ) report_writeTseriesErrorMsg(err, &Tseries[i]);
}
// --- validate hydrology objects
// (NOTE: order is important !!!!)
climate_validate();
lid_validate();
if ( Nobjects[SNOWMELT] == 0 ) IgnoreSnowmelt = TRUE;
if ( Nobjects[AQUIFER] == 0 ) IgnoreGwater = TRUE;
for ( i=0; i<Nobjects[GAGE]; i++ ) gage_validate(i);
for ( i=0; i<Nobjects[AQUIFER]; i++ ) gwater_validateAquifer(i);
for ( i=0; i<Nobjects[SUBCATCH]; i++ ) subcatch_validate(i);
for ( i=0; i<Nobjects[SNOWMELT]; i++ ) snow_validateSnowmelt(i);
// --- compute geometry tables for each shape curve
j = 0;
for ( i=0; i<Nobjects[CURVE]; i++ )
{
if ( Curve[i].curveType == SHAPE_CURVE )
{
Curve[i].refersTo = j;
Shape[j].curve = i;
if ( !shape_validate(&Shape[j], &Curve[i]) )
report_writeErrorMsg(ERR_CURVE_SEQUENCE, Curve[i].ID);
j++;
}
}
// --- validate links before nodes, since the latter can
// result in adjustment of node depths
for ( i=0; i<Nobjects[NODE]; i++) Node[i].oldDepth = Node[i].fullDepth;
for ( i=0; i<Nobjects[LINK]; i++) link_validate(i);
for ( i=0; i<Nobjects[NODE]; i++) node_validate(i);
// --- adjust time steps if necessary
if ( DryStep < WetStep )
{
report_writeWarningMsg(WARN06, "");
DryStep = WetStep;
}
if ( RouteStep > (double)WetStep )
{
report_writeWarningMsg(WARN07, "");
RouteStep = WetStep;
}
// --- adjust individual reporting flags to match global reporting flag
if ( RptFlags.subcatchments == ALL )
for (i=0; i<Nobjects[SUBCATCH]; i++) Subcatch[i].rptFlag = TRUE;
if ( RptFlags.nodes == ALL )
for (i=0; i<Nobjects[NODE]; i++) Node[i].rptFlag = TRUE;
if ( RptFlags.links == ALL )
for (i=0; i<Nobjects[LINK]; i++) Link[i].rptFlag = TRUE;
// --- adjust DYNWAVE options
if ( MinSurfArea == 0.0 ) MinSurfArea = DEFAULT_SURFAREA;
else MinSurfArea /= UCF(LENGTH) * UCF(LENGTH);
if ( HeadTol == 0.0 ) HeadTol = DEFAULT_HEADTOL;
else HeadTol /= UCF(LENGTH);
if ( MaxTrials == 0 ) MaxTrials = DEFAULT_MAXTRIALS;
}
void project_validate()
//
// Input: none
// Output: none
// Purpose: checks validity of project data.
//
{
int i;
int j;
int err;
// --- validate Curves and TimeSeries
for ( i=0; i<Nobjects[CURVE]; i++ )
{
err = table_validate(&Curve[i]);
if ( err ) report_writeErrorMsg(ERR_CURVE_SEQUENCE, Curve[i].ID);
}
for ( i=0; i<Nobjects[TSERIES]; i++ )
{
err = table_validate(&Tseries[i]);
if ( err ) report_writeTseriesErrorMsg(err, &Tseries[i]);
}
// --- validate hydrology objects
// (NOTE: order is important !!!!)
climate_validate();
lid_validate();
if ( Nobjects[SNOWMELT] == 0 ) IgnoreSnowmelt = TRUE;
if ( Nobjects[AQUIFER] == 0 ) IgnoreGwater = TRUE;
for ( i=0; i<Nobjects[GAGE]; i++ ) gage_validate(i);
for ( i=0; i<Nobjects[AQUIFER]; i++ ) gwater_validateAquifer(i);
for ( i=0; i<Nobjects[SUBCATCH]; i++ ) subcatch_validate(i);
for ( i=0; i<Nobjects[SNOWMELT]; i++ ) snow_validateSnowmelt(i);
// --- compute geometry tables for each shape curve
j = 0;
for ( i=0; i<Nobjects[CURVE]; i++ )
{
if ( Curve[i].curveType == SHAPE_CURVE )
{
Curve[i].refersTo = j;
Shape[j].curve = i;
if ( !shape_validate(&Shape[j], &Curve[i]) )
report_writeErrorMsg(ERR_CURVE_SEQUENCE, Curve[i].ID);
j++;
}
}
// --- validate links before nodes, since the latter can
// result in adjustment of node depths
for ( i=0; i<Nobjects[NODE]; i++) Node[i].oldDepth = Node[i].fullDepth;
for ( i=0; i<Nobjects[LINK]; i++) link_validate(i);
for ( i=0; i<Nobjects[NODE]; i++) node_validate(i);
// --- adjust time steps if necessary
if ( DryStep < WetStep )
{
report_writeWarningMsg(WARN06, "");
DryStep = WetStep;
}
if ( RouteStep > (double)WetStep )
{
report_writeWarningMsg(WARN07, "");
RouteStep = WetStep;
}
// --- adjust individual reporting flags to match global reporting flag
if ( RptFlags.subcatchments == ALL )
for (i=0; i<Nobjects[SUBCATCH]; i++) Subcatch[i].rptFlag = TRUE;
if ( RptFlags.nodes == ALL )
for (i=0; i<Nobjects[NODE]; i++) Node[i].rptFlag = TRUE;
if ( RptFlags.links == ALL )
for (i=0; i<Nobjects[LINK]; i++) Link[i].rptFlag = TRUE;
// --- validate dynamic wave options
if ( RouteModel == DW ) dynwave_validate(); //(5.1.008)
#pragma omp parallel //(5.1.008)
{
if ( NumThreads == 0 ) NumThreads = omp_get_num_threads(); //(5.1.008)
else NumThreads = MIN(NumThreads, omp_get_num_threads()); //(5.1.008)
}
if ( Nobjects[LINK] < 4 * NumThreads ) NumThreads = 1; //(5.1.008)
}
