CHANGES IN GMT from v3.0 - 3.3
--------------------------------------

		@(#)CHANGES	3.144  07/01/99

1. Overview

This file shows all details regarding specific enhancements and bug fixes
since the 1995 release (3.0).  Unfortunately, sometimes we
forget to add an entry so there may actually be more fixes than these.

2. Bug fixes.

The following is a list of specific bug fixes and enhancements since the
release of version 3.0 in August 1995.

Date:		Program:	Fix:
-----------------------------------------------------------------------------------------
===> 08/11/1995		tar 3.0 made and released
-----------------------------------------------------------------------------------------
08/14/1995	gmt_customio.c		_bin_read and _short_read allocated insufficient
					memory when subregion was specified.
		gmt_support.c		bicubic interpolation now returns node value
					when point is within SMALL of the node position.
		gmt.txt			Fixed spelling mistakes
08/18/1995	gmt_plot.c		Added seconds and year to Unix time stamp.
		gmt_support.c		Skips blank lines in .cpt files
08/23/1995	gmt.h			Added inverse hyperbolic trig macros for AUX
		gmt_init.c		1:xxxx scale info for Hammer, Eckert, Mollweide,
					Plate Carree, Winkel, Robinson, Sinusoidal was
					set incorrectly (scaled by radius of Earth...)
		grdimage.c		size of subset image scaled incorrectly for -Jx
		pscoast.c		Painted lakes incorrectly - now works.
08/25/1995	gmt_customio.c		write pointers for *write_grd_info was set to
					read for id>0
08/28/1995	gmt_plot.c		Map scale can now also be in nautical miles (n).
08/29/1995	psxy.c			-Se didnt work without a dummy size.
					-Sf with lines outside -R crashed with memory
					error.
		gmt_defaults.h		flattening for WGS-84 was off a bit.
08/30/1995	gmt_funcnames.h 	Had redeclaration of gmt_input and gmt_output.
08/31/1995	gmt_map.c		map_path needed to have some slop in its == test
09/05/1995	gmt_shore.h		Removed uchar typedef since IBMs complain.
09/07/1995	grdview.c		With -Qs and a continuous .cpt file grdview
					would sometimes pick a z_ave outside the range
					of a polygon.  Now checks.
09/13/1995	gmt_plot.c		Second anotation screwed up for negative
					longitude.
09/18/1995	psscale.c		-V did nothing; now reports range in cpt file.
09/28/1995	pstext.c		Now longer needs leading hyphen to give font
					names.
		gmt_init.c		fonts in .gmtdefaults and via gmtset can be now
					be font # or name.
		psxyz.c			Radius instead of diameter expected for circles.
10/17/1995	gmt_init.h		Forgot to add -JQ/q to list of unique options.
10/24/1995	gmt_defaults.h		RECT_GRATICULE did not include linear proj.
10/25/1995	psxy.c			greenwich set incorrectly.
11/22/1995	pscoast.c		variable key was not initialized to NULL.
11/30/1995	grdsample.c		Used incorrect dx to get new longitudes
12/07/1995	grdmath.c		Had memory error for DUP and EXCH operators.
12/27/1995	grdcontour.c		Allocated mem before increaseing n_alloc
		gmt_plot.c		Added support for POW gridlines
01/05/1996	grdcontour.c		Option -N to append units to contour anotations
01/31/1996	pshistogram.c		Added log10 option and restored -R as optional.
02/01/1996	grdproject.l		Did not explain -M at all.
02/14/1996	gmt_map.c		General cylindrical (-JY) did not give right
					map width.
02/19/1996	pstext.c		Accepts 2-char codes for justification, e.g.,
					LL for lower left, CM = centered middle, etc.
					Added -D to shift the PROJECTED x,y point by
					dx/dy.
		gmt_init.c		Fixed problem with exponential notation in -B
					option!
02/20/1996	pscontour.c		Added -L option to draw underlying triangular
					mesh.
		grdcontour.c		Optionally specify shift in the -Z option
02/27/1996	psscale.c		Mentions that a unit label can also be set with
					-B.
03/01/1996	makecpt.c		Added -I to reverse sense of color/gray-ramps.
		minmax.c		Added -C for simple column output.
03/05/1996	triangulate.l		Had no explanation for -G option.
		surface.l		Had no mention of -Q
03/07/1996	grdlandmask.c		Mask values not set correctly for -Nwet/dry
		multiple files		Allow for grdfile extension shorthand so that
					file.i2 is interpreted as file.i2=2/1/0/32767,
					etc.  These suffix abbreviations are kept in
					~/.gmt_io and GRIDFILE_SHORTHAND must be set to
					TRUE (this is a new parameter)
		gmt_grdio.c		Added byte and bit binary scanline formats.
		grdmath.c		Added cosd, sind, tand, r2d, d2r, gdist
		grdgradient.c		Options to find steepest slope direction and
					|grad z|.
03/08/1996	psmask.c		Added -M for map unit distances.
		grdgradient.c		Now checks if input grdfile exists!
		grdfft.c		--"--
		grdhisteq.c		--"--
		grdtrend.c		--"--
		psclip.c		Replaced -S with -C
		psmask.c		exchanged -S and -C
		grdmask.c		Added -S as in psmask
		nearneighbor.c		removed -M; added k suffix to -S instead
		psmask.c		--"--
03/12/1996	grdhisteq.c		Now takes <norm> with -N as an option.
03/13/1996	*.txt			Fixed formatting errors thruout
		grdvector.c		Now explains -I
04/02/1996	triangulate.c		Now initialize the header struct
04/03/1996	xyz2grd.c		Now averages multiple node entries.
04/10/1996	grd2xyz.c		Option -S to suppress output when a node == NaN
04/15/1996	gmt_map/init.c		-JM takes the alternate and more complete form
					-JM|m<clon/clat/width|scale> to specify a
					meridian and standard parallell other than 0,0.
04/16/1996	gmt_init.c		Writing of .gmtcommands did not check if an
					option started with '-' (which it must).
04/24/1996	pshistogram.c		Bug in usage message caused core dump.
04/25/1996	mapproject.c		Added option for nautical miles output.
05/08/1996	grdmath.c		Needed to allocate stack space when 2 constants
					were operated on to produce a grid
05/24/1996	gmt.h			ps_heximage declared double rather than BOOLEAN
05/28/1996	grdtrack.c		Added -R to work on subsets
07/05/1996	grdmask.c		Added missing -H option
07/10/1996	gmtselect.c		New program gmtselect.c for subset extraction
					based on spatial criteria such as closeness to
					points or lines, inside/outside polygons, etc.
07/15/1996	*.*			Renamed all shared functions for increased
					portability.
08/16/1996	gmt_map.c		Fixed bug in gmt_ialbers and failure to use
					sphere.  gmt_eckert must use fabs(delta) in
					iteration.  Equatorial stereographic view should
					now work again.
08/23/1996	binned_river*		Fixed a stray straight line connection between
					2 segments.
		binned_shore*		A lake crossing Greenwich had level 1 instead
					of 2
08/26/1996	gmt_support.c		Fixed the BoundingBox problem for landscape
					plots
		pslib.c			---"---
11/01/1996	pstext.c		Text on back side of Earth did not get clipped.
11/18/1996	psxy.c			-Sf option incorrectly connected fault endpoints
					automatically.
		sample1d.c		Now reports line # where time is not strictly
					monotonically increasing.  Can now also handle
					monotonically decreasing time.
11/26/1996	psxyz.c			Did not check z-range on data.
01/28/1997	pslib.c			Did not recognize user patterns.
02/07/1997	pshistogram.c		Now can do transpartent histogram using -L
					thanks to John L Lillibridge.
03/20/1997	gmt_init.c		Bug that messed up DPI for metric users have
					been killed.
03/28/1997	makecpt.c		Read problem with old-style tables fixed
					Can now output HSV values with or without -Z
					No more negative values (i.e. -1) for r,g,b
		psscale.c		Uses D_FORMAT when no decimal places are needed
		gmt_plot.c		-----------"-------------
		grdview.c		New option -L for bilinear rather than bicubic
					resampling. When resampling now checks for NaN
					return values.
		xyz2grd.c		Now can handle -R<w/s/e/n>r syntax
		psxy.c			Now deals correctly with combination -Sc and -Ey
		sample1d.c		Options -A -C -L now taken as -Fa, -Fl, -Fc
					(but is backwards compatible).
04/16/1997	grdvector.c		Added PI instead of 180 for Cartesian case...
04/30/1997	splitxyz.c		Now includes the last input point.
08/12/1997	psclip.c		No longer sets -O implicitly with -C
08/26/1997	grdtrack.c		Added -S to suppress NaN output
		*.c			Changed rint() to irint() where integer result
					is needed.  This should fix roundoff in log
					-anotations on HPs.
09/11/1997	gmt_cdf.c		Commented out old netCDF-2 function calls and
					rewrote using netCDF-3 with the new syntax
		gmt_unique.h		New include file
		gmt_keywords.h		----"----

09/14/1997	gmt_cdf.c		Finished updating function calls with netcdf 3.0
					(compiled but not tested)
		gmt_shore.c		Continue updating function calls with netcdf 3.0
					(not compiled yet)

09/22/1997	gmt_cdf.c		Tested and found errors; recompiled again
09/24/1997	grdview.c		Fixed bug that made pixel files go belly-up upon
					surface shading
09/28/1997	gmt_cdf.c		Added function calls check_nc_status; all
					nc_ functions will be checked for any errors.
09/29/1997	gmt_shore.c		Added function calls check_nc_status; all
					nc_ functions will be checked for any errors.
10/01/1997	gmt_cdf.c		Compiled, tested and found errors.
		gmt_shore.c		Compiled, tested and found errors.

10/09/1997	*.c			Updated all POSIX functions with the correct
					argument types; ex. free ((char *)) to
					free((void *)),fread, fwrite...

10/18/1997	makefile		Added -Xc (Ansi C) option to Sun compiler flags.

10/20/1997	grdimage.c		cast variable "bitimage" with "unsigned char".

10/21/1997	grdview.c		cast variable "bitimage" with "unsigned char".
		psscale.c

10/21/1997	pscontour.c		changed function "dump_contour" from type int
					to void; added function declaration in the
					beginning.

10/21/1997	pstext.c		erase "break" statement in the "default" switch.

10/22/1997	fitcircle.c		temporarily added #define M_PI_2 in the file;
					because the file won't compile even with math.h.

10/22/1997	gmt_shore.c		cast variable "c->node_level" with "unsigned
					char".

10/22/1997	grdraster.c		rename gmt functions because they are missing
					"GMT" in the beginning of the name of the
					functions; or "gmt" was in lower case.

11/04/1997	gmt_funcnames.h		added some of the function prototypes that are
					missing.  See file "functions" for detail.

11/08/1997	pscoast.c		added "&" in front of variable "ytmp[2]"
		gmt.h			moved around include files for convenience
02/19/1998	*.*			** All is now POSIX compliant and fully ANSI-C **
					Using void * pointers where old rule was char *
		gmt_plot.c		checker-board basemaps now start with correct
					shade of black or white
02/20/1998	gmt_map.c		Put the Robinson table into code; table file
					no longer used.
02/23/1998	*.*			Added GMT_IO structure to facilitate table i/o
		project.c		Support for multiple files, multisegment, and
					binary i/o
					Old -M renamed -Q and -M added for multisegment.
					Added -Z as in splitxyz for xy input only (no z)
030/1/1998	block*.c		Now quietly eats -M subheaders
		grdpaste		used -B which is reserved; now uses neither
					-A or -B but -Sa and -Sb.
03/03/1998	grdclip.c		Used -B also; now uses -Sa and -Sb.
03/04/1998	gmt*			Added -Jj -JJ Miller cylindrical projection
03/06/1998	gmt*			Added -Jf -JF Gnomonic azimuthal projection
		psmegaplot		Made part of supplemental archive
03/18/1998	psxy.c			Can now accept symbol code as part of input.

03/21/1998	*.c			Updated GMT_memory 2nd and 3rd argument to 
					(size_t), 1st argument to (void *)
		*.c			Changed buffer size 512 to BUFSIZ(1024)
		*.c			Updated memset 3rd argument to (size_t)
		*.c			Updated memcpy 3rd argument to (size_t)
		*.c			Updated fread 2nd and 3rd argument to (size_t)
		*.c			Updated fwrite 2nd and 3rd argument to (size_t)			
03/24/1998	*.c			-Gp now takes Sun rasterfiles instead of Sun
					icons and allow you to assign colors to 1-bit
					images.  8-bit & 24-bit images are also allowed.
		pscoast.c		Now may issue clip paths instead of painting
					polygons (makes pscoastclip obsolete).
03/26/1998	*.*			-F no longer used to set frame pen (use gmtset)
03/30/1998	psxy.c			Empty -G now means get -G string from subheaders
					in -M files (polygons only)
04/15/1998	gmt_support.c		ASCII read skips bad records (x and/or y is NaN)
04/17/1998	psscale.c		Now has -E to indicate back/forground colors +
					can use full -B option syntax
04/20/1998	gmt_grdio.*		Verify that -R and -I are compatible.
04/21/1998	xyz2grd.c		Now can deal with a variety of formats and
					ordering schemes, plus can byte swap.
04/22/1998	grd2xyz.c		Ditto.
04/27/1998	binned_shore_*		Fixed problems with 10x10 degree blocks (-Dl)
					getting wrong color.
04/28/1998	nearneighbor.c		Added option for boundary conditions
05/01/1998	*.*			Introduced gmt_nan.h which is created by the
					program gmt_nan_init and defines both the float
					and double NaNs patterns.
05/05/1998	pslib.c			Now write Boundingbox at end.
05/06/1998	gmt_support.c		Walter completed boundary functions
05/10/1998	grdimage.c		Fixed problem with global projections
		grdcontour.c		Fixed problem with -R...r setting
		grdimage.c		--"--
		grdvector.c		--"--
		shorelines		binned_shore_? renamed to binned_GSHHS_? to
					avoid confusion when using GMT 3.1
05/11/1998	grdview.c		Removed -D, merged -M and -W, and expanded -L
					to do BCs.
05/12/1998	grdinfo.c		Now uses D_FORMAT for output and has an output
					layout that is easier to parse with awk etc.
-----------------------------------------------------------------------------------------
===> 06/24/1998		tar 3.1b made and released to testers
-----------------------------------------------------------------------------------------
08/20/1998	gmt_map.c		Fixed clipping of data in polar projections.
08/20/1998	gmtselect.c		Always ate the first data record.  Now also uses
					Cartesian distances and units unless -J is set.
08/20/1998	*.c *.h			All GMT_boundcond_* functions now use pointers
					to struct.
08/21/1998	gmt_cdf.c		Use strlen(text)+1 so that the leading \0 is
					written to the grdfile.
		gmt_plot.c		Tried to free xings nc times in GMT_map_lontick
					(not good if nc > 1)
		pswiggle.c		Swapped filename and program-name
		gmt.h			Use POSIX-defined floor and ceil instead of
					buggy GMT macros.
		gmt_math_init.csh	Now also looks for an externally defined rint(x)
08/24/1998	psxy.c			Once again plots faults with -Sf.
					-SE instead of -Se expects azimuth and axes
					lengths in km.
08/25/1998	pslib_inc.h		Use POSIX-defined floor and ceil instead of
					buggy GMT macros.
08/26/1998	surface.c		converge_limit set before z_scale calculated.
					Thus, this limit was always 0.001 unless set
					with -C.
		nearneighbor.l		Makes it clear that nearneighbor only considers
					the nearest value inside each sector; not all
					values inside the radius.
		grdinfo.c		Option -C to format report on a single line for
					shell scripting convenience.
08/27/1998	*.*			FRAME_PEN, GRID_PEN, and TICK_PEN can now be
					generic GMT pens (e.g., 2/255/0/0ta), and new
					function GMT_setpen is used throughout GMT for
					setting pens.
		pscoast.c		New option -M to dump shorelines which makes the
					utility shoredump obsolete.
-----------------------------------------------------------------------------------------
===> 09/04/1998		tar 3.1brev2 made and released to testers
-----------------------------------------------------------------------------------------
09/05/1998	gmt_math_init.bash	Bash version for install under Cygwin32 or DJGPP
		job*.bash		All example scripts now come in csh and bash
					flavors.
		do_examples.bash	Driver script for all examples
09/07/1998	pslib_inc.h		Added 6 more patterns to make it an even 90.
		*.*			-Gp/-GP now takes dpi instead of size. 0
					still means current DPI setting.
09/09/1998	gmt_support.c		replaced pw_gecos with pw_name since the
					former is not guaranteed to be there by POSIX.
09/11/1998	minmax.c		Reported junk if -I was set and no files found.
09/16/1998	gmt_io.c		Added DOS portability for binary/text fopens.
		psscale.c		y-label now plotted using ANOT_FONT settings.
		psxy.c			Upper case symbol have equal size, lower case
					symbol act as in GMT 3.0.
		gmt*.c			Added DOS directory delimeter support and WIN32
					preprocessor flag.
		Appendix L		New info for using GMT on non-Unix machines.
09/17/1998	psxyz.c			Now enjoys the same 2-D symbols as psxy,
					and may read symbol types from data.
					data files may now be binary as in psxy
		psrose.c		Now can handle binary or ascii data.
		pshistogram.c		Now can do patterns.
09/23/1998	gmtdefaults		Added new parameter MAP_SCALE_HEIGHT for size of
					map scales on map (used to be TICK_LEN)
		gmt_math_win32.h	The gmt_math.h file for MSDOS
		gmt_win32.h		Extra include file for MSDOS compilation
09/26/1998	pshistogram.c		Fixed bug in -S.
		gmtselect.c		Fixed bug in memory allocations
		GMT			Added Examples 17-20 and updated cookbook
09/28/1998	grd*			Makes sure that grd files created by GMT have no
					more than 0.01% slop in dx/dy relative to data
					range.
					Likewise, require the same max slop upon reading
09/29/1998	grdedit.c		New option -A to adjust old grdfiles with too much
					slop in -R -I.
09/30/1998	spectrum1d.c		Allow for single window.
-----------------------------------------------------------------------------------------
===> 09/30/1998		tar 3.1brev3 made and released to testers
-----------------------------------------------------------------------------------------
10/08/1998	psxy.c			Will automatically look for -W and -G strings
					in multisegment headers and let those settings
					override the command line settings
		psxyz.c			Same as psxy.c above
10/09/1998	gmtdefaults.c		Needed GMT_stdout set separately since GMT_begin
					is not used in this program only.
		xyz2grd.c		Default for -Z should be ascii (-ZTLa)
		grd2xyz.c		Same.  Also updated man pages.
10/13/1998	grdview.c		Now checks for improper file sizes AFTER having
					the right information in hand...
		gmtenv.bat		DOS bat file to set GMT environment paths on
					Windows 95/1998.
		gmtinstall.bat		DOS bat file to compile GMT using Microsoft
					Visual C/C++ tools.
		letter2A4.sh		Supply shell script to convert documentation to A4
-----------------------------------------------------------------------------------------
===> 10/14/1998		tar 3.1brev4 made and released to testers
-----------------------------------------------------------------------------------------
10/15/1998	gmt_cdf.c		GMT_cdf_write_grd_info was not in define mode so
					changes to information was not allowed in grdedit.
		gmt.h			Fixed typos in ellipse parameters and added WGS-66
10/16/1998	blockm*.c		Incorrectly set n_expected_fields when -bi gave
					the number of fields.
10/18/1998	mapproject.c		Binary output did not happen
10/19/1998	grdtrack.c		Same problem for some cases of io.
		gmt_io.c		Ascii read function now chops off trailing blanks,
					commas, and tabs so # cols is set correctly.
		gmtselect.c		Simplified output format test.
		psclip.c		Now can handle empty multisegment segments.
		psxy.c			Now can handle empty multisegment segments.
		psxyz.c			Now can handle empty multisegment segments.
		spectrum1d.c		Now also may output using binary format with -bo[s]
		splitxyz.c		Now also may output using binary format with -bo[s]
		minmax.c		Must specify # of columns in -bi
		*.c			All programs that take -M (multisegment file input)
					now checks for empty segments.
		psmask.c		Added -F for pixel registrered masks.
		grdproject.c		Now sets measure unit to m with -Mm.
		grdcontour.c		-W+ will use the cpt colors when drawing contours
		pscontour.c		-W+ will use the cpt colors when drawing contours
-----------------------------------------------------------------------------------------
===> 10/20/1998		tar 3.1brev5 made and released to testers
-----------------------------------------------------------------------------------------
10/20/1998	sample1d.c		Fixed typos in usage message.
10/21/1998	block*.c		Finally, binary i/o works...
		*.c			Binary i/o only worked for 1st input file
					because GMT_io.status was never reset to 0.
10/22/1998	psxy[z].c		-L now ONLY means to close the polygons.
					Only -W[pen] sets symbol outline.
		job*			Edited all affected example scripts.
		pslib.c			If SI is set then A4 pagesize is assumed.
		letter2A4.sh		Adjust differently for landscape & portrait.
		Documentation		Fixes spellings and added more info.
		Scripts			Shrunk Ex 12 and 14 by 10%, added title to 17.
10/26/1998	pscoast.c		Now checks for all of the 5 coastline bases.
		grdlandmask.c		Same.
		*.c, *.txt		All -H explanations are now [-H[<nrec>]]
		pscoast.c		-M now accepts optional multisegment flag
		*.man			Cleaned up manpage typos and inconsistencies.
		Archives		Split crude, low, and intermediate coastlines
					and rivers/borders into separate _lib.tar.gz
					archive to make source archive small.
10/28/1998	GMT			Now under GNU General Public License.
		gmt_init.c		Updating of .gmtdefaults is now using
					advisory file locking.  This means that
					only one GMT program will be able to
					read and write to the .gmtcommands file
					at any time.  In GMT 3.0 and earlier
					versions it happened that a later GMT
					call in a series of pipes would try
					to read while another program was busy
					updating the file.
		*.c			Programs now update the .gmtcommands
					file immediately after error checking.
10/29/1998	blockmode.c		New program for block mode estimates.
		blockmedian.c		Now takes -Q to return (x,y) at median z
					value rather than median x, median y.
10/29/1998	psmask.c		When optimized the distance variable was
					not exactly zero when I thought it was,
					so a test failed that assumed it was zero.
10/29/1998	gmt_init.c		Permission problems with i/o to .gmtcommands
					are no longer fatal errors, just warnings.
		gmt_grdio.c		Now tells users to use grdedit -A upon
					finding incompatible -R -I in grdfiles.
10/30/1998	psxy[z].c		Needed to differentiate between number
					of required versus expected columns,
					since binary must read all columns.
		*.c			All programs that have -bi option are given
					reasonable default values for the number of
					columns that is expected.  This can always be
					overridden by -bi[s]<n>.
		job20.*			Assumed ./ in users' $path.
11/02/1998	grdvolume.c		Added -Z and -L options; also tighter error
					checking on arguments.
		grdcontour.c		Upper left node was not changed by -Z.
11/03/1998	project.c		Can handle any number of input columns beyond
					the required x, y.  Now behaves in the same way
					as mapproject and grdtrack in that text fields
					are passed through if both input and output
					modes are ASCII.
11/04/1998	psxyz.c			Used old, wrong variables in outside test.
11/05/1998	pslib.c			Only issue paint and pen changes when the new
					values differ from the current values.
		gmt_math_init.*		Now use isnanf = isnan[d]((double)(x)) instead
					of poor man's x != x macro.
11/08/1998	gmt.h			Special #ifdef __MACHTEN__ to fix bad machten
					implementation of ANSI/POSIX.
		pslib_inc.h		-------"-------
		gmthex2ras.c		-------"-------
		job09.csh		Made safe for broken csh versions that do not
					like if (..) action 1-liners.
		gmt_support.c		Added alternate GMT_delaunay function based on
					Jonathan Shewchuck's triangle routine.  The
					instructions on how to use this alternate code is
					placed in TRIANGLE.HOWTO in the src directory.
-----------------------------------------------------------------------------------------
===> 11/08/1998		tar 3.1 made and released via email
-----------------------------------------------------------------------------------------
11/09/1998	gmt_init.c		passed char rather than int to isdigit, isalpha
					and strchr.
		grdmath_main.c		passed char rather than int to isdigit
		pstext.c		Same.
		gmt_io.c		Passed char rather than int to strchr
		gmt_support.c		Arg to one free call not cast to void
					Several GMT_memory returns not cast to char *
		gmt_map.c		Arg to one free call not cast to void
		grdcontour.c		Several GMT_memory returns not cast to char *
		grdtrack.c		Same.
		mapproject.c		Same.  All these fixes are not bugs but bad style.
11/10/1998	Makefile		Replaced USER with PROGS to avoid conflicts
					with $USER in shells.  Added gmt_boundcond.h
					to list of include files in makefile.
11/13/1998	gmt_io.c		ascii_input now returns -1 if EOF instead of 0.
		gmtselect.c		do_project was not set if -J was not given but
					-R was used and -Jx1 was implied.
					check for n_fields >= 0 since -M files' multi-
					segment headers gives n_fields = 0
		grdtrack.c		Failed for -M files; now checks for n_fields >= 0
		mapproject.c		Same.
		pshistogram.c		Same.
		psmask.c		Same.
		spectrum1d.c		Same.
		trend?d.c		Same.
		psxy[z].c		-G -W in multiple segment headers now also apply
					to symbols.
11/16/1998	xyz2grd.c		Fixed Verbose output report to align names and
					values.
11/23/1998	psxy.c			Warns if error bars exceed allowable domain,
					i.e., goes negative for log plots.
		minmax.l		Said -K instead of -M (explain_-M.txt had typo)
		gmtconvert.l		Same.
		project.l		Same.
		psmask.l		Same.
		gmt_support.c		R. Krueger: Fix to yield the same pattern & line
					thickness regardless of SI or INCH.
		pslib.c			---"---
		pshistogram.c		Added -A option and new 'o' modifier for -I.
		pshistogram.l		Had no mention of -I.
		grdhisteq.c		Did not update grdfile history remarks.
-----------------------------------------------------------------------------------------
===> 11/23/1998		tar 3.1.1 made and released via email & EOS
-----------------------------------------------------------------------------------------
11/25/1998	project.c		Got hung if -W or -L was set.  Fixed many typos
					in both man page and usage message.
		psscale.l		Had old-style -B example.
		gmt.h			Added support for paper Media formats
		gmt_init.h		Same
		gmt_init.c		Same
		gmt_media_name.h	New include file for media names
		gmt_media_size.h	New include file for media sizes
		gmtmedia.d		User-editable file for custom paper media kept
					in $GMTHOME/lib
		gmt_keywords.h		Changed PAPER_WIDTH to PAPER_MEDIA
		gmtdefaults.l		Same, plus added list of all media
		pslib.c			Changed paper_width from double to integer array
					of length 2.  People usling pslib, beware!
		pslib.h			Same
		pslib_inc.h		Same
11/26/1998	grd2xyz.c		Now gives warning if -bo and -Z are selected.
		xyz2grd.c		Now gives warning if -bi and -Z are selected.
11/26/1998	pslib.c			Bug in ps_image caused 1-bit images to fail.
		pscoast.c		Now uses GMT_prep_polygons which is defined in
					gmt_shore.c as a new function.
		grdlandmask.c		Now uses GMT_prep_polygons instead of duplicate
					code from pscoast.  All progs needing to extract
					and use coastline polygons will use this routine.
11/28/1998	img2mercgrd.c		Supplemental program: Fixed bug in -T0 that gave
					0 data
11/29/1998	pslib.c			Read Sun rasters incorrectly on bigendian
					machines.  Now should work ok on all platforms.
		gmtselect.c		Now takes -N (and -A, D) to keep/skip points that
					are inside/outside any geographical features
					(coastlines).
11/30/1998	pslib.c			Added support for Run-Length Encoded rasterfiles.
					PS. this feature has not been extensively tested
					at this point since I haven't gotten many RLE
					files to test on.  Also fixed a bug for 24-bit
					images of odd width.
		gmthex2ras.c		Obsolete since Sun rasterfiles are defined at
					the byte level so new read-module in pslib works
					on all platforms.  The GMT patterns will now be
					distributed as binary rasterfiles.
12/01/1998	gmt_customio.c		Now reads/writes format # 3 Sun rasterfile header
					correctly on all platforms, like pslib.c
					If you ever used GMT to write Sun rasterfiles
					under Linux, you need to fix the raster headers.
12/02/1998	psclip.c		fclose outside file loop so only ~95 files could
					be openened before we run out of file pointers...
		gmtselect.c		Added friendly check to ensure -L and -F gets a
					multisegment file.
12/08/1998	grdimage.c		Had two problems: 1) Clip path did not follow
					the grid but the -R region. 2) images sometimes
					were shifted incorrectly.  Both problems fixed I hope.
		grdview.c		Clip path now based on grid, not -R.
					Pixel grids are now processed correctly.
