Offline nesting

[andres]

Hi,

I've read there are tools to create the boundary conditions for a small domain using a previous run of ROMS of a larger domain.

I've used nestgui from ROMSTOOLS to create the smaller domain and the initial and forcin files, but seems I can use it to extract the information from the his or avg file.

Thanks,

[mathieu]

This is a classic operation.
Everybody has different scripts. Mine are available from http://www.liga.ens.fr/~dutour/Roms_Interpolation/.
The idea is to interpolate from the history file (or from stations
if prefered) to the boundary points and then create a boundary forcing
file.

[andres]

This is a classic operation.
Everybody has different scripts. Mine are available from http://www.liga.ens.fr/~dutour/Roms_Interpolation/.
The idea is to interpolate from the history file (or from stations
if prefered) to the boundary points and then create a boundary forcing
file.

Thank you for your code. In my attempt to use it I got the following error

Using method 3
Reduce from nbNode=304515 to nbNodeRed=30744
nearest interpolation done
nbTrig=606816
ansatz success : 43599 cases
direct success : 0 cases
failures : 26581 cases
Sorry, but you are trying to interpolate on wrong place

But I am shure the nested file is a subdomain.

A couple of other questions:

What prefix do I use if I have only one his file?

I am using a climatological run and my fist element in scrum_time is 2592000 (seconds, i.e. 30 days)

What should I put in BeginTime EndTime?

Regards

[nganju]

I used Mathieu's toolbox for interpolation, but found a few bugs where the main function (e.g. NEST_CreateNestingTotalArray) calls other routines (e.g. InterpolGetSpMat...) with a different set of arguments than the routine is expecting. I think there were about 5 of these I had to correct. Once these are fixed, however, Mathieu's code for 2D and 3D fields performs as well as the method I was using before for 2D fields only, with over a 90% reduction in computational time.