L'interpolation verticale est donc, dans le cas général, un spline cubique.
Les routines utilisées sont issues du NUMERICAL RECIPES IN FORTRAN 773.1 avec une modification
dans la recherche des points. En effet, les grilles verticales utilisées ne varient pas et sont
monotones. Il est donc facile de calculer, au préalable, les points sources entourant le point cible.
Les conditions limites sont des splines naturels. Les dérivées secondes sont donc nulles.
La valeur des points cibles en surface est fixée à celle des points sources du premier niveaux.
Le spline est utilisé pour extrapoler les points en profondeur si nécessaire.
Nicolas Daget
2006-03-21