Another xy plot qustion

From: ugo merlini <ugomerlini_at_nyahnyahspammersnyahnyah>
Date: Sat Apr 23 2011 - 11:33:24 MDT


I'm able to read from opendap server where a I grab GFS model data. In particular I want to make an xy plot with 2m temperature on Y axis and time on X axis

here my script

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"


  url = ""
  filename = url + "gfs_master_00z"
  exists = isfilepresent(filename)
  if(.not.exists) then
    print("OPeNDAP isfilepresent test unsuccessful.")
    print("Either file doesn't exist, or NCL does not have OPeNDAP capabilities on this system")
    print("OPeNDAP isfilepresent test successful.")
    a = addfile(filename,"r")
    vars = getfilevarnames(a)
  end if

temp = a->tmp2m(:,48,9) ; temperature
temp = temp - 273.15
temp@units = "(C)"

;--- Open workstation
  wks_type = "png"
  wks = gsn_open_wks(wks_type,"2m Temperature")
;---Set some resources for the map.
  res = True
  res@cnMonoLineColor = True ; Turn off the drawing of contours lines in one color.
  res@tiMainString = "Temperature (C)" ; Create a title.
field1 = create "field1" coordArraysClass defaultapp
  "caYArray": (temp)
  "caXArray": (temp!0)
end create
; Create XyPlot object and assign data to it
xy_plot = create "xy_plot" xyPlotClass wks
    "xyCoordData": field1
end create

tmp2m data is a three dimensional array (time,lat,lon)

print(temp) give me the value so the file is readed correctly

I don't undestand how assign the Y axis and the X

I'm doing a 2 dimension plot?



ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
Received on Sat Apr 23 11:34:41 2011

This archive was generated by hypermail 2.1.8 : Tue May 03 2011 - 14:47:35 MDT