关于格点数据存成ascii文件的问题

gercibee

我的脚本是把b值都存起来了，正常的写ascii文件是需要把经纬度也存里面吗？这个需要怎么写进去呢？

begin
    fils = systemfunc ("ls /mnt/lustre/hanzq/ss/et/et*.h5")
    f = addfiles(fils, "r")

    do i = 0 ,dimsizes(fils)-1

        a = new((/12,360,720/), "float")
        a(0,:,:) = f ->$("ET_monthly_Jan_"+(i+2005))$
        a(1,:,:) = f ->$("ET_monthly_Feb_"+(i+2005))$
        a(2,:,:) = f ->$("ET_monthly_Mar_"+(i+2005))$
        a(3,:,:) = f ->$("ET_monthly_Apr_"+(i+2005))$
        a(4,:,:) = f ->$("ET_monthly_May_"+(i+2005))$
        a(5,:,:) = f ->$("ET_monthly_Jun_"+(i+2005))$
        a(6,:,:) = f ->$("ET_monthly_Jul_"+(i+2005))$
        a(7,:,:) = f ->$("ET_monthly_Aug_"+(i+2005))$
        a(8,:,:) = f ->$("ET_monthly_Sep_"+(i+2005))$
        a(9,:,:) = f ->$("ET_monthly_Oct_"+(i+2005))$
        a(10,:,:) = f ->$("ET_monthly_Nov_"+(i+2005))$
        a(11,:,:) = f ->$("ET_monthly_Dec_"+(i+2005))$

        b = dim_sum_n_Wrap(a(:,:,:),0)

；属性信息
        lat = fspan(90,-89.5,360)
        lon = fspan(-180,179.5,720)
        lat@long_name = "latitude"
        lat@units  = "degrees_north"
        lat!0      = "lat"
        lat&lat    =  lat
        lon@long_name = "longitude"
        lon@units  = "degrees_east"
        lon!0      = "lon"
        lon&lon    =  lon
        b!0 = "lat"
        b!1 = "lon"
        b&lat    =  lat                    ; create coordinate variable
        b&lon    =  lon
        b@_FillValue = -119988

        opt = True
        opt@fout = "./"+(i+2005)+".txt"

        write_matrix(b, "f10.1", opt)
        delete(a)
        delete(b)
    end do
end

longlivehj

http://www.ncl.ucar.edu/Document ... in/asciiwrite.shtml
example 4

gercibee

longlivehj 发表于 2015-3-12 10:30
http://www.ncl.ucar.edu/Document ... in/asciiwrite.shtml
example 4

这个几天不稳定，刚上来，多谢版主提醒，官网例子有一个小错误，修改下已经存完