wind meter idl routines

Last modified: Thu Mar 16 13:07:14 2023.


List of Routines


Routine Descriptions

WINDDIR - RETURN THE PATH TO THE WIND DIRECTORY

[Next Routine] [List of Routines]
NAME:
winddir - return the path to the wind directory

SYNTAX: dirname=winddir()
ARGS:
RETURNS:
	dirname: string directory name for wind data. It includes the trailing
                   /
DESCRIPTION:
	Return the directory name for the wind directory:
EXAMPLE:
   print,winddir()
   /share/megs2_u1/wind/

(See /pkg/rsi/local/libao/phil/wind/winddir.pro)


WINDGETMONTH - INPUT A MONTHS WORTH OF DATA

[Previous Routine] [Next Routine] [List of Routines]
NAME:
windgetmonth - input a months worth of data

SYNTAX: nrecs=windgetmonth(year,month,wd,smo=smo)
ARGS:
	year	: int 4 digit year
   month   : 1 through 12 month of year
KEYWORDS
		smo : int smooth and decimate the data to this many seconds
			      the default is 1 second resolution
RETURNS:
	wd[nrecs]:{windstr} data returned
   nrecs   : long  number or records found
			  -1 this months data not available
DESCRIPTION:
	Input a months worth of wind data. It is return in the array of
wind structures d. Each element of d contains:
** Structure WINDSTR, 3 tags, length=16:
   JULDAY DOUBLE   Julian day of measurement.AST noon rather than utc noon
					starts the day.
   VEL    FLOAT    wind velocity in mph
   DIR    FLOAT    direction in degree from where the wind is blowing.

(See /pkg/rsi/local/libao/phil/wind/windgetmonth.pro)


WINDINPRAW - INPUT THE RAW WIND MONITOR DATA.

[Previous Routine] [Next Routine] [List of Routines]
NAME:
windinpraw - input the raw wind monitor data.

SYNTAX: nrecs=windinpraw(yymmdd,b,daynum=daynum,year=year,$
						  append=append,recsToSkip=recsToSkip)
ARGS:
	yymmdd: long day to read in. -1 is current day. this can be overridden
				 with the daynum,year keywords.  
	b[]	  :{windstr} return the data here. If the append keyword is used
			     then the user should pass in the current buffer and the 
				 new data will be appended to it.
KEYWORDS:
   daynum: long daynumber of year for data to read.
   year  : long if daynumber is used then this is the corresponding year.
				 default is the current year.
 recsToSkip: long records to skip on input. Use to skip over bad data..
				
DESCRIPTION:
   This routine is called by windmon to input the next set of monitor
data. This reads the raw ascii data written by the wind monitor 
(regular or compressed files). 

 It returns number of recs read. returns -1 if it cannot open the file.
SEE ALSO:
 WIND
NOTES:
 	The jd time stamp in the log file is not accurate. It is coming from
the sbc clock which wanders by up to 1/2 hour. The hhmmss time stamp
is accurate since it comes from obsdisplay which runs ntpd

(See /pkg/rsi/local/libao/phil/wind/windinpraw.pro)


WINDMKMON - MAKE THE MONTH ARCHIVE

[Previous Routine] [Next Routine] [List of Routines]
NAME:
windmkmon - make the month archive

SYNTAX: nrecs=windmkmon(yy,mon,d)
ARGS:
	year	: int year to get
	mon	    : int month of year (1..12)
RETURNS:
	d[nrecs]:{windstr} data returned for this month
   nrecs   : long  number or records found
			  -1 this months data not available
DESCRIPTION:
	Input a months worth of data from the windmeter raw data archive.
Create the idl structure {windstr} and then write this out to the
idl save file ('wind_yymm.sav'). The output directory is in 
/share/megs2_u1/wind/

(See /pkg/rsi/local/libao/phil/wind/windmkmon.pro)


WINDPLOTDAY - PLOT A DAYS WORTH OF WIND DATA.

[Previous Routine] [Next Routine] [List of Routines]
NAME:
windplotday - plot a days worth of wind data.
SYNTAX: windplotday,d
ARGS:
 d[n]:{}	wind data input by windinpraw()
RETURNS:
		plots the wind data.
DESCRIPTION:
	Plot the wind data for a day. The top plot is the wind data versus hour 
of data (ast). The data has bin binned to 1 minute steps. The white (black)
plot is the average over each minute. The red plot is the peak hold over
each minute.
	The bottom plot is the median wind velocity for each minute versus 
