I have some portable weather stations that are using the CR10X datalogger, and we are updating to the CR1000 because we need some additional functionality.

First, a description of our weather stations. We are setting up portable stations similar to the RAWS systems, but ours are used for flight testing, often in remote locations where there is no other weather data available. Our systems are set up temporarily, monitored in real time (using RTMC), and then torn down at the end of the day with data parsed by SPLIT and saved for historical data.

We have the following sensors:
HMP45C (temp, rh)
CM3 (solar radiation)
CS105 (barometric pressure)
05103 wind sensors with 32500 electronic compass

We'll be updating the wind sensors soon, likely to RM Young 85000 series ultrasonic sensors.

We display measurements of wind speed/direction and solar radiation once every second. Every 60 seconds I display temp, RH, and pressure. Every five minutes I display gust data, prevailing wind direction, and standard deviation of wind direction for the past five minutes.

I have all this working on the CR10X and ported over to the CR1000 as well just fine, but I would like to make some refinements on the CR1000. List of things I want to change:

1. Add GPS for the following:
a. position data, only needs to be stored once per test setup.
b. time data, used to set station time to local time. This presents a couple issues (how does end user set UTC offset without editing .CR1 file? Is daylight saving time taken into consideration automatically, or need to be set by user?)
c. Use $GPRMC string to get Magnetic Variation, which would then allow true north to be displayed instead of magnetic north.

2. I'd like to change the five minute averages and standard deviation to be a rolling five minutes instead of 'chunk' data. I assume this is possible with the CR1000 but haven't started working on it yet.

Here's where I'm stuck. I have a test program working for the GPS data, using the samples in the GPS16X-HVS manual. There were some errors in the example that I fixed, but I am still having a problem with running the examples concurrently. It seems I can either store the GPS data in an array and then use the GPS command in the main program to set station time first example), or I can store the individual GPS strings in table format (second example) but if I use the GPS command in the main program, this breaks the tables and fills them with empty data.

I also noticed that there is a problem with the second example in that it is looking for the position of data in string, not for the comma delimiters. This means that if any data is absent, all of the subsequent pieces are incorrect. I'd rather not use this method as it appears to be unreliable.

So how do I write the arrays from the first example (page 8,9 in the manual) to a data table as in the second example (page 10,11)? I can't seem to get it to work.

[Report Inappropriate or Erroneous Content]

I've answered my last question in the above post with regard to how to get the data out of the array in example one. I still need answers to the other questions, but now I've got a place to start.

This is the working GPS program (adapted from example 1, with the fixes).

'CR1000 Series Datalogger

'GPS Variables
  Const LOCAL_TIME_OFFSET = -5         'Local Time offset relative to UTC time
  
  Dim nmea_sentence(2) As String * 90
  Public gps_data(15)
    Alias gps_data(1) = latitude_deg    'Degrees latitude (+ = North; - = South)
    Alias gps_data(2) = latitude_min    'Minutes latitude
    Alias gps_data(3) = longitude_deg   'Degress longitude (+ = East; - = West)
    Alias gps_data(4) = longitude_min   'Minutes longitude
    Alias gps_data(5) = speed           'Speed
    Alias gps_data(6) = course          'Course over ground
    Alias gps_data(7) = mag_var         'Magnetic variation from true north (+ = East; - = West)
    Alias gps_data(8) = fix_quality     'GPS fix quality: 0 = invalid, 1 = GPS, 2 = 'differential GPS, 6 = estimated
    Alias gps_data(9) = nmbr_sat        'Number of satellites used for fix
    Alias gps_data(10) = altitude       'Antenna altitude
    Alias gps_data(11) = pps            'Elapsed ms since last pulse per second (PPS) from GPS
    Alias gps_data(12) = dt_since_gprmc 'Time since last GPRMC string, normally less than '1 second
    Alias gps_data(13) = gps_ready      'Counts from 0 to 10, 10 = ready
    Alias gps_data(14) = max_clock_change 'Maximum value the clock was changed
    Alias gps_data(15) = nmbr_clock_change 'Number of times the clock was changed

    'Define Units to be used in data file header
    Units latitude_deg = degrees
    Units latitude_min = minutes
    Units longitude_deg = degrees
    Units longitude_min = minutes
    Units speed = m/s
    Units course = degrees
    Units mag_var = degrees '(+ = East; - = West)
    Units fix_quality = unitless
    Units nmbr_sat = unitless
    Units altitude = m
    Units pps = ms
    Units dt_since_gprmc = s
    Units gps_ready = unitless
    Units max_clock_change = ms
    Units nmbr_clock_change = samples


