# Quick program to fetch NWM raw output. # Logan Karsten # National Center for Atmospheric Research # Research Applications Laboratory import datetime import subprocess hoursBack = 36 hoursLag = 6 bDate = datetime.datetime.now() for step in range(hoursBack,hoursLag,-1): dCurrent = bDate - datetime.timedelta(hours=step) aacPath = "http://www.ftp.ncep.noaa.gov/data/nccf/com/nwm/prod/nwm." + \ dCurrent.strftime('%Y%m%d') + "/analysis_assim/nwm.t" + \ dCurrent.strftime('%H') + 'z.analysis_assim.channel_rt.tm00.conus.nc.gz' cmd = 'wget ' + aacPath subprocess.call(cmd,shell=True) aacPath = "http://www.ftp.ncep.noaa.gov/data/nccf/com/nwm/prod/nwm." + \ dCurrent.strftime('%Y%m%d') + "/analysis_assim/nwm.t" + \ dCurrent.strftime('%H') + 'z.analysis_assim.land.tm00.conus.nc.gz' cmd = 'wget ' + aacPath subprocess.call(cmd,shell=True) aacPath = "http://www.ftp.ncep.noaa.gov/data/nccf/com/nwm/prod/nwm." + \ dCurrent.strftime('%Y%m%d') + "/analysis_assim/nwm.t" + \ dCurrent.strftime('%H') + 'z.analysis_assim.terrain_rt.tm00.conus.nc.gz' cmd = 'wget ' + aacPath subprocess.call(cmd,shell=True)