void gage_initState(int j)
//
// Input: j = rain gage index
// Output: none
// Purpose: initializes state of rain gage.
//
{
// --- assume gage not used by any subcatchment
// (will be updated in subcatch_initState)
Gage[j].isUsed = FALSE;
Gage[j].rainfall = 0.0;
Gage[j].reportRainfall = 0.0;
if ( IgnoreRainfall ) return;
// --- for gage with file data:
if ( Gage[j].dataSource == RAIN_FILE )
{
// --- set current file position to start of period of record
Gage[j].currentFilePos = Gage[j].startFilePos;
// --- assign units conversion factor
// (rain depths on interface file are in inches)
if ( UnitSystem == SI ) Gage[j].unitsFactor = MMperINCH;
}
// --- get first & next rainfall values
if ( getFirstRainfall(j) )
{
// --- find date at end of starting rain interval
Gage[j].endDate = datetime_addSeconds(
Gage[j].startDate, Gage[j].rainInterval);
// --- if rainfall record begins after start of simulation,
if ( Gage[j].startDate > StartDateTime )
{
// --- make next rainfall date the start of the rain record
Gage[j].nextDate = Gage[j].startDate;
Gage[j].nextRainfall = Gage[j].rainfall;
// --- make start of current rain interval the simulation start
Gage[j].startDate = StartDateTime;
Gage[j].endDate = Gage[j].nextDate;
Gage[j].rainfall = 0.0;
}
// --- otherwise find next recorded rainfall
else if ( !getNextRainfall(j) ) Gage[j].nextDate = NO_DATE;
}
else Gage[j].startDate = NO_DATE;
}
void gage_setState(int j, DateTime t)
//
// Input: j = rain gage index
// t = a calendar date/time
// Output: none
// Purpose: updates state of rain gage for specified date.
//
{
// --- return if gage not used by any subcatchment
if ( Gage[j].isUsed == FALSE ) return;
// --- set rainfall to zero if disabled
if ( IgnoreRainfall )
{
Gage[j].rainfall = 0.0;
return;
}
// --- use rainfall from co-gage (gage with lower index that uses
// same rainfall time series or file) if it exists
if ( Gage[j].coGage >= 0)
{
Gage[j].rainfall = Gage[Gage[j].coGage].rainfall;
return;
}
// --- otherwise march through rainfall record until date t is bracketed
t += OneSecond;
for (;;)
{
// --- no rainfall if no interval start date
if ( Gage[j].startDate == NO_DATE )
{
Gage[j].rainfall = 0.0;
return;
}
// --- no rainfall if time is before interval start date
if ( t < Gage[j].startDate )
{
Gage[j].rainfall = 0.0;
return;
}
// --- use current rainfall if time is before interval end date
if ( t < Gage[j].endDate )
{
return;
}
// --- no rainfall if t >= interval end date & no next interval exists
if ( Gage[j].nextDate == NO_DATE )
{
Gage[j].rainfall = 0.0;
return;
}
// --- no rainfall if t > interval end date & < next interval date
if ( t < Gage[j].nextDate )
{
Gage[j].rainfall = 0.0;
return;
}
// --- otherwise update next rainfall interval date
Gage[j].startDate = Gage[j].nextDate;
Gage[j].endDate = datetime_addSeconds(Gage[j].startDate,
Gage[j].rainInterval);
Gage[j].rainfall = Gage[j].nextRainfall;
if ( !getNextRainfall(j) ) Gage[j].nextDate = NO_DATE;
}
}
void gage_initState(int j)
//
// Input: j = rain gage index
// Output: none
// Purpose: initializes state of rain gage.
//
{
int i, k;
// --- assume gage not used by any subcatchment
// (will be updated in subcatch_initState)
Gage[j].isUsed = FALSE;
Gage[j].rainfall = 0.0; //(5.0.010 - LR)
Gage[j].reportRainfall = 0.0; //(5.0.010 - LR)
if ( IgnoreRainfall ) return; //(5.0.010 - LR)
// --- for gage with file data:
if ( Gage[j].dataSource == RAIN_FILE )
{
// --- set current file position to start of period of record
Gage[j].currentFilePos = Gage[j].startFilePos;
// --- assign units conversion factor
// (rain depths on interface file are in inches)
if ( UnitSystem == SI ) Gage[j].unitsFactor = MMperINCH;
}
// --- for gage with time series data, determine if gage uses
// same time series as another gage with lower index
// (no check required for gages sharing same rain file)
if ( Gage[j].dataSource == RAIN_TSERIES )
{
k = Gage[j].tSeries;
for (i=0; i<j; i++)
{
if ( Gage[i].dataSource == RAIN_TSERIES &&
Gage[i].tSeries == k)
{
Gage[j].coGage = i;
return;
}
}
}
// --- get first & next rainfall values
if ( getFirstRainfall(j) )
{
// --- find date at end of starting rain interval
Gage[j].endDate = datetime_addSeconds(
Gage[j].startDate, Gage[j].rainInterval);
// --- if rainfall record begins after start of simulation,
if ( Gage[j].startDate > StartDateTime )
{
// --- make next rainfall date the start of the rain record
Gage[j].nextDate = Gage[j].startDate;
Gage[j].nextRainfall = Gage[j].rainfall;
// --- make start of current rain interval the simulation start
Gage[j].startDate = StartDateTime;
Gage[j].endDate = Gage[j].nextDate;
Gage[j].rainfall = 0.0;
}
// --- otherwise find next recorded rainfall
else if ( !getNextRainfall(j) ) Gage[j].nextDate = NO_DATE;
}
else Gage[j].startDate = NO_DATE;
}