DataTable (gpsdata,True,-1)
  DataInterval (0,1,Sec,10)
     Sample (1,latitude_deg,String)
     Sample (1,latitude_min,String)
     Sample (1,longitude_deg,String)
     Sample (1,longitude_min,String)
     Sample (1,speed,String)
     Sample (1,course,String)
     Sample (1,mag_var,String)
     Sample (1,fix_quality,String)
     Sample (1,nmbr_sat,String)
     Sample (1,altitude,String)
     Sample (1,pps,String)
     Sample (1,dt_since_gprmc,String)
     Sample (1,gps_ready,String)
     Sample (1,max_clock_change,String)
     Sample (1,nmbr_clock_change,String)
EndTable

'Main Program
BeginProg
  Scan (1,Sec,0,0)
     GPS (gps_data(),Com4,LOCAL_TIME_OFFSET*3600,0,nmea_sentence(1))
     CallTable gpsdata
  NextScan
EndProg

* Last updated by: MichaelM on 8/3/2010 @ 9:59 AM *

[Report Inappropriate or Erroneous Content]

position data, only needs to be stored once per test setup.
This will read the GPS once and store the Lat and Lon for you. If you want to update, just have something make GPSread false. I used Port4 to turn the GPS on and off when not in use, it draws ~70 mA if memory serves me correctly.

'Routine to Read GPS
If GPSread = False Then
If gps_ready <> 10 Then
PortSet (4,0 ) 'Turn on GPS
GPS (latitude_a,Com3,UTC_OFFSET,0,nmea_sentence(1)) ' Read GPS
ElseIf gps_ready = 10 Then
Lat = latitude_a + (latitude_b / 60) 'Calculate Decimal Latitude
Lon = longitude_a + (longitude_b / 60) 'Calculate Decimal Longitude
Alt = altitude
PortSet (4,1 ) 'Turn off GPS
GPSread = True 'Mark GPS as read
EndIf
EndIf

Change this:
Const LOCAL_TIME_OFFSET = -5
To This:
ConstTable
Const LOCAL_TIME_OFFSET = -5
EndConstTable

It will allow the customer to change the offset with a keyboard display, BUT BE CAREFUL, it will recompile the program so make sure to save data if required before you change.

Magnetic Variation:
Alias gps_data(7) = mag_var

I think that may answer the questions.....

Good Luck!

* Last updated by: IslandMan on 8/4/2010 @ 4:31 AM *

[Report Inappropriate or Erroneous Content]

Islandman, thank you. The routine for the GPS appears to work perfectly for my needs.

I can make the ConstTable thing work too, but it won't be pretty. I'll probably use AutoHotKey to write a script to see if I can emulate the terminal with minimal user interface. I think it can be done.

I've already got the Magnetic Variation working shortly after I posted my original question and response.

I'll keep this thread updated as I come across other issues/solutions.

[Report Inappropriate or Erroneous Content]

If you explain a little more about how you want to execute the 5 minute running average I may be able to pull some code for that as well.

[Report Inappropriate or Erroneous Content]

We are looking at 'gust' data (maximum wind speed kn knots) and prevailing wind direction with the standard deviation over five minute intervals.

With the CR10X, we were limited to chunks. In other words, our five minute averages were, for example, at 11:00, 11:05, 11:10, etc.

What I would like is to know the prevailing wind direction and maximum wind speed over the last five minutes at any given time.

So, once the station is powered up, it would take five minutes to compute the first averages (term used loosely to describe both the gust and the prevailing wind direction with standard deviation), then after that it would constantly compute this data, updated once per second, going back the last five minutes. Make sense?

I appreciate your help in this regard.

[Report Inappropriate or Erroneous Content]

Michael,