12/09/1998	*			All internal calculations in GMT are now done
					using inches.  All input lengths (e.g., map
					scales or widths, distance gaps, symbols sizes,
					etc) that previously were taken as inches or as
					either inch/cm are now scanned to see if they
					have the unit c, i, m, or p appended (c = cm,
					i = inch, m = meter, p = point).  If no unit is
					appended, GMT will default to the unit set by
					MEASURE_UNIT.  Thus, even if MEASURE_UNIT is cm
					one can specify a 6 inch map with -JM6i.
					All the man pages have also been updated.
		psmegaplot.c		Updated to deal with GMT 3.0 PS format.  No
					longer need PS files produced with -K at the end.
		job??.{csh,bash}	Now employ unit-independent widths, sizes and
					pens.  Should give identical results regardless
					of which MEASURE_UNIT is in effect.
12/14/1998	gmt_plot.c		FANCY frame for Lambert Conic had problems
					when south was not multiple of y_inc.
12/29/1998	gmt_inic.c		Now accepts -X -Y where a prepended a indicates
					absolute coordinates [Default is relative,
					prepending r has no effect]
		pslib.c			Added absolute mode according to above.
		gmt_support.c		Same
01/02/1999	grdmath.c		Added j0, j1, jn, y0, y1, yn, erf, erfc
					which is present in most math libraries,
					and bei, ber, kei, ker from gmt_math.c
					Also log1p which we already use elsewhere.
		gmtmath.c		New program with same operators as grdmath
					expect it works on tables instead of grids.
		gmt_math.c		Contains definitions of bei, ber, kei, ker.
		gmt_math_init.csh/bash	Looks for j0, j1, jn, y0, y1, yn, erf, erfc
					and puts dummy macros in gmt_math.h if not found.
					(erf and erfc also conditionally included
					in gmt_stat.c if not found locally).
