test11.c

Aller à la documentation de ce fichier.
00001 /*  This file is part of MED.
00002  *
00003  *  COPYRIGHT (C) 1999 - 2015  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_NAME_SIZE+1]="";
00076   char ipointstructmeshname[MED_NAME_SIZE+1]="";
00077   med_mesh_type type;
00078   med_sorting_type sort;
00079   med_int nstep=0;
00080   med_axis_type rep;
00081   med_int nsectionmeshcell;
00082   med_geometry_type sectiongeotype;
00083 
00084   if (argc != 2) {
00085     MESSAGE("Aucun nom de fichier precise, fichier test10.med utilise ");
00086     fichier = "test10.med";
00087   } else {
00088     fichier = argv[1];
00089   };
00090 
00091 
00092   /* Ouverture du fichier med */
00093   if ((fid = MEDfileOpen(fichier,MED_ACC_RDONLY)) < 0){
00094     MESSAGE("Erreur a l'ouverture du fichier : ");SSCRUTE(fichier);
00095     return -1;
00096   }
00097 
00098   ret = 0;
00099 
00100 
00101  /* Lecture des infos concernant le premier maillage */
00102   if ( MEDmeshInfo( fid, 1,  maa, &sdim, &mdim, &type, desc, dtunit, &sort,
00103                     &nstep,  &rep, nomcoo,unicoo) < 0 ) {
00104     MESSAGE("Erreur a la lecture des informations sur le maillage : ");SSCRUTE(maa);
00105     return -1;
00106   } else {
00107     printf("Maillage de nom : |%s| , de dimension : "IFORMAT" , et de type %d\n",maa,mdim,type);
00108     printf("\t -Dimension de l'espace : "IFORMAT"\n",sdim);
00109     printf("\t -Description du maillage : %s\n",desc);
00110     printf("\t -Noms des axes : |%s|\n",nomcoo);
00111     printf("\t -Unités des axes : |%s|\n",unicoo);
00112     printf("\t -Type de repère : %d\n",rep);
00113     printf("\t -Nombre d'étapes de calcul : "IFORMAT"\n",nstep);
00114     printf("\t -Unité des dates : |%s|\n",dtunit);
00115   }
00116 
00117 
00118   /* combien de champs dans le fichier */
00119   if ((ncha = MEDnField(fid)) < 0) {
00120     MESSAGE("Impossible de lire le nombre de champs : ");ISCRUTE(ncha);
00121     return ncha;
00122   }
00123 
00124   printf("Nombre de champs : "IFORMAT" \n",ncha);
00125 
00126   /* lecture de tous les champs  */
00127   for (i =0;i<ncha;i++) {
00128     lret = 0;
00129     printf("\nChamp numero : %d \n",i+1);
00130 
00131     /* Lecture du nombre de composantes */
00132     if ((ncomp = MEDfieldnComponent(fid,i+1)) < 0) {
00133       MESSAGE("Erreur a la lecture du nombre de composantes : "); ISCRUTE(ncomp);
00134       ret = -1; continue;
00135     }
00136 
00137     /* Lecture du type du champ, des noms des composantes et du nom de l'unite*/
00138     comp = (char*) malloc(ncomp*MED_SNAME_SIZE+1);
00139     EXIT_IF(comp == NULL,NULL,NULL);
00140     unit = (char*) malloc(ncomp*MED_SNAME_SIZE+1);
00141     EXIT_IF(unit == NULL,NULL,NULL);
00142 
00143     if ( MEDfieldInfo(fid,i+1,nomcha,_meshname,&_local,&typcha,comp,unit,_dtunit,&_ncstp) < 0 ) {
00144       MESSAGE("Erreur a la demande d'information sur les champs : ");
00145       ISCRUTE_int(i+1);SSCRUTE(nomcha);ISCRUTE_int(typcha);SSCRUTE(comp);SSCRUTE(unit);
00146       ISCRUTE(ncomp);
00147       ret = -1; continue;
00148     }
00149 
00150 
00151     printf("Nom du champ : |%s| de type %d\n",nomcha,typcha);
00152     printf("Nombre de composantes : |"IFORMAT"|\n",ncomp);
00153     printf("Nom des composantes : |%s|\n",comp);
00154     printf("Unites des composantes : |%s| \n",unit);
00155     printf("Unites des dates  : |%s| \n",_dtunit);
00156     printf("Le maillage associé est |%s|\n",_meshname);
00157     printf("Nombre de séquences de calcul |"IFORMAT"|\n",_ncstp);
00158 
00159       /* Le maillage reference est-il porte par un autre fichier */
00160     if ( !_local ) {
00161 
00162       if ( (lnsize=MEDlinkInfoByName(fid,_meshname) ) < 0 )  {
00163         MESSAGE("Erreur a la lecture de la taille du lien : ");
00164         SSCRUTE(_meshname);
00165         ret = -1;
00166       } else {
00167 
00168         lien = malloc((lnsize+1)*sizeof(char));
00169           EXIT_IF(lien == NULL,NULL,NULL);
00170 
00171           if ( MEDlinkRd(fid,_meshname, lien) < 0 )  {
00172             MESSAGE("Erreur a la lecture du lien : ");
00173             SSCRUTE(_meshname);SSCRUTE(lien);
00174             ret = -1;
00175           } else {
00176             printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",_meshname,lien);
00177           }
00178           free(lien);
00179         }
00180       }
00181     
00182     free(comp);
00183     free(unit);
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,&nsectionmeshcell,
00263                              &sectiongeotype) < 0)  {
00264       MESSAGE("Erreur a la demande d'information sur la localisation n° : "); ISCRUTE_int(i);
00265       ret = -1;continue;
00266     }
00267     printf("\t- Loc. n°%i de nom |%s| de dimension |"IFORMAT"| avec "IFORMAT" pts de GAUSS \n",i,locname,locsdim,ngauss);
00268     t1 = (type_geo%100)*(type_geo/100);
00269     t2 = ngauss*(type_geo/100);
00270     t3 = ngauss;
00271     refcoo = (med_float *) malloc(sizeof(med_float)*t1 );
00272     gscoo  = (med_float *) malloc(sizeof(med_float)*t2 );
00273     wg     = (med_float *) malloc(sizeof(med_float)*t3 );
00274 
00275     if ( MEDlocalizationRd(fid, locname, USER_INTERLACE, refcoo, gscoo, wg  ) < 0) {
00276       MESSAGE("Erreur a la lecture des valeurs de la localisation : ");
00277       SSCRUTE(locname);
00278       ret = -1;
00279     } else {
00280       printf("\t  Coordonnees de l'element de reference de type %i :\n\t\t",type_geo);
00281       for (j=0;j<t1;j++) printf(" %f ",*(refcoo+j));
00282       printf("\n");
00283       printf("\t  Localisation des points de GAUSS : \n\t\t");
00284       for (j=0;j<t2;j++) printf(" %f ",*(gscoo+j));
00285       printf("\n");
00286       printf("\t  Poids associes aux points de GAUSS :\n\t\t");
00287       for (j=0;j<t3;j++) printf(" %f ",*(wg+j));
00288       printf("\n\n");
00289     }
00290     free(refcoo);
00291     free(gscoo);
00292     free(wg);
00293   }
00294 
00295 
00296   /* fermeture du fichier */
00297   if ( MEDfileClose(fid) < 0) return -1;
00298 
00299   return ret;
00300 }
00301 
00302 med_err getFieldsOn(med_idt fid, char * nomcha, med_field_type typcha, med_int ncomp,
00303                     med_entity_type entite, med_switch_mode stockage, med_int ncstp) {
00304 
00305   int j,k,l,m,n,nb_geo=0;
00306   med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,ngroup,*vale=NULL,nval;
00307   med_int numdt=0,numo=0,_nprofile;
00308   med_float *valr=NULL,dt=0.0;
00309   med_err ret=0;
00310   char pflname [MED_NAME_SIZE+1]="";
00311   char locname [MED_NAME_SIZE+1]="";
00312   char * lien = NULL;
00313   char dt_unit [MED_SNAME_SIZE+1]="unknown";
00314 
00315 
00316   med_geometry_type * type_geo;
00317 
00318   const char * const * AFF;
00319   const char * const * AFF_ENT=MED_GET_ENTITY_TYPENAME+1;
00320   switch (entite) {
00321   case MED_NODE :
00322     type_geo = MED_GET_NODE_GEOMETRY_TYPE;
00323     nb_geo   = MED_N_NODE_FIXED_GEO;
00324     AFF      = MED_GET_NODE_GEOMETRY_TYPENAME;
00325     break;
00326   case  MED_CELL :
00327   case  MED_NODE_ELEMENT :
00328     type_geo = MED_GET_CELL_GEOMETRY_TYPE;
00329     nb_geo   = MED_N_CELL_FIXED_GEO;
00330     AFF      = MED_GET_CELL_GEOMETRY_TYPENAME;
00331     break;
00332   case  MED_DESCENDING_FACE :
00333     type_geo = MED_GET_FACE_GEOMETRY_TYPE;
00334     nb_geo   = MED_N_FACE_FIXED_GEO;
00335     AFF      = MED_GET_FACE_GEOMETRY_TYPENAME;
00336     break;
00337   case  MED_DESCENDING_EDGE :
00338     type_geo = MED_GET_EDGE_GEOMETRY_TYPE;
00339     nb_geo   = MED_N_EDGE_FIXED_GEO;
00340     AFF      = MED_GET_EDGE_GEOMETRY_TYPENAME;
00341     break;
00342   }
00343 
00344   for (k=1;k<=nb_geo;k++) {
00345 
00346     /* Combien de (PDT,NOR) a lire */
00347     nbpdtnor = ncstp;
00348     if (nbpdtnor < 1 ) continue;
00349 
00350     for (j=0;j<nbpdtnor;j++) {
00351 
00352       if ( MEDfieldComputingStepInfo(fid,nomcha,j+1, &numdt, &numo, &dt ) <0) {
00353         MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
00354         ISCRUTE(numdt); ISCRUTE(numo);
00355         ret = -1; continue;
00356       }
00357 
00358       if ( (_nprofile = MEDfieldnProfile(fid,nomcha,numdt,numo,entite,type_geo[k],
00359                                             pflname,locname   ) ) < 0 ) {
00360         MESSAGE("Erreur a la demande du nombre de profils referencés par le champ : ");
00361         SSCRUTE(nomcha);
00362         ISCRUTE(numdt); ISCRUTE(numo);
00363         ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
00364         SSCRUTE(AFF_ENT[(int)entite]);SSCRUTE(AFF[k]);
00365         ret = -1; continue;
00366       };
00367 
00368       for (l=0;l<_nprofile;l++) {
00369 
00370 
00371         if ( (nval = MEDfieldnValueWithProfile(fid, nomcha, numdt, numo, entite, type_geo[k],
00372                                                l+1,  USER_MODE, pflname,&pflsize,
00373                                                locname, &ngauss) ) < 0 ) {
00374           MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
00375           SSCRUTE(nomcha);ISCRUTE(numdt);ISCRUTE(numo);
00376           ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
00377           ISCRUTE_int(USER_MODE);
00378           ret = -1; continue;
00379         };
00380 
00381         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);
00382         printf("\t- Il y a "IFORMAT" entités qui portent des valeurs en mode %i. Chaque entite %s\
00383  de type geometrique %s associes au profile |%s| a "IFORMAT" valeurs associées \n",
00384                nval,USER_MODE,AFF_ENT[(int)entite],AFF[k],pflname,ngauss);
00385 
00386         /*Lecture des valeurs du champ */
00387         if (typcha == MED_FLOAT64) {
00388 
00389 
00390           valr = (med_float*) calloc(ncomp*nval*ngauss,sizeof(med_float));
00391           EXIT_IF(valr == NULL,NULL,NULL);
00392 
00393           if (MEDfieldValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],
00394                                          USER_MODE, pflname, stockage,MED_ALL_CONSTITUENT,
00395                                          (unsigned char*) valr) < 0 ) {
00396             MESSAGE("Erreur a la lecture des valeurs du champ : ");
00397             SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
00398             ISCRUTE(numdt);ISCRUTE(numo);
00399             ret = -1;
00400           }
00401 
00402         } else {
00403 
00404           vale = (med_int*) calloc(ncomp*nval*ngauss,sizeof(med_int));
00405           EXIT_IF(vale == NULL,NULL,NULL);
00406 
00407           if (MEDfieldValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],
00408                                          USER_MODE, pflname, stockage,MED_ALL_CONSTITUENT,
00409                                          (unsigned char*) vale) < 0 ) {
00410             MESSAGE("Erreur a la lecture des valeurs du champ : ");
00411             SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
00412             ISCRUTE(numdt);ISCRUTE(numo);
00413             ret = -1;
00414           };
00415 
00416 
00417         }
00418 
00419         if ( strlen(locname) )
00420           printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
00421 
00422         if (entite == MED_NODE_ELEMENT)
00423           ngroup = (type_geo[k] % 100);
00424         else
00425           ngroup = ngauss;
00426 
00427         switch (stockage) {
00428 
00429         case MED_FULL_INTERLACE :
00430           printf("\t- Valeurs :\n\t");
00431           for (m=0;m<(nval*ngauss)/ngroup;m++) {
00432             printf("|");
00433             for (n=0;n<ngroup*ncomp;n++)
00434               if (typcha == MED_FLOAT64)
00435                 printf(" %f ",*(valr+(m*ngroup*ncomp)+n));
00436               else
00437                 printf(" "IFORMAT" ",*(vale+(m*ngroup*ncomp)+n));
00438 
00439           }
00440           break;
00441 
00442           /*Affichage en fonction du profil à traiter*/
00443         case MED_NO_INTERLACE :
00444           printf("\t- Valeurs :\n\t");
00445           for (m=0;m<ncomp;m++) {
00446             printf("|");
00447             for (n=0;n<(nval*ngauss);n++)
00448               if (typcha == MED_FLOAT64)
00449                 printf(" %f ",*(valr+(m*nval)+n));
00450               else
00451                 printf(" "IFORMAT" ",*(vale+(m*nval)+n));
00452           }
00453           break;
00454         }
00455 
00456         printf("|\n");
00457         if (typcha == MED_FLOAT64) {
00458           if ( valr ) {free(valr);valr = NULL;}}
00459         else
00460           if (vale) { free(vale);vale = NULL; }
00461 
00462         /*Lecture du profil associe */
00463         if (strcmp(pflname,MED_NO_PROFILE) == 0 )
00464           printf("\t- Profil : MED_NO_PROFILE\n");
00465         else {
00466 
00467           if ( (pflsize = MEDprofileSizeByName(fid,pflname)) <0 )  {
00468             MESSAGE("Erreur a la lecture du nombre de valeurs du profil : ");
00469             SSCRUTE(pflname);
00470             ret = -1; continue;
00471           }
00472 
00473           printf("\t- Profil : |%s| de taille "IFORMAT"\n",pflname,pflsize);
00474 
00475           pflval = (med_int*) malloc(sizeof(med_int)*pflsize);
00476           EXIT_IF(pflval == NULL,NULL,NULL);
00477           if ( MEDprofileRd(fid,pflname,pflval) <0) {
00478             MESSAGE("Erreur a la lecture des valeurs du profil : ");
00479             SSCRUTE(pflname);
00480             ret = -1;
00481           }
00482           printf("\t");
00483           for (m=0;m<pflsize;m++) printf(" "IFORMAT" ",*(pflval+m));
00484           printf("\n");
00485           free(pflval);
00486 
00487         }
00488 
00489       }
00490     }
00491   } /* fin for sur les mailles*/
00492 
00493   return ret;
00494 }
00495 

Généré le Thu Oct 8 14:26:17 2015 pour MED fichier par  doxygen 1.6.1