Dear NCL users,
I am trying to plot a 1-dimensional variable: the results from a channel
model that has a spatial dimension of 1x50. I understand from other
emails in the NCL-Talk archive that the resources sfXArray and sfYArray
are used when plotting 1-dimensional data, but I have been unsuccessful
in producing a plot with my data. Thanks in advance for the help. My
script is attached.
-Dan
;================================================;
; hovmoller.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/contributed.ncl"
; ================================================;
begin
;=================================================;
; open file and read in data
;=================================================;
in = readAsciiTable("I.txt",52,"float",0)
I = in(:,1:50)
D = floattoint(readAsciiTable("D.txt",1,"float",0)/365)
;D = readAsciiTable("D.txt",1,"float",0)
;D = in2(:,1:50)
LAT = readAsciiTable("lat.txt",50,"float",0)
dimension = dimsizes(I)
I!0 = "time"
I!1 = "Latitude"
I&time = D(:,0)
I&Latitude = LAT(0,:)
printVarSummary(I)
;=================================================;
; PLOT
;=================================================;
yy = new((/50/),float)
do j = 0,49
yy(j) = 1
end do
print(yy)
do i = 0,20000,10
name = sprinti("%0.4i",i)
times = D(i,0)/365
wks = gsn_open_wks("ps",""+name+"") ; open a ncgm file
gsn_define_colormap(wks,"rainbow") ; choose a colormap
colors = gsn_retrieve_colormap(wks)
do i = 2,3,1
colors(i,:) = (/1,1,1/)
end do
;gsn_define_colormap(wks,colors)
res = True ; plot mods desired
;res@gsnDraw = False
;res@gsnFrame = False
res@cnFillOn = True ; turn on color
res@gsnSpreadColors = True ; use full range of colormap
res@gsnMaximize = True
;res@cnFillMode = "CellFill"
res@cnLinesOn = False
res@tiMainFont = 22
res@tmXBLabelFont = 22
res@tmYLLabelFont = 22
res@tiXAxisFont = 22
res@tiYAxisFont = 22
res@lbLabelFont = 22
res@gsnStringFont = 22
res@lbLabelStride = 2
res@lbLabelFontHeightF = 0.015
res@lbOrientation = "Vertical"
res@tiMainString = "Temperature (degrees celsius)"
res@tiMainJust = "BottomLeft"
res@tiMainPosition = "Left"
res@gsnLeftString = ""+times+" days"
res@gsnRightString = "Black line = ice edge"
res@tiYAxisString = "Time (years)"
res@tiXAxisString = "Latitude"
res@cnLevelSelectionMode = "ManualLevels"
res@cnMinLevelValF = -20
res@cnMaxLevelValF = 30
res@cnLevelSpacingF = 5
res@lbLabelStride = 2
res@gsnXAxisIrregular2Linear = True
res@sfXArray = LAT(0,:)
res@sfYArray = yy
plot1 = gsn_csm_contour(wks,I(i,:),res) ; contour the variable
;res1 = True
;res1@gsnDraw = False
;res1@gsnFrame = False
;res1@cnLinesOn = True
;res1@cnFillOn = False
;res1@cnLevelSelectionMode = "ExplicitLevels"
;res1@cnLevels = (/-10/)
;res1@cnInfoLabelOn = False
;res1@cnLineLabelsOn = False
;res1@cnLineThicknessF = 6.
;plot2 = gsn_csm_contour(wks,I(i,:),res1)
;overlay(plot1,plot2)
;draw(plot1)
;frame(wks)
end do
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Dec 16 13:26:07 2009
This archive was generated by hypermail 2.1.8 : Thu Dec 17 2009 - 17:15:52 MST