01/12/1999	grdmath.c		Added Legendre polynomical of degree L and order
					M (PLM), the fully normalized surface spherical
					harmonic YLM, and POP (delete top stack element).
		gmtmath.c		Added PLM and POP, see above.
		pstext.c		Now clips same way as psxy.
01/13/1999	gmt_map.c		Improved map-jump mechanism.
01/14/1999	pslib.c			-G with black failed in psbasemap.
01/15/1999	GMT			Added ZapfDingbats as font # 34
		psbasemap		Usage did not list -E
01/17/1999	pstext.c		-W can now take optional pen after 'o', and
					will not fill without a shade/color.
01/19/1999	Documentation		All rewritten in Latex and translated to
					both PostScript and HTML
01/20/1999	grd/gmtmath.c		Added modified Bessel I0 and modified
					Kelvin K0 functions.
01/21/1999	grdmath.c		Added 5 operators: DDX, DDY, D2DX2, D2DY2, CURV
		gmtmath.c		Added 2 operators: DDT and D2DT2
01/23/1999	psclip.c		The -M option wrote the multisegment header
					to the PS file without leading %.
		gmt_plot.c		Now can use grid crosses for linear projections
01/25/1999	psbasemap.c		Now works with -JX<width>p<power> when the
					<power> is negative.
		grdmath.c		Added modified Bessel functions I1, IN, K1, KN.
		gmtmath.c		Same.
