# 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)