use FRECVATM_model2_01.nc
shade FRECVATM[k=1,d=1]