# 2.3.6/test11.c

Aller à la documentation de ce fichier.

    00001 /*  This file is part of MED.
    00002  *
    00003  *  COPYRIGHT (C) 1999 - 2011  EDF R&D, CEA/DEN
    00004  *  MED is free software: you can redistribute it and/or modify
    00005  *  it under the terms of the GNU Lesser General Public License as published by
    00006  *  the Free Software Foundation, either version 3 of the License, or
    00007  *  (at your option) any later version.
    00008  *
    00009  *  MED is distributed in the hope that it will be useful,
    00010  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
    00011  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    00012  *  GNU Lesser General Public License for more details.
    00013  *
    00014  *  You should have received a copy of the GNU Lesser General Public License
    00015  *  along with MED.  If not, see <http://www.gnu.org/licenses/>.
    00016  */
    00017
    00018 /******************************************************************************
    00019  * - Nom du fichier : test11.c
    00020  *
    00021  * - Description : lecture de champs de resultats MED
    00022  *
    00023  *****************************************************************************/
    00024
    00025 #include <med.h>
    00026 #define MESGERR 1
    00027 #include "med_utils.h"
    00028 #include <string.h>
    00029
    00030 #ifdef DEF_LECT_ECR
    00031 #define MODE_ACCES MED_LECTURE_ECRITURE
    00032 #elif DEF_LECT_AJOUT
    00033 #define MODE_ACCES MED_LECTURE_AJOUT
    00034 #else
    00035 #define MODE_ACCES MED_CREATION
    00036 #endif
    00037
    00038 #ifndef USER_INTERLACE
    00039 #define USER_INTERLACE MED_FULL_INTERLACE
    00040 #endif
    00041
    00042 #define USER_MODE MED_COMPACT
    00043
    00044 med_err getFieldsOn(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp,
    00045                     med_entite_maillage entite, med_mode_switch stockage);
    00046
    00047 int main (int argc, char **argv)
    00048
    00049
    00050 {
    00051   med_err ret,lret;
    00052   med_idt fid;
    00053   char * fichier = NULL;
    00054   char maa[MED_TAILLE_NOM+1]="";
    00055   char desc[MED_TAILLE_DESC+1]="";
    00056   char pflname[MED_TAILLE_NOM+1]="",nomlien[MED_TAILLE_NOM+1]="";
    00057   char locname[MED_TAILLE_NOM+1]="";
    00058   char * lien = NULL;
    00059   char *comp= NULL, *unit= NULL;
    00060   char nomcha  [MED_TAILLE_NOM+1]="";
    00061   med_int mdim,ncomp,ncha,npro,nln,pflsize,*pflval,nval;
    00062   med_int ngauss,nloc;
    00063   int t1,t2,t3;
    00064   med_type_champ typcha;
    00065   med_maillage type;
    00066   med_geometrie_element type_geo;
    00067   med_float *refcoo, *gscoo, *wg;
    00068   int i,j;
    00069
    00070   if (argc != 2) {
    00071     MESSAGE("Aucun nom de fichier precise, fichier test10.med utilise ");
    00072     fichier = "test10.med";
    00073   } else {
    00074     fichier = argv[1];
    00075   };
    00076
    00077
    00078   /* Ouverture du fichier med */
    00079   if ((fid = MEDouvrir(fichier,MED_LECTURE)) < 0){
    00080     MESSAGE("Erreur a l'ouverture du fichier : ");SSCRUTE(fichier);
    00081     return -1;
    00082   }
    00083
    00084   ret = 0;
    00085
    00086   /* infos sur le premier maillage */
    00087   if (  MEDmaaInfo(fid,1,maa,&mdim,&type,desc) < 0 ) {
    00088     MESSAGE("Erreur a la lecture des informations sur le maillage : ");
    00089     SSCRUTE(maa);ISCRUTE(mdim);ISCRUTE_int(type);SSCRUTE(desc);
    00090     return -1;
    00091   }
    00092
    00093   printf("Maillage de nom |%s| et de dimension "IFORMAT" \n",maa,mdim);
    00094
    00095   /* combien de champs dans le fichier */
    00096   if ((ncha = MEDnChamp(fid,0)) < 0) {
    00097     MESSAGE("Impossible de lire le nombre de champs : ");ISCRUTE(ncha);
    00098     return ncha;
    00099   }
    00100
    00101   printf("Nombre de champs : "IFORMAT" \n",ncha);
    00102
    00103   /* lecture de tous les champs  */
    00104   for (i =0;i<ncha;i++) {
    00105     lret = 0;
    00106     printf("\nChamp numero : %d \n",i+1);
    00107
    00108     /* Lecture du nombre de composantes */
    00109     if ((ncomp = MEDnChamp(fid,i+1)) < 0) {
    00110       MESSAGE("Erreur a la lecture du nombre de composantes : "); ISCRUTE(ncomp);
    00111       ret = -1; continue;
    00112     }
    00113
    00114     /* Lecture du type du champ, des noms des composantes et du nom de l'unite*/
    00115     comp = (char*) malloc(ncomp*MED_TAILLE_PNOM+1);
    00116     EXIT_IF(comp == NULL,NULL,NULL);
    00117     unit = (char*) malloc(ncomp*MED_TAILLE_PNOM+1);
    00118     EXIT_IF(unit == NULL,NULL,NULL);
    00119
    00120     if ( MEDchampInfo(fid,i+1,nomcha,&typcha,comp,unit,ncomp) < 0 ) {
    00121       MESSAGE("Erreur a la demande d'information sur les champs : ");
    00122       ISCRUTE_int(i+1);SSCRUTE(nomcha);ISCRUTE_int(typcha);SSCRUTE(comp);SSCRUTE(unit);
    00123       ISCRUTE(ncomp);
    00124       ret = -1; continue;
    00125     }
    00126
    00127     printf("Nom du champ : |%s| de type %d\n",nomcha,typcha);
    00128     printf("Nom des composantes : |%s|\n",comp);
    00129     printf("Unites des composantes : |%s| \n",unit);
    00130
    00131     free(comp);
    00132     free(unit);
    00133
    00134
    00135     lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD, USER_INTERLACE );
    00136
    00137     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_MAILLE, USER_INTERLACE );
    00138     else { MESSAGE("Erreur a la lecture des champs aux noeuds "); ret = -1; continue;}
    00139
    00140     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_FACE,USER_INTERLACE);
    00141     else { MESSAGE("Erreur a la lecture des champs aux mailles "); ret = -1; continue;}
    00142
    00143     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_ARETE,USER_INTERLACE);
    00144     else {MESSAGE("Erreur a la lecture des champs aux faces "); ret = -1; continue;}
    00145
    00146     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD_MAILLE,USER_INTERLACE);
    00147     else {MESSAGE("Erreur a la lecture des champs aux aretes"); ret = -1; continue;}
    00148
    00149     if  (lret != 0) {MESSAGE("Erreur a la lecture des champs aux noeuds des mailles "); ret = -1;};
    00150   }
    00151
    00152
    00153   /* Interrogation des profils */
    00154   npro = MEDnProfil(fid);
    00155
    00156   printf("\nNombre de profils stockes : "IFORMAT"\n\n",npro);
    00157   for (i=1 ; i <= npro ; i++ ) {
    00158     if ( MEDprofilInfo(fid, i, pflname, &nval) < 0)  {
    00159       MESSAGE("Erreur a la demande d'information sur le profil n° : "); ISCRUTE_int(i);
    00160       ret = -1;continue;
    00161     }
    00162     printf("\t- Profil n°%i de nom |%s| et de taille "IFORMAT"\n",i,pflname,nval);
    00163     pflval = (med_int*) malloc(sizeof(med_int)*nval);
    00164     if ( MEDprofilLire(fid, pflval, pflname) < 0) {
    00165       MESSAGE("Erreur a la lecture des valeurs du profil : ");
    00166       SSCRUTE(pflname);
    00167       ret = -1;
    00168     } else {
    00169       printf("\t");
    00170       for (j=0;j<nval;j++) printf(" "IFORMAT" ",*(pflval+j));
    00171       printf("\n\n");
    00172     }
    00173     free(pflval);
    00174   }
    00175
    00176   /* Interrogation des liens */
    00177   nln = MEDnLien(fid);
    00178
    00179   printf("\nNombre de liens stockes : "IFORMAT"\n\n",nln);
    00180   for (i=1 ; i <= nln ; i++ ) {
    00181     if ( MEDlienInfo(fid, i, nomlien, &nval) < 0)  {
    00182       MESSAGE("Erreur a la demande d'information sur le lien n° : "); ISCRUTE_int(i);
    00183       ret = -1;continue;
    00184     }
    00185     printf("\t- Lien n°%i de nom |%s| et de taille "IFORMAT"\n",i,nomlien,nval);
    00186
    00187     lien = malloc((nval+1)*sizeof(char));
    00188     EXIT_IF(lien == NULL,NULL,NULL);
    00189
    00190     if ( MEDlienLire(fid, lien, nomlien) < 0 )  {
    00191       MESSAGE("Erreur a la lecture du lien : ");
    00192       SSCRUTE(nomlien);SSCRUTE(lien);
    00193       ret = -1;
    00194     } else {
    00195       lien[nval] = '\0';
    00196       printf("\t\t|%s|\n\n",lien);
    00197     }
    00198     free(lien);
    00199   }
    00200
    00201   /* Interrogation des localisations des points de GAUSS */
    00202   nloc = MEDnGauss(fid);
    00203
    00204   printf("\nNombre de localisations stockees : "IFORMAT"\n\n",nloc);
    00205   for (i=1 ; i <= nloc ; i++ ) {
    00206     if ( MEDgaussInfo(fid, i, locname, &type_geo, &ngauss) < 0)  {
    00207       MESSAGE("Erreur a la demande d'information sur la localisation n° : "); ISCRUTE_int(i);
    00208       ret = -1;continue;
    00209     }
    00210     printf("\t- Loc. n°%i de nom |%s| et nbr. de pts de GAUSS "IFORMAT"\n",i,locname,ngauss);
    00211     t1 = (type_geo%100)*(type_geo/100);
    00212     t2 = ngauss*(type_geo/100);
    00213     t3 = ngauss;
    00214     refcoo = (med_float *) malloc(sizeof(med_float)*t1 );
    00215     gscoo  = (med_float *) malloc(sizeof(med_float)*t2 );
    00216     wg     = (med_float *) malloc(sizeof(med_float)*t3 );
    00217
    00218     if ( MEDgaussLire(fid, refcoo, gscoo, wg, USER_INTERLACE, locname ) < 0) {
    00219       MESSAGE("Erreur a la lecture des valeurs de la localisation : ");
    00220       SSCRUTE(locname);
    00221       ret = -1;
    00222     } else {
    00223       printf("\t  Coordonnees de l'element de reference de type %i :\n\t\t",type_geo);
    00224       for (j=0;j<t1;j++) printf(" %f ",*(refcoo+j));
    00225       printf("\n");
    00226       printf("\t  Localisation des points de GAUSS : \n\t\t");
    00227       for (j=0;j<t2;j++) printf(" %f ",*(gscoo+j));
    00228       printf("\n");
    00229       printf("\t  Poids associes aux points de GAUSS :\n\t\t");
    00230       for (j=0;j<t3;j++) printf(" %f ",*(wg+j));
    00231       printf("\n\n");
    00232     }
    00233     free(refcoo);
    00234     free(gscoo);
    00235     free(wg);
    00236   }
    00237
    00238
    00239   /* fermeture du fichier */
    00240   if ( MEDfermer(fid) < 0) return -1;
    00241
    00242   return ret;
    00243 }
    00244
    00245 med_err getFieldsOn(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp,
    00246                     med_entite_maillage entite, med_mode_switch stockage) {
    00247
    00248   int j,k,l,m,n,nb_geo;
    00249   med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,ngroup,*vale=NULL,nval;
    00250   med_int numdt=0,numo=0,lnsize,nbrefmaa;
    00251   med_float *valr=NULL,dt=0.0;
    00252   med_err ret=0;
    00253   med_booleen local;
    00254   char pflname [MED_TAILLE_NOM+1]="";
    00255   char locname [MED_TAILLE_NOM+1]="";
    00256   char * lien = NULL;
    00257   char maa_ass [MED_TAILLE_NOM+1]="";
    00258   char dt_unit [MED_TAILLE_PNOM+1]="";
    00259
    00260
    00261   med_geometrie_element * type_geo;
    00262   med_geometrie_element typ_noeud[1] = { MED_NONE };
    00263   med_geometrie_element typmai[MED_NBR_GEOMETRIE_MAILLE+2] = {MED_POINT1, MED_SEG2, MED_SEG3, MED_TRIA3,
    00264                                                               MED_QUAD4, MED_TRIA6,MED_QUAD8, MED_TETRA4,
    00265                                                               MED_PYRA5, MED_PENTA6, MED_HEXA8, MED_TETRA10,
    00266                                                               MED_PYRA13, MED_PENTA15, MED_HEXA20,
    00267                                                               MED_POLYGONE, MED_POLYEDRE};
    00268   med_geometrie_element typfac[MED_NBR_GEOMETRIE_FACE+1] = {MED_TRIA3,MED_TRIA6,
    00269                                                             MED_QUAD4,MED_QUAD8,
    00270                                                             MED_POLYGONE};
    00271   med_geometrie_element typare[MED_NBR_GEOMETRIE_ARETE] = {MED_SEG2,MED_SEG3};
    00272
    00273   char ** AFF;
    00274
    00275   switch (entite) {
    00276   case MED_NOEUD :
    00277     type_geo = typ_noeud;
    00278     nb_geo   = 1;
    00279     AFF      = MED_GEOMETRIE_NOEUD_AFF;
    00280     break;
    00281   case  MED_MAILLE :
    00282   case  MED_NOEUD_MAILLE :
    00283     type_geo = typmai;
    00284     nb_geo   = MED_NBR_GEOMETRIE_MAILLE+2;
    00285     AFF      = MED_GEOMETRIE_MAILLE_AFF;
    00286     break;
    00287   case  MED_FACE :
    00288     type_geo = typfac;
    00289     nb_geo   = MED_NBR_GEOMETRIE_FACE+1;
    00290     AFF      = MED_GEOMETRIE_FACE_AFF;
    00291     break;
    00292   case  MED_ARETE :
    00293     type_geo = typare;
    00294     nb_geo   = MED_NBR_GEOMETRIE_ARETE;
    00295     AFF      = MED_GEOMETRIE_ARETE_AFF;
    00296     break;
    00297   }
    00298
    00299
    00300   for (k=0;k<nb_geo;k++) {
    00301
    00302     /* Combien de (PDT,NOR) a lire */
    00303     nbpdtnor = MEDnPasdetemps(fid,nomcha,entite,type_geo[k]);
    00304     if (nbpdtnor < 1 ) continue;
    00305
    00306     for (j=0;j<nbpdtnor;j++) {
    00307
    00308       if ( MEDpasdetempsInfo(fid,nomcha,entite,type_geo[k],
    00309                              j+1, &ngauss, &numdt, &numo, dt_unit,
    00310                              &dt, maa_ass, &local, &nbrefmaa) <0) {
    00311         MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
    00312         ISCRUTE(numdt); ISCRUTE(numo);
    00313         ret = -1; continue;
    00314       };
    00315
    00316       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);
    00317
    00318       printf("\tLe maillage par defaut est : |%s|, sur un total de : "IFORMAT" maillages associes\n",
    00319              maa_ass, nbrefmaa);
    00320
    00321       /* Le maillage reference est-il porte par un autre fichier */
    00322       if ( !local ) {
    00323
    00324         if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 )  {
    00325           MESSAGE("Erreur a la lecture de la taille du lien : ");
    00326           SSCRUTE(maa_ass);
    00327           ret = -1;
    00328         } else {
    00329
    00330           lien = malloc(lnsize*sizeof(char));
    00331           EXIT_IF(lien == NULL,NULL,NULL);
    00332
    00333           if ( MEDlienLire(fid, lien, maa_ass) < 0 )  {
    00334             MESSAGE("Erreur a la lecture du lien : ");
    00335             SSCRUTE(maa_ass);SSCRUTE(lien);
    00336             ret = -1;
    00337           } else {
    00338             printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
    00339           }
    00340           free(lien);
    00341         }
    00342       }
    00343
    00344       /* Combien de maillages lies aux (nomcha,ent,geo,numdt,numo)  */
    00345       /* Notons que cette information est egalement disponible a partir de MEDpasdetempsInfo */
    00346       if ( (nbrefmaa = MEDnChampRef(fid,nomcha,entite,type_geo[k],numdt,numo) ) < 0 ) {
    00347         MESSAGE("Erreur a la demande du nombre de maillages references par le champ : ");
    00348         SSCRUTE(nomcha);
    00349         ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
    00350         SSCRUTE(MED_ENTITE_MAILLAGE_AFF[(int)entite]);SSCRUTE(AFF[k]);
    00351         ISCRUTE(numdt); ISCRUTE(numo);
    00352         ret = -1; continue;
    00353       };
    00354
    00355       for (l=0;l<nbrefmaa;l++) {
    00356
    00357         if ( MEDchampRefInfo(fid,nomcha,entite,type_geo[k],
    00358                              l+1,numdt, numo, maa_ass, &local, &ngauss) <0 ) {
    00359           MESSAGE("Erreur a la demande d'information sur le maillage utilise par le champ n° : ");
    00360           SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
    00361           ISCRUTE_int(l+1);ISCRUTE(numdt); ISCRUTE(numo);SSCRUTE(maa_ass);
    00362           ret = -1; continue;
    00363         };
    00364
    00365         /* Prend en compte le nbre de pt de gauss automatiquement */
    00366         if ((nval = MEDnVal(fid,nomcha,entite,type_geo[k],numdt,numo,maa_ass,USER_MODE)) <= 0)   {
    00367           MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
    00368           SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
    00369           ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);ISCRUTE_int(USER_MODE);
    00370           ret = -1; continue;
    00371         };
    00372
    00373
    00374         printf("\t- Il y a "IFORMAT" valeurs en mode %i. Chaque entite %s\
    00375  de type geometrique %s associes au maillage |%s| a "IFORMAT" pts de gauss \n",
    00376                nval,USER_MODE,MED_ENTITE_MAILLAGE_AFF[(int)entite],AFF[k],maa_ass,ngauss);
    00377
    00378         /* Le maillage reference est-il porte par un autre fichier */
    00379         if ( !local ) {
    00380
    00381           if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 )  {
    00382             MESSAGE("Erreur a la lecture de la taille du lien : ");
    00383             SSCRUTE(maa_ass);
    00384             ret = -1;
    00385           } else {
    00386
    00387             lien = malloc(lnsize*sizeof(char) + 1);
    00388             EXIT_IF(lien == NULL,NULL,NULL);
    00389
    00390             if ( MEDlienLire(fid, lien, maa_ass) < 0 )  {
    00391               MESSAGE("Erreur a la lecture du lien : ");
    00392               SSCRUTE(maa_ass);SSCRUTE(lien);
    00393               ret = -1;
    00394             } else {
    00395               lien[lnsize] = '\0';
    00396               printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
    00397             }
    00398             free(lien);
    00399           }
    00400         }
    00401
    00402         /*Lecture des valeurs du champ */
    00403         if (typcha == MED_FLOAT64) {
    00404
    00405           valr = (med_float*) calloc(ncomp*nval,sizeof(med_float));
    00406           EXIT_IF(valr == NULL,NULL,NULL);
    00407           if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)valr,stockage,MED_ALL,locname,
    00408                             pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
    00409             MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
    00410             SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
    00411             ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
    00412             ret = -1;
    00413           };
    00414
    00415         } else {
    00416
    00417           vale = (med_int*) calloc(ncomp*nval,sizeof(med_int));
    00418           EXIT_IF(vale == NULL,NULL,NULL);
    00419           if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)vale,stockage,MED_ALL,locname,
    00420                             pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
    00421             MESSAGE("Erreur a la lecture des valeurs du champ : ");
    00422             SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
    00423             ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
    00424             ret = -1;
    00425           };
    00426
    00427         }
    00428
    00429         if ( ngauss > 1 )
    00430           printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
    00431
    00432         if (entite == MED_NOEUD_MAILLE)
    00433           ngroup = type_geo[k] % 100;
    00434         else
    00435           ngroup = ngauss;
    00436
    00437         switch (stockage) {
    00438
    00439         case MED_FULL_INTERLACE :
    00440           printf("\t- Valeurs :\n\t");
    00441           for (m=0;m<nval/ngroup;m++) {
    00442             printf("|");
    00443             for (n=0;n<ngroup*ncomp;n++)
    00444               if (typcha == MED_FLOAT64)
    00445                 printf(" %f ",*(valr+(m*ngroup*ncomp)+n));
    00446               else
    00447                 printf(" "IFORMAT" ",*(vale+(m*ngroup*ncomp)+n));
    00448
    00449           }
    00450           break;
    00451
    00452           /*Affichage en fonction du profil à traiter*/
    00453         case MED_NO_INTERLACE :
    00454           printf("\t- Valeurs :\n\t");
    00455           for (m=0;m<ncomp;m++) {
    00456             printf("|");
    00457             for (n=0;n<nval;n++)
    00458               if (typcha == MED_FLOAT64)
    00459                 printf(" %f ",*(valr+(m*nval)+n));
    00460               else
    00461                 printf(" "IFORMAT" ",*(vale+(m*nval)+n));
    00462           }
    00463           break;
    00464         }
    00465
    00466         printf("|\n");
    00467         if (typcha == MED_FLOAT64) {
    00468           if ( valr ) {free(valr);valr = NULL;}}
    00469         else
    00470           if (vale) { free(vale);vale = NULL; }
    00471
    00472         /*Lecture du profil associe */
    00473         if (strcmp(pflname,MED_NOPFL) == 0 )
    00474           printf("\t- Profil : MED_NOPFL\n");
    00475         else {
    00476
    00477           if ( (pflsize = MEDnValProfil(fid,pflname)) <0 )  {
    00478             MESSAGE("Erreur a la lecture du nombre de valeurs du profil : ");
    00479             SSCRUTE(pflname);
    00480             ret = -1; continue;
    00481           }
    00482
    00483           printf("\t- Profil : |%s| de taille "IFORMAT"\n",pflname,pflsize);
    00484
    00485           pflval = (med_int*) malloc(sizeof(med_int)*pflsize);
    00486           EXIT_IF(pflval == NULL,NULL,NULL);
    00487           if ( MEDprofilLire(fid,pflval,pflname) <0) {
    00488             MESSAGE("Erreur a la lecture des valeurs du profil : ");
    00489             SSCRUTE(pflname);
    00490             ret = -1;
    00491           }
    00492           printf("\t");
    00493           for (m=0;m<pflsize;m++) printf(" "IFORMAT" ",*(pflval+m));
    00494           printf("\n");
    00495           free(pflval);
    00496
    00497         }
    00498
    00499       }
    00500     }
    00501   } /* fin for sur les mailles*/
    00502
    00503   return ret;
    00504 }
    00505