How to plot boundary line

Hi Everyone,

I need to plot many fields in the same map but I can not get the boundary map and there is a problem of mismatching the plots.
Here is my code:

load"/home/NCARG/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load"/home/NCARG/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load"/home/NCARG/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
load"/home/NCARG/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin
  a = addfile("wrfhrly_d02_2000-06-01.nc", "r")
  wks = gsn_open_wks("x11","test")
  gsn_define_colormap(wks,"WhiteYellowOrangeRed")

  time = 1
  tc2 = wrf_user_getvar(a,"T2",time) ; T2 in Kelvin
  tc2 = tc2-273.16
  tc2@description = "Surface Temperature"

  u10 = wrf_user_getvar(a,"U10",time) ; u at 10 m
  v10 = wrf_user_getvar(a,"V10",time) ; v at 10 m
  psl = wrf_user_getvar(a,"PSFC",time)
  psl = (psl/100)

; TEMPERATURE 2 M
  rest = True
  rest@cnFillOn = True
  rest@UnitLabel = "C"
  rest@cnLevelSelectionMode = "ManualLevels"
  rest@cnMinLevelValF = -1.
  rest@cnMaxLevelValF = 30.
  rest@cnLevelSpacingF = 1.
  contour = wrf_contour(a,wks,tc2,rest)

; WIND VECTOR
  resw = True
  resw@UnitLabel = "m/s"
  resw@NumVectors = 30
  vector = wrf_vector(a,wks,u10,v10,resw)

; PRESSURE
  resp = True
  resp@cnLineColor = "White"
  resp@UnitLabel = "hPa"
  resp@InitTime = False
  resp@Footer = False
  resp@cnLevelSelectionMode = "ManualLevels"
  resp@cnMinLevelValF = 850.
  resp@cnMaxLevelValF = 1020.
  resp@cnLevelSpacingF = 20.
  contourp = wrf_contour(a,wks,psl,resp)

  wrf_overlay(wks,(/contour,contourp,vector/),True)

end

Thank you very much for you help,

Virginia

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon May 26 21:00:24 2014