把NCL算出来的值写入txt文件

黄小仙儿

算出来了各个时间段的最大值，想把对应的时间和最大值写入另一个文件，比如txt，然后用origin画时间序列图。    if(isfilevar(a,"QICE"))
      qi = wrf_user_getvar(a,"QICE",it)
      qi = qi*1000.
      qi@units = "g/kg"   
    end if
    max = max(qi)

但是官网上好像只有写入ascii和nc文件的例子。
asciiwrite(ascii_file,temp(7:3:2,0:4))
请问大家有没有把数据写入过txt之类的，我们可以打开并且可以看到数据的文件
请指点一二

longlivehj

ascii是一种编码类型，txt仅仅是一个文件后缀名，两者不是同一回事。txt可以有多种编码，ascii也是其中之一。所以，ncl官网上有关写ascii文件的例子都能满足你的要求。

longlivehj

黄小仙儿，你牛人啊，grads、ncl、meteoinfo、origin、python、fortran你同时来啊！o(∩∩)o...哈哈。

askira2011

asciiwrite("*.txt",data)不就行了，另外你可以试试ascii_table ,ascii_matrix

黄小仙儿

longlivehj 发表于 2014-4-23 13:25
黄小仙儿，你牛人啊，grads、ncl、meteoinfo、origin、python、fortran你同时来啊！o(∩∩)o...哈哈。

因为太菜了，什么都做不好啊~

黄小仙儿

askira2011 发表于 2014-4-23 13:46
asciiwrite("*.txt",data)不就行了，另外你可以试试ascii_table ,ascii_matrix

我试着改了下，但是好像循环有问题，出来的txt只有一个数值
do it =0, ntimes-1,1        ; TIME LOOP

    print("Working on time: " + times(it) )
    res@TimeLabel = times(it)   ; Set Valid time to use on plots

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; First get the variables we will need        
    if(isfilevar(a,"QICE"))
      qi = wrf_user_getvar(a,"QICE",it)
      qi = qi*1000.
      qi@units = "g/kg"   
    end if
    wmax = max(qi)
        print(wmax)
ascii_file = "data.txt"
asciiwrite(ascii_file, sprintf("%12.8f", wmax ))
end do

黄小仙儿

longlivehj 发表于 2014-4-23 13:25
黄小仙儿，你牛人啊，grads、ncl、meteoinfo、origin、python、fortran你同时来啊！o(∩∩)o...哈哈。

我试着改了下，但是好像循环有问题，出来的txt只有一个数值
do it =0, ntimes-1,1        ; TIME LOOP

    print("Working on time: " + times(it) )
    res@TimeLabel = times(it)   ; Set Valid time to use on plots

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; First get the variables we will need        
    if(isfilevar(a,"QICE"))
      qi = wrf_user_getvar(a,"QICE",it)
      qi = qi*1000.
      qi@units = "g/kg"   
    end if
    wmax = max(qi)
        print(wmax)
ascii_file = "data.txt"
asciiwrite(ascii_file, sprintf("%12.8f", wmax ))
end do

longlivehj

黄小仙儿 发表于 2014-4-23 16:06
我试着改了下，但是好像循环有问题，出来的txt只有一个数值
do it =0, ntimes-1,1        ; TIME LOOP
...

wmax = max(qi)就一个值么，每次都写这个值到同一个文件。
要输出什么？是要把每个时次QICE的最大值输出到一个文件中吗？如果是，代码应该如下：

ascii_file = "data.txt"
str = new(ntimes, "string")

do it =0, ntimes-1,1        ; TIME LOOP

    print("Working on time: " + times(it) )
    res@TimeLabel = times(it)   ; Set Valid time to use on plots

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; First get the variables we will need        
    if(isfilevar(a,"QICE"))
      qi = wrf_user_getvar(a,"QICE",it)
      qi = qi*1000.
      qi@units = "g/kg"   
    end if
    wmax = max(qi)
    str(it) = sprintf("%12.8f", wmax)
end do

asciiwrite(ascii_file, str)

黄小仙儿

longlivehj 发表于 2014-4-23 16:18
wmax = max(qi)就一个值么，每次都写这个值到同一个文件。
要输出什么？是要把每个时次QICE的最大值输出 ...

哇，谢谢大神~

longlivehj

黄小仙儿 发表于 2014-4-23 18:33
哇，谢谢大神~

不客气。
想跟你说，其实可以用wrf_user_getvar一次性把数据全部取出，然后用dim_max_n直接生成最大值系列，避免循环。