overlay two more variables

From: Jin Young Kim <Jin-Young.Kim_at_nyahnyahspammersnyahnyah>
Date: Mon, 11 May 2009 18:38:55 -0600

Could someone help me diagnose a problem I'm having when overlaying
wind vector and scalar contours plot? Here's a part of script what I'm
trying
to do.

  minlat = 0
  maxlat = 60
  minlon = 200
  maxlon = 320

  gfs_g500 = smth9(gfs_g500,0.5,-0.25,False)
  jetstm = sqrt(gfs_u200*gfs_u200+gfs_v200*gfs_v200)

; note that the gsn_csm_*map_ce templates automatically set
; res_at_mpLimitMode="LatLon" for you. If you are plotting a different
projection,
; you may have to set this resource.

  mres = True
  mres_at_mpLimitMode = "LatLon"
  mres_at_mpMinLatF = minlat ; range to zoom in on
  mres_at_mpMaxLatF = maxlat
  mres_at_mpMinLonF = minlon
  mres_at_mpProjection = "LambertEqualArea"
  mres_at_mpCenterLonF = (mres_at_mpMaxLonF + mres_at_mpMinLonF)/2.
  mres_at_gsnDraw = False
  mres_at_gsnFrame = False

  map = gsn_csm_map(wks,mres)

  res = True

  res_at_gsnFrame = False
; res_at_gsnMaximize = True

  res_at_cnFillOn = False
  res_at_cnLinesOn = True ; turn off contour lines
  res_at_cnLineLabelsOn = True ; turn off contour line
labels
  res_at_cnLevelSelectionMode = "ManualLevels"
  res_at_cnLevelSpacingF = 40 ; set contour spacing
  res_at_cnLineThicknessF = 2. ;line thickness

  g500 = gsn_csm_contour(wks,gfs_g500,res)
  overlay(map, g500)

  cnres= True
  cnres_at_cnFillOn = True
  cnres_at_gsnSpreadColors = True
  cnres_at_cnLinesOn = True
  cnres_at_cnLineLabelsOn = True
  cnres_at_cnLevelSelectionMode = "ManualLevels" ; set manual
contour levels
  cnres_at_cnMinLevelValF = 25. ; set min contour level

  jet200 = gsn_csm_contour(wks,jetstm,cnres)
  overlay(map, jet200)

  vcres = True

  vcres_at_vcRefMagnitudeF = 10. ; make vectors larger
  vcres_at_vcRefLengthF = 0.025 ; ref vec length
  vcres_at_vcGlyphStyle = "WindBarb"
  vcres_at_vcMinDistanceF = 0.015 ; thin out windbarbs
  vcres_at_vcRefAnnoOn = False ; turn off ref wind
barb
  vcres_at_vcWindBarbTickLengthF = 0.5 ; length of the
barb tick

  wb500 = gsn_csm_vector(wks,gfs_u500,gfs_v500,vcres)

 overlay(map, wb500)
  draw(map)
  frame(wks)
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon May 11 2009 - 18:38:55 MDT

This archive was generated by hypermail 2.2.0 : Wed May 13 2009 - 08:40:16 MDT