diff -uNr GMT4.1.1/man/manl/GMT.l GMT4.1.1_new/man/manl/GMT.l
--- GMT4.1.1/man/manl/GMT.l	2006-03-01 15:22:47.000000000 +0900
+++ GMT4.1.1_new/man/manl/GMT.l	2006-04-21 01:27:26.000000000 +0900
@@ -52,6 +52,8 @@
 .br
 \fBgrd2xyz\fP		Conversion from 2-D gridded file to table data
 .br
+\fBgrdbarb\fP		Plot wind barbs from grdfiles
+.br
 \fBgrdblend\fP		Blend several partially over-lapping grdfiles onto one grid
 .br
 \fBgrdclip\fP		Limit the z-range in gridded data
diff -uNr GMT4.1.1/man/manl/grdbarb.l GMT4.1.1_new/man/manl/grdbarb.l
--- GMT4.1.1/man/manl/grdbarb.l	1970-01-01 09:00:00.000000000 +0900
+++ GMT4.1.1_new/man/manl/grdbarb.l	2006-04-21 01:23:26.000000000 +0900
@@ -0,0 +1,217 @@
+.TH GRDBARB l "13 Mar 2005" "GMT4.0" "GMT4.0"
+.SH NAME
+grdbarb \- Plot wind barbs from grdfiles
+.SH SYNOPSIS
+\fBgrdbarb\fP \fIcompx.grd\fP \fIcompy.grd\fP \fB\-J\fP\fIparameters\fP [ \fB\-A\fP ]
+[ \fB\-B\fP\fItickinfo\fP ] [ \fB\-C\fP\fIcptfile\fP ] [ \fB\-E\fP ]
+[ \fB\-G\fP\fIfill\fP] [ \fB\-I\fP\fIx_inc\fP[\fBm|c\fP][/\fIy_inc\fP[\fBm|c\fP]] ] [ \fB\-K\fP ] [ \fB\-N\fP ] [ \fB\-O\fP ] [ \fB\-P\fP ]
+[ \fB\-Q\fP\fIparameters\fP ] [ \fB\-R\fP\fIwest/east/south/north\fP[\fBr\fP] ] [ \fB\-S\fP[\fBl\fP]\fIscale\fP ]
+[ \fB\-T\fP ] [ \fB\-U\fP[\fI/dx/dy/\fP][\fIlabel\fP] ] [ \fB\-V\fP ] [ \fB\-W\fP\fIcontourpen\fP ] [ \fB\-X\fP\fIx-shift\fP ]
+[ \fB\-Y\fP\fIy-shift\fP ] [ \fB\-Z\fP ] [ \fB\-c\fP\fIcopies\fP ]
+.SH DESCRIPTION
+\fBgrdbarb\fP reads two 2-D gridded files which represents the x- and y-components of a vector field and
+produces a wind field plot by
+drawing barbs with orientation and length according to the information in the files. Alternatively,
+polar coordinate components may be used (r, theta).
+\fBgrdbarb\fP is basically a short-hand for using 2 calls to \fBgrd2xyz\fP and pasting the output through
+\fBpsxy \-SU\fP.
+.TP
+\fIcompx.grd\fP
+Contains the x-component of the vector field.
+.TP
+\fIcompy.grd\fP
+Contains the y-component of the vector field.
+.TP
+.B \-J
+Selects the map projection. Scale is UNIT/degree, 1:xxxxx, or width in UNIT (upper case modifier).
+UNIT is cm, inch, or m, depending on the MEASURE_UNIT setting in \.gmtdefaults4, but this can be
+overridden on the command line by appending \fBc\fP, \fBi\fP, or \fBm\fP to the scale/width value.
+For map height, max dimension, or min dimension, append \fBh\fP, \fB+\fP, or \fB-\fP to the width,
+respectively.
+.br
+.sp
+\fBCYLINDRICAL PROJECTIONS:\fP
+.br
+.sp
+\fB\-Jc\fP\fIlon0/lat0/scale\fP (Cassini)
+.br
+\fB\-Jj\fP\fIlon0/scale\fP (Miller)
+.br
+\fB\-Jm\fP\fIscale\fP (Mercator - Greenwich and Equator as origin)
+.br
+\fB\-Jm\fP\fIlon0/lat0/scale\fP (Mercator - Give meridian and standard parallel)
+.br
+\fB\-Joa\fP\fIlon0/lat0/azimuth/scale\fP (Oblique Mercator - point and azimuth)
+.br
+\fB\-Job\fP\fIlon0/lat0/lon1/lat1/scale\fP (Oblique Mercator - two points)
+.br
+\fB\-Joc\fP\fIlon0/lat0/lonp/latp/scale\fP (Oblique Mercator - point and pole)
+.br
+\fB\-Jq\fP\fIlon0/scale\fP (Equidistant Cylindrical Projection (Plate Carree))
+.br
+\fB\-Jt\fP\fIlon0/scale\fP (TM - Transverse Mercator, with Equator as y = 0)
+.br
+\fB\-Jt\fP\fIlon0/lat0/scale\fP (TM - Transverse Mercator, set origin)
+.br
+\fB\-Ju\fP\fIzone/scale\fP (UTM - Universal Transverse Mercator)
+.br
+\fB\-Jy\fP\fIlon0/lats/scale\fP (Basic Cylindrical Projection)
+.br
+.sp
+\fBAZIMUTHAL PROJECTIONS:\fP
+.br
+.sp
+\fB\-Ja\fP\fIlon0/lat0/scale\fP (Lambert).
+.br
+\fB\-Je\fP\fIlon0/lat0/scale\fP (Equidistant).
+.br
+\fB\-Jf\fP\fIlon0/lat0/horizon/scale\fP (Gnomonic).
+.br
+\fB\-Jg\fP\fIlon0/lat0/scale\fP (Orthographic).
+.br
+\fB\-Js\fP\fIlon0/lat0/\fP[\fIslat/\fP]\fIscale\fP (General Stereographic)
+.br
+.sp
+\fBCONIC PROJECTIONS:\fP
+.br
+.sp
+\fB\-Jb\fP\fIlon0/lat0/lat1/lat2/scale\fP (Albers)
+.br
+\fB\-Jd\fP\fIlon0/lat0/lat1/lat2/scale\fP (Equidistant)
+.br
+\fB\-Jl\fP\fIlon0/lat0/lat1/lat2/scale\fP (Lambert)
+.br
+.sp
+\fBMISCELLANEOUS PROJECTIONS:\fP
+.br
+.sp
+\fB\-Jh\fP\fIlon0/scale\fP (Hammer)
+.br
+\fB\-Ji\fP\fIlon0/scale\fP (Sinusoidal)
+.br
+\fB\-Jk\fP[\fBf|s\fP]\fIlon0/scale\fP (Eckert IV (f) and VI (s))
+.br
+\fB\-Jn\fP\fIlon0/scale\fP (Robinson)
+.br
+\fB\-Jr\fP\fIlon0/scale\fP (Winkel Tripel)
+.br
+\fB\-Jv\fP\fIlon0/scale\fP (Van der Grinten)
+.br
+\fB\-Jw\fP\fIlon0/scale\fP (Mollweide)
+.br
+.sp
+\fBNON-GEOGRAPHICAL PROJECTIONS:\fP
+.br
+.sp
+\fB\-Jp\fP[\fBa\fP]\fIscale\fP[\fI/origin\fP] (polar (theta,r) coordinates, optional \fBa\fP for azimuths and offset theta [0])
+.br
+\fB\-Jx\fP\fIx-scale\fP[\fBl|p\fP\fIpow\fP][\fI/y-scale\fP[\fBl|p\fP\fIpow\fP]][\fBd\fP] (Linear, log, and power scaling)
+.br
+More details can be found in the \fBpsbasemap\fP man pages.
+.br
+.SH OPTIONS
+No space between the option flag and the associated arguments.
+.TP
+.B \-A
+Means grdfiles have polar (r, theta) components instead of Cartesian (x, y).
+.TP
+.B \-B
+Sets map boundary annotation and tickmark intervals; see the
+\fBpsbasemap\fP man page for all the details.
+.TP
+.B \-C
+Use \fIcptfile\fP to assign colors based on vector length.
+.TP
+.B \-E
+Center wind barbs on grid nodes [Default draws from grid node].
+.TP
+.B \-G
+Sets color or shade for barb interiors [Default is no fill].
+Specify the shade (0\-255) or color (r/g/b, each in 0\-255, h/s/v, with
+h in 0\-360, s, and v in 0\-1, and set \fBCOLOR_MODEL\fP = hsv, or
+c/m/y/k, each in 0\-100 %).
+Alternatively, specify \fB\-Gp\fP\fIdpi/pattern\fP, where \fIpattern\fP gives
+the number of the built-in pattern (1-90) \fIor\fP the name of a Sun 1-, 8-,
+or 24-bit raster file. The \fIdpi\fP sets the resolution of the image. For
+1-bit rasters: use \fB\-GP\fP for inverse video, or append
+\fB:F\fP\fIr/g/b\fP[\fBB\fP[\fIr/g/b\fP]] to specify fore- and background
+colors (use \fIr/g/b\fP = - for transparency).
+See \fBGMT\fP Cookbook & Technical Reference Appendix E for information
+on individual patterns.
+.TP
+.B \-I
+Only plot barbs at nodes every \fIx_inc, y_inc\fP apart (must be multiples of
+original grid spacing). Append \fBm\fP for minutes or \fBc\fP for seconds. [Default plots every node].
+.TP
+.B \-K
+More \fIPostScript\fP code will be appended later [Default terminates the plot system].
+.TP
+.B \-N
+Do NOT clip barbs at map boundaries [Default will clip].
+.TP
+.B \-O
+Selects Overlay plot mode [Default initializes a new plot system].
+.TP
+.B \-P
+Selects Portrait plotting mode [\fBGMT\fP Default is Landscape, see \fBgmtdefaults\fP to change this].
+.TP
+.B \-Q
+Select barb plot. Optionally, specify \fIparameters\fP which
+are \fIbarbwidth/barblength/barbangle/barbscale\fP [Default is 0.025\fBc\fP/0.051\fBc\fP/120/5 (or 0.1\fBi\fP/0.2\fBi\fP/120/5)]. 
+.TP
+.B \-R
+\fIxmin, xmax, ymin,\fP and \fIymax\fP specify the Region of interest. For geographic
+regions, these limits correspond to \fIwest, east, south,\fP and \fInorth\fP and you may specify them
+in decimal degrees or in [+-]dd:mm[:ss.xxx][W|E|S|N] format. Append \fBr\fP if lower left and upper right
+map coordinates are given instead of wesn. The two shorthands \fB\-Rg \-Rd\fP stand for global domain
+(0/360 or -180/+180 in longitude respectively, with -90/+90 in latitude).
+For calendar time coordinates you may either give relative
+time (relative to the selected TIME_EPOCH and in the selected TIME_UNIT; append \fBt\fP to \fB\-JX|x\fP),
+or absolute time of the form [\fIdate\fP]\fBT[\fP\fIclock\fP] (append \fBT\fP to \fB\-JX|x\fP). At least one of \fIdate\fP and \fIclock\fP
+must be present; the \fBT\fP is always required. The \fIdate\fP string must be of the form [-]yyyy[-mm[-dd]]
+(Gregorian calendar) or yyyy[-Www[-d]] (ISO week calendar), while the \fIclock\fP string must be of
+the form hh:mm:ss[.xxx]. The use of delimiters and their type and positions must be as indicated
+(however, input/output and plotting formats are flexible).
+Specify a subset of the grid.
+.TP
+.B \-T
+Means azimuth should be converted to angles based on the selected map projection.
+.TP
+.B \-U
+Draw Unix System time stamp on plot. User may specify where the lower left corner
+of the stamp should fall on the page relative to lower left corner of plot. Optionally,
+append a label, or \fBc\fP (which will plot the command string.). The \fBGMT\fP parameters
+UNIX_TIME and UNIX_TIME_POS can affect the appearance; see the \fBgmtdefaults\fP man
+page for details.
+.TP
+.B \-V
+Selects verbose mode, which will send progress reports to stderr [Default runs "silently"].
+.TP
+.B \-W
+Set pen attributes used for vector outlines [Default: width = 1, color = 0/0/0, texture = solid].
+.TP
+.B \-X \-Y
+Shift origin of plot by (\fIx-shift,y-shift\fP).
+Prepend \fBa\fP for absolute coordinates;
+the default (\fBr\fP) will reset plot origin.
+Give \fBc\fP to center plot using current page size.
+.TP
+.B \-Z
+Means the angles provided are azimuths rather than direction (requires \fB\-A\fP).
+.TP
+.B \-c
+Specifies the number of plot copies. [Default is 1].
+.SH EXAMPLES
+To draw the vector field given by the files r.grd and theta.grd on a linear plot
+with scale 5 cm per data unit,
+using vector rather than stick plot, and scale vector magnitudes so that 10 units
+equal 1 inch, run
+.br
+.sp
+\fBgrdvector\fP r.grd theta.grd \fB\-Jx\fP5\fBc\fP \fB\-A \-Q\fP \fB\-S\fP10\fBi\fP > gradient.ps
+.br
+.sp
+.SH "SEE ALSO"
+.IR gmt (l),
+.IR grdcontour (l),
+.IR psxy (l)
diff -uNr GMT4.1.1/man/manl/psxy.l GMT4.1.1_new/man/manl/psxy.l
--- GMT4.1.1/man/manl/psxy.l	2006-03-01 15:22:48.000000000 +0900
+++ GMT4.1.1_new/man/manl/psxy.l	2006-04-21 01:23:26.000000000 +0900
@@ -414,6 +414,14 @@
 .B \-St
 \fBt\fPriangle.  \fIsize\fP is diameter of circumscribing circle.
 .TP
+.B \-Su
+wind barb. Direction (in degrees counter-clockwise from horizontal) and length must be found in columns 3 and 4. \fIsize\fP,
+if present, will be interpreted as barbwidth/barblength/barbangle/barbscale [Default is 0.025\fBc\fP/0.051\fBc\fP/120/5 (or 0.1\fBi\fP/0.2\fBi\fP/120/5)].
+.B \-SU
+Same as \fB\-Su\fP, except azimuth (in degrees east of north) should be given instead of direction. The azimuth will
+be mapped into an angle based on the chosen map projection (\fB\-Su\fP leaves the directions
+unchanged.)
+.TP
 .B \-Sv
 \fBv\fPector.  Direction (in degrees counter-clockwise from horizontal) and length must be found in columns 3 and 4.  \fIsize\fP,
 if present, will be interpreted as arrowwidth/headlength/headwidth [Default is 0.075\fBc\fP/0.3\fBc\fP/0.25\fBc\fP (or 0.03\fBi\fP/0.12\fBi\fP/0.1\fBi\fP)].