02/04/1999	grdhisteq.c		Had error in testing filenames.
02/05/1999	gmtio.c			-M is allowed with -b and means look for records
					with all fields == NaN as to imply binary
					multiple segment headers.
		*.c			Pens are internally maintained in units of
					points to avoid DPI problems.
		mapproject.c		Added -D option to temporarily reset MEASURE_UNIT
02/08/1999	pslib.{c,man}		All symbols sizes are now given as the
					diameter of circumscribing circle.
		psxy[z].{c,man}		Same, except for upper case as before.
02/11/1999	psscale.c		Giving negative length now reverses scale
					Will read stdin if -C not given.
02/17/1999	pscoast.c		Now may take only -C and only plot lakes.
		gmt_map.c		Gives error if longitude range exceeds 360
03/02/1999	triangulate.c		Synopsis or man page did not mention -E option.
					Added -Dx|y option to take derivatives.
					-J was ignored when -G -I were set.
03/03/1999	grdcontour.c		-Z[fact[/shift]]p allows for 360-degree
					periodicity in z (i.e., phase data).
03/03/1999	psxyz.c			Did not let users specify -Su (3-D cubes)
03/08/1999	gmtmath.c		Added Heaviside step function STEP
		grdmath.c		Added Heaviside step functions STEPX and STEPY
