float sensorAccumGetAverage (WV_ACCUM_ID id)
{
float samples, sum;
sum = sensorAccumGetTotal(id);
samples = (float)radListGetNumberOfNodes(&id->samples);
if (samples > 0)
{
return (sum/samples);
}
else
{
return 0;
}
}
static void storeLoopPkt (WVIEWD_WORK *work, LOOP_PKT *dest, WMR918_DATA *src)
{
float tempfloat;
WMR918_IF_DATA* ifWorkData = (WMR918_IF_DATA*)work->stationData;
time_t nowTime = time(NULL);
// Clear optional data:
stationClearLoopData(work);
// WMR918 produces station pressure
if (10 < src->pressure && src->pressure < 50)
{
dest->stationPressure = src->pressure;
// Apply calibration here so the computed values reflect it:
dest->stationPressure *= work->calMPressure;
dest->stationPressure += work->calCPressure;
if (-150 < src->outTemp && src->outTemp < 150)
{
// compute sea-level pressure (BP)
dest->barometer = wvutilsConvertSPToSLP (dest->stationPressure,
src->outTemp,
(float)ifWorkData->elevation);
}
// calculate altimeter
dest->altimeter = wvutilsConvertSPToAltimeter (dest->stationPressure,
(float)ifWorkData->elevation);
}
if (-150 < src->outTemp && src->outTemp < 150)
dest->outTemp = src->outTemp;
if (0 <= src->outHumidity && src->outHumidity <= 100)
{
tempfloat = src->outHumidity;
tempfloat += 0.5;
dest->outHumidity = (USHORT)tempfloat;
}
if (0 <= src->windSpeed && src->windSpeed <= 250)
{
tempfloat = src->windSpeed;
tempfloat += 0.5;
dest->windSpeed = (USHORT)tempfloat;
}
if (0 <= src->windDir && src->windDir <= 360)
{
tempfloat = src->windDir;
tempfloat += 0.5;
dest->windDir = (USHORT)tempfloat;
}
if (0 <= src->maxWindSpeed && src->maxWindSpeed <= 250)
{
tempfloat = src->maxWindSpeed;
tempfloat += 0.5;
dest->windGust = (USHORT)tempfloat;
}
if (0 <= src->maxWindDir && src->maxWindDir <= 360)
{
tempfloat = src->maxWindDir;
tempfloat += 0.5;
dest->windGustDir = (USHORT)tempfloat;
}
if (0 <= src->rain)
{
if (!work->runningFlag)
{
// just starting, so start with whatever the station reports:
ifWorkData->totalRain = src->rain;
dest->sampleRain = 0;
}
else
{
// process the rain accumulator
if (src->rain - ifWorkData->totalRain >= 0)
{
dest->sampleRain = src->rain - ifWorkData->totalRain;
ifWorkData->totalRain = src->rain;
}
else
{
// we had a counter reset...
dest->sampleRain = src->rain;
ifWorkData->totalRain = src->rain;
}
}
if (dest->sampleRain > 2)
{
// Not possible, filter it out:
dest->sampleRain = 0;
}
// Compute rain rate - the WMR918 rain rate is poor...
// Update the rain accumulator:
sensorAccumAddSample (ifWorkData->rainRateAccumulator, nowTime, dest->sampleRain);
dest->rainRate = sensorAccumGetTotal (ifWorkData->rainRateAccumulator);
dest->rainRate *= (60/WMR918_RAIN_RATE_PERIOD);
}
else
{
dest->sampleRain = 0;
sensorAccumAddSample (ifWorkData->rainRateAccumulator, nowTime, dest->sampleRain);
dest->rainRate = sensorAccumGetTotal (ifWorkData->rainRateAccumulator);
dest->rainRate *= (60/WMR918_RAIN_RATE_PERIOD);
}
dest->inTemp = src->inTemp;
tempfloat = src->inHumidity;
tempfloat += 0.5;
dest->inHumidity = (USHORT)tempfloat;
dest->extraTemp1 = (float)src->extraTemp[0];
tempfloat = src->extraHumidity[0];
tempfloat += 0.5;
dest->extraHumid1 = (USHORT)tempfloat;
dest->extraTemp2 = (float)src->extraTemp[1];
tempfloat = src->extraHumidity[1];
tempfloat += 0.5;
dest->extraHumid2 = (USHORT)tempfloat;
dest->extraTemp3 = (float)src->extraTemp[2];
// WMR918 specific:
tempfloat = src->extraHumidity[2];
tempfloat += 0.5;
dest->wmr918Humid3 = (USHORT)tempfloat;
dest->wmr918Pool = src->pool;
dest->wmr918WindBatteryStatus = src->windBatteryStatus;
dest->wmr918RainBatteryStatus = src->rainBatteryStatus;
dest->wmr918OutTempBatteryStatus = src->outTempBatteryStatus;
dest->wmr918InTempBatteryStatus = src->inTempBatteryStatus;
dest->wmr918poolTempBatteryStatus = src->poolTempBatteryStatus;
dest->wmr918extra1BatteryStatus = src->extraBatteryStatus[0];
dest->wmr918extra2BatteryStatus = src->extraBatteryStatus[1];
dest->wmr918extra3BatteryStatus = src->extraBatteryStatus[2];
dest->wmr918Tendency = src->tendency;
return;
}
static void storeLoopPkt (WVIEWD_WORK* work, LOOP_PKT *dest, TWI_DATA *src)
{
float tempfloat;
static time_t lastTempUpdate;
time_t nowTime = time(NULL);
// Clear optional data:
stationClearLoopData(work);
// Pressure:
if ((10 < src->bp && src->bp < 50) &&
(-150 < src->outTemp && src->outTemp < 200))
{
if ((time(NULL) - lastTempUpdate) >= (twiWorkData.archiveInterval * 60))
{
// Add for the 12-hour temp average sample:
// mimic an archive interval:
sensorAccumAddSample(twi12HourTempAvg, time(NULL), src->outTemp);
lastTempUpdate = time(NULL);
}
// TWI produces sea level pressure (SLP):
dest->barometer = src->bp;
// Apply calibration here so the computed values reflect it:
dest->barometer *= work->calMBarometer;
dest->barometer += work->calCBarometer;
// calculate station pressure:
dest->stationPressure = twiConvertSLPToSP(dest->barometer,
sensorAccumGetAverage(twi12HourTempAvg),
(float)twiWorkData.elevation);
// now calculate altimeter:
dest->altimeter = wvutilsConvertSPToAltimeter(dest->stationPressure,
(float)twiWorkData.elevation);
}
if (-150 < src->inTemp && src->inTemp < 200)
{
dest->inTemp = src->inTemp;
}
if (-150 < src->outTemp && src->outTemp < 200)
{
dest->outTemp = src->outTemp;
}
if (0 <= src->humidity && src->humidity <= 100)
{
dest->outHumidity = (uint16_t)src->humidity;
dest->inHumidity = (uint16_t)src->humidity;
}
if (0 <= src->windSpeed && src->windSpeed <= 250)
{
tempfloat = src->windSpeed;
tempfloat += 0.5;
dest->windSpeed = (uint16_t)tempfloat;
dest->windGust = (uint16_t)tempfloat;
}
if (0 <= src->windDir && src->windDir < 360)
{
tempfloat = src->windDir;
tempfloat += 0.5;
dest->windDir = (uint16_t)tempfloat;
dest->windGustDir = (uint16_t)tempfloat;
}
dest->rainRate = src->rainrate;
// process the rain accumulator:
if (0 <= src->dayrain)
{
if (! work->runningFlag)
{
// just starting, so start with whatever the station reports:
twiWorkData.totalRain = src->dayrain;
dest->sampleRain = 0;
}
else
{
// process the rain accumulator
if (src->dayrain - twiWorkData.totalRain >= 0)
{
dest->sampleRain = src->dayrain - twiWorkData.totalRain;
twiWorkData.totalRain = src->dayrain;
}
else
{
// we had a counter reset...
dest->sampleRain = src->dayrain;
twiWorkData.totalRain = src->dayrain;
}
}
if (dest->sampleRain > 2)
{
// Not possible, filter it out:
dest->sampleRain = 0;
}
// Update the rain accumulator:
sensorAccumAddSample (twiWorkData.rainRateAccumulator, nowTime, dest->sampleRain);
dest->rainRate = sensorAccumGetTotal (twiWorkData.rainRateAccumulator);
dest->rainRate *= (60/TWI_RAIN_RATE_PERIOD);
}
else
{
dest->sampleRain = 0;
sensorAccumAddSample (twiWorkData.rainRateAccumulator, nowTime, dest->sampleRain);
dest->rainRate = sensorAccumGetTotal (twiWorkData.rainRateAccumulator);
dest->rainRate *= (60/TWI_RAIN_RATE_PERIOD);
}
return;
}
static void storeLoopPkt (WVIEWD_WORK *work, LOOP_PKT *dest, WH1080_DATA *src)
{
float tempfloat;
WH1080_IF_DATA* ifWorkData = (WH1080_IF_DATA*)work->stationData;
time_t nowTime = time(NULL);
// Clear optional data:
stationClearLoopData(work);
if ((10 < src->pressure && src->pressure < 50) &&
(-150 < src->outtemp && src->outtemp < 150))
{
// WH1080 produces station pressure
dest->stationPressure = src->pressure;
// Apply calibration here so the computed values reflect it:
dest->stationPressure *= work->calMPressure;
dest->stationPressure += work->calCPressure;
// compute sea-level pressure (BP)
tempfloat = wvutilsConvertSPToSLP(dest->stationPressure,
src->outtemp,
(float)ifWorkData->elevation);
dest->barometer = tempfloat;
// calculate altimeter
tempfloat = wvutilsConvertSPToAltimeter(dest->stationPressure,
(float)ifWorkData->elevation);
dest->altimeter = tempfloat;
}
if (-150 < src->outtemp && src->outtemp < 150)
{
dest->outTemp = src->outtemp;
}
if (0 <= src->outhumidity && src->outhumidity <= 100)
{
tempfloat = src->outhumidity;
tempfloat += 0.5;
dest->outHumidity = (uint16_t)tempfloat;
}
if (0 <= src->windAvgSpeed && src->windAvgSpeed <= 250)
{
tempfloat = src->windAvgSpeed;
dest->windSpeedF = tempfloat;
}
if (0 <= src->windDir && src->windDir <= 360)
{
tempfloat = src->windDir;
tempfloat += 0.5;
dest->windDir = (uint16_t)tempfloat;
dest->windGustDir = (uint16_t)tempfloat;
}
if (0 <= src->windGustSpeed && src->windGustSpeed <= 250)
{
tempfloat = src->windGustSpeed;
dest->windGustF = tempfloat;
}
if (0 <= src->rain)
{
if (!work->runningFlag)
{
// just starting, so start with whatever the station reports:
ifWorkData->totalRain = src->rain;
dest->sampleRain = 0;
}
else
{
// process the rain accumulator
if (src->rain - ifWorkData->totalRain >= 0)
{
dest->sampleRain = src->rain - ifWorkData->totalRain;
ifWorkData->totalRain = src->rain;
}
else
{
// we had a counter reset...
dest->sampleRain = src->rain;
dest->sampleRain += (WH1080_RAIN_MAX - ifWorkData->totalRain);
ifWorkData->totalRain = src->rain;
}
}
if (dest->sampleRain > 2)
{
// Not possible, filter it out:
dest->sampleRain = 0;
}
// Update the rain accumulator:
sensorAccumAddSample (ifWorkData->rainRateAccumulator, nowTime, dest->sampleRain);
dest->rainRate = sensorAccumGetTotal (ifWorkData->rainRateAccumulator);
dest->rainRate *= (60/WH1080_RAIN_RATE_PERIOD);
}
else
{
dest->sampleRain = 0;
sensorAccumAddSample (ifWorkData->rainRateAccumulator, nowTime, dest->sampleRain);
dest->rainRate = sensorAccumGetTotal (ifWorkData->rainRateAccumulator);
dest->rainRate *= (60/WH1080_RAIN_RATE_PERIOD);
}
dest->inTemp = src->intemp;
tempfloat = src->inhumidity;
tempfloat += 0.5;
dest->inHumidity = (uint16_t)tempfloat;
return;
}
