Zooming in on native lambert conformal data

Hi,

I'm having troubles plotting native lambert conformal data from my netcdf wrf output...

A simple script is below, basically, I want to plot zoomed views of my domain. The first plot from the script (the entire domain) works fine but if I specify new "corners" and plot, I get a zoomed background map but the data is at its original scale.

I would prefer to keep everything Lambert Conformal but am prepared to change if neccessary!

Any help much appreciated.

Thanks very much,

Tom

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin

  f = addfile("wrfout_d01_2005-11-23_03:00:00.nc", "r")

  lat2d = f->XLAT(0,:,:)
  lon2d = f->XLONG(0,:,:)
  dimll = dimsizes(lat2d)
  nlat = dimll(0)
  mlon = dimll(1)

  x = f->HGT(0,:,:) ; surface height

  wks = gsn_open_wks("x11" ,"WRF_test") ; ps,pdf,x11,ncgm,eps
  res = True ; plot mods desired

  res_at_gsnAddCyclic = False ; regional data: not cyclic
  res_at_tfDoNDCOverlay = True ; set True for native projection

  res_at_gsnMaximize = True
  res_at_mpFillOn = True
  res_at_pmTickMarkDisplayMode = "Always"

  res_at_mpProjection = "LambertConformal"
  res_at_mpLambertParallel1F = f_at_TRUELAT1
  res_at_mpLambertParallel2F = f_at_TRUELAT2
  res_at_mpLambertMeridianF = f_at_CEN_LON

  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)

  plot = gsn_csm_contour_map(wks, x, res) ; This works fine

  res_at_mpLimitMode = "Corners"
  res_at_mpLeftCornerLatF = -45
  res_at_mpLeftCornerLonF = 170
  res_at_mpRightCornerLatF = -40
  res_at_mpRightCornerLonF = 175

  plot = gsn_csm_contour_map(wks, x, res) ; The map is zoomed but not the data (x)

end
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 06 2005 - 15:24:14 MST