I will let Dave answer your running average question (since he's probably got it half written by now...), but also take a look at the custom menu capability in the datalogger, and in particular the MenuRecompile instruction. It will let you add a pick list for a Constant Table item, which may allow for a little easier UI for the end user.

Dana W.

[Report Inappropriate or Erroneous Content]

Michael,

Give this a try. You have to add the Wind sensor measurements.

'CR1000 Series Datalogger

'Declare Public Variables
Public Wspd
Public WspdGust(300)
Public GustCount
Public Wdir
Public SinWdir
Public CosWdir
Public AvgSinWdir
Public AvgCosWdir
Public e
Public AvgWspd
Public AvgWdir
Public Sigma
Public Valid
Public MaxGust(2)
Alias MaxGust(1) = Gust

'define all angles in degrees not radians
AngleDegrees

'Main Program
BeginProg

'Init Variables
Valid = 0
GustCount = 1

Scan (1,Sec,0,0)

' Do not output Averages or Sigma until all 300 samples are acquired
If Valid < 300 Then
AvgWspd = 0
AvgWdir = 0
Sigma =0
Gust = 0
Valid = Valid + 1
EndIf

'Measure Wind Speed an Direction Sensor

'Store Wind Speed readings in array for '"Gust" calculation
If Gustcount < 300 Then
WspdGust(Gustcount) = Wspd
Gustcount = Gustcount + 1
Else
Gustcount = 1
EndIf

SinWdir = SIN(Wdir) ' 10M WD Sine
CosWdir = COS(Wdir) ' 10M WD Cosine

AvgRun (AvgWspd,1,Wspd,300) 'Average Wind Speed
AvgRun (AvgSinWdir,1,SinWdir,300) 'Average Wdir Sine
AvgRun (AvgCosWdir,1,CosWdir,300) 'Average Wdir Cosine

'Calculate 5 Minute Running averages for Wind Direction
AvgWdir = ATN(AvgSinWdir/AvgCosWdir)
If AvgCosWdir < 0 Then
AvgWdir = 180 + AvgWdir
ElseIf AvgSinWdir < 0 Then
AvgWdir = 360 + AvgWdir
EndIf

'Calculate 5 Minute Running Sigmas for Wind Direction (yamartino method)
e = SQR (1-((AvgSinWdir^2)+(AvgCosWdir^2)))
Sigma = ASIN (e*(1+(0.1547*(e^3))))

'find Max Wind Speed (GUST) over 300 samples
MaxSpa (Gust,300,WspdGust())

NextScan
EndProg

[Report Inappropriate or Erroneous Content]

Dave,

Thanks again for sharing your expertise!

[Report Inappropriate or Erroneous Content]

I'm bumping an old thread with a new (related) question.

After some testing in the field and some toying with data, I'd like to read the GPS more than once, something on the order of once every 15 minutes or so. The trick is, I don't want to start recording any data for any of the other sensors unless and until there is an initial GPS lock, then start recording from there.

Essentially (leaving out the non-relevant parts):

'Declares
Public GPSread = False
Public GPSInit = False

'Tables
DataTable (gpsdata,True,-1)
'No DataInterval since calling table conditionally
Sample (1,lat,string)
Sample (1,lon,string)
Sample (etc., etc. etc.)
EndTable

DataTable (1_Sec_Data,True,-1)
DataInterval (0,1,Sec,10)
Sample
Minimum
etc.
EndTable

DataTable (1_Min_Data,True,-1)
DataInterval (0,1,Min,10)
Sample
Minimum
etc.
EndTable

'Main Program
BeginProg
Scan (500,mSec,10,0)
'Routine to Read GPS
If GPSread = False Then
If gps_ready <> 10 Then
PortSet (3,0) 'Turn on GPS
GPS (gps_data(),Com4,UTC_OFFSET,0,nmea_sentence(1)) ' Read GPS
ElseIf gps_ready = 10 Then
Lat = latitude_deg + (latitude_min / 60) 'Calculate Decimal Latitude
Lon = longitude_deg + (longitude_min / 60) 'Calculate Decimal Longitude
Alt = altitude
PortSet (3,1) 'Turn off GPS
CallTable gpsdata
GPSread = True 'Mark GPS as read
GPSInit = True 'Start recording all other sensors
EndIf
EndIf

'Measure One Second data here
'...

'Store Data in 1_Sec_Data table
If GPSInit = True Then CallTable 1_Sec_Data

'Measure One Minute data here
'...

'Store Data in 1_Min_Data table
If GPSInit = True Then CallTable 1_Min_Data

'Here's where I run into problems
If IfTime (14,15,Min) Then GPSread = False

NextScan
EndProg

The GPS data in the Public table (where it is initially stored) is never updated after the first run, so all records in the DataTable gpsdata are copies of the original reading.

I'm sure I'm missing something simple here, but I can't figure where I'm going wrong. Everything else in my code works just as I expect it to, only the subsequent GPS reads are currently not working.

* Last updated by: MichaelM on 2/4/2011 @ 12:17 PM *

[Report Inappropriate or Erroneous Content]