NCL读取micpas4的网格数据读入

王长豫

NCL读取micpas4的网格数据读入，从气象家园看了一些代码修改的，但总是不对，请问该如何读入数据，并格式摆放正确，谢谢。
filename = ".\23110120.000"
  data = asciiread(filename,-1,"string")
  nlat = 1041；用的是micpas4的EC预报数据，但不知道这个个数是多少个
  nlon = 1176
  lon = fspan(30,180,nlon);fspan(开始，结束，个数)，等间距扩展一列数（前面备注都是别人代码里备注的，不知道如何正确填写）
  lat = fspan(0,90,nlat)
  lon2d = onedtond(lon,(/nlat,nlon/));将lon从一维扩展到多维（nlat,nlon）
  lat2d = onedtond(lat,(/nlat,nlon/));
  a = onedtond(data(23::1),(/nlat,nlon/));读数据从第四行读到最后,数据先排一列再排一行（看了micaps4的数据说明文件，但是也不太清楚数据格式摆放，和经纬度网格数）

  a!0="lat"
  a!1="lon"
  a&lat = lat
  a&lon = lon
  a&lat@units="degrees_north"
  a&lon@units="degrees_east"
  a@_FillValue = -9999.