Multi-threaded optimal interpolation module for octave

[Alexander Barth]

Hi all,

I wrote a Fortran 90 module performing a n-dimensional local optimal 
interpolation. The optimal interpolation allows to interpolate 
arbitrarily located observations to a regular grid using a background 
field as first guess. It can be used as an alternative to griddata. The 
optimal interpolated module can be called from Octave using a C++ wrapper.

I was able to create a working oct-file with gfortran, g95, Intel 
Fortran Compiler (free as in beer) and PGI Compiler (commercial) . 
Unfortunately, gfortran version 4.0 doesn't work but version 4.1 works fine.

It is also possible to use multiple CPUs during the optimal 
interpolation in octave if the Fortran code is compiled with OpenMP 
support enabled. This works with the Intel and PGI compilers. OpenMP 
support for GCC is planned for version 4.2.

Instructions and the source code can be found at:
http://ocgmod1.marine.usf.edu/OI

Suggestions are very welcome.

Cheers
Alex

--
_______________________________________________________________

 Alexander Barth

 Ocean Circulation Group
 University of South Florida
 College of Marine Science
 140 Seventh Avenue South 
 St. Petersburg, Florida  33701
 USA

 Phone: +1-727-553-3508     FAX:   +1-727-553-1189

_______________________________________________________________