03/11/1999	makecpt.c		Rewritten to take a master cptfile as starting
					point and -T option to set z-steps.  About a
					dozen master cpt files distributed with GMT.
		grd2cpt.c		Same.  Takes -C<mastercpt> and -Tz0/z1/dz as
					makecpt.c
03/16/1999	gmt_services.html	Added sub-page for supplmental man pages
03/17/1999	manpages		Fixed misc. typos and inconsistencies.
03/18/1999	splitxyz.c		z-output not passed through correctly.
03/22/1999	grdedit.c		Using -R now changes the dx/dy automatically.
03/23/1999	gmt_math_init.csh	More robust test for which awk to use
		do_examples.csh		Same.
-----------------------------------------------------------------------------------------
===> 03/25/1999		tar 3.2 made and released via email
-----------------------------------------------------------------------------------------
03/26/1999	gmtmath.c		Forgot to implement -H, now done
03/26/1999	pslib.c			Added ps_imagemask for transparent bitmaps
		gmt_support.c		Understands that fore- or background color in
					-Gp can be - which means transparency
		psimage.c		-Gf- or -Gb- means transparency
		grdimage.c		New -G[f|b]<rgb> option to make transparent
					color masks from otherwise black/white 1-bit
					images
04/02/1999	grdview.c		-Qsm drew contours even though -W was not set.
		pscontour.c		Usage did not explain -W+.  Added -A- to
					force deactivation of all contour labels
		grdcontour.c		Same.
		sample1d.c		Had memory allocation errors with -M.
					Now also works with -N file that spans the
					range of several -M segments.
		gmt_plot.c		GMT_polar_adjust sensitive to round-off which
					could affect orientation of basemap anotations
		gmt_math_init.csh/bash	Had wrong #define for the dummy yn(n,x) function
