# 2.3v3.0.old/test11.c

Aller à la documentation de ce fichier.

    00001 /*  This file is part of MED.
    00002  *
    00003  *  COPYRIGHT (C) 1999 - 2009  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_ACC_RDWR
    00032 #elif DEF_LECT_AJOUT
    00033 #define MODE_ACCES MED_ACC_RDEXT
    00034 #else
    00035 #define MODE_ACCES MED_ACC_CREAT
    00036 #endif
    00037
    00038 #ifndef USER_INTERLACE
    00039 #define USER_INTERLACE MED_FULL_INTERLACE
    00040 #endif
    00041
    00042 #define USER_MODE MED_COMPACT_PFLMODE
    00043
    00044 med_err getFieldsOn(med_idt fid, char * nomcha, med_field_type typcha, med_int ncomp,
    00045                     med_entity_type entite, med_switch_mode stockage, med_int ncstp);
    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_NAME_SIZE+1]="";
    00055   char desc[MED_COMMENT_SIZE+1]="";
    00056   char pflname[MED_NAME_SIZE+1]="",nomlien[MED_NAME_SIZE+1]="";
    00057   char _meshname [MED_NAME_SIZE+1]="";
    00058   char _dtunit [MED_SNAME_SIZE+1]="";
    00059   char locname[MED_NAME_SIZE+1]="";
    00060   char * lien = NULL;
    00061   char *comp= NULL, *unit= NULL;
    00062   char nomcha  [MED_NAME_SIZE+1]="";
    00063   med_int mdim=0,sdim=0,ncomp,ncha,npro,nln,pflsize,*pflval,nval;
    00064   med_int _ncstp=0,ngauss=0,nloc=0,locsdim=0,lnsize=0;
    00065   int t1,t2,t3;
    00066   med_field_type    typcha;
    00067   med_geometry_type type_geo;
    00068   med_float *refcoo, *gscoo, *wg;
    00069   int i,j;
    00070   med_bool _local;
    00071
    00072   char dtunit[MED_SNAME_SIZE+1]="";
    00073   char nomcoo[3*MED_SNAME_SIZE+1]="";
    00074   char unicoo[3*MED_SNAME_SIZE+1]="";
    00075   char geointerpname[MED_SNAME_SIZE+1]="";
    00076   char ipointstructmeshname[MED_SNAME_SIZE+1]="";
    00077   med_mesh_type type;
    00078   med_sorting_type sort;
    00079   med_int nstep=0;
    00080   med_axis_type rep;
    00081   if (argc != 2) {
    00082     MESSAGE("Aucun nom de fichier precise, fichier test10.med utilise ");
    00083     fichier = "test10.med";
    00084   } else {
    00085     fichier = argv[1];
    00086   };
    00087
    00088
    00089   /* Ouverture du fichier med */
    00090   if ((fid = MEDfileOpen(fichier,MED_ACC_RDONLY)) < 0){
    00091     MESSAGE("Erreur a l'ouverture du fichier : ");SSCRUTE(fichier);
    00092     return -1;
    00093   }
    00094
    00095   ret = 0;
    00096
    00097
    00098  /* Lecture des infos concernant le premier maillage */
    00099   if ( MEDmeshInfo( fid, 1,  maa, &sdim, &mdim, &type, desc, dtunit, &sort,
    00100                     &nstep,  &rep, nomcoo,unicoo) < 0 ) {
    00101     MESSAGE("Erreur a la lecture des informations sur le maillage : ");SSCRUTE(maa);
    00102     return -1;
    00103   } else {
    00104     printf("Maillage de nom : |%s| , de dimension : %d , et de type %d\n",maa,mdim,type);
    00105     printf("\t -Dimension de l'espace : %d\n",sdim);
    00106     printf("\t -Description du maillage : %s\n",desc);
    00107     printf("\t -Noms des axes : |%s|\n",nomcoo);
    00108     printf("\t -Unités des axes : |%s|\n",unicoo);
    00109     printf("\t -Type de repère : %d\n",rep);
    00110     printf("\t -Nombre d'étapes de calcul : %d\n",nstep);
    00111     printf("\t -Unité des dates : |%s|\n",dtunit);
    00112   }
    00113
    00114
    00115   /* combien de champs dans le fichier */
    00116   if ((ncha = MEDnField(fid)) < 0) {
    00117     MESSAGE("Impossible de lire le nombre de champs : ");ISCRUTE(ncha);
    00118     return ncha;
    00119   }
    00120
    00121   printf("Nombre de champs : "IFORMAT" \n",ncha);
    00122
    00123   /* lecture de tous les champs  */
    00124   for (i =0;i<ncha;i++) {
    00125     lret = 0;
    00126     printf("\nChamp numero : %d \n",i+1);
    00127
    00128     /* Lecture du nombre de composantes */
    00129     if ((ncomp = MEDfieldnComponent(fid,i+1)) < 0) {
    00130       MESSAGE("Erreur a la lecture du nombre de composantes : "); ISCRUTE(ncomp);
    00131       ret = -1; continue;
    00132     }
    00133
    00134     /* Lecture du type du champ, des noms des composantes et du nom de l'unite*/
    00135     comp = (char*) malloc(ncomp*MED_SNAME_SIZE+1);
    00136     EXIT_IF(comp == NULL,NULL,NULL);
    00137     unit = (char*) malloc(ncomp*MED_SNAME_SIZE+1);
    00138     EXIT_IF(unit == NULL,NULL,NULL);
    00139
    00140     if ( MEDfieldInfo(fid,i+1,nomcha,_meshname,&_local,&typcha,comp,unit,_dtunit,&_ncstp) < 0 ) {
    00141       MESSAGE("Erreur a la demande d'information sur les champs : ");
    00142       ISCRUTE_int(i+1);SSCRUTE(nomcha);ISCRUTE_int(typcha);SSCRUTE(comp);SSCRUTE(unit);
    00143       ISCRUTE(ncomp);
    00144       ret = -1; continue;
    00145     }
    00146
    00147
    00148     printf("Nom du champ : |%s| de type %d\n",nomcha,typcha);
    00149     printf("Nom des composantes : |%s|\n",comp);
    00150     printf("Unites des composantes : |%s| \n",unit);
    00151     printf("Unites des dates  : |%s| \n",_dtunit);
    00152     printf("Le maillage associé est |%s|\n",_meshname);
    00153     printf("Nombre de séquences de calcul |%d|\n",_ncstp);
    00154
    00155       /* Le maillage reference est-il porte par un autre fichier */
    00156     if ( !_local ) {
    00157
    00158       if ( (lnsize=MEDlinkInfoByName(fid,_meshname) ) < 0 )  {
    00159         MESSAGE("Erreur a la lecture de la taille du lien : ");
    00160         SSCRUTE(_meshname);
    00161         ret = -1;
    00162       } else {
    00163
    00164           lien = malloc(lnsize*sizeof(char));
    00165           EXIT_IF(lien == NULL,NULL,NULL);
    00166
    00167           if ( MEDlinkRd(fid,_meshname, lien) < 0 )  {
    00168             MESSAGE("Erreur a la lecture du lien : ");
    00169             SSCRUTE(_meshname);SSCRUTE(lien);
    00170             ret = -1;
    00171           } else {
    00172             printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",_meshname,lien);
    00173           }
    00174           free(lien);
    00175         }
    00176       }
    00177
    00178     free(comp);
    00179     free(unit);
    00180
    00181     /*TODO : Créer les API30 spécifiques pour la lecture de champs multi-maillages*/
    00182     if (strcmp(nomcha,"champ entier")) {
    00183       MESSAGE("There is no API yest for reading field on multiple meshes"); continue;
    00184     }
    00185
    00186     lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NODE, USER_INTERLACE,_ncstp );
    00187
    00188     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_CELL, USER_INTERLACE,_ncstp );
    00189     else { MESSAGE("Erreur a la lecture des champs aux noeuds "); ret = -1; continue;}
    00190
    00191     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_DESCENDING_FACE,USER_INTERLACE,_ncstp);
    00192     else { MESSAGE("Erreur a la lecture des champs aux mailles "); ret = -1; continue;}
    00193
    00194     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_DESCENDING_EDGE,USER_INTERLACE,_ncstp);
    00195     else {MESSAGE("Erreur a la lecture des champs aux faces "); ret = -1; continue;}
    00196
    00197     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NODE_ELEMENT,USER_INTERLACE,_ncstp);
    00198     else {MESSAGE("Erreur a la lecture des champs aux aretes"); ret = -1; continue;}
    00199
    00200     /*TODO */
    00201 /*     if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_STRUCT_ELEMENT,USER_INTERLACE,_ncstp); */
    00202 /*     else {MESSAGE("Erreur a la lecture des champs aux éléments de structure"); ret = -1; continue;} */
    00203
    00204     if  (lret != 0) {MESSAGE("Erreur a la lecture des champs aux noeuds des mailles "); ret = -1;};
    00205   }
    00206
    00207
    00208   /* Interrogation des profils */
    00209   npro = MEDnProfile(fid);
    00210
    00211   printf("\nNombre de profils stockes : "IFORMAT"\n\n",npro);
    00212   for (i=1 ; i <= npro ; i++ ) {
    00213     if ( MEDprofileInfo(fid, i, pflname, &nval) < 0)  {
    00214       MESSAGE("Erreur a la demande d'information sur le profil n° : "); ISCRUTE_int(i);
    00215       ret = -1;continue;
    00216     }
    00217     printf("\t- Profil n°%i de nom |%s| et de taille "IFORMAT"\n",i,pflname,nval);
    00218     pflval = (med_int*) malloc(sizeof(med_int)*nval);
    00219     if ( MEDprofileRd(fid, pflname, pflval) < 0) {
    00220       MESSAGE("Erreur a la lecture des valeurs du profil : ");
    00221       SSCRUTE(pflname);
    00222       ret = -1;
    00223     } else {
    00224       printf("\t");
    00225       for (j=0;j<nval;j++) printf(" "IFORMAT" ",*(pflval+j));
    00226       printf("\n\n");
    00227     }
    00228     free(pflval);
    00229   }
    00230
    00231   /* Interrogation des liens */
    00232   nln = MEDnLink(fid);
    00233
    00234   printf("\nNombre de liens stockes : "IFORMAT"\n\n",nln);
    00235   for (i=1 ; i <= nln ; i++ ) {
    00236     if ( MEDlinkInfo(fid, i, nomlien, &nval) < 0)  {
    00237       MESSAGE("Erreur a la demande d'information sur le lien n° : "); ISCRUTE_int(i);
    00238       ret = -1;continue;
    00239     }
    00240     printf("\t- Lien n°%i de nom |%s| et de taille "IFORMAT"\n",i,nomlien,nval);
    00241
    00242     lien = malloc((nval+1)*sizeof(char));
    00243     EXIT_IF(lien == NULL,NULL,NULL);
    00244
    00245     if ( MEDlinkRd(fid, nomlien, lien ) < 0 )  {
    00246       MESSAGE("Erreur a la lecture du lien : ");
    00247       SSCRUTE(nomlien);SSCRUTE(lien);
    00248       ret = -1;
    00249     } else {
    00250       lien[nval] = '\0';
    00251       printf("\t\t|%s|\n\n",lien);
    00252     }
    00253     free(lien);
    00254   }
    00255
    00256   /* Interrogation des localisations des points de GAUSS */
    00257   nloc = MEDnLocalization(fid);
    00258
    00259   printf("\nNombre de localisations stockees : "IFORMAT"\n\n",nloc);
    00260   for (i=1 ; i <= nloc ; i++ ) {
    00261     if ( MEDlocalizationInfo(fid, i, locname, &type_geo, &locsdim,&ngauss,
    00262                              geointerpname, ipointstructmeshname) < 0)  {
    00263       MESSAGE("Erreur a la demande d'information sur la localisation n° : "); ISCRUTE_int(i);
    00264       ret = -1;continue;
    00265     }
    00266     printf("\t- Loc. n°%i de nom |%s| de dimension %i avec "IFORMAT" pts de GAUSS \n",i,locname,locsdim,ngauss);
    00267     t1 = (type_geo%100)*(type_geo/100);
    00268     t2 = ngauss*(type_geo/100);
    00269     t3 = ngauss;
    00270     refcoo = (med_float *) malloc(sizeof(med_float)*t1 );
    00271     gscoo  = (med_float *) malloc(sizeof(med_float)*t2 );
    00272     wg     = (med_float *) malloc(sizeof(med_float)*t3 );
    00273
    00274     if ( MEDlocalizationRd(fid, locname, USER_INTERLACE, refcoo, gscoo, wg  ) < 0) {
    00275       MESSAGE("Erreur a la lecture des valeurs de la localisation : ");
    00276       SSCRUTE(locname);
    00277       ret = -1;
    00278     } else {
    00279       printf("\t  Coordonnees de l'element de reference de type %i :\n\t\t",type_geo);
    00280       for (j=0;j<t1;j++) printf(" %f ",*(refcoo+j));
    00281       printf("\n");
    00282       printf("\t  Localisation des points de GAUSS : \n\t\t");
    00283       for (j=0;j<t2;j++) printf(" %f ",*(gscoo+j));
    00284       printf("\n");
    00285       printf("\t  Poids associes aux points de GAUSS :\n\t\t");
    00286       for (j=0;j<t3;j++) printf(" %f ",*(wg+j));
    00287       printf("\n\n");
    00288     }
    00289     free(refcoo);
    00290     free(gscoo);
    00291     free(wg);
    00292   }
    00293
    00294
    00295   /* fermeture du fichier */
    00296   if ( MEDfileClose(fid) < 0) return -1;
    00297
    00298   return ret;
    00299 }
    00300
    00301 med_err getFieldsOn(med_idt fid, char * nomcha, med_field_type typcha, med_int ncomp,
    00302                     med_entity_type entite, med_switch_mode stockage, med_int ncstp) {
    00303
    00304   int j,k,l,m,n,nb_geo=0;
    00305   med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,ngroup,*vale=NULL,nval;
    00306   med_int numdt=0,numo=0,_nbofprofile;
    00307   med_float *valr=NULL,dt=0.0;
    00308   med_err ret=0;
    00309   char pflname [MED_NAME_SIZE+1]="";
    00310   char locname [MED_NAME_SIZE+1]="";
    00311   char * lien = NULL;
    00312   char dt_unit [MED_SNAME_SIZE+1]="unknown";
    00313
    00314
    00315   med_geometry_type * type_geo;
    00316 /*   med_geometry_type typ_noeud[1] = { MED_NONE }; */
    00317 /*   med_geometry_type typmai[MED_NBR_GEOMETRIE_MAILLE+2] = {MED_POINT1, MED_SEG2, MED_SEG3, MED_TRIA3, */
    00318 /*                                                            MED_QUAD4, MED_TRIA6,MED_QUAD8, MED_TETRA4, */
    00319 /*                                                            MED_PYRA5, MED_PENTA6, MED_HEXA8, MED_TETRA10, */
    00320 /*                                                            MED_PYRA13, MED_PENTA15, MED_HEXA20, */
    00321 /*                                                            MED_POLYGON, MED_POLYHEDRON}; */
    00322 /*   med_geometry_type typfac[MED_NBR_GEOMETRIE_FACE+1] = {MED_TRIA3,MED_TRIA6, */
    00323 /*                                                          MED_QUAD4,MED_QUAD8, */
    00324 /*                                                          MED_POLYGON}; */
    00325 /*   med_geometry_type typare[MED_NBR_GEOMETRIE_ARETE] = {MED_SEG2,MED_SEG3}; */
    00326
    00327   const char * const * AFF;
    00328   const char * const * AFF_ENT=MED_GET_ENTITY_TYPENAME+1;
    00329   switch (entite) {
    00330   case MED_NODE :
    00331     type_geo = MED_GET_NODE_GEOMETRY_TYPE+1;
    00332     nb_geo   = 1;
    00333     AFF      = MED_GET_NODE_GEOMETRY_TYPENAME+1;
    00334     break;
    00335   case  MED_CELL :
    00336   case  MED_NODE_ELEMENT :
    00337     type_geo = MED_GET_CELL_GEOMETRY_TYPE+1;
    00338     nb_geo   = MED_NOF_CELL_TYPES;
    00339     AFF      = MED_GET_CELL_GEOMETRY_TYPENAME+1;
    00340     break;
    00341   case  MED_DESCENDING_FACE :
    00342     type_geo = MED_GET_FACE_GEOMETRY_TYPE+1;
    00343     nb_geo   = MED_NOF_FACE_TYPES;
    00344     AFF      = MED_GET_FACE_GEOMETRY_TYPENAME+1;
    00345     break;
    00346   case  MED_DESCENDING_EDGE :
    00347     type_geo = MED_GET_EDGE_GEOMETRY_TYPE+1;
    00348     nb_geo   = MED_NOF_EDGE_TYPES;
    00349     AFF      = MED_GET_EDGE_GEOMETRY_TYPENAME+1;
    00350     break;
    00351   }
    00352
    00353   for (k=0;k<nb_geo;k++) {
    00354
    00355     /* Combien de (PDT,NOR) a lire */
    00356     nbpdtnor = ncstp;
    00357     if (nbpdtnor < 1 ) continue;
    00358
    00359     for (j=0;j<nbpdtnor;j++) {
    00360
    00361       if ( MEDfieldComputingStepInfo(fid,nomcha,j+1, &numdt, &numo, &dt ) <0) {
    00362         MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
    00363         ISCRUTE(numdt); ISCRUTE(numo);
    00364         ret = -1; continue;
    00365       }
    00366
    00367       if ( (_nbofprofile = MEDfieldnProfile(fid,nomcha,numdt,numo,entite,type_geo[k],
    00368                                        &pflname,&locname   ) ) < 0 ) {
    00369         MESSAGE("Erreur a la demande du nombre de profils referencés par le champ : ");
    00370         SSCRUTE(nomcha);
    00371         ISCRUTE(numdt); ISCRUTE(numo);
    00372         ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
    00373         SSCRUTE(AFF_ENT[(int)entite]);SSCRUTE(AFF[k]);
    00374         ret = -1; continue;
    00375       };
    00376
    00377       for (l=0;l<_nbofprofile;l++) {
    00378
    00379
    00380         if ( (nval = MEDfieldnValueWithProfile(fid, nomcha, numdt, numo, entite, type_geo[k],
    00381                                                l+1,  USER_MODE, pflname,&pflsize,
    00382                                                locname, &ngauss) ) < 0 ) {
    00383           MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
    00384           SSCRUTE(nomcha);ISCRUTE(numdt);ISCRUTE(numo);
    00385           ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
    00386           ISCRUTE_int(USER_MODE);
    00387           ret = -1; continue;
    00388         };
    00389
    00390         printf("\n  +Pas de Temps n."IFORMAT" (%f) [%s], n. d'ordre "IFORMAT", avec "IFORMAT" valeur(s) par entité.\n",numdt,dt,dt_unit,numo,ngauss);
    00391         printf("\t- Il y a "IFORMAT" entités qui portent des valeurs en mode %i. Chaque entite %s\
    00392  de type geometrique %s associes au profile |%s| a "IFORMAT" valeurs associées \n",
    00393                nval,USER_MODE,AFF_ENT[(int)entite],AFF[k],pflname,ngauss);
    00394
    00395         /*Lecture des valeurs du champ */
    00396         if (typcha == MED_FLOAT64) {
    00397
    00398
    00399           valr = (med_float*) calloc(ncomp*nval*ngauss,sizeof(med_float));
    00400           EXIT_IF(valr == NULL,NULL,NULL);
    00401
    00402           if (MEDfieldValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],
    00403                                          USER_MODE, pflname, stockage,MED_ALL_CONSTITUENT,
    00404                                          (unsigned char*) valr) < 0 ) {
    00405             MESSAGE("Erreur a la lecture des valeurs du champ : ");
    00406             SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
    00407             ISCRUTE(numdt);ISCRUTE(numo);
    00408             ret = -1;
    00409           }
    00410
    00411         } else {
    00412
    00413           vale = (med_int*) calloc(ncomp*nval*ngauss,sizeof(med_int));
    00414           EXIT_IF(vale == NULL,NULL,NULL);
    00415
    00416           if (MEDfieldValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],
    00417                                          USER_MODE, pflname, stockage,MED_ALL_CONSTITUENT,
    00418                                          (unsigned char*) vale) < 0 ) {
    00419             MESSAGE("Erreur a la lecture des valeurs du champ : ");
    00420             SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
    00421             ISCRUTE(numdt);ISCRUTE(numo);
    00422             ret = -1;
    00423           };
    00424
    00425
    00426         }
    00427
    00428         if ( strlen(locname) )
    00429           printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
    00430
    00431         if (entite == MED_NODE_ELEMENT)
    00432           ngroup = (type_geo[k] % 100);
    00433         else
    00434           ngroup = ngauss;
    00435
    00436         switch (stockage) {
    00437
    00438         case MED_FULL_INTERLACE :
    00439           printf("\t- Valeurs :\n\t");
    00440           for (m=0;m<(nval*ngauss)/ngroup;m++) {
    00441             printf("|");
    00442             for (n=0;n<ngroup*ncomp;n++)
    00443               if (typcha == MED_FLOAT64)
    00444                 printf(" %f ",*(valr+(m*ngroup*ncomp)+n));
    00445               else
    00446                 printf(" "IFORMAT" ",*(vale+(m*ngroup*ncomp)+n));
    00447
    00448           }
    00449           break;
    00450
    00451           /*Affichage en fonction du profil à traiter*/
    00452         case MED_NO_INTERLACE :
    00453           printf("\t- Valeurs :\n\t");
    00454           for (m=0;m<ncomp;m++) {
    00455             printf("|");
    00456             for (n=0;n<(nval*ngauss);n++)
    00457               if (typcha == MED_FLOAT64)
    00458                 printf(" %f ",*(valr+(m*nval)+n));
    00459               else
    00460                 printf(" "IFORMAT" ",*(vale+(m*nval)+n));
    00461           }
    00462           break;
    00463         }
    00464
    00465         printf("|\n");
    00466         if (typcha == MED_FLOAT64) {
    00467           if ( valr ) {free(valr);valr = NULL;}}
    00468         else
    00469           if (vale) { free(vale);vale = NULL; }
    00470
    00471         /*Lecture du profil associe */
    00472         if (strcmp(pflname,MED_NO_PROFILE) == 0 )
    00473           printf("\t- Profil : MED_NO_PROFILE\n");
    00474         else {
    00475
    00476           if ( (pflsize = MEDprofileSizeByName(fid,pflname)) <0 )  {
    00477             MESSAGE("Erreur a la lecture du nombre de valeurs du profil : ");
    00478             SSCRUTE(pflname);
    00479             ret = -1; continue;
    00480           }
    00481
    00482           printf("\t- Profil : |%s| de taille "IFORMAT"\n",pflname,pflsize);
    00483
    00484           pflval = (med_int*) malloc(sizeof(med_int)*pflsize);
    00485           EXIT_IF(pflval == NULL,NULL,NULL);
    00486           if ( MEDprofileRd(fid,pflname,pflval) <0) {
    00487             MESSAGE("Erreur a la lecture des valeurs du profil : ");
    00488             SSCRUTE(pflname);
    00489             ret = -1;
    00490           }
    00491           printf("\t");
    00492           for (m=0;m<pflsize;m++) printf(" "IFORMAT" ",*(pflval+m));
    00493           printf("\n");
    00494           free(pflval);
    00495
    00496         }
    00497
    00498       }
    00499     }
    00500   } /* fin for sur les mailles*/
    00501
    00502   return ret;
    00503 }
    00504