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