Draw regional figures using global data

Hi:

I have a global dataset (see z(17,73,144)) and I'd like to draw figures
only for the North Pacific and the North America. The figure I get covers
the global with contour lines outside 0-80N and 120-300E were set as
missing. The contours west of 180E are also separated from the contours
east of 180E (see the attachment for example). Can anybody help me with
this problem? The codes I used to draw the figure are listed below,

  wks = gsn_open_wks ("pdf", "test" ) ; define the name of the plot
  plot = new(3, graphic) ; create a plot array
  res1 = True
  res2 = True
  res3 = True

  res1_at_tiMainString = "El Nino years"
  res2_at_tiMainString = "La Nina years"
  res3_at_tiMainString = "Warm-Cold years"

  res1_at_gsnDraw = False ; don't draw
  res1_at_gsnFrame = False ; don't advance frame
  res1_at_gsnAddCyclic = False
  res1_at_gsnContourNegLineDashPattern = 14
  res1_at_gsnContourZeroLineThicknessF = 3
  res2_at_gsnDraw = False ; don't draw
  res2_at_gsnFrame = False ; don't advance frame
  res2_at_gsnAddCyclic = False
  res2_at_gsnContourNegLineDashPattern = 14
  res2_at_gsnContourZeroLineThicknessF = 3
  res3_at_gsnDraw = False ; don't draw
  res3_at_gsnFrame = False ; don't advance frame
  res3_at_gsnAddCyclic = False
  res3_at_gsnContourNegLineDashPattern = 14
  res3_at_gsnContourZeroLineThicknessF = 3

  plot(0) = gsn_csm_contour_map(wks,zw(2,{0:80},{120:300}), res1)
  plot(1) = gsn_csm_contour_map(wks,zc(2,{0:80},{120:300}), res2)
  plot(2) = gsn_csm_contour_map(wks,ztest(2,{0:80},{120:300}), res3)
 gsn_panel(wks,plot,(/3,1/),False) ; now draw as one plot

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk