- Sat Apr 01, 2017 2:04 pm
#194226
I have two questions / issues
1. I am able to get the wind speed and wind direction, but I never get wind gust / direction readings...
2. I have the wind settings setup for 8 directions, but understand that there are 16 available directions. I am unable to find a adc chart to match up with normal compass readings. So I am missing the degrees 22.5, 67.5, 112.5, 157.5, 202.5, 247.5, 292.5, 337.5 .
Here is the code I am using, pieced together from here and there
1. I am able to get the wind speed and wind direction, but I never get wind gust / direction readings...
2. I have the wind settings setup for 8 directions, but understand that there are 16 available directions. I am unable to find a adc chart to match up with normal compass readings. So I am missing the degrees 22.5, 67.5, 112.5, 157.5, 202.5, 247.5, 292.5, 337.5 .
Here is the code I am using, pieced together from here and there
Code: Select all
/******************************************************************************
SparkFun_Photon_Weather_Wunderground.ino
SparkFun Photon Weather Shield basic example
Joel Bartlett @ SparkFun Electronics
Original Creation Date: May 18, 2015
Updated August 21, 2015
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_Photon_Weather_Shield_Particle_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/PWS_-_Upload_Protocol
Sign up at http://www.wunderground.com/personal-weather-station/signup.asp
*******************************************************************************/
#define RBG_NOTIFICATIONS_OFF //Turn off LED
void ledOff() {
RGB.control(true);
RGB.color(0,0,0);
RGB.brightness(0);
}
STARTUP(Time.setTime(0));
#if defined(RBG_NOTIFICATIONS_OFF)
STARTUP(ledOff();)
#endif
STARTUP(WiFi.selectAntenna(ANT_EXTERNAL));
#include "SparkFun_Photon_Weather_Shield_Library.h"
#include "math.h" //For Dew Point Calculation
int WDIR = A0;
int RAIN = D2;
int WSPEED = D3; //march2017
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] march2017
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;
//Wunderground Vars
char SERVER[] = "rtupdate.wunderground.com"; //Rapidfire update server - for multiple sends per minute
//char SERVER [] = "/rtupdate.wunderground.com/weatherstation/updateweatherstation.php?"; //Standard server - for sends once per minute or less
char WEBPAGE [] = "GET /weatherstation/updateweatherstation.php?";
//Station Identification
char ID [] = "KTNALAMO3"; //Your station ID here
char PASSWORD [] = "passwordhere"; //your Weather Underground password here
TCPClient client;
// volatiles are subject to modification by IRQs
volatile long lastWindIRQ = 0; //march 2017
volatile byte windClicks = 0;
volatile unsigned long raintime, rainlast, raininterval, rain; // march2017
//Create Instance of HTU21D or SI7021 temp and humidity sensor and MPL3115A2 barometric sensor
Weather sensor;
//Interrupt routines (these are called by the hardware interrupts, not by the main code)
//-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=march2017
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.
}
}
//march2017
//---------------------------------------------------------------
void setup()
{
Particle.connect();
Serial.begin(9600); // open serial over USB at 9600 baud
//Initialize the I2C sensors and ping them
sensor.begin();
pinMode(WSPEED, INPUT_PULLUP); // input from wind meters windspeed sensor
pinMode(RAIN, INPUT_PULLUP); // input from wind meters rain gauge sensor march 2017
/*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
seconds = 0;
lastSecond = millis();
// attach external interrupt pins to IRQ functions
attachInterrupt(RAIN, rainIRQ, FALLING);
attachInterrupt(WSPEED, wspeedIRQ, FALLING);
// turn on interrupts
interrupts();
}
//---------------------------------------------------------------
void loop()
{
//Get readings from all sensors
getWeather();
//Print to console //console off during wireless
//printInfo();
//Send data to Weather Underground
sendToWU();
//Power down between sends to save power, measured in seconds.
//System.sleep(SLEEP_MODE_DEEP,25); //for Particle Photon
//Spark.sleep(SLEEP_MODE_DEEP,300); //for Spark Core
WiFi.off();
delay(25000);
WiFi.connect();
}
//---------------------------------------------------------------
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/pressureConversion.pdf
//If in altitude mode, print with these lines
//Serial.print("Altitude:");
//Serial.print(altf);
//Serial.println("ft.");
//This function prints the weather data out to the default Serial Port
Serial.print("Wind_Dir:");
switch (winddir)
{
case 0:
Serial.print("North");
break;
case 45:
Serial.print("NE");
break;
case 90:
Serial.print("East");
break;
case 135:
Serial.print("SE");
break;
case 180:
Serial.print("South");
break;
case 225:
Serial.print("SW");
break;
case 270:
Serial.print("West");
break;
case 315:
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., ");
}
//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 (45);//NE
if(adc > 3890 && adc < 3999) return (90);//East
if(adc > 3780 && adc < 3850) return (135);//SE
if(adc > 3570 && adc < 3650) return (180);//South
if(adc > 2790 && adc < 2850) return (225);//SW
if(adc > 1580 && adc < 1610) return (270);//West
if(adc > 1930 && adc < 1950) return (315);//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);
}
//---------------------------------------------------------------
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-5);
client.print("&dewptf=");
client.print(dewptF);
client.print("&humidity=");
client.print(humidity);
client.print("&baromin=");
client.print(inches);
client.print("&rainin=");
client.print(rainin);
client.print("&winddir="); //- [0-360 instantaneous wind direction]
client.print(winddir);
client.print("&windspeedmph="); //- [mph instantaneous wind speed]
client.print(windspeedmph);
client.print("&windgustmph="); //- [mph current wind gust, using software specific time period]
client.print(windgustmph);
client.print("&windgustdir="); //- [0-360 using software specific time period]
client.print(windgustdir);
client.print("&windspdmph_avg2m="); //- [mph 2 minute average wind speed mph]
client.print(windspdmph_avg2m);
client.print("&winddir_avg2m="); //- [0-360 2 minute average wind direction]
client.print(winddir_avg2m);
client.print("&windgustmph_10m="); // - [mph past 10 minutes wind gust mph ]
client.print(windgustmph_10m);
client.print("&windgustdir_10m=");// - [0-360 past 10 minutes wind gust direction]
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=30"); //Rapid Fire update rate - for sending multiple times per minute, specify frequency in seconds
client.println(" HTTP 1.1"); // note the space before HTTP
client.println("Host: rtupdate.wunderground.com");
client.println("Connection: close");
client.println();
Serial.println("Upload complete");
delay(1000); //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[i];
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[i];
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[i] > windgustmph_10m)
{
windgustmph_10m = windgust_10m[i];
windgustdir_10m = windgustdirection_10m[i];
}
}
//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[i];
}
//---------------------------------------------------------------
// dewPoint function from NOAA
// reference (1) : http://wahiduddin.net/calc/density_algorithms.htm
// reference (2) : http://www.colorado.edu/geography/weather_station/Geog_site/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);
// Delay for 10 sec before next reading // added by JDM
//delay(10000);
}
//---------------------------------------------------------------