[ncl-talk] contour plot with 2d x,y locations ...

[Ufuk Utku Turuncoglu (BE)]

Hi,

I am trying to create a contour plot from ocean float dataset but i 
could not create correct plot. To map data to x,y locations, i tried to 
use "TriangularMesh" resource but in this case gsn_csm_contour does not 
produce anything (just empty plot with weird axis). If i remove 
"TriangularMesh" definition, then NCL creates the plot but it does not 
create the axises and gives following warnings.

warning:PlotManagerSetValues: TickMark annotation cannot be added after 
NhlCreate
warning:tmXBLabelFontHeightF is not a valid resource in 
plot_glpr_2d_contour.PlotManager at this time
warning:tmYLLabelFontHeightF is not a valid resource in 
plot_glpr_2d_contour.PlotManager at this time

So, what is the correct way to create this kind of float plots.

PS: i also attached the data file to mail.

Thanks,
Regards,

--ufuk

;-----------------------------------------------------------
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
   ;--- parameters ---
   ifile = "GL_PR_PF_1901200.nc"

   ;--- open file ---
   nc1 = addfile(ifile, "r")

   ;--- read variables ---
   time = nc1->TIME
   pres = nc1->PRES
   date = cd_calendar(time, 4)
   plat = nc1->LATITUDE
   plon = nc1->LONGITUDE
   data = nc1->TEMP

   dims = dimsizes(pres)
   ntime = plat(0)
   nlev = plat(1)

   ;--- calculate depth from pressure ---
   ;--- UNESCO Technical Papers in Marine Science No. 44 ---
   ;--- http://www.seabird.com/document/an69-conversion-pressure-depth ---
   x = (sin(plat/57.29578))^2
   x2d = conform(pres, x, (/ 0 /))
   g = 9.780318*(1.0+(5.2788e-3+2.36e-5*x2d)*x2d)+1.092e-6*pres
   depth = pres
   depth = 
((((-1.82e-15*pres+2.279e-10)*pres-2.2512e-5)*pres+9.72659)*pres)/g

   ;--- plot data ---
   wks  = gsn_open_wks ("newpdf", "plot_glpr_2d")
   gsn_define_colormap (wks, "matlab_jet")
   i = NhlNewColor(wks,0.8,0.8,0.8)

   res = True
   res at vpHeightF = 0.15
   res at vpWidthF = 0.75
   res at tiYAxisString = "Depth (m)"
   res at tiXAxisString = "Time"
   res at gsnLeftString = ""
   res at gsnRightString = ""
   res at gsnSpreadColors = True
   res at gsnSpreadColorEnd = -3
   res at cnFillOn = True
   res at cnFillMode = "RasterFill"
   res at cnInfoLabelOn = False
   res at cnLinesOn = False
   res at cnLineLabelsOn = False
   res at lbLabelBarOn = False
   res at pmTickMarkDisplayMode = "Always"
   res at tiXAxisFontHeightF = 0.008
   res at tiYAxisFontHeightF = 0.008
   res at tmXBLabelFontHeightF = 0.008
   res at tmYLLabelFontHeightF = 0.008

   xa = conform(pres, date, (/ 0 /))
   copy_VarCoords(pres, xa)
   xa = mask(xa, ismissing(depth), False)
   res at sfXArray = xa(DEPTH|:,TIME|:)
   res at sfYArray = depth(DEPTH|:,TIME|:)
   res at trGridType = "TriangularMesh"
   res at trYReverse = True

   plot1 = gsn_csm_contour(wks, data(DEPTH|:,TIME|:), res)
end
-------------- next part --------------
A non-text attachment was scrubbed...
Name: GL_TS_PF_1901200.nc
Type: application/octet-stream
Size: 724436 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141117/783e0e74/attachment-0001.obj