Hi
I am trying to overlay and increment variable with colour filled contour
and the background variable with line contour using
gsn_csm_contour_map_overlay.
However the cnLevelSpacingF for the line contour has a warning saying
that its not valid at time of plotting.
I did set my resource for the line contour, hence not exactly sure how
to rectify this.
Attached below is the code segment for the plotting
Thanks
Agnes
====================================================
wks = gsn_open_wks("x11"
,"2008-05-22_00_00_GHT_anal_incre_"+i) ; open a ps file
res = True ; plot mods desired
gsn_define_colormap(wks,"BlWhRe") ; select color map
res@gsnSpreadColors = True
res@mpLandFillColor = -1 ; set land to be transperant
res@mpOutlineOn = True
res@mpOutlineBoundarySets = "AllBoundaries"
res@mpProjection = "LambertConformal"
res@mpEllipticalBoundary = False ; If you want
elliptical
res@cnFillOn = True
res@cnLinesOn = False
res@lbBoxLinesOn = True
res@mpMinLatF = 7 ; range to zoom in on
res@mpMaxLatF = 56
res@mpMinLonF = -170
res@mpMaxLonF = -53
res@tmXTOn = False ; Turn off top tickmarks and labels
res@tmYROn = False ; Turn off left tickmarks and labels
res@tmXBLabelFontHeightF = 0.01 ; resize tick labels
res@tmYLLabelFontHeightF = 0.01
res@pmTickMarkDisplayMode = "Always"
res@lbLabelAutoStride = True
res@cnLevelSelectionMode = "AutomaticLevels" ; set manual
contour levels
res@mpLambertParallel1F = wrf_bkg_stdQC@TRUELAT1
res@mpLambertParallel2F = wrf_bkg_stdQC@TRUELAT2
res@mpLambertMeridianF = wrf_bkg_stdQC@CEN_LON ; or f@CEN_LON
res@gsnMaskLambertConformal = True
res@lbLabelBarOn = False
res@cnLevelSelectionMode = "ManualLevels"
res@cnMinLevelValF = -60
res@cnMaxLevelValF = 60
res@cnLevelSpacingF = 5
res@gsnAddCyclic = False
res@gsnDraw=False
res@gsnFrame=False
res1 = True
res1@cnFillOn = False
res1@cnLinesOn = True
res1@cnLevelSelectionMode = "ManualLevels"
res1@cnMinLevelValF = 5000
res1@cnMaxLevelValF = 6000
res1@cnLevelSpacingiF = 100
plot(0)=gsn_csm_contour_map_overlay(wks, H_increment_stdQC(i,:,:),
H_bkg_stdQC(i,:,:), res, res1)
plot(1)=gsn_csm_contour_map_overlay(wks, H_increment_extQC(i,:,:),
H_bkg_extQC(i,:,:), res, res1)
plot(2)=gsn_csm_contour_map_overlay(wks, H_increment_extQC_1(i,:,:),
H_bkg_extQC_1(i,:,:), res, res1)
resP=True
resP@gsnPanelFigureStings=titles
resP@amJust="TopLeft"
resPgsnPanelFigureStngsPerimOn=False
resP@txString ="Analysis Increment of Geopotential Height @ Model
Level="+i+" for 2008-05-22 00:00:00"
resP@gsnPanelLabelBar = True
resP@lbLabelFontHeightF=0.015
resP@lbLabelStride = 2
gsn_panel(wks, plot, (/3,1/), resP)
gsnMaximize = True
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon May 16 14:57:43 2011
This archive was generated by hypermail 2.1.8 : Tue May 17 2011 - 15:11:27 MDT