求助NCL数组的问题

段段

各位大侠：我有一个txt数据，数据格式是：站号，经纬度，高度，温度，风速等要素，总共有11列，25个时次，每个时次是600个站，我想把这里面的要素按站号和时次单独输出到一个txt文件中，
即：
stations  lat  lon  hgt  t  td  u
52181   
..........
.........
52181
...........
.........
总共25个时次

提取其中单个要素输出成如下格式：
station  lat  lon  1  2  3  4  5  6  ..........25
52181 ............................................



文件按以下方式读出：
begin
;************ read station domain 1 datas ***********
; diri          = "./"
filename    = "wrf-to-stations_d01.txt"
    nrows    = 600
    ncols    = 11
    ntimes    = 25
    station_sta   =(/ntimes,nrows,1/)

    do i= 0 , ntimes-1
       data         = asciiread(filename, (/ntimes,nrows,ncols/), "float")
       data@_FillValue = -99999.00
       station_sta = data(i,:,0)
       lat_sta      = data(i,:,1)
       lon_sta      = data(i,:,2)
       height_sta   = data(i,:,3)
       height_wrf   = data(i,:,4)
       t2m_sta      = data(i,:,5)
       t2m_wrf      = data(i,:,6)
       uspd_sta     = data(i,:,7)
       uspd_wrf     = data(i,:,8)
       udir_sta     = data(i,:,9)
       udir_wrf     = data(i,:,10)
   end do

请问我该如何输出数据呢？

段段

果断转用fortran解决问题了！NCL在处理数组和输出数据，特别是输出ascii文件的时候很不方便！

yetsyu

NCL输出ASCII文件确实不是很方便，但是基本和fortran是一致的，可以参考一下NCL主站，仔细看看还是蛮easy的。

longlivehj

段段 发表于 2015-4-9 20:39
果断转用fortran解决问题了！NCL在处理数组和输出数据，特别是输出ascii文件的时候很不方便！

先恭喜楼主，不管用什么工具，解决问题就行！
同时为ncl鸣冤，数组处理和ascii输出对ncl来说还是很容易的。

niuda

ncl处理数组应该挺方便的，都有许多现成的函数只是速度不如fortran

lideqin

覃建明