
sys.exit()

# soil moisture data (10-minute)
obs   = nc.Dataset('observations/cesar_soil_water_lb1_t10_v1.1_200607.nc')
sm_03 = obs.variables['TH03'][tidx]
sm_08 = obs.variables['TH08'][tidx]
sm_20 = obs.variables['TH20'][tidx]

# soil moisture data (daily average from prior day)
# Julian day for July 2, 2006 is 183, so day prior is 182
jdi   = 181