1 /*************************************************************************
2 * COPYRIGHT (C) 1999 - 2003 EDF R&D
3 * THIS LIBRARY IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
4 * IT UNDER THE TERMS OF THE GNU LESSER GENERAL PUBLIC LICENSE
5 * AS PUBLISHED BY THE FREE SOFTWARE FOUNDATION;
6 * EITHER VERSION 2.1 OF THE LICENSE, OR (AT YOUR OPTION) ANY LATER VERSION.
7 *
8 * THIS LIBRARY IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
9 * WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
10 * MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
11 * LESSER GENERAL PUBLIC LICENSE FOR MORE DETAILS.
12 *
13 * YOU SHOULD HAVE RECEIVED A COPY OF THE GNU LESSER GENERAL PUBLIC LICENSE
14 * ALONG WITH THIS LIBRARY; IF NOT, WRITE TO THE FREE SOFTWARE FOUNDATION,
15 * INC., 59 TEMPLE PLACE, SUITE 330, BOSTON, MA 02111-1307 USA
16 *
17 *************************************************************************/
18
19
20 /******************************************************************************
21 * - Nom du fichier : test22.c
22 *
23 * - Description : lecture de valeurs scalaires numeriques crees dans test21.
24 *
25 *****************************************************************************/
26
27 #include <med.h>
28 #include <med_utils.h>
29 #include <stdio.h>
30
31 int main ()
32 {
33 med_idt fid;
34 char nom_scalaire[MED_TAILLE_NOM+1];
35 char description[MED_TAILLE_DESC+1];
36 med_int vali;
37 med_float valr;
38 med_int i,n,npdt,j;
39 med_type_champ type;
40 med_int numdt,numo;
41 med_float dt;
42 char dt_unit[MED_TAILLE_PNOM+1];
43
44 /* Ouverture du fichier test21.med en lecture seule */
45 if ((fid = MEDouvrir("test21.med",MED_LECTURE)) < 0) {
46 MESSAGE("Erreur a l'ouverture du fichier test21.med");
47 return -1;
48 }
49
50 /* Lecture du nombre de variable scalaire */
51 n = MEDnScalaire(fid);
52 if (n < 0) {
53 MESSAGE("Erreur a la lecture du nombre de variable scalaire");
54 return -1;
55 }
56 printf("Nombre de variables scalaires dans tes21.med = %d\n",n);
57
58 /* Lecture des infos sur les variables (type,description) */
59 for (i=1;i<=n;i++) {
60
61 if (MEDscalaireInfo(fid,i,nom_scalaire,&type,description) < 0) {
62 MESSAGE("Erreur a la lecture des infos sur la variable d'indice : ");
63 ISCRUTE(i);
64 return -1;
65 }
66 printf("- Scalaire n°%d de nom %s \n",i,nom_scalaire);
67 if (type == MED_FLOAT64)
68 printf(" Type flottant. \n");
69 else
70 printf(" Type entier. \n");
71 printf(" Description associee : [%s] \n",description);
72
73 /* Pour chaque scalaire on regarde les valeurs associees
74 eventuellement a des pas de temps et des numeros d'ordre */
75 npdt = MEDnScalairePasdetemps(fid,nom_scalaire);
76 if (npdt < 0) {
77 MESSAGE("Erreur a la lecture du nombre de pas de temps");
78 return -1;
79 }
80 printf(" Nombre de valeurs stockees : %d \n",npdt);
81
82 for (j=1;j<=npdt;j++) {
83
84 if (MEDscalairePasdetempsInfo(fid,nom_scalaire,j,&numdt,dt_unit,&dt,&numo) < 0) {
85 MESSAGE("Erreur a la lecture des informations sur le pas de temps d'indice :");
86 ISCRUTE(j);
87 return -1;
88 }
89
90 printf(" Valeur n°%d : \n",j);
91 if (numdt == MED_NOPDT)
92 printf(" - Aucun de pas de temps \n");
93 else
94 printf(" - Pas de de temps de numero %d de valeur %f [%s] \n",numdt,dt,dt_unit);
95 if (numo == MED_NONOR)
96 printf(" - Aucun numero d'ordre \n");
97 else
98 printf(" - Numero d'ordre : %d \n",numo);
99
100 /* Lecture de la valeur flottante associee au pas de temps */
101 if (type == MED_FLOAT64) {
102 if (MEDscalaireFlottantLire(fid,nom_scalaire,&valr,numdt,numo) < 0) {
103 MESSAGE("Erreur a la lecture de la valeur flottante : ");
104 SSCRUTE(nom_scalaire); ISCRUTE(numdt);ISCRUTE(numo);
105 return -1;
106 }
107 printf(" - Valeur : %f \n",valr);
108 } else {
109 /* Lecture de la valeur scalaire associee au pas de temps */
110 if (MEDscalaireEntierLire(fid,nom_scalaire,&vali,numdt,numo) < 0) {
111 MESSAGE("Erreur a la lecture de la valeur entiere : ");
112 SSCRUTE(nom_scalaire); ISCRUTE(numdt);ISCRUTE(numo);
113 return -1;
114 }
115 printf(" - Valeur : %d \n",vali);
116 }
117 }
118
119 }
120
121 /* Fermeture du fichier */
122 if (MEDfermer(fid) < 0) {
123 MESSAGE("Erreur a la fermeture du fichier ");
124 return -1;
125 }
126
127 return 0;
128 }