The ``testinterp'' test-case can be run to test the interpolation of different source fields corresponding to analytical functions and to evaluate the error between the interpolated fields and the same analytical functions calculated on the different target grids.
All files needed to run this test-case can be found in oasis3/examples/testinterp/input and oasis3/examples/testinterp/restart.
To run ``testinterp'', OASIS3 first has to be compiled (see section 8.1.1) in interpolator-only mode NONE, i.e. by putting CHAN = NONE in the TopMakefileOasis3 header file. Then the programs that will calculate the interpolation error, i.e. gen_error.f90 and gen_error_vector.f90 (for vector fields) in directory oasis3/examples/testinterp/error have to be compiled (see script sc_comp_error).
Then, one has to adapt and execute the running script oasis3/examples/testinterp/
sc_run_testinterp.
With TIME=ONE, the configuration file oasis3/examples/testinterp/
input/namcouple_ONE, the input files
flda1.nc, flda2.nc, flda3.nc, fldb1.nc,
fldo1.nc
and fldz1.nc from oasis3/examples/testinterp/restart and the
input files aaIin.nc and caJin.nc from oasis3/examples/testinterp/restart/vector are used. This example also
shows one vector interpolation (field components a_at42_I and
c_at42_J). The test-case automatically writes the error
fields in error_*.nc files and error statistics in log_*
files.
To run the example into which OASIS3 interpolates many time
occurrences from one input file, put TIME=MANY in sc_run_testinterp. The configuration file oasis3/examples/testinterp/
input/namcouple_MANY and the input
file fldin.nc in oasis3/examples/testinterp/
restart are then
used.
The results obtained after running the testinterp test-case should
match the ones in oasis3/examples/
testinterp/outdata.
(XXX this directory is not updated yet).