use FSENDOCN_model1_02.nc
shade FSENDOCN[k=3,d=1]