This directory contains the files which were used to do the inversion and make the figures for the paper "Approximation and inversion of a complex meteorological system via local linear filters" by Schoenberg, Berk, Fovell, Li, Lu, and Weiss (2000).

The files are designed for the program R, which is a free unix version of S-plus. The files will thus also work in Splus.

The file getdata.s will not be useful for most users. It just documents commands we used to get the data, and describes the datasets briefly.

The files fig1to2.s and fig3to7.s were used to make figures 1 through 2 and 3 through 7, respectively.

Then the file findh.s was used to construct the matrix H. Its inverse is computed using invert.s . The inverse is applied to obtain input estmates and construct figures 8 and 9 in the file applyinverse.s .

Any questions can be addressed to Rick Schoenberg at frederic@stat.ucla.edu.