WY=2019;
matdates=datenum([WY-1 10 1]):datenum([WY 9 30]);

fill_and_run_modis(tiles,r0dates,matdates,...
    hdfbasedir,topofile,mask,Ffile,shade,grain_thresh,dust_thresh,tolval,...
    outloc,nameprefix);

smoothSPIREScube(nameprefix,outloc,matdates,...
    windowSize,windowThresh,mingrainradius,maxgrainradius,mindust,maxdust,...
    mask,topofile,el_cutoff,fsca_thresh,cc,fice,b_R,dust_rg_thresh);