04/05/1999	grdreformat.c		Now skips comments in gmtformats.d file
		grdraster.c		Forgot to deal with DOS directory slash
					and drive letters
		gmt_mgg.c		Forgot to deal with DOS directory slash
04/06/1999	gmt_map.c		Hopefully fixed a bug in finding map crossings
					for linear projection.
04/07/1999	pslib.c			Now knows about DOS drive letters in any
					pattern filenames
		gmtmath.c		Added -Q shorthand option for scalar math
					Same as giving -Ca -N2/0 -T0/0/1
04/08/1999	pslib.c			Missed final linefeed when -K and absolute
					positioning (-Xa? -Ya?) were used.
		grdimage.c		Shit! Somehow, MAPPING was replaced by
					!RECT_GRATICULE on 12/8/98, i.e. after 3.1.1
					Now fixed (again). All grdimage plots made since
					3.2 came out have problems if they use a map
					projection that had rectangular graticules, like
					Mercator and many cylindrical projections.
04/09/1999	pslib.txt		ps_command and ps_colortiles are now documented
		pslib.c			Added ps_flush for buffer flushing.
		pscoast.c		Added -bo option for use with -M
04/12/1999	grdraster.c		Did not open files as binary (which is required
					under Win32)
		dat2gmt.c		Same.
		gmt2dat.c		Same.
		gmt2bin.c		Same.
		gmtinfo.c		Same.
		gmtselect.c		Sped up test for proximity to points by first
					checking for inside a square (fast), then the
					circle (slow).  Suggested by Joaquim Luis.
					Also added -Cf for using Flat earth (fast)
					rather than great circle distance calculations.
					Also, use original x,y values for the output
		psscale.c		Did not respond to any -X -Y options.
		gmt_init.c		Complained if one of -X, -Y was set as absolute
					and the other was implicitly 0.
		grdcut.c		Did not check that output -R was ok
		grdfilter.c		Did not check that output -R, -I were ok
