Robert Sadourny et Le Van Phu
Laboratoire de Météorologie Dynamique
École Normale Supérieure
Mars 2001
La fonction d'étirement est basée sur des raccordements de tangentes hyperboliques, ce qui permet d'obtenir des résolutions relativement uniformes dans le zoom et à l'extérieur.
On notera que
et
ne sont pas indépendants,
car il ne faut pas dégrader exagérément la grille
à l'extérieur du zoom : cette dégradation (ou rapport
de la maille avec zoom à la maille sans zoom) est de l'ordre de
en
longitude, ou
en latitude.
Une limite raisonnable consiste à ne pas accepter plus de la moitié
des points à l'intérieur du zoom :
,
(1)
soit une dégradation maximale de l'ordre de
à l'extérieur. La condition (1) n'est suffisante que
pour une raideur infinie. On verra sur des exemples qu'il faut être
plus contraignant si l'on use de grossissements forts et de raideurs faibles.
On notera aussi que la transition entre le zoom et
l'extérieur est d'autant plus rapide que
est grand. Une raideur moyenne est de l'ordre de 1. Les effets
de distorsion dus aux variations de la résolution à travers
les bords de la fenêtre zoomée seront d'autant plus marqués
que la raideur choisie sera grande.
Pour ne pas dégrader l'ordre des approximations, la structure du zoom sera définie par des applications analytiques. Pour obtenir des résolutions quasi-uniformes à l'intérieur et à l'extérieur de la fenêtre du zoom (en dehors des zones de transition), le grossissement sera construit par des recollements de fonctions tangente hyperbolique, l'intérieur et l'extérieur du zoom correspondant aux asymptotes où cette fonction varie peu. Nous définirons
(
),
,
. (2)
Pour faciliter l'écriture, nous userons aussi
de la fonction de Heavyside
(définie par
pour
,
et par
pour
)
pour formaliser la façon dont on recadre des intervalles ou dont
on réordonne des indices.
Nous choisissons pour valuation de la longitude
de référence le domaine
:
l'origine des longitudes est ainsi située au centre de la grille.
Commençons par définir une longitude
recentrée au centre du zoom :
,
.
Elle sera liée à la longitude initiale par
,
(3)
.
(4)
On se donne alors la dérivée (ou grossissement
local) de la coordonnée de travail
,
variant elle aussi sur (
) et
centrée au centre du zoom :
,
(5)
avec
,
(6)
pour
.
La condition que, sur cet intervalle,
lui-même doit varier de 0 à
entraîne
,
(7)
d'où l'expression de
:
.
(8)
De
on déduit
de 0 à
sur une grille de travail régulière
suffisamment fine (de l'ordre de 30 000 points) pour minimiser l'erreur
dans les calculs suivants. On antisymétrise ensuite la fonction
pour traiter l'intervalle
:
ou
. (9)
La coordonnée de travail
varie ainsi de
à
comme
.
Enfin, on se donne la grille régulière du modèle dans l'espace de travail centré sur le zoom :
,
, (10)
et l'on résout
(11)
pour obtenir les longitudes des points de la grille zoomée :
,
(12)
(13)
Il reste à réordonner ces points dans
l'intervalle
:
on définira donc
tel que
,
(14)
qui nous amène au champ réordonné
,
(15)
.
(16)
Le facteur d'échelle associé à la grille sera alors
.
(17)
Le problème est ici légèrement
différent à cause de l'asymétrie qui résulte
d'un zoom non centré sur l'équateur, et aussi de l'absence
de valuation multiple pour la latitude. Désignons par
la latitude du centre du zoom. La formulation du zoom latitudinal va dépendre
de la position du centre par rapport à l'équateur. Nous le
définirons comme précédemment, mais sur un domaine
,
symétrique par rapport au centre du zoom, s'étendant donc
au delà du pôle de l'hémisphère où est
centré le zoom.
On définira alors le grossissement local par l'analogue de (5)
,
(18)
avec
,
, (19)
,
. (20)
La condition à imposer est que
varie de
à
quand
varie de
à
:
,
(21)
d'où l'expression de
:
.
(22)
De
on déduit
sur une
grille de travail régulière
suffisamment fine (de l'ordre de 30 000 points) pour apporter le moins
d'erreur possible dans les calculs qui suivent.
On se donne alors la grille régulière du modèle dans l'espace de travail, allant du pôle sud au pôle nord inclus :
,
, (23)
où
correspond à la résolution méridienne du modèle
en nombre d'intervalles, et l'on résout
(24)
pour obtenir les latitudes
des points de la grille zoomée.
Le facteur d'échelle associé sera alors
.
(25)
.
(26)
Si cette condition n'est pas satisfaite, un message d'erreur apparaît :
Votre choix de paramètres n'est pas acceptable. Vous devez
Ou tolérer plus de dégradation hors
du zoom (diminuer
)