shift in wind field

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