Description of data files:

The coefficients in the tables are all time rate of change in normalized Stokes 
coefficients in units of 10^-11 per years.  The normalization convention is as 
in James and Ivins [1998] or Ivins et al., [1993].  In other words, these are 
the same conventions as used by Lambeck or Chao.

Each data file has at the very top:

viscosity factors for lower mantle x 10^21 Pa S
    1.0000      1.5000      2.5000      3.5000      5.0000      7.0000
   10.0000     15.0000     25.0000     40.0000     70.0000    100.0000

which gives the 12 values (in increasing order) for the lower mantle viscosity 
variation (i.e., the "7th" value is 10.0000x10^21 Pa S and so on). After this 
top set of numbers comes the coefficients:

L  =       2
 cosine terms

   m
   0      0.7957      1.3200      1.9559      2.2809      2.6025      2.9788
          3.4701      3.9693      4.0960      3.6082      2.6577      2.0719
   1      0.0595      0.1085      0.1737      0.2062      0.2329      0.2627
          0.3080      0.3641      0.3925      0.3548      0.2649      0.2070
   2     -0.0733     -0.1460     -0.2500     -0.3019     -0.3396     -0.3784
         -0.4444     -0.5384     -0.6006     -0.5534     -0.4174     -0.3267

 sine terms
   m
   0      0.0000      0.0000      0.0000      0.0000      0.0000      0.0000
          0.0000      0.0000      0.0000      0.0000      0.0000      0.0000
   1     -0.1097     -0.2617     -0.4884     -0.5931     -0.6523     -0.7132
         -0.8511     -1.0736     -1.2321     -1.1385     -0.8538     -0.6662
   2     -0.0102     -0.0124     -0.0091     -0.0057     -0.0055     -0.0094
         -0.0130     -0.0090      0.0054      0.0151      0.0162      0.0135

and then L = 3... and so on.  

So, for example,

   m
   0      0.7957      1.3200   #### #### ...

corresponds to the first two viscosity calculations, eta = 1x10^21 and
eta = 1.5x10^21, and so on, twelve values in all (So there are actually
12 coefficient sets here!).