ast hour of day. The direction is where the wind is coming from.

EXAMPLE:
	Get the data and then plot it.

 yymmdd=070412
 n=windinpraw(yymmdd,d)
 windplotday,d

(See /pkg/rsi/local/libao/phil/wind/windplotday.pro)


WINDPLOTDEAN - PLOT A DAYS WORTH OF WIND DATA.

[Previous Routine] [Next Routine] [List of Routines]
NAME:
windplotdean - plot a days worth of wind data.
SYNTAX: windplotday,d
ARGS:
 d[n]:{}	wind data input by windinpraw()
RETURNS:
		plots the wind data.
DESCRIPTION:
	Plot the wind data for a day. The top plot is the wind data versus hour 
of data (ast). The data has bin binned to 1 minute steps. The white (black)
plot is the average over each minute. The red plot is the peak hold over
each minute.
	The bottom plot is the median wind velocity for each minute versus 
ast hour of day. The direction is where the wind is coming from.

EXAMPLE:
	Get the data and then plot it.

 yymmdd=070412
 n=windinpraw(yymmdd,d)
 windplotday,d

(See /pkg/rsi/local/libao/phil/wind/windplotdean.pro)


WINDPLOTEARL - PLOT A DAYS WORTH OF WIND DATA.

[Previous Routine] [Next Routine] [List of Routines]
NAME:
windplotearl - plot a days worth of wind data.
SYNTAX: windplotearl,d
ARGS:
 d[n]:{}	wind data input by windinpraw()
RETURNS:
		plots the wind data.
DESCRIPTION:
	Plot the wind data for a day. The top plot is the wind data versus hour 
of data (ast). The data has bin binned to 1 minute steps. The white (black)
plot is the average over each minute. The red plot is the peak hold over
each minute.
	The bottom plot is the median wind velocity for each minute versus 
ast hour of day. The direction is where the wind is coming from.

EXAMPLE:
	Get the data and then plot it.

 yymmdd=070412
 n=windinpraw(yymmdd,d)
 windplotday,d

(See /pkg/rsi/local/libao/phil/wind/windplotearl.pro)


WINDPLOTSTORM - PLOT A DAYS WORTH OF WIND DATA.

[Previous Routine] [Next Routine] [List of Routines]
NAME:
windplotstorm - plot a days worth of wind data.
SYNTAX: windplotstorm,d
ARGS:
 d[n]:{}	wind data input by windinpraw()
RETURNS:
		plots the wind data.
DESCRIPTION:
	Plot the wind data for a day. The top plot is the wind data versus hour 
of data (ast). The data has bin binned to 1 minute steps. The white (black)
plot is the average over each minute. The red plot is the peak hold over
each minute.
	The bottom plot is the median wind velocity for each minute versus 
ast hour of day. The direction is where the wind is coming from.

EXAMPLE:
	Get the data and then plot it.

 yymmdd=070412
 n=windinpraw(yymmdd,d)
 windplotday,d

(See /pkg/rsi/local/libao/phil/wind/windplotstorm.pro)


WINDPLOTYR - MAKE THE WINDMETER YEARLY SUMMARY PLOTS

[Previous Routine] [Next Routine] [List of Routines]
NAME:
windplotyr - make the windmeter yearly summary plots
SYNTAX: windplotyr,year
ARGS:
 year: int	year to plot
DESCRIPTION:
	Make the yearly summary wind plots for the specified year.

(See /pkg/rsi/local/libao/phil/wind/windplotyr.pro)


WINDPLTDATE - PLOT WIND DATA VERSUS DATE

[Previous Routine] [List of Routines]
NAME:
windpltdate - plot wind data versus date
SYNTAX: windplot,wd,h1=h1,h2=h2,datef=datef,ast=ast
ARGS:
	wd[n]	: {windstr} wind data to plot
KEYWORDS:
    h1: double yymmdd.fract date for left edge of plot. default is the start of
               the data.
    h2: double yymmdd.fract date for right edge of plot. default is the end of
               the data.
	AST:	    if set the subtract 4./24. so the date is ast rather the gmt.
datef : string date format to use when plotting. The default is 
               '%D%M%Z' (see label_date() routine).

DESCRIPTION:
	Plot the window velocity vs date for the wind data structure. This can
be for multiple days.

(See /pkg/rsi/local/libao/phil/wind/windpltdate.pro)