use FRECVOCN_model1_03.nc
shade FRECVOCN[k=1,d=1]