04/27/1999	gmt_defaults.h		Without .gmtdefaults, paper size was 0x0
		install_gmt		Now figures out if binary files written on Suns
					(e.g., data files like world_grav etc) need to be
					swapped on the user's machine.  This is used to
					compile in swapping in some of the supplemental
					programs. Not used by main GMT programs.
		img2mercgrd.c		Compiles in swapping if install_gmt says so.
04/28/1999	xyz2grd.c		Added -S[file] option to simply swap the
					byte order and write out new table.
					Added H and I for unsigned short and ints
		grd2xyz.c		Added H and I for unsigned short and ints
		gmt_io.c		Fixed problems with swapping of signed
					short and ints.  Added i/o for unsigned
					short and int.
04/29/1999	blockm*.c		Inadvertently turned -W on if data had 4 cols.
		gmtdefaults		Added case-insensitive checking for most
					GMT parameters unless noted.  Also enforcing
					stricter error checking on given values
		gmt_map.c		Sped up several projections by eliminating
					repeated trig calls with same argument.
		gmt_support.c		GMT_LUT structure now holds the quantity
					i_dz = 1 / (high_z - low_z) which got
					recalculated every time in GMT_get_rgb24.
05/04/1999	pstext.c		Added -M mode for paragraph typesetting.
					-Dx/y[v[<pen>]] will draw line from point to
					the offset location (-M only).
		pslib.c			Added ps_words function for paragraph
					typesetting.  Will include file PSL_text.ps
					from lib until it is stable. ps_words knows about
					3 additional escape sequences:
					@_ toggles between normal and underlined text
					@:size: changes font size, @:: resets it.
					@;r/g/b; changes font color, @;; resets it.
05/06/1999	grdmath.c		Added function DILOG.
		gmtmath.c		Same.
