[gradsusr] Extract time series from gridded data

[Lydia Rill]

Hi all,

I am new to using GrADS. Currently I have a script that uses the function
gr2stn to pull the nearest grid point from the latitude and longitude I
specify from the gridded data in a large binary file.
Because I want a time series, I have t set to vary, while x and y are set
to 1.
I grab the values from the result using sublin and subwrd. I do this for
each of the 4 variables that I want separately, and I have this set in a
loop to get data for every day over multiple years.

This script works, but it is very slow as grabbing the daily data for 35
years takes almost a half hour!

Is there a more efficient way to pull a time series of a point from gridded
data in a binary file?

Other information:
The gridded data is NLDAS2 forcing A which I aggregated from hourly to
daily time step.
The daily gridded data is in a 22G binary file.
I am running GrADS 2.0.2 on a Mac OSX, and I call the attached script from
a bash script.

Thank you for your help!

Lydia

-- 
Lydia D. Rill
M.S., Michigan State University 2016
B.S., Valparaiso University 2014
(224) 406-5130
Lydia.D.Rill at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160628/0b73485d/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: station.gs
Type: application/octet-stream
Size: 2673 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20160628/0b73485d/attachment.obj