plotting wind between two values

From: ahmed lasheen <ahmed4kernel_at_nyahnyahspammersnyahnyah>
Date: Wed Apr 18 2012 - 10:12:56 MDT

Hello
I am trying to plot wind between two values using the resources
vcres@vcLevelSelectionMode ="ManualLevels"
vcres@vcMinLevelValF = 0
vcres@vcMaxLevelValF =5
but it doesnot work
any suggestion?
the script , I download it from the NCL website and i just add the previous
three line to it .
----------------------------------------------------------------------

; vector_3.ncl
;
; Concepts illustrated:
; - Drawing a black-and-white vector plot over a cylindrical equidistant
map
; - Adding a time stamp to a plot
; - Thinning vectors using a minimum distance resource
; - Moving the vector reference annotation to the top right of the plot
;----------------------------------------------------------------------
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
  a = addfile("uv300.nc","r")
;************************************************
; read in zonal [u] and meridional [v] winds (July)
;************************************************
  u = a->U(1,:,:)
  v = a->V(1,:,:)
;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("x11","vector") ; open a ps file

  vcres = True ; plot mods desired
  vcres@gsnFrame = False ; so we can draw time
stamp
  vcres@vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector up
  vcres@vcRefMagnitudeF = 10.0 ; define vector ref mag
  vcres@vcRefLengthF = 0.045 ; define length of vec
ref
  vcres@vcGlyphStyle = "CurlyVector" ; turn on curly vectors
  vcres@vcMinDistanceF = 0.017

vcres@vcLevelSelectionMode ="ManualLevels"
vcres@vcMinLevelValF = 0 ; m/s
vcres@vcMaxLevelValF =5; m/s

  vector = gsn_csm_vector_map_ce(wks,u,v,vcres)
;************************************************
; add time stamp to plot
;************************************************
  infoTimeStamp(wks,0.025,"uv300.nc") ; arguments are workstation
                                           ; size of text desired, and
                                           ; right string.
  frame (wks)
end

thanks in advance

-- 
===============
Ahmed Lasheen
Junior researcher at Cairo Numerical Weather Prediction Center (CNWPC)
Egyptian Meteorological Authority(EMA)
Cairo,Egypt
===============

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Apr 18 10:13:50 2012

This archive was generated by hypermail 2.1.8 : Fri Apr 20 2012 - 16:21:18 MDT