R/fill_grd.R
    fill_grd.RdAdd date, year, month, and day columns to the interpolation grids using dates from dat_in. The interp.surface function is used after splitting the interpolation matrix by month to fill missing values.  Function is used in wrtds.
fill_grd(grd_in, dat_in, interp = FALSE)