RE: [ncl-talk] add plots in ncl

Hello, Debasish,

Your information is very helpful, and the error is gone now:)
However, I got only one line instead of 3 as expected in my new plot. Here is the simple case I am trying:
-------------------------------------
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
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/shea_util.ncl"

begin
 path = "/Volumes/Data/Models/cam/AOD/aeronet/"
 fileHandle5= path + "test.txt"
 DUSTAOD5 = asciiread(fileHandle5, (/4,4/), "float")
 DUSTAOD5@_FillValue = -9999

 aod5 = new((/3,4/),"float", DUSTAOD5@_FillValue)
 day5 = new((/1,4/),"float","No_FillValue")
 day5(0,:) = DUSTAOD5(:,0)
 aod5(0,:) = DUSTAOD5(:,1)
 aod5(1,:) = DUSTAOD5(:,2)
 aod5(2,:) = DUSTAOD5(:,3)
 
 delete(path)

 print("plotting")
 wks = gsn_open_wks("pdf", "test")
 plot = new(1,graphic)

 res1 = True
 res1_at_gsnDraw = False ;don't draw
 res1_at_gsnFrame = False
 res1_at_xyMarkLineModes = "Markers" ;1st line
 res1_at_xyMarkers = 3
 res1_at_xyMarkerSizeF = 0.02
 res1_at_xyLineThicknesses = (/2.0,2.0,2.0/)
 res1_at_xyLineColors = (/"black","blue","red"/)
 res1_at_pmLegendDisplayMode = "Always" ; turn on legend
 res1_at_pmLegendSide = "Top" ; Change location of
 res1_at_pmLegendParallelPosF = .6 ; move units right
 res1_at_pmLegendOrthogonalPosF = -0.3 ; move units down
 res1_at_pmLegendWidthF = 0.15 ; Change width and
 res1_at_pmLegendHeightF = 0.09 ; height of legend.
 res1_at_lgPerimOn = False ; turn off box around
 res1_at_lgTitleFontHeightF = .018
 res1_at_lgTitleString = "Beijing(39.98N 116.38E)" ; create legend title
 res1_at_lgLabelFontHeightF = .015 ; label font height
 res1_at_xyExplicitLegendLabels = (/"AERONET","10-m wind","ustar"/)
 res1_at_tiXAxisString = "days in April 2001"
 res1_at_tiYAxisString = "AOD[500nm] "
 res1_at_tmLabelAutoStride = True
 res1_at_trXLog = False
 res1_at_trYLog = False
 res1_at_trXMinF = 0.0
 res1_at_trXMaxF = 30
 res1_at_trYMinF = 0.0
 res1_at_trYMaxF = 6.0

 plot(0) = gsn_csm_xy(wks, day5, aod5, res1)

end
---------------------------

And here is the data (test.txt with 1st column Julian days):

------------------------------
1.1 0.3 0.1 0.2
1.2 0.4 -9999 -9999
1.3 0.5 -9999 -9999
2.1 0,6 0.2 0.3
----------------------------------

I appreciate any input to correct my wrong scripts!

Thank you,
-Lin

