[ncl-talk] NCL code problem

snow snowflake_822 at yahoo.com
Mon Jun 13 23:52:23 MDT 2016


 

  To whom it may concern,
I have written a simple NCL code to extract and take the average of my data and write the table of my data. but in the write part the output is only one data instead of 4.
I wrote this email. to find any kind of help.
I appreciate your time and kindness.
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/20160614/ce2b5d21/attachment.html 


More information about the ncl-talk mailing list