1 /*************************************************************************
2 * COPYRIGHT (C) 1999 - 2007 EDF R&D, CEA/DEN
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 * - Nom du fichier : test11.c
21 *
22 * - Description : lecture de champs de resultats MED
23 *
24 *****************************************************************************/
25
26 #include <med.h>
27 #define MESGERR
28 #include <med_utils.h>
29
30 #ifdef DEF_LECT_ECR
31 #define MODE_ACCES MED_LECTURE_ECRITURE
32 #elif DEF_LECT_AJOUT
33 #define MODE_ACCES MED_LECTURE_AJOUT
34 #else
35 #define MODE_ACCES MED_CREATION
36 #endif
37
38 #ifndef USER_INTERLACE
39 #define USER_INTERLACE MED_FULL_INTERLACE
40 #endif
41
42 #define USER_MODE MED_COMPACT
43
44 med_err getFieldsOn(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp,
45 med_entite_maillage entite, med_mode_switch stockage);
46
47 int main (int argc, char **argv)
48
49
50 {
51 med_err ret,lret;
52 med_idt fid;
53 char * fichier = NULL;
54 char maa[MED_TAILLE_NOM+1]="";
55 char desc[MED_TAILLE_DESC+1]="";
56 char pflname[MED_TAILLE_NOM+1]="",nomlien[MED_TAILLE_NOM+1]="";
57 char locname[MED_TAILLE_NOM+1]="";
58 char * lien = NULL;
59 char *comp= NULL, *unit= NULL;
60 char nomcha [MED_TAILLE_NOM+1]="";
61 med_int mdim,ncomp,ncha,npro,nln,pflsize,*pflval,nval;
62 med_int ngauss,nloc;
63 int t1,t2,t3;
64 med_type_champ typcha;
65 med_maillage type;
66 med_geometrie_element type_geo;
67 med_float *refcoo, *gscoo, *wg;
68 int i,j;
69
70 if (argc != 2) {
71 MESSAGE("Aucun nom de fichier precise, fichier test10.med utilise ");
72 fichier = "test10.med";
73 } else {
74 fichier = argv[1];
75 };
76
77
78 /* Ouverture du fichier med */
79 if ((fid = MEDouvrir(fichier,MED_LECTURE)) < 0){
80 MESSAGE("Erreur a l'ouverture du fichier : ");SSCRUTE(fichier);
81 return -1;
82 }
83
84 ret = 0;
85
86 /* infos sur le premier maillage */
87 if ( MEDmaaInfo(fid,1,maa,&mdim,&type,desc) < 0 ) {
88 MESSAGE("Erreur a la lecture des informations sur le maillage : ");
89 SSCRUTE(maa);ISCRUTE(mdim);ISCRUTE(type);SSCRUTE(desc);
90 return -1;
91 }
92
93 printf("Maillage de nom |%s| et de dimension %d \n",maa,mdim);
94
95 /* combien de champs dans le fichier */
96 if ((ncha = MEDnChamp(fid,0)) < 0) {
97 MESSAGE("Impossible de lire le nombre de champs : ");ISCRUTE(ncha);
98 return ncha;
99 }
100
101 printf("Nombre de champs : %d \n",ncha);
102
103 /* lecture de tous les champs */
104 for (i =0;i<ncha;i++) {
105 lret = 0;
106 printf("\nChamp numero : %d \n",i+1);
107
108 /* Lecture du nombre de composantes */
109 if ((ncomp = MEDnChamp(fid,i+1)) < 0) {
110 MESSAGE("Erreur a la lecture du nombre de composantes : "); ISCRUTE(ncomp);
111 ret = -1; continue;
112 }
113
114 /* Lecture du type du champ, des noms des composantes et du nom de l'unite*/
115 comp = (char*) malloc(ncomp*MED_TAILLE_PNOM+1);
116 EXIT_IF(comp == NULL,NULL,NULL);
117 unit = (char*) malloc(ncomp*MED_TAILLE_PNOM+1);
118 EXIT_IF(unit == NULL,NULL,NULL);
119
120 if ( MEDchampInfo(fid,i+1,nomcha,&typcha,comp,unit,ncomp) < 0 ) {
121 MESSAGE("Erreur a la demande d'information sur les champs : ");
122 ISCRUTE(i+1);SSCRUTE(nomcha);ISCRUTE(typcha);SSCRUTE(comp);SSCRUTE(unit);
123 ISCRUTE(ncomp);
124 ret = -1; continue;
125 }
126
127 printf("Nom du champ : |%s| de type %d\n",nomcha,typcha);
128 printf("Nom des composantes : |%s|\n",comp);
129 printf("Unites des composantes : |%s| \n",unit);
130
131 free(comp);
132 free(unit);
133
134
135 lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD, USER_INTERLACE );
136
137 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_MAILLE, USER_INTERLACE );
138 else { MESSAGE("Erreur a la lecture des champs aux noeuds "); ret = -1; continue;}
139
140 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_FACE,USER_INTERLACE);
141 else { MESSAGE("Erreur a la lecture des champs aux mailles "); ret = -1; continue;}
142
143 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_ARETE,USER_INTERLACE);
144 else {MESSAGE("Erreur a la lecture des champs aux faces "); ret = -1; continue;}
145
146 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD_MAILLE,USER_INTERLACE);
147 else {MESSAGE("Erreur a la lecture des champs aux aretes"); ret = -1; continue;}
148
149 if (lret != 0) {MESSAGE("Erreur a la lecture des champs aux noeuds des mailles "); ret = -1;};
150 }
151
152
153 /* Interrogation des profils */
154 npro = MEDnProfil(fid);
155
156 printf("\nNombre de profils stockes : %i\n\n",npro);
157 for (i=1 ; i <= npro ; i++ ) {
158 if ( MEDprofilInfo(fid, i, pflname, &nval) < 0) {
159 MESSAGE("Erreur a la demande d'information sur le profil n° : "); ISCRUTE(i);
160 ret = -1;continue;
161 }
162 printf("\t- Profil n°%i de nom |%s| et de taille %i\n",i,pflname,nval);
163 pflval = (med_int*) malloc(sizeof(med_int)*nval);
164 if ( MEDprofilLire(fid, pflval, pflname) < 0) {
165 MESSAGE("Erreur a la lecture des valeurs du profil : ");
166 SSCRUTE(pflname);
167 ret = -1;
168 } else {
169 printf("\t");
170 for (j=0;j<nval;j++) printf(" %i ",*(pflval+j));
171 printf("\n\n");
172 }
173 free(pflval);
174 }
175
176 /* Interrogation des liens */
177 nln = MEDnLien(fid);
178
179 printf("\nNombre de liens stockes : %i\n\n",nln);
180 for (i=1 ; i <= nln ; i++ ) {
181 if ( MEDlienInfo(fid, i, nomlien, &nval) < 0) {
182 MESSAGE("Erreur a la demande d'information sur le lien n° : "); ISCRUTE(i);
183 ret = -1;continue;
184 }
185 printf("\t- Lien n°%i de nom |%s| et de taille %i\n",i,nomlien,nval);
186
187 lien = malloc((nval+1)*sizeof(char));
188 EXIT_IF(lien == NULL,NULL,NULL);
189
190 if ( MEDlienLire(fid, lien, nomlien) < 0 ) {
191 MESSAGE("Erreur a la lecture du lien : ");
192 SSCRUTE(nomlien);SSCRUTE(lien);
193 ret = -1;
194 } else {
195 lien[nval] = '\0';
196 printf("\t\t|%s|\n\n",lien);
197 }
198 free(lien);
199 }
200
201 /* Interrogation des localisations des points de GAUSS */
202 nloc = MEDnGauss(fid);
203
204 printf("\nNombre de localisations stockees : %i\n\n",nloc);
205 for (i=1 ; i <= nloc ; i++ ) {
206 if ( MEDgaussInfo(fid, i, locname, &type_geo, &ngauss) < 0) {
207 MESSAGE("Erreur a la demande d'information sur la localisation n° : "); ISCRUTE(i);
208 ret = -1;continue;
209 }
210 printf("\t- Loc. n°%i de nom |%s| et nbr. de pts de GAUSS %i\n",i,locname,ngauss);
211 t1 = (type_geo%100)*(type_geo/100);
212 t2 = ngauss*(type_geo/100);
213 t3 = ngauss;
214 refcoo = (med_float *) malloc(sizeof(med_float)*t1 );
215 gscoo = (med_float *) malloc(sizeof(med_float)*t2 );
216 wg = (med_float *) malloc(sizeof(med_float)*t3 );
217
218 if ( MEDgaussLire(fid, refcoo, gscoo, wg, USER_INTERLACE, locname ) < 0) {
219 MESSAGE("Erreur a la lecture des valeurs de la localisation : ");
220 SSCRUTE(locname);
221 ret = -1;
222 } else {
223 printf("\t Coordonnees de l'element de reference de type %i :\n\t\t",type_geo);
224 for (j=0;j<t1;j++) printf(" %f ",*(refcoo+j));
225 printf("\n");
226 printf("\t Localisation des points de GAUSS : \n\t\t");
227 for (j=0;j<t2;j++) printf(" %f ",*(gscoo+j));
228 printf("\n");
229 printf("\t Poids associes aux points de GAUSS :\n\t\t");
230 for (j=0;j<t3;j++) printf(" %f ",*(wg+j));
231 printf("\n\n");
232 }
233 free(refcoo);
234 free(gscoo);
235 free(wg);
236 }
237
238
239 /* fermeture du fichier */
240 if ( MEDfermer(fid) < 0) return -1;
241
242 return ret;
243 }
244
245 med_err getFieldsOn(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp,
246 med_entite_maillage entite, med_mode_switch stockage) {
247
248 int j,k,l,m,n,nb_geo;
249 med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,ngroup,*vale=NULL,nval;
250 med_int numdt=0,numo=0,lnsize,nbrefmaa;
251 med_float *valr=NULL,dt=0.0;
252 med_err ret=0;
253 med_booleen local;
254 char pflname [MED_TAILLE_NOM+1]="";
255 char locname [MED_TAILLE_NOM+1]="";
256 char * lien = NULL;
257 char maa_ass [MED_TAILLE_NOM+1]="";
258 char dt_unit [MED_TAILLE_PNOM+1]="";
259
260
261 med_geometrie_element * type_geo;
262 med_geometrie_element typ_noeud[1] = { MED_NONE };
263 med_geometrie_element typmai[MED_NBR_GEOMETRIE_MAILLE+2] = {MED_POINT1, MED_SEG2, MED_SEG3, MED_TRIA3,
264 MED_QUAD4, MED_TRIA6,MED_QUAD8, MED_TETRA4,
265 MED_PYRA5, MED_PENTA6, MED_HEXA8, MED_TETRA10,
266 MED_PYRA13, MED_PENTA15, MED_HEXA20,
267 MED_POLYGONE, MED_POLYEDRE};
268 med_geometrie_element typfac[MED_NBR_GEOMETRIE_FACE+1] = {MED_TRIA3,MED_TRIA6,
269 MED_QUAD4,MED_QUAD8,
270 MED_POLYGONE};
271 med_geometrie_element typare[MED_NBR_GEOMETRIE_ARETE] = {MED_SEG2,MED_SEG3};
272
273 char ** AFF;
274
275 switch (entite) {
276 case MED_NOEUD :
277 type_geo = typ_noeud;
278 nb_geo = 1;
279 AFF = MED_GEOMETRIE_NOEUD_AFF;
280 break;
281 case MED_MAILLE :
282 case MED_NOEUD_MAILLE :
283 type_geo = typmai;
284 nb_geo = MED_NBR_GEOMETRIE_MAILLE+2;
285 AFF = MED_GEOMETRIE_MAILLE_AFF;
286 break;
287 case MED_FACE :
288 type_geo = typfac;
289 nb_geo = MED_NBR_GEOMETRIE_FACE+1;
290 AFF = MED_GEOMETRIE_FACE_AFF;
291 break;
292 case MED_ARETE :
293 type_geo = typare;
294 nb_geo = MED_NBR_GEOMETRIE_ARETE;
295 AFF = MED_GEOMETRIE_ARETE_AFF;
296 break;
297 }
298
299
300 for (k=0;k<nb_geo;k++) {
301
302 /* Combien de (PDT,NOR) a lire */
303 nbpdtnor = MEDnPasdetemps(fid,nomcha,entite,type_geo[k]);
304 if (nbpdtnor < 1 ) continue;
305
306 for (j=0;j<nbpdtnor;j++) {
307
308 if ( MEDpasdetempsInfo(fid,nomcha,entite,type_geo[k],
309 j+1, &ngauss, &numdt, &numo, dt_unit,
310 &dt, maa_ass, &local, &nbrefmaa) <0) {
311 MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
312 ISCRUTE(numdt); ISCRUTE(numo);
313 ret = -1; continue;
314 };
315
316 printf("\n +Pas de Temps n.%i (%f) [%s], n. d'ordre %i, avec %i pts de gauss sur le maillage par defaut.\n",numdt,dt,dt_unit,numo,ngauss);
317
318 printf("\tLe maillage par defaut est : |%s|, sur un total de : %i maillages associes\n",
319 maa_ass, nbrefmaa);
320
321 /* Le maillage reference est-il porte par un autre fichier */
322 if ( !local ) {
323
324 if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 ) {
325 MESSAGE("Erreur a la lecture de la taille du lien : ");
326 SSCRUTE(maa_ass);
327 ret = -1;
328 } else {
329
330 lien = malloc(lnsize*sizeof(char));
331 EXIT_IF(lien == NULL,NULL,NULL);
332
333 if ( MEDlienLire(fid, lien, maa_ass) < 0 ) {
334 MESSAGE("Erreur a la lecture du lien : ");
335 SSCRUTE(maa_ass);SSCRUTE(lien);
336 ret = -1;
337 } else {
338 printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
339 }
340 free(lien);
341 }
342 }
343
344 /* Combien de maillages lies aux (nomcha,ent,geo,numdt,numo) */
345 /* Notons que cette information est egalement disponible a partir de MEDpasdetempsInfo */
346 if ( (nbrefmaa = MEDnChampRef(fid,nomcha,entite,type_geo[k],numdt,numo) ) < 0 ) {
347 MESSAGE("Erreur a la demande du nombre de maillages references par le champ : ");
348 SSCRUTE(nomcha);
349 ISCRUTE(entite);ISCRUTE(type_geo[k]);
350 SSCRUTE(MED_ENTITE_MAILLAGE_AFF[(int)entite]);SSCRUTE(AFF[k]);
351 ISCRUTE(numdt); ISCRUTE(numo);
352 ret = -1; continue;
353 };
354
355 for (l=0;l<nbrefmaa;l++) {
356
357 if ( MEDchampRefInfo(fid,nomcha,entite,type_geo[k],
358 l+1,numdt, numo, maa_ass, &local, &ngauss) <0 ) {
359 MESSAGE("Erreur a la demande d'information sur le maillage utilise par le champ n° : ");
360 SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
361 ISCRUTE(l+1);ISCRUTE(numdt); ISCRUTE(numo);SSCRUTE(maa_ass);
362 ret = -1; continue;
363 };
364
365 /* Prend en compte le nbre de pt de gauss automatiquement */
366 if ((nval = MEDnVal(fid,nomcha,entite,type_geo[k],numdt,numo,maa_ass,USER_MODE)) <= 0) {
367 MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
368 SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
369 ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);ISCRUTE(USER_MODE);
370 ret = -1; continue;
371 };
372
373
374 printf("\t- Il y a %d valeurs en mode %i. Chaque entite %s\
375 de type geometrique %s associes au maillage |%s| a %i pts de gauss \n",
376 nval,USER_MODE,MED_ENTITE_MAILLAGE_AFF[(int)entite],AFF[k],maa_ass,ngauss);
377
378 /* Le maillage reference est-il porte par un autre fichier */
379 if ( !local ) {
380
381 if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 ) {
382 MESSAGE("Erreur a la lecture de la taille du lien : ");
383 SSCRUTE(maa_ass);
384 ret = -1;
385 } else {
386
387 lien = malloc(lnsize*sizeof(char) + 1);
388 EXIT_IF(lien == NULL,NULL,NULL);
389
390 if ( MEDlienLire(fid, lien, maa_ass) < 0 ) {
391 MESSAGE("Erreur a la lecture du lien : ");
392 SSCRUTE(maa_ass);SSCRUTE(lien);
393 ret = -1;
394 } else {
395 lien[lnsize] = '\0';
396 printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
397 }
398 free(lien);
399 }
400 }
401
402 /*Lecture des valeurs du champ */
403 if (typcha == MED_FLOAT64) {
404
405 valr = (med_float*) calloc(ncomp*nval,sizeof(med_float));
406 EXIT_IF(valr == NULL,NULL,NULL);
407 if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)valr,stockage,MED_ALL,locname,
408 pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
409 MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
410 SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
411 ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
412 ret = -1;
413 };
414
415 } else {
416
417 vale = (med_int*) calloc(ncomp*nval,sizeof(med_int));
418 EXIT_IF(vale == NULL,NULL,NULL);
419 if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)vale,stockage,MED_ALL,locname,
420 pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
421 MESSAGE("Erreur a la lecture des valeurs du champ : ");
422 SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
423 ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
424 ret = -1;
425 };
426
427 }
428
429 if ( ngauss > 1 )
430 printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
431
432 if (entite == MED_NOEUD_MAILLE)
433 ngroup = type_geo[k] % 100;
434 else
435 ngroup = ngauss;
436
437 switch (stockage) {
438
439 case MED_FULL_INTERLACE :
440 printf("\t- Valeurs :\n\t");
441 for (m=0;m<nval/ngroup;m++) {
442 printf("|");
443 for (n=0;n<ngroup*ncomp;n++)
444 if (typcha == MED_FLOAT64)
445 printf(" %f ",*(valr+(m*ngroup*ncomp)+n));
446 else
447 printf(" %d ",*(vale+(m*ngroup*ncomp)+n));
448
449 }
450 break;
451
452 /*Affichage en fonction du profil à traiter*/
453 case MED_NO_INTERLACE :
454 printf("\t- Valeurs :\n\t");
455 for (m=0;m<ncomp;m++) {
456 printf("|");
457 for (n=0;n<nval;n++)
458 if (typcha == MED_FLOAT64)
459 printf(" %f ",*(valr+(m*nval)+n));
460 else
461 printf(" %d ",*(vale+(m*nval)+n));
462 }
463 break;
464 }
465
466 printf("|\n");
467 if (typcha == MED_FLOAT64) {
468 if ( valr ) {free(valr);valr = NULL;}}
469 else
470 if (vale) { free(vale);vale = NULL; }
471
472 /*Lecture du profil associe */
473 if (strcmp(pflname,MED_NOPFL) == 0 )
474 printf("\t- Profil : MED_NOPFL\n");
475 else {
476
477 if ( (pflsize = MEDnValProfil(fid,pflname)) <0 ) {
478 MESSAGE("Erreur a la lecture du nombre de valeurs du profil : ");
479 SSCRUTE(pflname);
480 ret = -1; continue;
481 }
482
483 printf("\t- Profil : |%s| de taille %i\n",pflname,pflsize);
484
485 pflval = (med_int*) malloc(sizeof(med_int)*pflsize);
486 EXIT_IF(pflval == NULL,NULL,NULL);
487 if ( MEDprofilLire(fid,pflval,pflname) <0) {
488 MESSAGE("Erreur a la lecture des valeurs du profil : ");
489 SSCRUTE(pflname);
490 ret = -1;
491 }
492 printf("\t");
493 for (m=0;m<pflsize;m++) printf(" %i ",*(pflval+m));
494 printf("\n");
495 free(pflval);
496
497 }
498
499 }
500 }
501 } /* fin for sur les mailles*/
502
503 return ret;
504 }
505