return curdate
# today's date is default.
if len(sys.argv) > 1:
YYYYMMDD = sys.argv[1]
else:
YYYYMMDD = datetime.datetime.today().strftime('%Y%m%d')
YYYYMM = YYYYMMDD[0:6]
# set OpenDAP server URL.
HH='09'
URLbase="http://nomad3.ncep.noaa.gov:9090/dods/sref/sref"
URL=URLbase+YYYYMMDD+"/sref_eta_ctl1_"+HH+"z"
print URL+'\n'
try:
data = client.open(URL)
except:
msg = """
opendap server not providing the requested data.
Try another date by providing YYYYMM on command line."""
raise IOError, msg
# read levels, lats,lons,times.
print data.keys()
levels = data['lev']
latitudes = data['lat']
longitudes = data['lon']
fcsttimes = data['time']
times = fcsttimes[:]
YYYYMMDDHH2 = '1993031700'
YYYY = YYYYMMDDHH1[0:4]
if YYYY != YYYYMMDDHH2[0:4]:
raise ValueError,'dates must be in same year'
# set OpenDAP server URL.
URLbase="http://nomad3.ncep.noaa.gov:9090/dods/reanalyses/reanalysis-2/6hr/pgb/"
URL=URLbase+'pres'
URLu=URLbase+'wind'
URLv=URLbase+'wind'
print URL
print URLu
print URLv
try:
data = client.open(URL)
datau = client.open(URLu)
datav = client.open(URLv)
except:
raise IOError, 'opendap server not providing the requested data'
# read lats,lons,times.
print data.keys()
print datau.keys()
print datav.keys()
latitudes = data['lat'][:]
longitudes = data['lon'][:].tolist()
times = data['time'][:]
# put times in YYYYMMDDHH format.
dates=[]
for t in times:
YYYYMMDDHH2 = '1993031700'
YYYY = YYYYMMDDHH1[0:4]
if YYYY != YYYYMMDDHH2[0:4]:
raise ValueError, 'dates must be in same year'
# set OpenDAP server URL.
URLbase = "http://nomad3.ncep.noaa.gov:9090/dods/reanalyses/reanalysis-2/6hr/pgb/"
URL = URLbase + 'pres'
URLu = URLbase + 'wind'
URLv = URLbase + 'wind'
print URL
print URLu
print URLv
try:
data = client.open(URL)
datau = client.open(URLu)
datav = client.open(URLv)
except:
raise IOError, 'opendap server not providing the requested data'
# read lats,lons,times.
print data.keys()
print datau.keys()
print datav.keys()
latitudes = data['lat'][:]
longitudes = data['lon'][:].tolist()
times = data['time'][:]
# put times in YYYYMMDDHH format.
dates = []
for t in times:
