33 #define MODE_ACCES MED_ACC_RDWR
35 #define MODE_ACCES MED_ACC_RDEXT
37 #define MODE_ACCES MED_ACC_CREAT
40 #ifndef USER_INTERLACE
41 #define USER_INTERLACE MED_FULL_INTERLACE
44 #define USER_MODE MED_COMPACT_STMODE
49 int main (
int argc,
char **argv)
55 char * fichier = NULL;
63 char *comp= NULL, *unit= NULL;
65 med_int mdim=0,sdim=0,ncomp,ncha,npro,nln,pflsize,*pflval,nval;
66 med_int _ncstp=0,ngauss=0,nloc=0,locsdim=0,lnsize=0;
87 MESSAGE(
"Aucun nom de fichier precise, fichier test10.med utilise ");
88 fichier =
"test10.med";
104 if (
MEDmeshInfo( fid, 1, maa, &sdim, &mdim, &type, desc, dtunit, &sort,
105 &nstep, &rep, nomcoo,unicoo) < 0 ) {
106 MESSAGE(
"Erreur a la lecture des informations sur le maillage : ");
SSCRUTE(maa);
109 printf(
"Maillage de nom : |%s| , de dimension : "IFORMAT" , et de type %d\n",maa,mdim,type);
110 printf(
"\t -Dimension de l'espace : "IFORMAT"\n",sdim);
111 printf(
"\t -Description du maillage : %s\n",desc);
112 printf(
"\t -Noms des axes : |%s|\n",nomcoo);
113 printf(
"\t -Unités des axes : |%s|\n",unicoo);
114 printf(
"\t -Type de repère : %d\n",rep);
115 printf(
"\t -Nombre d'étapes de calcul : "IFORMAT"\n",nstep);
116 printf(
"\t -Unité des dates : |%s|\n",dtunit);
122 MESSAGE(
"Impossible de lire le nombre de champs : ");
ISCRUTE(ncha);
126 printf(
"Nombre de champs : "IFORMAT" \n",ncha);
129 for (i =0;i<ncha;i++) {
131 printf(
"\nChamp numero : %d \n",i+1);
135 MESSAGE(
"Erreur a la lecture du nombre de composantes : ");
ISCRUTE(ncomp);
141 EXIT_IF(comp == NULL,NULL,NULL);
143 EXIT_IF(unit == NULL,NULL,NULL);
145 if (
MEDfieldInfo(fid,i+1,nomcha,_meshname,&_local,&typcha,comp,unit,_dtunit,&_ncstp) < 0 ) {
146 MESSAGE(
"Erreur a la demande d'information sur les champs : ");
153 printf(
"Nom du champ : |%s| de type %d\n",nomcha,typcha);
154 printf(
"Nombre de composantes : |"IFORMAT"|\n",ncomp);
155 printf(
"Nom des composantes : |%s|\n",comp);
156 printf(
"Unites des composantes : |%s| \n",unit);
157 printf(
"Unites des dates : |%s| \n",_dtunit);
158 printf(
"Le maillage associé est |%s|\n",_meshname);
159 printf(
"Nombre d'étapes de calcul |"IFORMAT"|\n",_ncstp);
165 MESSAGE(
"Erreur a la lecture de la taille du lien : ");
170 lien = malloc((lnsize+1)*
sizeof(
char));
171 EXIT_IF(lien == NULL,NULL,NULL);
173 if (
MEDlinkRd(fid,_meshname, lien) < 0 ) {
174 MESSAGE(
"Erreur a la lecture du lien : ");
178 printf(
"\tLe maillage |%s| est porte par un fichier distant |%s|\n",_meshname,lien);
191 else {
MESSAGE(
"Erreur a la lecture des champs aux noeuds "); ret = -1;
continue;}
194 else {
MESSAGE(
"Erreur a la lecture des champs aux mailles "); ret = -1;
continue;}
197 else {
MESSAGE(
"Erreur a la lecture des champs aux faces "); ret = -1;
continue;}
200 else {
MESSAGE(
"Erreur a la lecture des champs aux aretes"); ret = -1;
continue;}
206 if (lret != 0) {
MESSAGE(
"Erreur a la lecture des champs aux noeuds des mailles "); ret = -1;};
213 printf(
"\nNombre de profils stockes : "IFORMAT"\n\n",npro);
214 for (i=1 ; i <= npro ; i++ ) {
219 printf(
"\t- Profil n°%i de nom |%s| et de taille "IFORMAT"\n",i,pflname,nval);
222 MESSAGE(
"Erreur a la lecture des valeurs du profil : ");
227 for (j=0;j<nval;j++) printf(
" "IFORMAT" ",*(pflval+j));
236 printf(
"\nNombre de liens stockes : "IFORMAT"\n\n",nln);
237 for (i=1 ; i <= nln ; i++ ) {
242 printf(
"\t- Lien n°%i de nom |%s| et de taille "IFORMAT"\n",i,nomlien,nval);
244 lien = malloc((nval+1)*
sizeof(
char));
245 EXIT_IF(lien == NULL,NULL,NULL);
247 if (
MEDlinkRd(fid, nomlien, lien ) < 0 ) {
248 MESSAGE(
"Erreur a la lecture du lien : ");
253 printf(
"\t\t|%s|\n\n",lien);
261 printf(
"\nNombre de localisations stockees : "IFORMAT"\n\n",nloc);
262 for (i=1 ; i <= nloc ; i++ ) {
264 geointerpname, ipointstructmeshname,&nsectionmeshcell,
265 §iongeotype) < 0) {
266 MESSAGE(
"Erreur a la demande d'information sur la localisation n° : ");
ISCRUTE_int(i);
269 printf(
"\t- Loc. n°%i de nom |%s| de dimension |"IFORMAT"| avec "IFORMAT" pts de GAUSS \n",i,locname,locsdim,ngauss);
270 t1 = (type_geo%100)*(type_geo/100);
271 t2 = ngauss*(type_geo/100);
278 MESSAGE(
"Erreur a la lecture des valeurs de la localisation : ");
282 printf(
"\t Coordonnees de l'element de reference de type %i :\n\t\t",type_geo);
283 for (j=0;j<t1;j++) printf(
" %f ",*(refcoo+j));
285 printf(
"\t Localisation des points de GAUSS : \n\t\t");
286 for (j=0;j<t2;j++) printf(
" %f ",*(gscoo+j));
288 printf(
"\t Poids associes aux points de GAUSS :\n\t\t");
289 for (j=0;j<t3;j++) printf(
" %f ",*(wg+j));
307 int j,k,l,m,n,nb_geo=0;
308 med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,ngroup,nval;
309 med_int numdt=0,numo=0,_nprofile;
312 unsigned char *val = NULL;
321 const char *
const * AFF;
347 for (k=1;k<=nb_geo;k++) {
351 if (nbpdtnor < 1 )
continue;
353 for (j=0;j<nbpdtnor;j++) {
356 MESSAGE(
"Erreur a la demande d'information sur (pdt,nor) : ");
361 if ( (_nprofile =
MEDfieldnProfile(fid,nomcha,numdt,numo,entite,type_geo[k],
362 pflname,locname ) ) < 0 ) {
363 MESSAGE(
"Erreur a la demande du nombre de profils referencés par le champ : ");
371 for (l=0;l<_nprofile;l++) {
376 locname, &ngauss) ) < 0 ) {
377 MESSAGE(
"Erreur a la lecture du nombre de valeurs du champ : ");
384 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);
385 printf(
"\t- Il y a "IFORMAT" entités qui portent des valeurs en mode %i. Chaque entite %s\
386 de type geometrique %s associes au profile |%s| a "IFORMAT" valeurs associées \n",
387 nval,
USER_MODE,AFF_ENT[(
int)entite],AFF[k],pflname,ngauss);
398 MESSAGE(
"Erreur a la lecture du type de champ : ");
400 EXIT_IF(NULL == NULL,NULL,NULL);
403 val = (
unsigned char*) calloc(ncomp*nval*ngauss,medtype_size);
404 EXIT_IF(val == NULL,NULL,NULL);
408 MESSAGE(
"Erreur a la lecture des valeurs du champ : ");
414 if ( strlen(locname) )
415 printf(
"\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
418 ngroup = (type_geo[k] % 100);
425 printf(
"\t- Valeurs :\n\t");
426 for (m=0;m<(nval*ngauss)/ngroup;m++) {
428 for (n=0;n<ngroup*ncomp;n++)
431 printf(
" %f ",*(((
med_double*)val)+(m*ngroup*ncomp)+n ) );
436 printf(
" %f ",*(((
med_float32*)val)+((m*ngroup*ncomp)+n)));
439 printf(
" %d ",*(((
med_int32*)val)+(m*ngroup*ncomp)+n));
442 printf(
" %lld ",*(((
med_int64*)val)+(m*ngroup*ncomp)+n));
454 printf(
"\t- Valeurs :\n\t");
455 for (m=0;m<ncomp;m++) {
457 for (n=0;n<(nval*ngauss);n++)
460 printf(
" %f ",*(((
med_double*)val)+(m*nval*ngauss)+n ) );
465 printf(
" %f ",*(((
med_float32*)val)+((m*nval*ngauss)+n)));
468 printf(
" %d ",*(((
med_int32*)val)+(m*nval*ngauss)+n));
471 printf(
" %lld ",*(((
med_int64*)val)+(m*nval*ngauss)+n));
484 if ( val ) {free(val);val = NULL;}
488 printf(
"\t- Profil : MED_NO_PROFILE\n");
492 MESSAGE(
"Erreur a la lecture du nombre de valeurs du profil : ");
497 printf(
"\t- Profil : |%s| de taille "IFORMAT"\n",pflname,pflsize);
500 EXIT_IF(pflval == NULL,NULL,NULL);
502 MESSAGE(
"Erreur a la lecture des valeurs du profil : ");
507 for (m=0;m<pflsize;m++) printf(
" "IFORMAT" ",*(pflval+m));