funtion [val]=pintrp(fname,lon,lat,z,t)

