Fix URLs in examples using broken addcyclic

Author: molinav

examples/fcstmaps.py

@@ -3,27 +3,25 @@
 from __future__ import unicode_literals
 # this example reads today's numerical weather forecasts
 # from the NOAA OpenDAP servers and makes a multi-panel plot.
+import datetime as dt
 import numpy as np
 import matplotlib.pyplot as plt
 import sys
 import numpy.ma as ma
-import datetime
 from mpl_toolkits.basemap import Basemap, addcyclic
 from netCDF4 import Dataset as NetCDFFile, num2date
 
 
 # today's date is default.
 if len(sys.argv) > 1:
-    YYYYMMDD = sys.argv[1]
+    date = dt.datetime.strptime(sys.argv[1], "%Y%m%d")
 else:
-    YYYYMMDD = datetime.datetime.today().strftime('%Y%m%d')
+    date = dt.datetime.today()
 
 # set OpenDAP server URL.
 try:
-    URLbase="http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs"
-    URL=URLbase+YYYYMMDD+'/gfs_00z'
-    print(URL)
-    data = NetCDFFile(URL)
+    urlbase = "http://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs%Y%m%d/gfs_0p25_00z"
+    data = NetCDFFile(date.strftime(urlbase))
 except:
     msg = """
 opendap server not providing the requested data.

examples/fcstmaps_axesgrid.py

@@ -4,35 +4,32 @@
 # this example reads today's numerical weather forecasts
 # from the NOAA OpenDAP servers and makes a multi-panel plot.
 # This version demonstrates the use of the AxesGrid toolkit.
+import datetime as dt
 import numpy as np
 import matplotlib.pyplot as plt
 import sys
 import numpy.ma as ma
-import datetime
 from mpl_toolkits.basemap import Basemap, addcyclic
 from mpl_toolkits.axes_grid1 import AxesGrid
 from netCDF4 import Dataset as NetCDFFile, num2date
 
 
 # today's date is default.
 if len(sys.argv) > 1:
-    YYYYMMDD = sys.argv[1]
+    date = dt.datetime.strptime(sys.argv[1], "%Y%m%d")
 else:
-    YYYYMMDD = datetime.datetime.today().strftime('%Y%m%d')
+    date = dt.datetime.today()
 
 # set OpenDAP server URL.
 try:
-    URLbase="http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs"
-    URL=URLbase+YYYYMMDD+'/gfs_00z'
-    print(URL)
-    data = NetCDFFile(URL)
+    urlbase = "http://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs%Y%m%d/gfs_0p25_00z"
+    data = NetCDFFile(date.strftime(urlbase))
 except:
     msg = """
 opendap server not providing the requested data.
 Try another date by providing YYYYMMDD on command line."""
     raise IOError(msg)
 
-
 # read lats,lons,times.
 
 print(data.variables.keys())

examples/plothighsandlows.py

@@ -4,9 +4,9 @@
 plot H's and L's on a sea-level pressure map
 (uses scipy.ndimage.filters and netcdf4-python)
 """
+import datetime as dt
 import numpy as np
 import matplotlib.pyplot as plt
-from datetime import datetime
 from mpl_toolkits.basemap import Basemap, addcyclic
 from scipy.ndimage.filters import minimum_filter, maximum_filter
 from netCDF4 import Dataset
@@ -22,15 +22,11 @@ def extrema(mat,mode='wrap',window=10):
     return np.nonzero(mat == mn), np.nonzero(mat == mx)
 
 # plot 00 UTC today.
-date = datetime.now().strftime('%Y%m%d')+'00'
+urlbase = "http://nomads.ncep.noaa.gov/dods/gfs_0p25/gfs%Y%m%d/gfs_0p25_00z"
+date = dt.datetime.now()
 
 # open OpenDAP dataset.
-#data=Dataset("http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs/%s/gfs_%sz_anl" %\
-#        (date[0:8],date[8:10]))
-data=Dataset("http://nomads.ncep.noaa.gov:9090/dods/gfs_hd/gfs_hd%s/gfs_hd_%sz"%\
-        (date[0:8],date[8:10]))
-
-
+data = Dataset(date.strftime(urlbase))
 
 # read lats,lons.
 lats = data.variables['lat'][:]