NCL 如何读取站点数据的问题

xhr2xx

刚接触ncl，接到任务，拿到105年美国2823个郡的粮食产量的数据，有站点的经纬度，以及value，全都是csv格式，由于发现每一年的有值的conty 不一样多，所以整理成一共出现过的2823个conty，105年的数据，站点没有值的年份用-9999作为缺测值，之后我不知道是否能用ncl读取csv，所以我用fortran转成了grd，一个grd存放经纬度，一个存放value,想这样读取到ncl 中生成nc文件绘制等值线图，但是遇到一个问题，我在读取时，不知道怎么样才能对应一个经纬度一个值，创建一个（105,2823,2823）的数组的时候，系统报错超出系统上限字节，所以不知道这么办 求各位大大教教如何读取站点，或如何插值？农作物的数据用什么插值法比较好?

xhr2xx

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"  
; ================================================;
begin
  

  year=ispan(1910,2014,1)
  nyear = dimsizes(year)
  lat=new(2823,float)
  lon=new(2823,float)
  do i = 1,2823,1  
   lat(i-1) = fbindirread("C:\Users\xhr2xx\Desktop\USDA data\corn\corn_grain_acres_harvested\jwd.grd",2*i-2,1,"float")
   lon(i-1) = fbindirread("C:\Users\xhr2xx\Desktop\USDA data\corn\corn_grain_acres_harvested\jwd.grd",2*i-1,1,"float")
   
  end do

  lat@units = "degrees_north"
  lon@units = "degrees_east"

  var = new((/nyear,2823,2823/),"float",-9999)
  var!0   = "year"
  var!1   = "lat"
  var!2   = "lon"
  var&year = ispan(1,365,1)
  var&lat = lat
  var&lon = lon

var = 0.0        
do i = 1,28823,1  
   do iy=1,105,1
   var(iy-1,lat(i-1),lon(i-1)) = fbindirread("C:\Users\xhr2xx\Desktop\USDA data\corn\corn_grain_acres_harvested\zzsj.grd",105*(i-1)+iy-1,1,"float")
   var@_FillValue = -9999
   end do
end do

xhr2xx

附我的脚本，还没写完 但是这样就已经无法继续了

xhr2xx

只知道谁能帮忙一下

尽头的尽头

ncl可以直接读取csv格式的数据，至于插值官网有教程，论坛也有教程