use error.nc shade error[k=1,d=1] shade/lev=(0,1,0.01) error