Version 3
Gérard Pellerin - CMC
Bernard Dugas - RPN
FÉVRIER 1994
Environnement Canada
Les travaux de B.Reid et de R.Benoit pour l'extraction des diagnostics zonaux forment la base d'une revision de code qui a été rendue nécessaire suite à la conversion au nouvel ordinateur NEC. La routine d'extraction des variables n'a pas change de nom, ni d'arguments. Elle a cependant été réécrite en entier et est maintenant indépendante de l'extraction des diagnostics ponctuels.
Le sous-programme d'extraction des diagnostics (MZONXST) n'a pas changé de nom, ni d'arguments, par contre il a été réécrit en entier pour permettre une plus grande flexibilité dans la spécification des zones (pas seulement des bandes de latitude). Une variable qui est dans la liste d'extraction, sera moyennée ou échantillonnée au cours de l'intégration (selon une pondération qui dépend des zones) et stockée dans des vecteurs d'accumulation. Ces vecteurs sont écrits dans le fichier NOUTZON à chaque pas de temps, dans l'éventualité d'une reprise, mais ne sont sauvés que si le pas de temps correspond à la période spécifiée par NPTEMPS, i.e. sauvegarde aux 24 heures. Ce fichier contient aussi toutes les informations pertinentes au modèle (longueur du pas de temps, dimensions de la grille, nombre de niveaux, etc). C'est de cette table qu'est régie l'extraction des diagnostics. La définition des zones (bins) et les pondérations corespondantes (poids) seront aussi écrites au préalable dans ce fichier. De même, pour un modèle où il y a une rotation par rapport au méridien de Greenwhich, comme dans le modèle régional, les angles sont calculés pour l'extraction correcte des composantes de vent (MVZNXST). Aussi, une section de code à été introduite (MZONOPR) pour assurer une séquence correcte des extractions.
Le problème principal pour les diagnostics zonaux dans le modèle régional, est l'approximation des cercles de latitude par des aires rectangulaires de la grille. Etant donnée la configuration variable de la grille, il faut choisir la largeur des bandes de facon à ce que le nombre de noeuds par bande soit suffisant pour obtenir une moyenne représentative du champ. Cette sélection est effectuée dans la routine LATCAL qui prend comme largeur de bande la plus grande distance entre deux noeuds. Les poids et rangs sont calculés dans ZEFRPRT. Une attention spéciale doit être portée lors de l'extraction des variables vectorielles. Pour le modèle spectral, chaque point de grille sur une même latitude (bande) recoit le même poids. Aussi, étant donné que ce modèle emploie les coordonnées géographiques, la rotation des vecteurs n'est pas un problème. Par contre pour le modèle regional, on tient compte de la rotation de la grille en définissant des angles de rotation pour chaque points de grille. Ce sont le sinus et le cosinus des angles qui sont écrits pour accélérer les extractions. L'initialisation des calculs pour les diagnostics se fait dans ZONINIT. C'est dans ce sous-programme qu'on appelle LATCAL et ZEFRPRT. Il y a aussi un équivalent de ZONINIT dans le modèle global, c'est dans ZONEPDS que sont calculés les poids et rangs de chaque point de grille.
Les enregistrements qui contiennent plusieurs niveaux peuvent être visualisés avec un utilitaire d'affichage comme REC. Par ailleurs, le programme GRAPHZN à été concu spécialement pour l'affichage des diagnostics zonaux. Le progiciel ZONTOOL offre une variété de menu permettant de mettre en évidence différents aspects d'une intégration.
dbzono -ISZONAL noutzon -OMZONAL zonfile -ECHOZ oui
La clef ECHOZ activée permet d'imprimer plus d'information, entre autre la moyenne des zones par niveaux et la moyenne totale. Le programme GRAPHZN requiet des niveaux du modèle pour l'affichage graphique, c'est pourquoi cette information est véhiculée dans les fichiers depuis le début. Le fichier NOUTZON n'est plus requis après la création de ZONFILE.
Ces deux clefs sont suivies d'une liste de variables ayant quatre caractères servant à définir les noms. La structure est la suivante:
Plusieurs appels à mzonxst sont insérés dans le code des programmathèques physique et dynamique. La liste des noms de variable est continuellement sollicité pour voir si les moyennes sont extraites ou non. Quand NPTEMPS=0, une clef de contrôle est activée au départ, cette clef impose une sortie sans sauvegarde. Cette clef sert aussi pour l'extraction ou non des variables vectorielles.
En exemple, dans la liste des variables de surface, '+VE.' indique que la variable VE sera accumulée sur NPTEMPS. Dans la liste des variables de profil, UU devra être immédiatement suivie de la variable VV, de même que dans l'extraction via un appel à MVZNXST. Le caractère de réserve pourrait servir éventuellement à spécifier qu'une variable a une dépendance à la coordonnée verticale ou autre, selon le besoin des usagers.
La première clef (+1) indique que le fichier doit être ouvert pour lire la table d'information, la liste des variables ou autres et allouer l'espace de travail requis pour le calcul des moyennes. La deuxième (+2) indique qu'il s'agit du début d'un pas de temps (lecture ou non des diagnostics dans le cas d'une reprise). La troisième (+3) indique qu'il s'agit du début d'une tranche verticale et la quatrième (+4) permet de spécifier le nombre de niveaux à extraire. Les mêmes clefs doivent être appelées dans l'ordre inverse, selon qu'il s'agisse de la fin d'une tranche (-3); les sommes sont effectuées pour celle-ci, de la fin d'un pas de temps (-2); les sommes et la table d'information sont écrites ou non, ou de la fin d'une intégration (-1); les mémoires sont désallouées et le fichier NOUTZON est fermé. Le numéro du fichier, le pas de temps et la tranche sont transmis via le deuxième argument de la routine avec le numéro approprié.
Le modèle peut aussi être intégré par tranches qui peuvent être plus courtes que la période spécifiée pour le moyennage. La table d'information contient un paramètre de contrôle qui en tient compte.
Tout d'abord, examinons le contenu du fichier obtenu pendant l'intégration. Le tableau I montre le contenu du fichier NOUTZON. Il y a 16 enregistrements dans le fichier,Le premier contient la table d'informations (T/) comprenant les variables pertinentes au système d'extraction. Les deux suivants contiennent la liste (V/) et la position respective (P/) de chaque variable dans le vecteur d'accumulation, les poids relatifs(W/) de chaque point de grille, les rangs (B/) son appartenance à chaque bins, le sinus (S/) et le cosinus (C/) de chaque point, qui correspondent à une rotation ou non de la grille. Les niveaux du modèle (S^) et la position des zones le long d'un méridien (L^) nécessaires au programme de graphisme. Notons que ces deux derniers portent l'étiquette de l'expérience tandis que les autres portent l'étiquette CONTROLE. Le fichier contient aussi les vecteurs d'accumulation (1/) par période de 12 heures et portent l'étiquette MZONXST.
Le travail du déballeur n'est pas très compliqué. Il s'agit de lire la table d'informations pour y extraire les dimensions pertinentes aux diagnostics, de lire la liste des variables, les vecteurs d'accumulations et de réécrire l'information par variables pour chaque période extraite du modèle dans le fichier ZONFILE. Le Tableau II montre le contenu du fichier pour une expérience du modele SEF. Il contient les enregistrements par variable pour des périodes de 12 heures. Chaque enregistrement porte l'étiquette de l'expérience SEFEXP33 et c'est ce fichier qui a servi à produire les graphiques de l'Appendice 1.
Examinons maintenant le contenu de la table d'informations pour une expérience du modèle régional obtenue de l'intégration d'un cas d'été débutant à 12Z le 23 juin 1993. Cette information sera imprimée lors de l'exécution du déballeur.
LITZON - TABLE DE CONTROLE
==========================
NDELTAT= 30
DELTAT= 720
MODE= 3
NI= 225
NJ= 175
NK= 25
NBIN= 15
SOMNK= 670
LATMIN= 0
ROT= 21
Il s'agit d'une intégration du modèle EFR, dont les dimensions de la grille sont 225 par 175 points avec un pas de temps de 720 secondes. Il y a 15 bandes de latitudes partant de l'équateur qui ont été extraites à tous les 30 pas, i.e. à tous les 6 heures. Les moyennes et le carré des moyennes ont été sauvés (MODE=3). Le programme lit ensuite la liste des variables qui est imprimée (Tableau III) ainsi que leurs positions respectives dans les vecteurs d'accummulations. Le chiffre dans la deuxième colonne indique s'il s'agit d'une variable moyennée (0) ou instantannée (1). La position relative dans les vecteurs d'accumulations est précalculée à l'entrée du modèle. Ce tableau est aussi imprimé lors de l'exécution de DBZONO.
Etant donné que ce fichier est du type standard, il est possible de visualiser certains champs à partir de l'application REC. Les Figures 1 et 2 suivantes montrent les champs S/ et C/ auquels sont superposées les bandes de latitudes, ces deux champs sont complémentaires. La Figure 3 montre le champ W/; le poid relatif de chaque point de grille. Notons que la grille est irrégulière d'où l'aspect bizarre de ce dernier. Le fichier noutzon contient aussi (Tableau 1V) les champs 1/ et 2/ pour 2 pas de temps différents;1, 30 i.e. 720*30/3600=6 heures. Il n'y a qu'un champ portant l'étiquette FJUI23F il contient les niveaux du modèle. Il ne sert pas pour l'extraction comme tel, mais il permet d'identifier le fichier car c'est le seul portant l'étiquette de l'expérience. Les variables échantillonnées au début le sont vraiment au premier pas de temps pour des raisons numériques.
Dans le fichier ZONFILE; chaque enregistrements contient une variable de surface ou de profil contenant les moyennes pour chaque bandes pour toutes les périodes spécifiées par NPTEMPS. L'heure zéro pour les variables (en réalité le premier pas) est nécessaire pour le calcul des tendances. Notons que les niveaux du modèle ont été recopiés du fichier NOUTZON, ces niveaux sont nécessaires pour l'affichage des variables de profil. Les caractéristiques du modèle apparaissent dans la définition des clefs IG1 à IG4. La compaction est à R24, i.e. il y a 24 bits stockées pour chaque valeur de nombre réel. Ce fichier est conforme à celui antérieurement obtenu de FESERI dans l'ancien système d'extraction.
Chartier Y.(1992): An introduction ot RPN standard files. Section Informatique RPN
Reid, B.(1989): Diagnostiques zonaux dans les modeles EFR et SEF. Document interne RPN
Benoit, R., Reid B, Bilodeau B. (1990): Diagnostics zonaux. Document interne CMC-RPN
Bourassa, Y. (1994) GRAPHZN - Guide de l'usager. Document interne RPN
1
********************************************************************************************
* *
* VOIR 3.4 *
* *
* *
* Fri Feb 25 14:25:36 1994 *
* *
* BEGIN EXECUTION *
* *
********************************************************************************************
UNIT = 10 EST OUVERT RANDOM
1 FILE=noutzon TYPE=RANDOM Fri Feb 25 1994 14:25:36 PAGE 1
KEY# ID IP1 IP2 IP3 NI NJ NK ETIQ. DATE ORIG. DEET NPAS GR IG1 IG2 IG3 IG4 DTY SWA LNG
0-T/ + 0 0 0 14 1 1 CONTROLE 2021279121 1800 48 G 0 0 0 0 I32 1531 14
1-V/ + 0 0 0 196 1 1 CONTROLE 2021279121 1800 0 G 0 0 0 0 C08 1561 49
2-P/ + 0 0 0 48 1 1 CONTROLE 2021279121 1800 0 G 0 0 0 0 I32 1621 48
3-W/ + 0 0 0 128 64 1 CONTROLE 2021279121 1800 0 G 0 0 0 0 R24 1681 6148
4-S/ + 0 0 0 128 64 1 CONTROLE 2021279121 1800 0 G 0 0 0 0 R24 7831 6148
5-C/ + 0 0 0 128 64 1 CONTROLE 2021279121 1800 0 G 0 0 0 0 R24 13981 6148
6-B/ + 0 0 0 128 64 1 CONTROLE 2021279121 1800 0 G 0 0 0 0 I24 20131 6144
7-S^ + 0 0 0 21 1 1 SEFEXP33 2021279121 1800 0 G 0 0 0 0 R24 26281 20
8-L^ + 0 0 0 64 1 1 SEFEXP33 2021279121 1800 0 G 0 0 0 0 R24 26311 52
9-1/ + 0 0 0 64 567 1 MZONXST 2021279121 1800 0 G 0 0 0 0 E32 26371 36288
10-2/ + 0 0 0 64 567 1 MZONXST 2021279121 1800 0 G 0 0 0 0 E32 62671 36288
11-1/ + 0 1 0 64 567 1 MZONXST 2021279121 1800 1 G 0 0 0 0 E32 98971 36288
12-2/ + 0 1 0 64 567 1 MZONXST 2021279121 1800 1 G 0 0 0 0 E32 135271 36288
13-1/ + 0 24 0 64 567 1 MZONXST 2021279121 1800 24 G 0 0 0 0 E32 171571 36288
14-2/ + 0 24 0 64 567 1 MZONXST 2021279121 1800 24 G 0 0 0 0 E32 207871 36288
15-1/ + 0 48 0 64 567 1 MZONXST 2021279121 1800 48 G 0 0 0 0 E32 244171 36288
16-2/ + 0 48 0 64 567 1 MZONXST 2021279121 1800 48 G 0 0 0 0 E32 280471 36288
STATISTIQUES
DIMENSION DU DIRECTEUR DISQUE 100
NOMBRE D ENTREES UTILISEES 17
LONGUEUR DU FICHIER 316770 MOTS
NOMBRE D ECRITURES 21
NOMBRE DE RE-ECRITURES 4
NOMBRE D EFFACAGES 4
NOMBRE D EXTENSIONS 0
NOMBRE DE CORRECTIONS 0
*****************************************
UNITE FORTRAN IUN= 10 EST FERME
********************************************************************************************
* *
* VOIR O.K. *
* *
* Fri Feb 25 14:25:36 1994 *
* *
* END EXECUTION *
* *
* CP SECS = 0.010 *
* *
********************************************************************************************
1
********************************************************************************************
* *
* VOIR 3.4 *
* *
* *
* Wed Feb 23 09:03:38 1994 *
* *
* BEGIN EXECUTION *
* *
********************************************************************************************
UNIT = 10 EST OUVERT RANDOM
1 FILE=zonfile TYPE=RANDOM Wed Feb 23 1994 09:03:38 PAGE 1
KEY# ID IP1 IP2 IP3 NI NJ NK ETIQ. DATE ORIG. DEET NPAS GR IG1 IG2 IG3 IG4 DTY SWA LNG
0-S^ Z 0 0 0 21 1 1 SEFEXP33 2021279121 1800 0 G 0 0 0 0 R24 57391 20
1-SH Z 0 0 0 21 1 1 SEFEXP33 2021279121 1800 0 G 0 0 0 0 R24 57541 20
2-L> Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 0 G 0 0 0 0 R24 57691 52
3-TS Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 57991 52
4-VE Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 58291 52
5-SE Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 58591 52
6-P0 Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 60811 52
7-FS Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 64951 52
8-FI Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 69091 52
9-FW Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 73231 52
10-FC Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 77371 52
11-FV Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 81511 52
12-FL Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2191 52
13-FQ Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2251 52
14-EI Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2311 52
15-EV Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2371 52
16-SI Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2431 52
17-AP Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2491 52
18-RR Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2551 52
19-PR Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2611 52
20-SN Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2671 52
21-SR Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2731 52
22-NT Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2791 52
23-RN Z 0 0 0 64 1 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2851 52
24-UU Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 2911 1012
25-VV Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 3931 1012
26-SS Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 4951 1012
27-TT Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 5971 1012
28-TF Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 6991 1012
29-TA Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 8011 1012
30-TI Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 9031 1012
31-T2 Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 10051 1012
32-HU Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 11071 1012
33-HR Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 12091 1012
34-NU Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 13111 1012
35-QF Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 14131 1012
36-QA Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 15151 1012
37-EN Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 16171 1012
38-KM Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 17191 1012
39-LM Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 18211 1012
40-TU Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 19231 1012
41-TV Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 20251 1012
42-GU Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 21271 1012
43-GV Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 22291 1012
44-TK Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 23311 1012
45-QK Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 24331 1012
46-NS Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 25351 1012
47-NC Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 26371 1012
1 FILE=zonfile TYPE=RANDOM Wed Feb 23 1994 09:03:38 PAGE 2
KEY# ID IP1 IP2 IP3 NI NJ NK ETIQ. DATE ORIG. DEET NPAS GR IG1 IG2 IG3 IG4 DTY SWA LNG
48-QL Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 27391 1012
49-EP Z 1 0 0 64 21 1 SEFEXP33 2021279121 1800 1 G 1 0 0 0 R24 28411 1012
50-TS Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29431 52
51-VE Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29491 52
52-SE Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29551 52
53-P0 Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29611 52
54-FS Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29671 52
55-FI Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29731 52
56-FW Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29791 52
57-FC Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29851 52
58-FV Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29911 52
59-FL Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 29971 52
60-FQ Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30031 52
61-EI Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30091 52
62-EV Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30151 52
63-SI Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30211 52
64-AP Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30271 52
65-RR Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30331 52
66-PR Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30391 52
67-SN Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30451 52
68-SR Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30511 52
69-NT Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30571 52
70-RN Z 0 12 24 64 1 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30631 52
71-UU Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 30691 1012
72-VV Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 31711 1012
73-SS Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 32731 1012
74-TT Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 33751 1012
75-TF Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 34771 1012
76-TA Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 35791 1012
77-TI Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 36811 1012
78-T2 Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 37831 1012
79-HU Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 38851 1012
80-HR Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 39871 1012
81-NU Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 40891 1012
82-QF Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 41911 1012
83-QA Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 42931 1012
84-EN Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 43951 1012
85-KM Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 44971 1012
86-LM Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 45991 1012
87-TU Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 47011 1012
88-TV Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 48031 1012
89-GU Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 49051 1012
90-GV Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 50071 1012
91-TK Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 51091 1012
92-QK Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 52111 1012
93-NS Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 53131 1012
94-NC Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 54151 1012
95-QL Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 55171 1012
1 FILE=zonfile TYPE=RANDOM Wed Feb 23 1994 09:03:38 PAGE 3
KEY# ID IP1 IP2 IP3 NI NJ NK ETIQ. DATE ORIG. DEET NPAS GR IG1 IG2 IG3 IG4 DTY SWA LNG
96-EP Z 1 12 24 64 21 1 SEFEXP33 2021279121 1800 24 G 1 0 0 0 R24 56191 1012
97-TS Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57211 52
98-VE Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57271 52
99-SE Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57331 52
100-P0 Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57421 52
101-FS Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57481 52
102-FI Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57571 52
103-FW Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57631 52
104-FC Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57751 52
105-FV Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57811 52
106-FL Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57871 52
107-FQ Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 57931 52
108-EI Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58051 52
109-EV Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58111 52
110-SI Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58171 52
111-AP Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58231 52
112-RR Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58351 52
113-PR Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58411 52
114-SN Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58471 52
115-SR Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58531 52
116-NT Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58651 52
117-RN Z 0 24 48 64 1 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58711 52
118-UU Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 58771 1012
119-VV Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 59791 1012
120-SS Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 60871 1012
121-TT Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 61891 1012
122-TF Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 62911 1012
123-TA Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 63931 1012
124-TI Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 65011 1012
125-T2 Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 66031 1012
126-HU Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 67051 1012
127-HR Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 68071 1012
128-NU Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 69151 1012
129-QF Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 70171 1012
130-QA Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 71191 1012
131-EN Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 72211 1012
132-KM Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 73291 1012
133-LM Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 74311 1012
134-TU Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 75331 1012
135-TV Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 76351 1012
136-GU Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 77431 1012
137-GV Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 78451 1012
138-TK Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 79471 1012
139-QK Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 80491 1012
140-NS Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 81571 1012
141-NC Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 82591 1012
142-QL Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 83611 1012
143-EP Z 1 24 48 64 21 1 SEFEXP33 2021279121 1800 48 G 1 0 0 0 R24 84631 1012
STATISTIQUES
DIMENSION DU DIRECTEUR DISQUE 144
NOMBRE D ENTREES UTILISEES 144
LONGUEUR DU FICHIER 85650 MOTS
NOMBRE D ECRITURES 144
NOMBRE DE RE-ECRITURES 0
NOMBRE D EFFACAGES 0
NOMBRE D EXTENSIONS 12
NOMBRE DE CORRECTIONS 0
*****************************************
UNITE FORTRAN IUN= 10 EST FERME
********************************************************************************************
* *
* VOIR O.K. *
* *
* Wed Feb 23 09:03:38 1994 *
* *
* END EXECUTION *
* *
* CP SECS = 0.130 *
* *
********************************************************************************************
LITZON - Liste des variables ============================= .TS. 1 2 +SE. 0 3 +P0. 0 4 +FS. 0 5 +FI. 0 6 +FW. 0 7 +FC. 0 8 +FV. 0 9 +FL. 0 10 +FQ. 0 11 +EI. 0 12 +EV. 0 13 +SI. 0 14 +AP. 0 15 +RR. 0 16 +PR. 0 17 +SN. 0 18 +RC. 0 19 +NT. 0 20 +RN. 0 21 .UU> 1 46 .VV> 1 71 .SS. 1 96 .TT. 1 121 .HU. 1 146 +TA. 0 171 +QA. 0 196 +EP. 0 221 +TK. 0 246 +QK. 0 271 .HR. 1 296 +NU. 0 321 +TI. 0 346 +T2. 0 371 +NS. 0 396 +NC. 0 421 +K2. 0 446 +GU> 0 471 +GV> 0 496 +KM. 0 521 +TU> 0 546 +TV> 0 571 +KT. 0 596 +KQ. 0 621 +QF. 0 646 +TF. 0 671
1
********************************************************************************************
* *
* VOIR 3.4 *
* *
* *
* Tue Mar 8 10:53:23 1994 *
* *
* BEGIN EXECUTION *
* *
********************************************************************************************
UNIT = 10 EST OUVERT RANDOM
1 FILE=noutzon TYPE=RANDOM Tue Mar 8 1994 10:53:23 PAGE 1
KEY# ID IP1 IP2 IP3 NI NJ NK ETIQ. DATE ORIG. DEET NPAS GR IG1 IG2 IG3 IG4 DTY SWA LNG
0-T/ + 0 0 0 14 1 1 CONTROLE 4062393121 1200 36 X 0 0 0 0 I32 1531 14
1-V/ + 0 0 0 192 1 1 CONTROLE 4062393121 1200 0 X 0 0 0 0 C08 1561 48
2-P/ + 0 0 0 47 1 1 CONTROLE 4062393121 1200 0 X 0 0 0 0 I32 1621 47
3-W/ + 0 0 0 125 101 1 CONTROLE 4062393121 1200 0 X 0 0 0 0 R24 1681 9473
4-S/ + 0 0 0 125 101 1 CONTROLE 4062393121 1200 0 X 0 0 0 0 R24 11161 9473
5-C/ + 0 0 0 125 101 1 CONTROLE 4062393121 1200 0 X 0 0 0 0 R24 20641 9473
6-B/ + 0 0 0 125 101 1 CONTROLE 4062393121 1200 0 X 0 0 0 0 I24 30121 9469
7-S^ + 0 0 0 25 1 1 FNECD07 4062393121 1200 0 X 0 0 0 0 R24 39601 23
8-1/ + 0 1 0 14 670 1 MZONXST 4062393121 1200 1 X 0 0 0 0 E32 39631 9380
9-2/ + 0 1 0 14 670 1 MZONXST 4062393121 1200 1 X 0 0 0 0 E32 49021 9380
10-1/ + 0 12 0 14 670 1 MZONXST 4062393121 1200 12 X 0 0 0 0 E32 58411 9380
11-2/ + 0 12 0 14 670 1 MZONXST 4062393121 1200 12 X 0 0 0 0 E32 67801 9380
12-1/ + 0 24 0 14 670 1 MZONXST 4062393121 1200 24 X 0 0 0 0 E32 77191 9380
13-2/ + 0 24 0 14 670 1 MZONXST 4062393121 1200 24 X 0 0 0 0 E32 86581 9380
14-1/ + 0 36 0 14 670 1 MZONXST 4062393121 1200 36 X 0 0 0 0 E32 95971 9380
15-2/ + 0 36 0 14 670 1 MZONXST 4062393121 1200 36 X 0 0 0 0 E32 105361 9380
STATISTIQUES
DIMENSION DU DIRECTEUR DISQUE 100
NOMBRE D ENTREES UTILISEES 16
LONGUEUR DU FICHIER 114750 MOTS
NOMBRE D ECRITURES 20
NOMBRE DE RE-ECRITURES 4
NOMBRE D EFFACAGES 4
NOMBRE D EXTENSIONS 0
NOMBRE DE CORRECTIONS 0
*****************************************
UNITE FORTRAN IUN= 10 EST FERME
********************************************************************************************
* *
* VOIR O.K. *
* *
* Tue Mar 8 10:53:23 1994 *
* *
* END EXECUTION *
* *
* CP SECS = 0.020 *
* *
********************************************************************************************
DBZONO(1) UTILITAIRE RPN DBZONO(1)
NOM
dbzono - Deballeur de diagnostics zonaux. Les vecteurs
d'accumulation ecrits durant l'integration sont
separes par variables pour les bandes specifiees.
USAGE
dbzono -ISZONAL noutzon -OMZONAL zonfile -ECHOZ non -DATE
DESCRIPTION
dbzono est un programme d'interface qui lit un fichier standard
NOUTZON et produit le fichier ZONFILE. Celui-ci est requis pour
le tracage de variables par bandes de latitude (GRAPHZN).
Les arguments:
-ISZONAL nom du fichier standard d'entree
-OMZONAL nom du fichier standard de sortie
-ECHOZ des diagnostics sont imprimes ou non durant l'execution
-DATE oprun/oui/non (facultatif)
Le fichier standard noutzon n'est plus requis apres l'execution
de dbzono. Toutes les informations requises pour les moyennes par
variables et bandes de latitude sont ecrites dans zonfile.
Quoiqu'il puisse traiter d'autre genre de moyenne(i.e. par regions)
dbzono a ete concu specifiquement pour GRAPHZN.
CONTENU DE NOUTZON
Une integration avec diagnostics zonaux aura produit un fichier
contenant des enregistrements de controle et des enregistrements
d'extraction. A chaque NPTEMPS, une table d'information est remise
a jour. Cette table contient les dimensions de la grille, le pas
de temps, le nombre de bandes, le nombre de pas pour l'experience
et une cle de controle requise pour une integration par clone.
Le fichier contient aussi la liste des variables a extraire,
les poids relatifs de chaque points de grille par rapport a la
surface du domaine et son rang (sa position dans le vecteur
d'accumulation). Cette information est ecrite une fois seulement
au debut de l'integration. De plus, des enregistrements contiennent
les niveaux du modele et la position relative des bandes. Ces
informations constituent les axes pour le logiciel de graphisme.
Les accumulateurs dans lesquels les variables du modele sont
moyennees ou echantilonnees sont ecrits dans un enregistrement dit
d'extraction. Le carre des moyennes peut aussi etre calcule.
Ces enregistrements seront sauves en meme temps.
CONTENU DE ZONFILE
La moyenne pour chaque variable se fait sur des periodes specifiees
et il peut y en avoir plusieurs selon la duree de l'integration.
Par exemple, pour une integration de 48 heures on peut vouloir
sauver les moyennes a toutes les 12 heures ou 24 heures.
Le fichier de sortie contient donc ces moyennes pour chaques
periodes. Les enregistrements qui contiennent plusieurs niveaux
peuvent etre visualiser avec REC. Ce fichier est conforme a ce qui
est requis pour GRAPHZN, le trace des variables par bandes.
INFORMATION ADDITIONELLE
Le deballeur peut manipuler un fichier provenant d'une integration
sur MIPS (32 bits) ou sur NEC (64bits/32bits)
Le contenu du fichier NOUTZON est ecrit de sorte qu'il puisse etre
lu a partir d'un serveur. La compaction est a 32 bits et la table
d'information contient seulement des entiers. La liste des
variables est ecrite en format R4A. Les vecteurs d'accumulation
sont ecrits en format IEEE.
Les routines requises: BSORT,DEBALZN,LITZON,OUTFLD,ZONFIN,WRITLZN
DIFUVD8,QQEXIT,R4STRG
EXEMPLES
Lors de l'execution de DBZONO, le premier enregistrement qui est
lu est la table d'information. Cette table apparait comme suit:
LITZON - table de controle
==========================
NDELTAT= 24
DELTAT= 1800
MODE= 3
NI= 128
NJ= 64
NK= 21
NBIN= 64
SOMNK= 107
LATMIN= -90
==========================
Ces parametres ont ete prealablement definis avant l'integration.
et ecrit dans un enregistrement de controle.
On a sauve les diagnostics a tous les 24 pas de 1800 secondes, soit
aux 12 heures. Les moyennes et le carre des moyennes sont sauves
(MODE=3). Il y a 64 bandes de latitude et la position finale de la
derniere variable extraite est a 107 dans le vecteur d'accumulation
Il y 128 par 64 points de grille dans l'horizontale et 21 niveauxi
dans la verticale.
La lecture des deux enregistrements suivant nous retourne les
niveaux et la position relative des bandes de latitude.
Les niveaux intermediaires seront calcules s'il n'ont pas ete
ecrit dans le fichier.
Suite a la lecture de l'enregistrement de controle des variables,
la liste des variables extraites sera ecrite de meme que les
proprietes de cette variable, i.e. s'il s'agit d'une variable
echantillonee ou moyennee et sa position relative dans le vecteur
d'accumulation.
LITZON - liste des variables
============================
.TS. 1 2
+VE. 0 3
.UU> 1 24
.VV> 1 45
.SS. 1 66
.TT. 1 87
.TT. 1 108
============================
ETIKEX= SEFEXP27
Meme si la variable comporte 4 caracteres, seul les trois premiers
sont utilises. La deuxieme colonne indique s'il s'agit d'une
variable echantillonee(1) ou moyennee(0). Le caractere . ou
+ dans le nom de la variable indique la meme chose. La troisieme
colonne indique s'il s'agit d'une tranche verticale ou non, c'est
l'indicateur de position des niveaux verticaux correspondant a
cette variable dans la tranche d'extraction.
Il y a deux variables de plus dans la liste ou l'etiquette de
l'experience a ete encodee. Dans l'exemple, l'etiquette indique
qu'il s'agit d'une integration du modele globale.
Avec ces informations, il est possible de lire les accumulateurs et
de les scinder par variables pour chaque periode sauvegardee.
EXPLICATION ADDITIONELLE
voir MZONOPR(1),MZONXST(2),MVZNXST(3) pour l'extraction des
variables du modele.
INZONO(1) UTILITAIRES RPN INZONO(1)
NOM
inzono - Routine d'emballage des diagnostics zonaux.
Les enregistrements de controle sont precalcules
APPEL
INZONO (POIDS, RANG, THETA, NDELTAT, DELTAT, MODE,
DZNSRF, ZSURFAC, DZNPRF, ZPROFIL, LATMIN,
ROT, IUN, S, ETIKX, IDAYO, NI, NJ, NK)
*ARGUMENTS
* POIDS - poids relatifs des points de grille pour extraction
* RANG - le numero de la bande pour position des accumulateurs
* THETA - angle de rotation de la grille p/r Greenwich
* NDELTAT - nombre de pas de temps d'accumulations
* DELTAT - nombre de seconde entre chaque pas de temps
* MODE - sauve les moyennes, la somme des moyennes ou les deux
* DZNSRF - nombre de variables de surface
* ZSURFAC - variables de surface demandees
* DZNPRF - nombre de variables de surface
* ZPROFIL - variables de profil demandees
* LATMIN - plus grand cercle de latitude inscrit dans la grille
* ROT - l'angle que fait l'axe des X avec le meridien Greenwhich
* IUN - le numero du fichier standard ou l'on ecrit l'information
* S - les niveaux du modele
* ETIKX - l'etiket de l'experience
* IDAYO - le datestamp de l'experience
* NI - dimension horizontale de la grille du modele
* NJ - deuxieme dimension horizontale de la grille
* NK - nombre de niveaux du modele
NOTES
Dans le modele global la tranche correspond a une bande de
latitude complete. Dans le modele regional, la tranche correspond
a la deuxieme dimension du plan horizontale de la grille, i.e. la
direction Y, la troixieme etant la direction Z, i.e. les niveaux du
modele. Une tranche verticale correspond au plan verticale dans la
direction X de la grille.
Les rangs assignes a chaque point de grille, assure l'appartenance
a une region d'interet. Par exemple une bande de latitude pour les
diagnostics zonaux.
Les poids assignes a chaque point de grille accorde une importance
relative du point dans la region definie. Ces vecteurs sont
precalcules dans la routine ZEFRPRT pour le modele regional et
ZONEPDS pour le modele global.
INFORMATIONS REQUISES POUR L'EXTRACTION
Les enregistrements de CONTROLE contiennent des definitions propres
a une extraction des diagnostics. Ces informations sont contenues
dans des enregistrements predefinis avant l'integration d'un modele
Les cles suivantes definissent ces enregistrements:
T/ la table d'information: le nombre d'extractions, les dimensions
de la grille, le pas de temps, le nombre de bandes de latitude,
l'angle de rotation si tel est le cas
V/ la liste des variables a echantillonner ou a moyenner
P/ la position relative des variables dans les accumulateurs
W/ les poids respectifs de chaque points de grille
B/ le rang i.e l'appartenance a une bande de latitude
S/ les sinus de l'angle de rotation i.e. angle que fait l'axe des
horizontale de la grille avec le meridien de Greenwhich
C/ les cosinus de l'angle de rotation. Ces deux derniers sont requis
pour l'extration des composantes vectorielles i.e. UU et VV
L^ la position des poins de grille le long d'un meridien
pour le modele SEF. Necessaire pour le logiciel de graphisme.
S^ les niveaux sigmas. Necessaire pour le logiciel de graphisme.
DOCUMENTATION ADDITIONELLE
voir MZONOPR, MZONXST, MVZNXST, DBZONO.
MVZNXST(3) UTILITAIRES RPN MVZNXST(3)
NOM
mvznxst - Routine d'extraction des composantes vectorielles des
diagnostics zonaux.
APPEL
MVZNXST ( X,Y,XNOM,YNOM,J,NDIM,FACF,ORD,xxxx)
ARGUMENTS
X,Y - Composantes dans la direction X et Y du plan horizontal
XNOM - Chaine de deux caracteres. Variables a extraire ou non.
YNOM - Idem. Peut aussi etre une tendance vectorielle.
J - Tranche a extraire
NDIM - Premiere dimension des composantes du vecteur X ou Y
FACF - Facteur multiplicatif a appliquer avant transfert
dans les accumulateurs
ORD - Cle specifiant ce qui doit etre sauve
(-1) tout le vecteur a sauve
(0) extraction du niveau NK seulement
(K) extraction du K ieme niveau du modele
xxxx - ne sert plus
NOTES
Dans le modele global la tranche correspond a une bande complete de
latitude. Dans le modele regional, la tranche correspond
a la deuxieme dimension du plan horizontal de la grille, i.e. la
direction Y, la troixieme etant la direction Z, i.e. les niveaux du
modele. Une tranche verticale correspond au plan verticale dans la
direction X de la grille.
Les rangs assignes a chaque point de grille, assure l'appartenance
a une region d'interet, par exemple une bande de latitude pour les
diagnostics zonaux. Le poid assigne a chaque point de grille accorde
une importance relative au point dans la region definie.
Ces vecteurs sont precalcules dans la routine ZEFRPRT pour le modele
regional et ZONEPDS pour le modele global.
Les enregistrements de controle sont definis dans la routine INZONO.
L'extraction se fait une tranche a la fois. L'extraction correcte
d'un vecteur a 3 dimensions est faite par la routine XST3D, qui elle
meme appelle MZONXST.
Une conversion d'unite peut etre effectuer via FAFC au moment de
l'extraction des composantes.
Pour minimiser le cout de l'extraction des composantes, il est
important que les deux variables soient une a la suite de l'autre
dans la liste des variables a extraire.
L'extraction se fait une tranche a la fois.
Une conversion d'unite peut etre effectuer via FAFC au moment de
l'extraction des composantes vectorielles.
Le mode multitache est limite a quatre processeurs( maximum 4 )
DOCUMENTATION ADDITIONELLE
voir MZONOPR, MZONXST, DBZONO.
MZONOPR(1) UTILITAIRES RPN MZONOPR(1)
NOM
mzonopr - Routine de controle pour l'extraction des
diagnostics zonaux.
APPEL
MZONOPR (OPR,VAL)
ARGUMENTS
OPR - clef de changement de transition (+:debut ou -:fin)
(1) ouverture et fermeture du fichier noutzon
(2) debut et fin d'un pas de temps
(3) debut et fin d'une tranche a extraire
(4) specification du nombre de niveau a extraire
VAL - Valeur associee avec la transition d'etat
si |opr| = 1 l'unite d'entree/sortie (IUN)
si |opr| = 2 le pas de temps en cours (KOUNT)
si |opr| = 3 le numero de la tranche (J)
si |opr| = 4 le nombre de niveau a extraire pour une
tranche verticale (KA)
NOTES
La specification des cles est indispensable au bon fonctionnement
de l'extraction des diagnostiques zonaux. Une variable d'etat
est maintenue a jour tout au long de l'integration et un test de
validite assure que l'ordre de transition est respectee.
Par exemple, il importe d'ouvrir le fichier avant tout autre
operation d'extraction. Une bonne marche requiert les transitions
suivantes: +1 ouverture du fichier (IUN) et allocation memoire des
vecteurs d'accumulations et autres vecteurs de travail; +2 Debut
d'une integration avec le passage d'un pas de temps approprie en
deuxieme argument. +3 Debut d'une extraction avec le numero de la
tranche en deuxieme argument; -3 indique la fin de l'extraction
pour cette tranche. Toutes les tranches doivent etre fermees avant
la transition suivante. -2 indique que le pas de temps courant est
termine. -1 indique que l'integration est terminee. Il y aura
ecriture des vecteurs d'accumulations, desallocation memoire des
champs et fermeture du fichier (IUN).
REMARQUES
Il y a une certaine rigidite dans la procedure d'extraction qui
est necessaire pour economiser le temps de calcul (CPU).
Quand le systeme d'extraction rencontre un etat innaproprie,
l'integration du modele avorte immediatement avec un message
d'erreur approprie. Les parametres de controle sont predefinis
et ecris dans la table d'information qui est mise a jour a chaque
ecriture des accumulateurs. De meme, les variables que l'on veut
extraire sont passees a la liste de comparaison qui est stockee
dans un enregistrement de controle lu a l'ouverture du fichier.
Une routine d'initialisation ZONINIT a ete ecrite pour les deux
modeles. Tous les champs qui sont requis sont precalcules et stores
a l'entre du modele.
L'extraction se fait une tranche a la fois. L'extraction correcte
d'un vecteur a trois dimensions est fait par la routine XST3D qui
elle meme appelle la routine MZONXST.
Le mode multitache est limite a quatre processeurs( maximum 4 )
DOCUMENTATION ADDITIONELLE
voir MZONXST, MVZNXST, DBZONO.
MZONXST(2) UTILITAIRES RPN MZONXST(2)
NOM
mzonxst - Routine d'extraction des diagnostics zonaux.
APPEL
MZONXST (Z,NOM,J,NDIM,xxxx,FACF,ORD,xxxx)
ARGUMENTS
Z - vecteur a extraire provenant du modele
NOM - Chaine de deux caracteres a extraire ou non.
J - Tranche a extraire
NDIM - Premiere dimension du vecteur Z
heure - heure de l'extraction, ne sert plus
FACF - Facteur multiplicatif a appliquer avant transfert
dans les accumulateurs
ORD - Cle specifiant ce qui doit etre sauve
(-1) tout le vecteur a sauve
(0) un niveau seulement a sauve
(K) extraction de KA niveaux du modele
xxxx - ne sert plus
NOTES
Dans le modele global la tranche correspond a une bande complete de
latitude. Dans le modele regional, la tranche correspond
a la deuxieme dimension du plan horizontal de la grille, i.e. la
direction Y, la troixieme etant la direction Z, i.e. les niveaux du
modele. Une tranche verticale correspond au plan verticale dans la
direction X de la grille.
Les rangs assignes a chaque point de grille, assure l'appartenance
a une region d'interet, par exemple une bande de latitude pour les
diagnostics zonaux. Le poid assigne a chaque point de grille accorde
une importance relative au point dans la region definie.
Ces vecteurs sont precalcules dans la routine ZEFRPRT pour le modele
regional et ZONEPDS pour le modele global.
Les enregistrements de controle sont definis dans la routine INZONO.
L'extraction se fait une tranche a la fois. L'extraction correcte
d'un vecteur a 3 dimensions est faite par la routine XST3D, qui elle
meme appelle MZONXST.
Une conversion d'unite peut etre effectuer via FAFC au moment de
l'extraction des composantes.
Le mode multitache est limite a quatre processeurs( maximum 4 )
DOCUMENTATION ADDITIONELLE
voir MZONOPR, MVZNXST, DBZONO.