[ncl-install] Question
Najib Yusuf
najibgal at yahoo.com
Tue May 12 09:09:08 MDT 2015
Hello,I am a new ncl learner, and i am trying to make a plot from an ascii data file but i am getting error messages, kindly guide me, this is how it looks like;
najib at UbuntuVaio:~$ ncl
Copyright (C) 1995-2012 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.1.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
ncl 0> dir = "/home/najib/Desktop/"
ncl 1> filen = "TRODAN.dat"
ncl 2> z1 = asciiread(dir+filen,(/3326/),"string")
ncl 3> printVarSummary(z1)
Variable: z1
Type: string
Total Size: 26608 bytes
3326 values
Number of Dimensions: 1
Dimensions and sizes: [3326]
Coordinates:
Number Of Attributes: 1
_FillValue : missing
ncl 4> new_str = str_split_csv(z1, ",", 1)
ncl 5> printVarSummary(new_str)
Variable: new_str
Type: string
Total Size: 266080 bytes
33260 values
Number of Dimensions: 2
Dimensions and sizes: [3326] x [10]
Coordinates:
Number Of Attributes: 1
_FillValue : missing
ncl 6> _FillValue = 999
ncl 7> variables = new_str(0,0:9)
ncl 8> print(variables)
Variable: variables
Type: string
Total Size: 80 bytes
10 values
Number of Dimensions: 1
Dimensions and sizes: [10]
Coordinates:
Number Of Attributes: 1
_FillValue : missing
(0) Year
(1) Month
(2) Day
(3) Hour
(4) Minute
(5) Second
(6) UT_sec
(7) Ozone_ppbv
(8) Temperature_K
(9) Pressure_Torr
ncl 9> year = stringtoint(new_str(1:3325,0))
ncl 10> month = stringtoint(new_str(1:3325,1))
ncl 11> day = stringtoint(new_str(1:3325,2))
ncl 12> hour = stringtoint(new_str(1:3325,3))
ncl 13> minute = stringtoint(new_str(1:3325,4))
ncl 14> second = stringtoint(new_str(1:3325,5))
ncl 15> UT_sec = stringtofloat(new_str(1:3325,6))
ncl 16> ozone_ppbv = stringtofloat(new_str(1:3325,7))
ncl 17> Temp_K = stringtofloat(new_str(1:3325,8))
ncl 18> Press_Pa = stringtofloat(new_str(1:3325,9))
ncl 19> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
ncl 20> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
ncl 21> begin
ncl 22> f = addfile ("$NCARG_ROOT/lib/ncarg/data/ascii/TRODAN.dat.asc","r")
ncl 23> u = f->U
ncl 24> wks = gsn_open_wks ("ps","xy")
ncl 25> res = True
ncl 26> rest at tiMainString = "Basic XY plot"
ncl 27> plot = gsn_csm_xy (wks,u&lat,u(1,{3325}),res)
ncl 28> end
fatal:["FileSupport.c":2758]:_NclFindFileExt: Requested file <$NCARG_ROOT/lib/ncarg/data/ascii/TRODAN.dat.asc> or <$NCARG_ROOT/lib/ncarg/data/ascii/TRODAN.dat> does not exist
fatal:["FileSupport.c":3052]:($NCARG_ROOT/lib/ncarg/data/ascii/TRODAN.dat.asc) has no file extension, can't determine type of file to open
fatal:file (f) isn't defined
fatal:["Execute.c":7741]:Execute: Error occurred at or near line 23
Thank you, Najib
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-install/attachments/20150512/9eeed8f5/attachment-0001.html
More information about the ncl-install
mailing list