I am sorry to have to give you so much "stuff" (pic and scripts), but I do have a question about suppressing or diminishing the "_Missing Values" for U / V plots.... I happen to have some ocean current (netcdf) files for the Northern Gulf (at sealevel) and the large arrows that you can see on the vector plot (below) are the "Missing_Values" data that I cannot seem to find anything (parameters, attributes, etc.) in the resources (list of ncl attributes) to stop these large ("Missing_Values") from being displayed. Please let me know if you might be able to give me any ideas on this ...
(Example Plot is below and also, I have included my entire ncl-script just below the plot)..
Thanks,
John McCormick
========================================================
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"
begin
; dir = ncargpath("data")
; uf = addfile(dir+"/cdf/Ustorm.cdf","r") ; Open two netCDF files.
; vf = addfile(dir+"/cdf/Vstorm.cdf","r")
uf = addfile("/amber/home/jpu/GulfBreeze/VectorPlots/u3d_2009061100_nc_ascii.nc","r") ; Open two netCDF files.
vf = addfile("/amber/home/jpu/GulfBreeze/VectorPlots/v3d_2009061100_nc_ascii.nc","r")
u = uf->uu(0,0,:,:) ; Get time, dep, lat, lon, data.
v = vf->vv(0,0,:,:)
lat = vf->v&lat
lon = vf->v&lon
; wks = gsn_open_wks("x11","color") ; Open a workstation.
wks = gsn_open_wks("png","color") ; Open a workstation.
gsn_define_colormap(wks, "rainbow")
nlon = dimsizes(lon)
nlat = dimsizes(lat)
; #####################################################################
; #####################################################################
; ###
; ### Map Base for U-V plot -- ("mapbase")
; ###
; #####################################################################
; #####################################################################
levels = (/30,35,40,45,50,55,60,65,70,75,80,85,90,95,100/)
labels = (/"30","35","40","45","50","55","60","65","70","75","80", "85",\
"90","95","100"/)
nlevels = dimsizes(levels)
mapbase = gsn_csm_map(wks,res1) ; Draw map.
; draw(mapbase)
; frame(wks)
res = True
vplot = gsn_vector(wks,u,v,res)
draw(vplot)
frame(wks)
; vplot = gsn_vector_map(wks,u,v,res) ; Draw a vector plot of u and v.
overlay(mapbase,vplot)
draw(mapbase)
frame(wks)
; vmap = gsn_vector_scalar(wks,u,v,False,res) ; Draw a vector plot of u and v.
; plot = gsn_vector_scalar(wks,u,v,res) ; Draw a vector plot of
; map = gsn_vector_map(wks,u,v,res) ; Draw a vector plot of
; u and v over a map.
end
John W. McCormick, Systems Programmer / DBA
High Performance Computing & Scientific Visualization
Lockheed Martin Information Technology, Supporting the EPA
Research Triangle Park, NC
919-541-0890
Federal Infrastructure Contact - Ravi Nair
919-541-5467 -
nair.ravi@epa.govFederal Infrastructure Contact - Joe Retzer, Ph.D.
919-541-4190 -
retzer.joseph@epa.gov############################################
### Work Information ###
############################################
John W. McCormick, Systems Programmer, MD: N127-01
Lockheed Martin, Contractor to US-EPA
US-EPA, Bldg NCC, Rm N135B
Durham, NC 27711
Office: 919-541-0890
Fax : 919-541-0056
Work Email:
mccormick.johnw@epa.gov