00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
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
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
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
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
00127 for (i =0;i<ncha;i++) {
00128 lret = 0;
00129 printf("\nChamp numero : %d \n",i+1);
00130
00131
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
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("Nom des composantes : |%s|\n",comp);
00153 printf("Unites des composantes : |%s| \n",unit);
00154 printf("Unites des dates : |%s| \n",_dtunit);
00155 printf("Le maillage associé est |%s|\n",_meshname);
00156 printf("Nombre de séquences de calcul |"IFORMAT"|\n",_ncstp);
00157
00158
00159 if ( !_local ) {
00160
00161 if ( (lnsize=MEDlinkInfoByName(fid,_meshname) ) < 0 ) {
00162 MESSAGE("Erreur a la lecture de la taille du lien : ");
00163 SSCRUTE(_meshname);
00164 ret = -1;
00165 } else {
00166
00167 lien = malloc((lnsize+1)*sizeof(char));
00168 EXIT_IF(lien == NULL,NULL,NULL);
00169
00170 if ( MEDlinkRd(fid,_meshname, lien) < 0 ) {
00171 MESSAGE("Erreur a la lecture du lien : ");
00172 SSCRUTE(_meshname);SSCRUTE(lien);
00173 ret = -1;
00174 } else {
00175 printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",_meshname,lien);
00176 }
00177 free(lien);
00178 }
00179 }
00180
00181 free(comp);
00182 free(unit);
00183
00184
00185 lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NODE, USER_INTERLACE,_ncstp );
00186
00187 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_CELL, USER_INTERLACE,_ncstp );
00188 else { MESSAGE("Erreur a la lecture des champs aux noeuds "); ret = -1; continue;}
00189
00190 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_DESCENDING_FACE,USER_INTERLACE,_ncstp);
00191 else { MESSAGE("Erreur a la lecture des champs aux mailles "); ret = -1; continue;}
00192
00193 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_DESCENDING_EDGE,USER_INTERLACE,_ncstp);
00194 else {MESSAGE("Erreur a la lecture des champs aux faces "); ret = -1; continue;}
00195
00196 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NODE_ELEMENT,USER_INTERLACE,_ncstp);
00197 else {MESSAGE("Erreur a la lecture des champs aux aretes"); ret = -1; continue;}
00198
00199
00200
00201
00202
00203 if (lret != 0) {MESSAGE("Erreur a la lecture des champs aux noeuds des mailles "); ret = -1;};
00204 }
00205
00206
00207
00208 npro = MEDnProfile(fid);
00209
00210 printf("\nNombre de profils stockes : "IFORMAT"\n\n",npro);
00211 for (i=1 ; i <= npro ; i++ ) {
00212 if ( MEDprofileInfo(fid, i, pflname, &nval) < 0) {
00213 MESSAGE("Erreur a la demande d'information sur le profil n° : "); ISCRUTE_int(i);
00214 ret = -1;continue;
00215 }
00216 printf("\t- Profil n°%i de nom |%s| et de taille "IFORMAT"\n",i,pflname,nval);
00217 pflval = (med_int*) malloc(sizeof(med_int)*nval);
00218 if ( MEDprofileRd(fid, pflname, pflval) < 0) {
00219 MESSAGE("Erreur a la lecture des valeurs du profil : ");
00220 SSCRUTE(pflname);
00221 ret = -1;
00222 } else {
00223 printf("\t");
00224 for (j=0;j<nval;j++) printf(" "IFORMAT" ",*(pflval+j));
00225 printf("\n\n");
00226 }
00227 free(pflval);
00228 }
00229
00230
00231 nln = MEDnLink(fid);
00232
00233 printf("\nNombre de liens stockes : "IFORMAT"\n\n",nln);
00234 for (i=1 ; i <= nln ; i++ ) {
00235 if ( MEDlinkInfo(fid, i, nomlien, &nval) < 0) {
00236 MESSAGE("Erreur a la demande d'information sur le lien n° : "); ISCRUTE_int(i);
00237 ret = -1;continue;
00238 }
00239 printf("\t- Lien n°%i de nom |%s| et de taille "IFORMAT"\n",i,nomlien,nval);
00240
00241 lien = (char * ) malloc((nval+1)*sizeof(char));
00242 EXIT_IF(lien == NULL,NULL,NULL);
00243
00244 if ( MEDlinkRd(fid, nomlien, lien ) < 0 ) {
00245 MESSAGE("Erreur a la lecture du lien : ");
00246 SSCRUTE(nomlien);SSCRUTE(lien);
00247 ret = -1;
00248 } else {
00249 lien[nval] = '\0';
00250 printf("\t\t|%s|\n\n",lien);
00251 }
00252 free(lien);
00253 }
00254
00255
00256 nloc = MEDnLocalization(fid);
00257
00258 printf("\nNombre de localisations stockees : "IFORMAT"\n\n",nloc);
00259 for (i=1 ; i <= nloc ; i++ ) {
00260 if ( MEDlocalizationInfo(fid, i, locname, &type_geo, &locsdim,&ngauss,
00261 geointerpname, ipointstructmeshname,&nsectionmeshcell,
00262 §iongeotype) < 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 "IFORMAT" 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
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 i,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,_nprofile;
00307 med_int meshnumdt=0,meshnumit=0;
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 meshname [MED_NAME_SIZE+1]="";
00313 char * lien = NULL;
00314 char dt_unit [MED_SNAME_SIZE+1]="unknown";
00315 med_bool localmesh;
00316 med_int nmesh=0;
00317
00318 med_geometry_type * type_geo;
00319
00320 const char * const * AFF;
00321 const char * const * AFF_ENT=MED_GET_ENTITY_TYPENAME+1;
00322 switch (entite) {
00323 case MED_NODE :
00324 type_geo = MED_GET_NODE_GEOMETRY_TYPE;
00325 nb_geo = MED_N_NODE_FIXED_GEO;
00326 AFF = MED_GET_NODE_GEOMETRY_TYPENAME;
00327 break;
00328 case MED_CELL :
00329 case MED_NODE_ELEMENT :
00330 type_geo = MED_GET_CELL_GEOMETRY_TYPE;
00331 nb_geo = MED_N_CELL_FIXED_GEO;
00332 AFF = MED_GET_CELL_GEOMETRY_TYPENAME;
00333 break;
00334 case MED_DESCENDING_FACE :
00335 type_geo = MED_GET_FACE_GEOMETRY_TYPE;
00336 nb_geo = MED_N_FACE_FIXED_GEO;
00337 AFF = MED_GET_FACE_GEOMETRY_TYPENAME;
00338 break;
00339 case MED_DESCENDING_EDGE :
00340 type_geo = MED_GET_EDGE_GEOMETRY_TYPE;
00341 nb_geo = MED_N_EDGE_FIXED_GEO;
00342 AFF = MED_GET_EDGE_GEOMETRY_TYPENAME;
00343 break;
00344 }
00345
00346 for (k=1;k<=nb_geo;k++) {
00347
00348
00349 nbpdtnor = ncstp;
00350 if (nbpdtnor < 1 ) continue;
00351
00352 for (j=0;j<nbpdtnor;j++) {
00353
00354 if ( MEDfield23ComputingStepMeshInfo(fid,nomcha,j+1, &numdt, &numo, &dt,
00355 &nmesh, meshname,&localmesh, &meshnumdt, &meshnumit ) <0) {
00356 MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
00357 ISCRUTE(numdt); ISCRUTE(numo);ISCRUTE(nmesh);SSCRUTE(meshname);ISCRUTE_int(localmesh);
00358 ISCRUTE(meshnumdt);ISCRUTE(meshnumit);
00359 ret = -1; continue;
00360 }
00361
00362 for (i=0;i< nmesh;++i) {
00363
00364 if ( (_nprofile = MEDfield23nProfile(fid,nomcha,numdt,numo,entite,type_geo[k],i+1,meshname,
00365 pflname,locname ) ) < 0 ) {
00366 MESSAGE("Erreur a la demande du nombre de profils referencés par le champ : ");
00367 SSCRUTE(nomcha); ISCRUTE(numdt); ISCRUTE(numo);SSCRUTE(meshname);
00368 ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);SSCRUTE(pflname);SSCRUTE(locname);
00369 SSCRUTE(AFF_ENT[(int)entite]);SSCRUTE(AFF[k]);
00370 ret = -1; continue;
00371 };
00372
00373 for (l=0;l<_nprofile;l++) {
00374
00375
00376 if ( (nval = MEDfield23nValueWithProfile(fid, nomcha, numdt, numo, entite, type_geo[k],meshname,
00377 l+1, USER_MODE, pflname,&pflsize,
00378 locname, &ngauss) ) < 0 ) {
00379 MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
00380 SSCRUTE(nomcha);ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(meshname);
00381 ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
00382 ISCRUTE_int(USER_MODE);
00383 ret = -1; continue;
00384 };
00385
00386 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);
00387 printf("\t- Il y a "IFORMAT" entités qui portent des valeurs en mode %i. Chaque entite %s\
00388 de type geometrique %s associes au profile |%s| a "IFORMAT" valeurs associées \n",
00389 nval,USER_MODE,AFF_ENT[(int)entite],AFF[k],pflname,ngauss);
00390 printf("\t- Le maillage associé est |%s|\n",meshname);
00391
00392
00393 if (typcha == MED_FLOAT64) {
00394
00395 valr = (med_float*) calloc(ncomp*nval*ngauss,sizeof(med_float));
00396 EXIT_IF(valr == NULL,NULL,NULL);
00397
00398 if (MEDfield23ValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],meshname,
00399 USER_MODE, pflname, stockage,MED_ALL_CONSTITUENT,
00400 (unsigned char*) valr) < 0 ) {
00401 MESSAGE("Erreur a la lecture des valeurs du champ : ");
00402 SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
00403 ISCRUTE(numdt);ISCRUTE(numo);
00404 ret = -1;
00405 }
00406 } else {
00407
00408 vale = (med_int*) calloc(ncomp*nval*ngauss,sizeof(med_int));
00409 EXIT_IF(vale == NULL,NULL,NULL);
00410
00411 if (MEDfield23ValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],meshname,
00412 USER_MODE, pflname, stockage,MED_ALL_CONSTITUENT,
00413 (unsigned char*) vale) < 0 ) {
00414 MESSAGE("Erreur a la lecture des valeurs du champ : ");
00415 SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
00416 ISCRUTE(numdt);ISCRUTE(numo);
00417 ret = -1;
00418 };
00419 }
00420
00421 if ( strlen(locname) )
00422 printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
00423
00424 if (entite == MED_NODE_ELEMENT)
00425 ngroup = (type_geo[k] % 100);
00426 else
00427 ngroup = ngauss;
00428
00429 switch (stockage) {
00430
00431 case MED_FULL_INTERLACE :
00432 printf("\t- Valeurs :\n\t");
00433 for (m=0;m<(nval*ngauss)/ngroup;m++) {
00434 printf("|");
00435 for (n=0;n<ngroup*ncomp;n++)
00436 if (typcha == MED_FLOAT64)
00437 printf(" %f ",*(valr+(m*ngroup*ncomp)+n));
00438 else
00439 printf(" "IFORMAT" ",*(vale+(m*ngroup*ncomp)+n));
00440
00441 }
00442 break;
00443
00444
00445 case MED_NO_INTERLACE :
00446 printf("\t- Valeurs :\n\t");
00447 for (m=0;m<ncomp;m++) {
00448 printf("|");
00449 for (n=0;n<(nval*ngauss);n++)
00450 if (typcha == MED_FLOAT64)
00451 printf(" %f ",*(valr+(m*nval)+n));
00452 else
00453 printf(" "IFORMAT" ",*(vale+(m*nval)+n));
00454 }
00455 break;
00456 }
00457
00458 printf("|\n");
00459 if (typcha == MED_FLOAT64) {
00460 if ( valr ) {free(valr);valr = NULL;}}
00461 else
00462 if (vale) { free(vale);vale = NULL; }
00463
00464
00465 if (strcmp(pflname,MED_NO_PROFILE) == 0 )
00466 printf("\t- Profil : MED_NO_PROFILE\n");
00467 else {
00468 if ( (pflsize = MEDprofileSizeByName(fid,pflname)) <0 ) {
00469 MESSAGE("Erreur a la lecture du nombre de valeurs du profil : ");
00470 SSCRUTE(pflname);
00471 ret = -1; continue;
00472 }
00473
00474 printf("\t- Profil : |%s| de taille "IFORMAT"\n",pflname,pflsize);
00475
00476 pflval = (med_int*) malloc(sizeof(med_int)*pflsize);
00477 EXIT_IF(pflval == NULL,NULL,NULL);
00478 if ( MEDprofileRd(fid,pflname,pflval) <0) {
00479 MESSAGE("Erreur a la lecture des valeurs du profil : ");
00480 SSCRUTE(pflname);
00481 ret = -1;
00482 }
00483 printf("\t");
00484 for (m=0;m<pflsize;m++) printf(" "IFORMAT" ",*(pflval+m));
00485 printf("\n");
00486 free(pflval);
00487 }
00488 }
00489 }
00490 }
00491 }
00492
00493 return ret;
00494 }
00495