05/11/1999	pslib.c			No longer issues %!PS statements more than once.
		gmt_map.c		Added spherical code for the projections Mercator,
					TM, UTM, Stereographic, Lambert Conic, Albers
					Conic, and Cassini.  Used when flattening == 0.
		psxy[z].c		Did not process -Sbb<base> correctly (when width
					is expected as one of the data columns).
05/12/1999	pslib.c			%%Orientation was not set correctly, and is now
					only issued for PS files, not EPS.  %%Page is now
					used for PS files.
		gmt_map.c		scale_factor was inadvertently set to 1 for the
					Stereographic projection.
05/18/1999	gmt_support.c		GMT_shift_grd had error for gridline-registered
					grids.  This affected grdedit -S operations.
		grdimage.c		Needlessly interpolated grid when -JX..d was used.
		xyz2grd.c		-L and xyz input will now enforce periodicity.
05/20/1999	GMT			Added -Jv (van der Grinten) projection
					Added -Jd (Equidistant conic) projection
					Added -Jkf Eckert IV (f for 4).  Also allow
					-Jks for Eckert VI (s for 6) but -Jk without
					modifier will default to Eckert VI.
05/25/1999	gmt_plot.c		Now FANCY basemap works for all Conic projections
05/26/1999	gmt_plot.c		Implemented global TM, and allow optional way to
					specify TM with latitude of origin (default 0).
05/27/1999	xyz2grd.c		Was not robust enough when calculating the node
					positions.  Warns about pixel/gridline confusion.
05/28/1999	grdproject.c		-F now means toggle node registration,
					and default is same as input file.
05/31/1999	xyz2grd.c		Test for -N nodata values mixed floats
					and doubles and would sometimes fail.
		gmt_init.c		More tolerant of blank lines in .gmtdefaults
					file and comments in ellipse files.
		grdraster.c		Now tolerates blank lines in grdraster.info file
06/04/1999	gmt_plot.c		Round-off problems would lead to some anotations
					being ommitted at xmax or ymax
		*.c			#ifdef SINCOS added to take advantage of sincos
					function where available.
06/07/1999	gmt_init.c		Make sure -b options are processed before file
					names are encountered since fopen mode must be
					set correctly. Should only affect WIN32 version
		examples		Now include DOS *.bat versions for all examples
					(except 03).  Added scripts to view all plots.
06/08/1999	gmt_init.c		flush .gmtcommands file before unlocking/closing
					and use same open file pointer for read/write.
-----------------------------------------------------------------------------------------
===> 06/09/1999		tar 3.3 made and released via email
-----------------------------------------------------------------------------------------
06/11/1999	xyz2grd.c		When -S is set the fopen modes must be
					rb and wb or it will fail under WIN32.
06/14/1999	psxy.c			Needed a "stroke" operator after drawing
					error bars.
		psxy.txt		Clarify that -N does not affect the PostScript
					BoundingBox; you may have to post-process in
					order to get a correct BB.
		psxyz.txt		Same.
		pstext.txt		Same.
		gmt_support.c		CPT files with gray-value columns only did not
					get processed correctly.
06/15/1999	gmt_map.c		Mindless insistance on |y| <= 90 applied even
					for non-geographic data in gridlines...
		grdraster.c		Added optional final column that can be set to L
					or B to indicate the endianness of the data file.
06/17/1999	gmt_init.c		Use GMT_DEFAULT_PATH instead of /usr/local/gmt.
					This is set in makegmt.macros
		pslib.c			Same.
06/18/1999	gmt_cdf.c		Reduced repeated code and simplified
		gmt_customio.c		Same.
		gmt_grdio.c		Same.  Combined, these 3 shrank in size by 30%
06/20/1999	psxy[z].c		-SV[]n<norm> did not work, now fixed.
		psxy.c			-E with no symbols gave SEGV.
					Now, -E with no xy means -Exy (or -Eyx, with -:)
06/21/1999	pstext.c		Better treatment of European font reencodeing.
		*.c *.h			Added MSC support for DLL libraries (WIN32)
		-B			Can also take :,"unit string": which will be
					appended to each anotation value.
06/21/1999	gmt_stat.c		Speed up calculation of kei, ker, bei, ber and
					extended domain to |x| > 8.
					Fixed m-loop bug in plm.
06/24/1999	fourt.c			Erroneously included sunmath.h instead of giving
					prototype.
		mgd77togmt.c		Now understands the new, Y2K-compliant MGD77
					format and also has a kludge Y2K fix for older
					(i.e., 2-digit year) files.
		spotter			New supplemental package for plate tectonics
		x2sys			New supplemental package for finding trackline
					intersections (crossovers)
		meca			Four new tools (pscoupe, psmeca, pspolar, and
					psvelo) replaces old psvelomeca and adds more
					functionality.
		Makefile		Now supports shared libraries on most systems
06/30/1999	psxy[z].c		-SE would scale axes incorrectly if MEASURE_UNIT
					did not happen to be INCH.
		gmt_map.c		GMT_left_albers, GMT_right_albers replaced with
					GMT_left_conic, GMT_right_conic that actually does
					the work (no more report to gmt@soest message).
					Same for Lambert and equidistant conic.
		gmttrack.c		Messed up measurement units if not INCH, now ok.
		grdview.c		-Qs did not process F/B colors correctly.
07/01/1999	PSL_text.ps		Typo PLS_top instead of PSL_top caused error
					for 1-line paragraphs.
-----------------------------------------------------------------------------------------
===> 07/01/1999		tar 3.3.1 made and released via email
-----------------------------------------------------------------------------------------

3. Known Problems that have not been fixed yet.

	* Still no explicit implementation of time-axes. Scheduled for 3.4
	* Filling of polygons that covers 360 degrees of longitude may not always work
	  in psxy.  Split them into pieces with < 180 degree longitude span.
