Changes between Version 2 and Version 3 of IPSLCM6/Coupling/BugRunOff
- Timestamp:
- 07/18/17 11:47:24 (7 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
IPSLCM6/Coupling/BugRunOff
v2 v3 19 19 La construction des poids pour le run-off se fait en plusieurs étapes. 20 20 21 1. Limitation à la côte 21 ==== 1. Limitation à la côte ==== 22 22 23 23 On détermine tout d’abord les points côtiers. Pour l’atmosphère, ce sont les points qui ont une sous-surface océanique non nulle. Pour l’océan, les points qui ont un voisin masqué. … … 25 25 Les liens qui relient des points non côtiers sont supprimés. 26 26 27 2. Normalisation 27 ==== 2. Normalisation ==== 28 28 29 29 On procède ensuite à une normalisation des poids restant pour que tout le flux d’une maille atmosphère arrive bien dans une maille d’océan. 30 30 31 3. Extension 31 ==== 3. Extension ==== 32 32 33 33 ORCHIDEE construit son schéma de routage à partir d’un fichier de rivière à résolution élevé. Ce fichier contient implicitement une côte, et ORCHIDEE va router l’eau jusqu’à cette côte. Dans LMDZ-ORCHIDEE, le trait de côte est celui de ORCA. Certains estuaires profonds comme l’Amazone ou l’Amour ne sont pas représentés. Dans ces estuaires, ORCHIDEE ne route pas l’eau jusqu’à un point côtier. Il faut donc traiter ces soucis. Lors-qu’un point LMDZ est proche de la côte (paramètre dist_atm=400 km), on cherche le point côtier le plus proche, et on crée des liens d’interpolation qui vont de cette maille vers les points océans qui sont sous cette maille côtière. Une phase de normalisation permet la conservation. 34 34 35 La conservation globale 35 == La conservation globale == 36 36 37 37 Il reste des points de l’atmosphère qui sont loin de la côte, et dont le run-off n’est routé nul part. Et il y a sur le continent des points de run-off endoréiques comme le Lac Tchad par exemple. Il faut fermer le bilan en eau du modèle, et empêcher ainsi une dérive du niveau marin. On utilise pour cela une fonctionnalité d’OASIS : la conservation globale d’un champ. OASIS calcule l’intégrale du flux sur la grille d’origine (atmosphère) et sur la grille de destination (océan). Il calcule la différence, et corrige ensuite le flux sur la grille de des-tination. La correction est soit additive (ajout d’une valeur uniforme), soit multiplicative (multiplication du champ par une constante). IPSLCM6 utilise une correction additive. … … 39 39 Cette correction globale empêche de diagnostiquer de façon simple les erreurs sur les poids : la conservation est toujours respectées … 40 40 41 Les différentes versions, et leurs bugs 41 == Les différentes versions, et leurs bugs == 42 42 43 Version 1 43 === Version 1 === 44 44 45 45 IPSLCM5 utilisait le coupleur OASIS 3. Cette version du coupleur ne gère pas de façon native l’orientation sud-nord de la grille LMDZ. Il y a en conséquence des manips dans MOSAIC pour retourner la grille, et générer des fichiers de poids acceptés par OASIS. … … 47 47 IPSLCM6 utilise la version OASIS-MCT qui est plus versatile. Il a fallu revenir en arrière. Un paramètre supplémentaire dans MOSAIC permet de gérer le besoin de retourner la grille. On utilise une valeur différente pour générer des poids pour le ‘vieux’ IPSLCM5 utili-sé par les paléoclimatologues et pour IPSLCM6. 48 48 49 Le bug 49 ==== Le bug ==== 50 50 51 51 Il y avait une erreur dans MOSAIC, et un masque atmosphère était resté dans la mauvaise orientation sud-nord. On appliquait un faux masque terre-océan (inversé nord-sud). Cela masquait certaines rivières, dont l’Amazone. 52 52 53 Version 2 53 === Version 2 === 54 54 55 55 Pour la version 2, on corrige le sens du masque. 56 56 57 Le bug 57 ==== Le bug ==== 58 58 59 59 Pour calculer les poids du run-off, il faut commencer par relire les poids des flux atmos-phère vers océan. Il y avait un bug de lecture dans cette version du code : les poids étaient à zéro. Cependant la phase d’extension fonctionne normalement. On arrive à ce paradoxe : les points à l’intérieur de la terre sont routés vers l’océan, mais pas les points côtiers (Figure 2). L’Amazone est donc correctement routée, mais pas le Nil par exemple !. … … 61 61 Figure 2 : points atmosphère dont le run-off est routé directement à la mer. A gauche pour la ver-sion 1, à droite pour la version 2. Les versions 3 et 4 ont les mêmes points routés que la version 1. 62 62 63 Version 3 63 === Version 3 === 64 64 65 65 La version 3 n’a, pour l’instant, pas de bug connu. Les points atmosphères routés vers la mer sont les mêmes que pour la version 1.