1 /* This file is part of MED.
2 *
3 * COPYRIGHT (C) 1999 - 2019 EDF R&D, CEA/DEN
4 * MED is free software: you can redistribute it and/or modify
5 * it under the terms of the GNU Lesser General Public License as published by
6 * the Free Software Foundation, either version 3 of the License, or
7 * (at your option) any later version.
8 *
9 * MED is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 * GNU Lesser General Public License for more details.
13 *
14 * You should have received a copy of the GNU Lesser General Public License
15 * along with MED. If not, see <http://www.gnu.org/licenses/>.
16 */
17
18 /******************************************************************************
19 * - Nom du fichier : test11.c
20 *
21 * - Description : lecture de champs de resultats MED
22 *
23 *****************************************************************************/
24
25 #include <med.h>
26 #define MESGERR 1
27 #include "med_utils.h"
28 #include <string.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
main(int argc,char ** argv)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_int(type);SSCRUTE(desc);
90 return -1;
91 }
92
93 printf("Maillage de nom |%s| et de dimension "IFORMAT" \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 : "IFORMAT" \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_int(i+1);SSCRUTE(nomcha);ISCRUTE_int(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 : "IFORMAT"\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_int(i);
160 ret = -1;continue;
161 }
162 printf("\t- Profil n°%i de nom |%s| et de taille "IFORMAT"\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(" "IFORMAT" ",*(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 : "IFORMAT"\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_int(i);
183 ret = -1;continue;
184 }
185 printf("\t- Lien n°%i de nom |%s| et de taille "IFORMAT"\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 : "IFORMAT"\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_int(i);
208 ret = -1;continue;
209 }
210 printf("\t- Loc. n°%i de nom |%s| et nbr. de pts de GAUSS "IFORMAT"\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
getFieldsOn(med_idt fid,char * nomcha,med_type_champ typcha,med_int ncomp,med_entite_maillage entite,med_mode_switch stockage)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."IFORMAT" (%f) [%s], n. d'ordre "IFORMAT", avec "IFORMAT" 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 : "IFORMAT" 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)+1);
331 EXIT_IF(lien == NULL,NULL,NULL);
332 if ( MEDlienLire(fid, lien, maa_ass) < 0 ) {
333 MESSAGE("Erreur a la lecture du lien : ");
334 SSCRUTE(maa_ass);SSCRUTE(lien);
335 ret = -1;
336 } else {
337 printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
338 }
339 free(lien);
340 }
341 }
342
343 /* Combien de maillages lies aux (nomcha,ent,geo,numdt,numo) */
344 /* Notons que cette information est egalement disponible a partir de MEDpasdetempsInfo */
345 if ( (nbrefmaa = MEDnChampRef(fid,nomcha,entite,type_geo[k],numdt,numo) ) < 0 ) {
346 MESSAGE("Erreur a la demande du nombre de maillages references par le champ : ");
347 SSCRUTE(nomcha);
348 ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
349 SSCRUTE(MED_ENTITE_MAILLAGE_AFF[(int)entite]);SSCRUTE(AFF[k]);
350 ISCRUTE(numdt); ISCRUTE(numo);
351 ret = -1; continue;
352 };
353
354 for (l=0;l<nbrefmaa;l++) {
355
356 if ( MEDchampRefInfo(fid,nomcha,entite,type_geo[k],
357 l+1,numdt, numo, maa_ass, &local, &ngauss) <0 ) {
358 MESSAGE("Erreur a la demande d'information sur le maillage utilise par le champ n° : ");
359 SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
360 ISCRUTE_int(l+1);ISCRUTE(numdt); ISCRUTE(numo);SSCRUTE(maa_ass);
361 ret = -1; continue;
362 };
363
364 /* Prend en compte le nbre de pt de gauss automatiquement */
365 if ((nval = MEDnVal(fid,nomcha,entite,type_geo[k],numdt,numo,maa_ass,USER_MODE)) <= 0) {
366 MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
367 SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
368 ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);ISCRUTE_int(USER_MODE);
369 ret = -1; continue;
370 };
371
372
373 printf("\t- Il y a "IFORMAT" valeurs en mode %i. Chaque entite %s\
374 de type geometrique %s associes au maillage |%s| a "IFORMAT" pts de gauss \n",
375 nval,USER_MODE,MED_ENTITE_MAILLAGE_AFF[(int)entite],AFF[k],maa_ass,ngauss);
376
377 /* Le maillage reference est-il porte par un autre fichier */
378 if ( !local ) {
379
380 if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 ) {
381 MESSAGE("Erreur a la lecture de la taille du lien : ");
382 SSCRUTE(maa_ass);
383 ret = -1;
384 } else {
385
386 lien = malloc(lnsize*sizeof(char) + 1);
387 EXIT_IF(lien == NULL,NULL,NULL);
388
389 if ( MEDlienLire(fid, lien, maa_ass) < 0 ) {
390 MESSAGE("Erreur a la lecture du lien : ");
391 SSCRUTE(maa_ass);SSCRUTE(lien);
392 ret = -1;
393 } else {
394 lien[lnsize] = '\0';
395 printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
396 }
397 free(lien);
398 }
399 }
400
401 /*Lecture des valeurs du champ */
402 if (typcha == MED_FLOAT64) {
403
404 valr = (med_float*) calloc(ncomp*nval,sizeof(med_float));
405 EXIT_IF(valr == NULL,NULL,NULL);
406 if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)valr,stockage,MED_ALL,locname,
407 pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
408 MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
409 SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
410 ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
411 ret = -1;
412 };
413
414 } else {
415
416 vale = (med_int*) calloc(ncomp*nval,sizeof(med_int));
417 EXIT_IF(vale == NULL,NULL,NULL);
418 if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)vale,stockage,MED_ALL,locname,
419 pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
420 MESSAGE("Erreur a la lecture des valeurs du champ : ");
421 SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
422 ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
423 ret = -1;
424 };
425
426 }
427
428 if ( ngauss > 1 )
429 printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
430
431 if (entite == MED_NOEUD_MAILLE)
432 ngroup = type_geo[k] % 100;
433 else
434 ngroup = ngauss;
435
436 switch (stockage) {
437
438 case MED_FULL_INTERLACE :
439 printf("\t- Valeurs :\n\t");
440 for (m=0;m<nval/ngroup;m++) {
441 printf("|");
442 for (n=0;n<ngroup*ncomp;n++)
443 if (typcha == MED_FLOAT64)
444 printf(" %f ",*(valr+(m*ngroup*ncomp)+n));
445 else
446 printf(" "IFORMAT" ",*(vale+(m*ngroup*ncomp)+n));
447
448 }
449 break;
450
451 /*Affichage en fonction du profil à traiter*/
452 case MED_NO_INTERLACE :
453 printf("\t- Valeurs :\n\t");
454 for (m=0;m<ncomp;m++) {
455 printf("|");
456 for (n=0;n<nval;n++)
457 if (typcha == MED_FLOAT64)
458 printf(" %f ",*(valr+(m*nval)+n));
459 else
460 printf(" "IFORMAT" ",*(vale+(m*nval)+n));
461 }
462 break;
463 }
464
465 printf("|\n");
466 if (typcha == MED_FLOAT64) {
467 if ( valr ) {free(valr);valr = NULL;}}
468 else
469 if (vale) { free(vale);vale = NULL; }
470
471 /*Lecture du profil associe */
472 if (strcmp(pflname,MED_NOPFL) == 0 )
473 printf("\t- Profil : MED_NOPFL\n");
474 else {
475
476 if ( (pflsize = MEDnValProfil(fid,pflname)) <0 ) {
477 MESSAGE("Erreur a la lecture du nombre de valeurs du profil : ");
478 SSCRUTE(pflname);
479 ret = -1; continue;
480 }
481
482 printf("\t- Profil : |%s| de taille "IFORMAT"\n",pflname,pflsize);
483
484 pflval = (med_int*) malloc(sizeof(med_int)*pflsize);
485 EXIT_IF(pflval == NULL,NULL,NULL);
486 if ( MEDprofilLire(fid,pflval,pflname) <0) {
487 MESSAGE("Erreur a la lecture des valeurs du profil : ");
488 SSCRUTE(pflname);
489 ret = -1;
490 }
491 printf("\t");
492 for (m=0;m<pflsize;m++) printf(" "IFORMAT" ",*(pflval+m));
493 printf("\n");
494 free(pflval);
495
496 }
497
498 }
499 }
500 } /* fin for sur les mailles*/
501
502 return ret;
503 }
504
505