[ncl-talk] Write table problem

Wei Huang whuang at univ-wea.com
Thu Jun 16 07:33:15 MDT 2016


Snow,

Please check the docs at: https://www.ncl.ucar.edu/Document/Functions/Built-in/write_table.shtml

You need to change line:

write_table("list.txt", "a", alist, "%5.5f")

To

write_table("list.txt", "a", alist, "%10.5f%10.5f%10.5f%10.5f")

You may add a separator, such as “,” as:

write_table("list.txt", "a", alist, "%10.5f,%10.5f,%10.5f,%10.5f")

Try this script:

     R5=(/189,193,191,192,192,193,192,188,198,196,192,187,197,194,190,187,197,196,195,195,196/)
     R6=(/285,286,287,288,289,290,291,292,293,294,295,296,296,297,298,299,298,299,300,301,302/)
     R7=(/180,174,172,181,179,173,184,175,185,174,184,173,183,193,182,191,179,189,178,188,176/)
     R8=(/264,266,267,267,268,269,270,270,271,271,272,273,273,273,274,273,275,274,276,275,277/)

     ca=tofloat(R5)
     cb=tofloat(R6)
     cc=tofloat(R7)
     cd=tofloat(R8)

    alist=[/ca,cb,cc,cd/]
write_table("list.txt", "w", alist, "%10.5f,%10.5f,%10.5f,%10.5f")


Regards,


Wei Huang

From: ncl-talk-bounces at ucar.edu [mailto:ncl-talk-bounces at ucar.edu] On Behalf Of snow
Sent: Wednesday, June 15, 2016 8:05 PM
To: Ncl Talk
Subject: [ncl-talk] Write table problem


To whom it may concern,

My below ncl code has aproblem in writting the table with 4 data in one column. It only writes the first data and doesnt continue to write the other data.

does anyone know how to fix this problem?

I appreciate your time and consideration.

Best regards,

Snow Sho


this is the code:
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
;----------------------------------------------------------------------
; WRF data translate NetCDF file to  ascii file
;----------------------------------------------------------------------

    infile=("wrfout_d03_2007-07-24_00_00_00")
    fn=dimsizes(infile)
    fin=addfile(infile,"r")
    R=fin->RAINNC(0,:,:)
     a=(/189,193,191,192,192,193,192,188,198,196,192,187,197,194,190,187,197,196,195,195,196/)
     b=(/285,286,287,288,289,290,291,292,293,294,295,296,296,297,298,299,298,299,300,301,302/)
     c=(/180,174,172,181,179,173,184,175,185,174,184,173,183,193,182,191,179,189,178,188,176/)
     d=(/264,266,267,267,268,269,270,270,271,271,272,273,273,273,274,273,275,274,276,275,277/)
     e=(/202,205,206,206,198,197,206,205,204,202,200,200,200,200,210,209,205,202,212,208,204,203/)
     f=(/270,271,272,273,274,276,277,278,279,280,281,282,283,284,284,284,285,286,287,288,289,290/)
     g=(/208,208,209,210,211,210,207,205,215,212,213,215,216,207,209,210,211,214,215,209,211,215,214/)
     h=(/263,264,265,266,267,267,268,269,269,270,271,273,274,276,277,278,279,279,280,282,283,284,287/)
        R1=new((/21/),"float")
        R2=new((/21/),"float")
        R3=new((/22/),"float")
        R4=new((/23/),"float")

        do k=0,20
        R1(k)=R(a(k),b(k))
        end do
     ;  print(R1)
        do k=0,20
        R2(k)=R(c(k),d(k))
        end do
        do k=0,21
        R3(k)=R(e(k),f(k))
        end do
        do k=0,22
        R4(k)=R(g(k),h(k))
        end do
        R5=avg(R1(:))
        R6=avg(R2(:))
        R7=avg(R3(:))
        R8=avg(R4(:))

        ca=new((/1/),"float")
        ca=(/R5/)
        cb=new((/1/),"float")
        cb=(/R6/)
        cc=new((/1/),"float")
        cc=(/R7/)
        cd=new((/1/),"float")
        cd=(/R8/)

;   alist=[/R5,R6,R7,R8/]
   alist=[/ca,cb,cc,cd/]
  write_table("list.txt", "a", alist, "%5.5f")
  end
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160616/4f8f38ec/attachment.html 


More information about the ncl-talk mailing list