Changes in gstat 1.9j, relative to 1.9g:

+  rewrote manual mostly
+  rewrote gstat home page, added gif's to examples
+  allow common parameters across variables (e.g. standardized cokriging),
   see the command ``merge''
+  added stratified universal kriging with user-specified base functions.
+  allows a more free variogram syntax, e.g.
    -1nug() - 1sph(30); or: - 1 nug() -1sph(30); etc.
+  calculate, report or plot kriging weights
+  reorganized nonlinear fitting procedure
+  error on absence of Linear Model of Coregionalization, except when
   `set nocheck = 1' is set
+  calculate block average coordinate polynomials if needed
+  added writing of debug information to log file 
+  added command line options -p plot_file (e.g. kriging weights plotting),
   and -l logfile; removed command line options -v, -b and -n
+  user interface:
   allow changing individual variogram parameters
   'p' prints to gnuplot plot file, term externally controlled (GNUTERM).
+  support GIF variogram plotting through gnuplot
+  added `make maputils', map utility/conversion programs
+  added map2gif.c (needs gd 1.2 gif library, from http://www.boutell.com/gd/)
+  replaced gram_schmidt orthogonalization with condition number estimation;
   now, if set, the condition number is checked against a maximum to detect
   singularities. In case of a singular covariance matrix, the prediction
   location is skipped.
+  data.c: allow one or more empty lines at end of column or GeoEAS file.

Bugs removed:
+  removed bug in standardized coordinate polynomial calculation 
+  removed bug if regression models were defined for multiple variables
+  several minor bugs