one dimensional data

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