HI,
I’m hoping someone can help me out.
I have a Photon, with a sparkfun weathershield and weather meter kit (wind sp, direction, precipitation), and have it uploading to personal weather station on weather underground (http://www.wunderground.com/personal-we … X3#history). I used the tutorial (sparfun weather/imp/phant) and have made it a ways. The problem I am running into is that the weather meter does not appear to accurately work and upload to wunderground. (Note: I also have not connected the photon to a serial term to check actual readings, since the code puts the photon into a deep sleep after upload).
I’m hoping someone can point me in the right direction.
1.) Wind speed is uploaded to wunderground in spurts of 7.2 no matter what the actual speed
2.) no perceptiona has been recorded wunderground even though I have had alot of rain
3.) Always reports wind for “N” (north) to underground even though other directions have been observed.
Any help would be appreciated: code from build.io below
// This #include statement was automatically added by the Particle IDE.
/******************************************************************************
SparkFun_Photon_Weather_Wunderground.ino
SparkFun Photon Weather Shield basic example
Joel Bartlett @ SparkFun Electronics
Original Creation Date: May 18, 2015
Upadted Spetmeber 14 2015
Modified by Clayton Wrisley
This sketch prints the temperature, humidity, and barometric pressure OR
altitude to the Serial port.
The library used in this example can be found here:
https://github.com/sparkfun/SparkFun_Ph … le_Library
Hardware Connections:
This sketch was written specifically for the Photon Weather Shield,
which connects the HTU21D and MPL3115A2 to the I2C bus by default.
If you have an HTU21D and/or an MPL3115A2 breakout, use the following
hardware setup:
HTU21D ------------- Photon
(-) ------------------- GND
(+) ------------------- 3.3V (VCC)
CL ------------------- D1/SCL
DA ------------------- D0/SDA
MPL3115A2 ------------- Photon
GND ------------------- GND
VCC ------------------- 3.3V (VCC)
SCL ------------------ D1/SCL
SDA ------------------ D0/SDA
Development environment specifics:
IDE: Particle Dev
Hardware Platform: Particle Photon
Particle Core
This code is beerware; if you see me (or any other SparkFun
employee) at the local, and you’ve found our code helpful,
please buy us a round!
Distributed as-is; no warranty is given.
//---------------------------------------------------------------
Weather Underground Upload sections: Dan Fein @ Weather Underground
Weather Underground Upload Protocol:
http://wiki.wunderground.com/index.php/ … d_Protocol
Sign up at http://www.wunderground.com/personal-we … signup.asp
*******************************************************************************/
#include “SparkFun_Photon_Weather_Shield_Library/SparkFun_Photon_Weather_Shield_Library.h”
#include “math.h” //For Dew Point Calculation
float humidity = 0;
float humTempF = 0; //humidity sensor temp reading, fahrenheit
float humTempC = 0; //humidity sensor temp reading, celsius
float baroTempF = 0; //barometer sensor temp reading, fahrenheit
float baroTempC = 0; //barometer sensor temp reading, celsius
float tempF = 0; //Average of the sensors temperature readings, fahrenheit
float tempC = 0; //Average of the sensors temperature readings, celsius
float dewptF = 0;
float dewptC = 0;
float pascals = 0;
float inches = 0;
int WDIR = A0;
int RAIN = D2;
int WSPEED = D3;
//Global Variables
//-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
long lastSecond; //The millis counter to see when a second rolls by
byte seconds; //When it hits 60, increase the current minute
byte seconds_2m; //Keeps track of the “wind speed/dir avg” over last 2 minutes array of data
byte minutes; //Keeps track of where we are in various arrays of data
byte minutes_10m; //Keeps track of where we are in wind gust/dir over last 10 minutes array of data
//We need to keep track of the following variables:
//Wind speed/dir each update (no storage)
//Wind gust/dir over the day (no storage)
//Wind speed/dir, avg over 2 minutes (store 1 per second)
//Wind gust/dir over last 10 minutes (store 1 per minute)
//Rain over the past hour (store 1 per minute)
//Total rain over date (store one per day)
byte windspdavg[120]; //120 bytes to keep track of 2 minute average
int winddiravg[120]; //120 ints to keep track of 2 minute average
float windgust_10m[10]; //10 floats to keep track of 10 minute max
int windgustdirection_10m[10]; //10 ints to keep track of 10 minute max
volatile float rainHour[60]; //60 floating numbers to keep track of 60 minutes of rain
//These are all the weather values that wunderground expects:
int winddir = 0; // [0-360 instantaneous wind direction]
float windspeedmph = 0; // [mph instantaneous wind speed]
float windgustmph = 0; // [mph current wind gust, using software specific time period]
int windgustdir = 0; // [0-360 using software specific time period]
float windspdmph_avg2m = 0; // [mph 2 minute average wind speed mph]
int winddir_avg2m = 0; // [0-360 2 minute average wind direction]
float windgustmph_10m = 0; // [mph past 10 minutes wind gust mph ]
int windgustdir_10m = 0; // [0-360 past 10 minutes wind gust direction]
float rainin = 0; // [rain inches over the past hour)] – the accumulated rainfall in the past 60 min
long lastWindCheck = 0;
volatile float dailyrainin = 0; // [rain inches so far today in local time]
int count = 0;
// volatiles are subject to modification by IRQs
volatile long lastWindIRQ = 0;
volatile byte windClicks = 0;
volatile unsigned long raintime, rainlast, raininterval, rain;
void rainIRQ()
// Count rain gauge bucket tips as they occur
// Activated by the magnet and reed switch in the rain gauge, attached to input D2
{
raintime = millis(); // grab current time
raininterval = raintime - rainlast; // calculate interval between this and last event
if (raininterval > 10) // ignore switch-bounce glitches less than 10mS after initial edge
{
dailyrainin += 0.011; //Each dump is 0.011" of water
rainHour[minutes] += 0.011; //Increase this minute’s amount of rain
rainlast = raintime; // set up for next event
}
}
void wspeedIRQ()
// Activated by the magnet in the anemometer (2 ticks per rotation), attached to input D3
{
if (millis() - lastWindIRQ > 10) // Ignore switch-bounce glitches less than 10ms (142MPH max reading) after the reed switch closes
{
lastWindIRQ = millis(); //Grab the current time
windClicks++; //There is 1.492MPH for each click per second.
}
}
//Wunderground Vars
//char SERVER = “rtupdate.wunderground.com”; //Rapidfire update server - for multiple sends per minute
char SERVER = “weatherstation.wunderground.com”; //Standard server - for sends once per minute or less
char WEBPAGE = “GET /weatherstation/updateweatherstation.php?”;
//Station Identification
char ID = “KNYESSEX3”; //Your station ID here
char PASSWORD = “XXRemovedXX”; //your Weather Underground password here
TCPClient client;
//Create Instance of HTU21D or SI7021 temp and humidity sensor and MPL3115A2 barometric sensor
Weather sensor;
//---------------------------------------------------------------
void setup()
{
pinMode(WSPEED, INPUT_PULLUP); // input from wind meters windspeed sensor
pinMode(RAIN, INPUT_PULLUP); // input from wind meters rain gauge sensor
Serial.begin(9600); // open serial over USB at 9600 baud
//Initialize the I2C sensors and ping them
sensor.begin();
/*You can only receive acurate barrometric readings or acurate altitiude
readings at a given time, not both at the same time. The following two lines
tell the sensor what mode to use. You could easily write a function that
takes a reading in one made and then switches to the other mode to grab that
reading, resulting in data that contains both acurate altitude and barrometric
readings. For this example, we will only be using the barometer mode. Be sure
to only uncomment one line at a time. */
sensor.setModeBarometer();//Set to Barometer Mode
//baro.setModeAltimeter();//Set to altimeter Mode
//These are additional MPL3115A2 functions the MUST be called for the sensor to work.
sensor.setOversampleRate(7); // Set Oversample rate
//Call with a rate from 0 to 7. See page 33 for table of ratios.
//Sets the over sample rate. Datasheet calls for 128 but you can set it
//from 1 to 128 samples. The higher the oversample rate the greater
//the time between data samples.
sensor.enableEventFlags(); //Necessary register calls to enble temp, baro ansd alt
// attach external interrupt pins to IRQ functions
attachInterrupt(RAIN, rainIRQ, FALLING);
attachInterrupt(WSPEED, wspeedIRQ, FALLING);
// turn on interrupts
interrupts();
}
//---------------------------------------------------------------
void loop()
{
//Keep track of which minute it is
if(millis() - lastSecond >= 1000)
{
lastSecond += 1000;
//Take a speed and direction reading every second for 2 minute average
if(++seconds_2m > 119) seconds_2m = 0;
//Calc the wind speed and direction every second for 120 second to get 2 minute average
float currentSpeed = get_wind_speed();
//float currentSpeed = random(5); //For testing
int currentDirection = get_wind_direction();
windspdavg[seconds_2m] = (int)currentSpeed;
winddiravg[seconds_2m] = currentDirection;
//if(seconds_2m % 10 == 0) displayArrays(); //For testing
//Check to see if this is a gust for the minute
if(currentSpeed > windgust_10m[minutes_10m])
{
windgust_10m[minutes_10m] = currentSpeed;
windgustdirection_10m[minutes_10m] = currentDirection;
}
//Check to see if this is a gust for the day
if(currentSpeed > windgustmph)
{
windgustmph = currentSpeed;
windgustdir = currentDirection;
}
if(++seconds > 59)
{
seconds = 0;
if(++minutes > 59) minutes = 0;
if(++minutes_10m > 9) minutes_10m = 0;
rainHour[minutes] = 0; //Zero out this minute’s rainfall amount
windgust_10m[minutes_10m] = 0; //Zero out this minute’s gust
}
//Get readings from all sensors
getWeather();
//Print to console
printInfo();
//Send data to Weather Underground
sendToWU();
//Power down between sends to save power, measured in seconds.
System.sleep(SLEEP_MODE_DEEP,300); //for Particle Photon
//Spark.sleep(SLEEP_MODE_DEEP,300); //for Spark Core
//delay(200)
}
}
//---------------------------------------------------------------
void printInfo()
{
//This function prints the weather data out to the default Serial Port
Serial.print(“Temp:”);
Serial.print(tempF);
Serial.print("F, ");
Serial.print(“Humidity:”);
Serial.print(humidity);
Serial.print("%, ");
Serial.print(“Baro_Temp:”);
Serial.print(baroTempF);
Serial.print("F, ");
Serial.print(“Humid_Temp:”);
Serial.print(humTempF);
Serial.print("F, ");
Serial.print(“Pressure:”);
Serial.print(pascals/100);
Serial.print("hPa, ");
Serial.print(inches);
Serial.println(“in.Hg”);
//The MPL3115A2 outputs the pressure in Pascals. However, most weather stations
//report pressure in hectopascals or millibars. Divide by 100 to get a reading
//more closely resembling what online weather reports may say in hPa or mb.
//Another common unit for pressure is Inches of Mercury (in.Hg). To convert
//from mb to in.Hg, use the following formula. P(inHg) = 0.0295300 * P(mb)
//More info on conversion can be found here:
//www.srh.noaa.gov/images/epz/wxcalc/pres … ersion.pdf
//If in altitude mode, print with these lines
//Serial.print(“Altitude:”);
//Serial.print(altf);
//Serial.println(“ft.”);
Serial.print(“Winddir:”);
switch (winddir)
{
case 0:
Serial.print(“North”);
break;
case 1:
Serial.print(“NE”);
break;
case 2:
Serial.print(“East”);
break;
case 3:
Serial.print(“SE”);
break;
case 4:
Serial.print(“South”);
break;
case 5:
Serial.print(“SW”);
break;
case 6:
Serial.print(“West”);
break;
case 7:
Serial.print(“NW”);
break;
default:
Serial.print(“No Wind”);
// if nothing else matches, do the
// default (which is optional)
}
Serial.print(" windspeedmph:");
Serial.print(windspeedmph, 1);
Serial.print("mph, ");
Serial.print(“Rainin:”);
Serial.print(rainin, 2);
Serial.print("in., ");
}
//---------------------------------------------------------------
void sendToWU()
{
Serial.println(“connecting…”);
if (client.connect(SERVER, 80)) {
Serial.println(“Connected”);
client.print(WEBPAGE);
client.print(“ID=”);
client.print(ID);
client.print(“&PASSWORD=”);
client.print(PASSWORD);
client.print(“&dateutc=now”); //can use ‘now’ instead of time if sending in real time
client.print(“&tempf=”);
client.print(tempF);
client.print(“&dewptf=”);
client.print(dewptF);
client.print(“&humidity=”);
client.print(humidity);
client.print(“&baromin=”);
client.print(inches);
client.print(“&winddir=”);
client.print(winddir);
client.print(“&windspeedmph=”);
client.print(windspeedmph);
client.print(“&rainin=”);
client.print(rainin);
client.print(“&dailyrainin=”);
client.print(dailyrainin);
client.print(“&windgustmph=”);
client.print(windgustmph);
client.print(“&wwindgustdir=”);
client.print(windgustdir);
client.print(“&windspdmph_avg2m=”);
client.print(windspdmph_avg2m);
client.print(“&winddir_avg2m=”);
client.print(winddir_avg2m);
client.print(“&windgustmph_10m=”);
client.print(windgustmph_10m);
client.print(“&windgustdir_10m=”);
client.print(windgustdir_10m);
client.print(“&action=updateraw”); //Standard update rate - for sending once a minute or less
//client.print(“&softwaretype=Particle-Photon&action=updateraw&realtime=1&rtfreq=5”); //Rapid Fire update rate - for sending multiple times per minute, specify frequency in seconds
client.println();
Serial.println(“Upload complete”);
delay(300); //Without the delay it goes to sleep too fast and the send is unreliable
}else{
Serial.println(F(“Connection failed”));
return;
}
}
//---------------------------------------------------------------
void getWeather()
{
// Measure Relative Humidity from the HTU21D or Si7021
humidity = sensor.getRH();
// Measure Temperature from the HTU21D or Si7021
humTempC = sensor.getTemp();
humTempF = (humTempC * 9)/5 + 32;
// Temperature is measured every time RH is requested.
// It is faster, therefore, to read it from previous RH
// measurement with getTemp() instead with readTemp()
//Measure the Barometer temperature in F from the MPL3115A2
baroTempC = sensor.readBaroTemp();
baroTempF = (baroTempC * 9)/5 + 32; //convert the temperature to F
//Measure Pressure from the MPL3115A2
pascals = sensor.readPressure();
inches = pascals * 0.0002953; // Calc for converting Pa to inHg (Wunderground expects inHg)
//If in altitude mode, you can get a reading in feet with this line:
//float altf = sensor.readAltitudeFt();
//Average the temperature reading from both sensors
tempC=((humTempC+baroTempC)/2);
tempF=((humTempF+baroTempF)/2);
//Calculate Dew Point
dewptC = dewPoint(tempC, humidity);
dewptF = (dewptC * 9.0)/ 5.0 + 32.0;
//Calc winddir
winddir = get_wind_direction();
//Calc windspeed
windspeedmph = get_wind_speed();
//Calc windgustmph
//Calc windgustdir
//Report the largest windgust today
windgustmph = 0;
windgustdir = 0;
//Calc windspdmph_avg2m
float temp = 0;
for(int i = 0 ; i < 120 ; i++)
temp += windspdavg*;*
temp /= 120.0;
windspdmph_avg2m = temp;
//Calc winddir_avg2m
temp = 0; //Can’t use winddir_avg2m because it’s an int
for(int i = 0 ; i < 120 ; i++)
temp += winddiravg;
temp /= 120;
winddir_avg2m = temp;
//Calc windgustmph_10m
//Calc windgustdir_10m
//Find the largest windgust in the last 10 minutes
windgustmph_10m = 0;
windgustdir_10m = 0;
//Step through the 10 minutes
for(int i = 0; i < 10 ; i++)
{
if(windgust_10m > windgustmph_10m)
{
windgustmph_10m = windgust_10m;
windgustdir_10m = windgustdirection_10m;
}
}
//Total rainfall for the day is calculated within the interrupt
//Calculate amount of rainfall for the last 60 minutes
rainin = 0;
for(int i = 0 ; i < 60 ; i++)
rainin += rainHour;
}
//---------------------------------------------------------------
// dewPoint function from NOAA
// reference (1) : http://wahiduddin.net/calc/density_algorithms.htm
// reference (2) : http://www.colorado.edu/geography/weath … /about.htm
//---------------------------------------------------------------
double dewPoint(double celsius, double humidity)
{
* // (1) Saturation Vapor Pressure = ESGG(T)*
* double RATIO = 373.15 / (273.15 + celsius);*
_ double RHS = -7.90298 * (RATIO - 1);
RHS += 5.02808 * log10(RATIO);
RHS += -1.3816e-7 * (pow(10, (11.344 * (1 - 1/RATIO ))) - 1) ;
RHS += 8.1328e-3 * (pow(10, (-3.49149 * (RATIO - 1))) - 1) ;
* RHS += log10(1013.246);
// factor -3 is to adjust units - Vapor Pressure SVP * humidity_
_ double VP = pow(10, RHS - 3) * humidity;
// (2) DEWPOINT = F(Vapor Pressure)
double T = log(VP/0.61078); // temp var*
return (241.88 * T) / (17.558 - T);
}
//---------------------------------------------------------------
//Read the wind direction sensor, return heading in degrees
int get_wind_direction()
{
unsigned int adc;
adc = analogRead(WDIR); // get the current reading from the sensor
// The following table is ADC readings for the wind direction sensor output, sorted from low to high.
// Each threshold is the midpoint between adjacent headings. The output is degrees for that ADC reading.
// Note that these are not in compass degree order! See Weather Meters datasheet for more information.
//Wind Vains may vary in the values they return. To get exact wind direction,
//it is recomended that you AnalogRead the Wind Vain to make sure the values
//your wind vain output fall within the values listed below.
if(adc > 2270 && adc < 2290) return (0);//North
if(adc > 3220 && adc < 3299) return (1);//NE
if(adc > 3890 && adc < 3999) return (2);//East
if(adc > 3780 && adc < 3850) return (3);//SE
if(adc > 3570 && adc < 3650) return (4);//South
if(adc > 2790 && adc < 2850) return (5);//SW
if(adc > 1580 && adc < 1610) return (6);//West
if(adc > 1930 && adc < 1950) return (7);//NW
return (-1); // error, disconnected?
}
//---------------------------------------------------------------
//Returns the instataneous wind speed
float get_wind_speed()
{
float deltaTime = millis() - lastWindCheck; //750ms
deltaTime /= 1000.0; //Covert to seconds
float windSpeed = (float)windClicks / deltaTime; //3 / 0.750s = 4
windClicks = 0; //Reset and start watching for new wind
lastWindCheck = millis();
windSpeed *= 1.492; //4 * 1.492 = 5.968MPH
/* Serial.println();
Serial.print(“Windspeed:”);
Serial.println(windSpeed);*/
return(windSpeed);
}
//---------------------------------------------------------------_