UN ZOOM DE DEUXIÈME GÉNÉRATION POUR LMD-Z

Robert Sadourny et Le Van Phu

Laboratoire de Météorologie Dynamique

École Normale Supérieure

Mars 2001


  1. Introduction
Le zoom (ou étirement de coordonnée) du modèle LMD-Z est défini séparément dans chacune des deux directions, la longitude  et la latitude ,    par  quatre  paramètres :
  1.  L'abscisse du centre du zoom ou  (en radians) ;
  2.  Le grossissement au centre du zoom (par rapport à une grille régulière de même nombre de points) ;
  3.  La largeur de la fenêtre du zoom (en radians) ;
  4.  La raideur de la transition de l'intérieur à l'extérieur du zoom
.

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.
 

  1. Le zoom en longitude

  2. 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)
     
     

  3. Le zoom en latitude

  4. 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)
     

  5. Contrôle des paramètres
Dans la pratique, on se fixe une limite considérée a priori comme acceptable pour la dégradation de la grille hors de la fenêtre :  ("grossissement" compris entre 0 et 1, typiquement de l'ordre de 0,2 ou 0,1). C'est un cinquième paramètre du programme. Le logiciel de construction du zoom teste la condition

.   (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 modérer le grossissement (diminuer ), Ou restreindre la fenêtre du zoom (diminuer ), Ou accentuer la raideur (augmenter )

Ou tolérer plus de dégradation hors du zoom (diminuer )
 

                                          Voir  Graphiques ?                   oui  !
                                                                                         non !