A question for you experts:
I'm plotting the wind speed field of a regional model, and find the field has
a southeastward shift relative to the map of the region.
The model projection is "polar_stereographic", and I plotted the field
as "LambertConformal". The plotting part of the file is copied here.
---------------------
w = ........
lat2d = f->lat
lon2d = f->lon
dimlc = dimsizes(lat2d)
nlat = dimlc(0)
mlon = dimlc(1)
wks = gsn_open_wks ("eps", "wind")
res = True
res_at_cnFillOn = True
res_at_cnLinesOn = False
res_at_gsnSpreadColors = True
res_at_mpGridAndLimbOn = True
res_at_pmTickMarkDisplayMode = "Always"
res_at_tmXTOn = False
res_at_gsnAddCyclic = False
res_at_tfDoNDCOverlay = True
res_at_sfXArray = lon2d
res_at_sfYArray = lat2d
res_at_mpLimitMode = "Corners"
res_at_mpLeftCornerLatF = lat2d(0,0)
res_at_mpLeftCornerLonF = lon2d(0,0)
res_at_mpRightCornerLatF = lat2d(nlat-1,mlon-1)
res_at_mpRightCornerLonF = lon2d(nlat-1,mlon-1)
res_at_mpProjection = "LambertConformal"
res_at_mpLambertParallel1F = 30.
res_at_mpLambertParallel2F = 60.
res_at_mpLambertMeridianF = -97.
res_at_mpOutlineDrawOrder = "PostDraw"
res_at_mpOutlineBoundarySets = "GeophysicalAndUSStates"
res_at_cnLevelSelectionMode = "ManualLevels"
res_at_cnMinLevelValF = 0.0
res_at_cnMaxLevelValF = 4.0
res_at_cnLevelSpacingF = 0.5
plot = gsn_csm_contour_map(wks,w,res)
end
------------------------------------
Thanks for help.
Eric
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Mar 17 2009 - 14:49:26 MDT
This archive was generated by hypermail 2.2.0 : Wed Mar 18 2009 - 14:50:21 MDT