How do I get the top x-axis label and the right y-axis label to disappear?
I don't see anything in my script that put them there in the first place.
;***************************************************************************
; narr_ncl_template.ncl
;
; Read NARR NCDC-NOMADS data
;
; Brittany Hailey 04/04/2013
; - Plotting NARR data
;
; Executea: ncl narr_ncl_template.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" ; plot.
Include before
begin
;*******************************************
; Open file and read in data
;*******************************************
dir = "/narrr-data/data/" ; input
directory
f = addfile (dir + "narr-a_221_20011231_0000_000.grb", "r")
;*******************************************
; Read variables and assign names
;*******************************************
; These are 2-D arrays
lat = f->gridlat_221
lon = f->gridlon_221
lev = f->lv_ISBL2 ; [0:19]
[150,200,300,500,550,600,650,700,725,750,775,800,825,850,875,900,925,950,975,1000]
nlevels =dimsizes(lev)
; print(nlevels)
; printVarSummary(lev)
; print(lev)
;*******************************************
; Remove ";" for required variable
;*******************************************
u_wind = f->U_GRD_221_ISBL(19,:,:) ; (lev,lat,lon)
u-wind_component
v_wind = f->V_GRD_221_ISBL(19,:,:)
u_wind@lat2d = lat
u_wind@lon2d = lon
v_wind@lat2d = lat
v_wind@lon2d = lon
speed = sqrt(u_wind^2+v_wind^2) ; Calculate windspeed
; copy_VarCoords(u_wind,speed) ; copy coord vars to speed
;*************create plots*******************
wks = gsn_open_wks ("ps", "winds_scalar") ; Create output
filename and type
gsn_define_colormap (wks,"gui_default") ; choose color map
res = True ; plot mods desired for
original grid
; res@cnFillOn = True ; turn on color
;************************************************************
; Define NARR Map
;
; Define how map is to be displayed [gridlines, boundaries]
res@mpLimitMode = "Corners" ; choose range of map
res@mpLeftCornerLatF = 24
res@mpLeftCornerLonF = -125
res@mpRightCornerLatF = 45
res@mpRightCornerLonF = -65
res@tfDoNDCOverlay = False
res@mpProjection = "LambertConformal"
res@mpLambertParallel1F = 50
res@mpLambertParallel2F = 50
res@mpLambertMeridianF = 253
res@pmTickMarkDisplayMode = "Always"
res@mpFillOn = False ; turn off map fill
res@mpDataBaseVersion = "MediumRes" ; Allow plotting of counties
res@mpOutlineDrawOrder = "PostDraw" ; draw continental outline
last
res@mpProvincialLineColor = "Blue" ; Border color
res@mpProvincialLineThicknessF = 5 ; Border line thickness
; res@mpOutlineBoundarySets = "AllBoundaries" ; state/county boundaries
res@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries
;******************************************************************
; Labels
res@tiMainString = "Vector CONUS at "+lev(19)+ "mb "
; res@tiXAxisString = "Longitude"
; res@tiYAxisString = "Latitude"
res@vcRefMagnitudeF = 4.0 ; add a reference vector
res@vcRefLengthF = 0.02 ; what the ref length is
res@gsnSpreadColors = True ; use full colormap
res@lbOrientation = "Horizontal" ; horizontal label bar
; res@pmLabelBarOrthogonalPosF = -0.01 ; move label bar closer
; res@lbLabelStride = 4
;*****************************************************************
; Plot Map
res@vcGlyphStyle = "CurlyVector" ; turn on curley vectors
res@vcMinDistanceF = 0.017 ; thin out vectors
res@gsnMaximize = True
plot = gsn_csm_vector_scalar_map(wks,u_wind,v_wind,speed,res)
end
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 4 11:15:22 2013
This archive was generated by hypermail 2.1.8 : Sun Apr 07 2013 - 21:13:30 MDT