---- Original message ----
>Date: Thu, 2 Aug 2007 13:33:26 -0800
>From: "Debasish Pai Mazumder" <debasish_at_gi.alaska.edu>
>Subject: RE: [ncl-talk] add plots in ncl
>To: "'Lin Su'" <Lin.Su_at_Colorado.EDU>
>
>Dear Lin,
>
>I guess that you have missing values in your some columns because of that
>you get dimension mismatch error.
>
>One way I can fill with -9999 for skipping data or
>
>Check the last example in
>
> http://www.ncl.ucar.edu/Applications/ascii_multicol.shtml
>
>and specify the numder of rows where you have values and rest of them assume
>-9999 i.e. missing values.
>
>Good luck
>
>Debasish
>
>
>-----Original Message-----
>From: Lin Su [mailto:Lin.Su_at_Colorado.EDU]
>Sent: Thursday, August 02, 2007 12:56 PM
>To: Debasish Pai Mazumder
>Subject: RE: [ncl-talk] add plots in ncl
>
>BTW, I have more than 400 lines data for each colummn, but sounds NCL only
>read 76 data of them? -Lin
>
>---- Original message ----
>>Date: Thu, 2 Aug 2007 12:44:44 -0800
>>From: "Debasish Pai Mazumder" <debasish_at_gi.alaska.edu>
>>Subject: RE: [ncl-talk] add plots in ncl
>>To: "'Lin Su'" <Lin.Su_at_Colorado.EDU>, <ncl-talk_at_ucar.edu>
>>
>>Dear Lin,
>>
>>You can declare
>>
>>data=new((/3,4/),float) ; dimension= column, time
>>
>>data(0,:)= column2(:)
>>data(1,:)= column3(:)
>>data(2,:)= column4(:)
>>
>>Then in plotting statement:
>>
>>plot= gsn_csm_xy (wks,day,data(:,:),res) ; create plot
>>
>>check the example2 in
>>
>>http://www.ncl.ucar.edu/Applications/xy.shtml
>>
>>Good luck
>>
>>Debasish
>>
>>-----Original Message-----
>>From: Lin Su [mailto:Lin.Su_at_Colorado.EDU]
>>Sent: Thursday, August 02, 2007 12:19 PM
>>To: Debasish Pai Mazumder; 'Lin Su'; ncl-talk_at_ucar.edu
>>Subject: RE: [ncl-talk] add plots in ncl
>>
>>Hello, Debasish,
>>
>>I appreciate your kindly reply:)
>>My case is as follows, and sorry that I didn't give more clearly
>expressions
>>in my last email:
>>
>>Now I have 4 columns data which looks like this:
>>----------------------------
>>1 0.1 0.2 0.1
>>2 0.2
>>3 0.2
>>4 0.4 0.3 0.5
>>---------------------------
>>The first column is the time (day), the second column is the data with the
>>whole column. However, the other two columns with only half data of column
>>2. How can I plot these 3 lines of time series in one plot? I tried to
>>seperate these data into 3 files and used gsn_csm_xy3 to plot it. However,
>I
>>got 3 y axes which I didn't expect.
>>
>>Thanks again,
>>-Lin
>>
>>
>>---- Original message ----
>>>Date: Thu, 2 Aug 2007 11:25:47 -0800
>>>From: "Debasish Pai Mazumder" <debasish_at_gi.alaska.edu>
>>>Subject: RE: [ncl-talk] add plots in ncl
>>>To: "'Lin Su'" <Lin.Su_at_Colorado.EDU>, <ncl-talk_at_ucar.edu>
>>>
>>>Dear Lin,
>>>
>>>You can check the "overlay" in
>>>
>>>http://www.ncl.ucar.edu/Document/Functions/Built-in/overlay.shtml
>>>
>>>or the example2 in
>>>
>>>http://www.ncl.ucar.edu/Applications/xy.shtml
>>>
>>>This example produces multiple line plots with varying colors,
>thicknesses,
>>>and markers.
>>>
>>>Thanks
>>>
>>>Debasish
>>>
>>>
>>>-----Original Message-----
>>>From: ncl-talk-bounces_at_ucar.edu [mailto:ncl-talk-bounces_at_ucar.edu] On
>>Behalf
>>>Of Lin Su
>>>Sent: Thursday, August 02, 2007 11:10 AM
>>>To: ncl-talk_at_ucar.edu
>>>Subject: add plots in ncl
>>>
>>>Hello, all,
>>>
>>>Could anybody tell me how to add another two lines in the xy plot? I tried
>>>gsn_polyline and gsn_csm_xy3, and which are not what I need. Actually, I
>am
>>>looking for some function which likes IDL oplot.
>>>
>>>Thanks a lot,
>>>-Lin
>>>_______________________________________________
>>>ncl-talk mailing list
>>>ncl-talk_at_ucar.edu
>>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>>
>
>
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Aug 03 2007 - 11:51:27 MDT