function interpp1(X,Y,V)
    itp = interpolate((X,), Y, Gridded(Linear()))
    sz = size(V)
    y=zeros(sz[1],sz[2])
    for i=1:sz[1]
        for j=1:sz[2]
            y[i,j]=itp(V[i,j])
        end
    end
    return y
end
