1 /************************************************************************* 2 * COPYRIGHT (C) 1999 - 2007 EDF R&D, CEA/DEN 3 * THIS LIBRARY IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY 4 * IT UNDER THE TERMS OF THE GNU LESSER GENERAL PUBLIC LICENSE 5 * AS PUBLISHED BY THE FREE SOFTWARE FOUNDATION; 6 * EITHER VERSION 2.1 OF THE LICENSE, OR (AT YOUR OPTION) ANY LATER VERSION. 7 * 8 * THIS LIBRARY IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT 9 * WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF 10 * MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU 11 * LESSER GENERAL PUBLIC LICENSE FOR MORE DETAILS. 12 * 13 * YOU SHOULD HAVE RECEIVED A COPY OF THE GNU LESSER GENERAL PUBLIC LICENSE 14 * ALONG WITH THIS LIBRARY; IF NOT, WRITE TO THE FREE SOFTWARE FOUNDATION, 15 * INC., 59 TEMPLE PLACE, SUITE 330, BOSTON, MA 02111-1307 USA 16 * 17 *************************************************************************/ 18 19 /****************************************************************************** 20 * - Nom du fichier : test11.c 21 * 22 * - Description : lecture de champs de resultats MED 23 * 24 *****************************************************************************/ 25 26 #include <med.h> 27 #define MESGERR 28 #include <med_utils.h> 29 30 #ifdef DEF_LECT_ECR 31 #define MODE_ACCES MED_LECTURE_ECRITURE 32 #elif DEF_LECT_AJOUT 33 #define MODE_ACCES MED_LECTURE_AJOUT 34 #else 35 #define MODE_ACCES MED_CREATION 36 #endif 37 38 #ifndef USER_INTERLACE 39 #define USER_INTERLACE MED_FULL_INTERLACE 40 #endif 41 42 #define USER_MODE MED_COMPACT 43 44 med_err getFieldsOn(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp, 45 med_entite_maillage entite, med_mode_switch stockage); 46 47 int main (int argc, char **argv) 48 49 50 { 51 med_err ret,lret; 52 med_idt fid; 53 char * fichier = NULL; 54 char maa[MED_TAILLE_NOM+1]=""; 55 char desc[MED_TAILLE_DESC+1]=""; 56 char pflname[MED_TAILLE_NOM+1]="",nomlien[MED_TAILLE_NOM+1]=""; 57 char locname[MED_TAILLE_NOM+1]=""; 58 char * lien = NULL; 59 char *comp= NULL, *unit= NULL; 60 char nomcha [MED_TAILLE_NOM+1]=""; 61 med_int mdim,ncomp,ncha,npro,nln,pflsize,*pflval,nval; 62 med_int ngauss,nloc; 63 int t1,t2,t3; 64 med_type_champ typcha; 65 med_maillage type; 66 med_geometrie_element type_geo; 67 med_float *refcoo, *gscoo, *wg; 68 int i,j; 69 70 if (argc != 2) { 71 MESSAGE("Aucun nom de fichier precise, fichier test10.med utilise "); 72 fichier = "test10.med"; 73 } else { 74 fichier = argv[1]; 75 }; 76 77 78 /* Ouverture du fichier med */ 79 if ((fid = MEDouvrir(fichier,MED_LECTURE)) < 0){ 80 MESSAGE("Erreur a l'ouverture du fichier : ");SSCRUTE(fichier); 81 return -1; 82 } 83 84 ret = 0; 85 86 /* infos sur le premier maillage */ 87 if ( MEDmaaInfo(fid,1,maa,&mdim,&type,desc) < 0 ) { 88 MESSAGE("Erreur a la lecture des informations sur le maillage : "); 89 SSCRUTE(maa);ISCRUTE(mdim);ISCRUTE(type);SSCRUTE(desc); 90 return -1; 91 } 92 93 printf("Maillage de nom |%s| et de dimension %d \n",maa,mdim); 94 95 /* combien de champs dans le fichier */ 96 if ((ncha = MEDnChamp(fid,0)) < 0) { 97 MESSAGE("Impossible de lire le nombre de champs : ");ISCRUTE(ncha); 98 return ncha; 99 } 100 101 printf("Nombre de champs : %d \n",ncha); 102 103 /* lecture de tous les champs */ 104 for (i =0;i<ncha;i++) { 105 lret = 0; 106 printf("\nChamp numero : %d \n",i+1); 107 108 /* Lecture du nombre de composantes */ 109 if ((ncomp = MEDnChamp(fid,i+1)) < 0) { 110 MESSAGE("Erreur a la lecture du nombre de composantes : "); ISCRUTE(ncomp); 111 ret = -1; continue; 112 } 113 114 /* Lecture du type du champ, des noms des composantes et du nom de l'unite*/ 115 comp = (char*) malloc(ncomp*MED_TAILLE_PNOM+1); 116 EXIT_IF(comp == NULL,NULL,NULL); 117 unit = (char*) malloc(ncomp*MED_TAILLE_PNOM+1); 118 EXIT_IF(unit == NULL,NULL,NULL); 119 120 if ( MEDchampInfo(fid,i+1,nomcha,&typcha,comp,unit,ncomp) < 0 ) { 121 MESSAGE("Erreur a la demande d'information sur les champs : "); 122 ISCRUTE(i+1);SSCRUTE(nomcha);ISCRUTE(typcha);SSCRUTE(comp);SSCRUTE(unit); 123 ISCRUTE(ncomp); 124 ret = -1; continue; 125 } 126 127 printf("Nom du champ : |%s| de type %d\n",nomcha,typcha); 128 printf("Nom des composantes : |%s|\n",comp); 129 printf("Unites des composantes : |%s| \n",unit); 130 131 free(comp); 132 free(unit); 133 134 135 lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD, USER_INTERLACE ); 136 137 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_MAILLE, USER_INTERLACE ); 138 else { MESSAGE("Erreur a la lecture des champs aux noeuds "); ret = -1; continue;} 139 140 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_FACE,USER_INTERLACE); 141 else { MESSAGE("Erreur a la lecture des champs aux mailles "); ret = -1; continue;} 142 143 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_ARETE,USER_INTERLACE); 144 else {MESSAGE("Erreur a la lecture des champs aux faces "); ret = -1; continue;} 145 146 if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD_MAILLE,USER_INTERLACE); 147 else {MESSAGE("Erreur a la lecture des champs aux aretes"); ret = -1; continue;} 148 149 if (lret != 0) {MESSAGE("Erreur a la lecture des champs aux noeuds des mailles "); ret = -1;}; 150 } 151 152 153 /* Interrogation des profils */ 154 npro = MEDnProfil(fid); 155 156 printf("\nNombre de profils stockes : %i\n\n",npro); 157 for (i=1 ; i <= npro ; i++ ) { 158 if ( MEDprofilInfo(fid, i, pflname, &nval) < 0) { 159 MESSAGE("Erreur a la demande d'information sur le profil n° : "); ISCRUTE(i); 160 ret = -1;continue; 161 } 162 printf("\t- Profil n°%i de nom |%s| et de taille %i\n",i,pflname,nval); 163 pflval = (med_int*) malloc(sizeof(med_int)*nval); 164 if ( MEDprofilLire(fid, pflval, pflname) < 0) { 165 MESSAGE("Erreur a la lecture des valeurs du profil : "); 166 SSCRUTE(pflname); 167 ret = -1; 168 } else { 169 printf("\t"); 170 for (j=0;j<nval;j++) printf(" %i ",*(pflval+j)); 171 printf("\n\n"); 172 } 173 free(pflval); 174 } 175 176 /* Interrogation des liens */ 177 nln = MEDnLien(fid); 178 179 printf("\nNombre de liens stockes : %i\n\n",nln); 180 for (i=1 ; i <= nln ; i++ ) { 181 if ( MEDlienInfo(fid, i, nomlien, &nval) < 0) { 182 MESSAGE("Erreur a la demande d'information sur le lien n° : "); ISCRUTE(i); 183 ret = -1;continue; 184 } 185 printf("\t- Lien n°%i de nom |%s| et de taille %i\n",i,nomlien,nval); 186 187 lien = malloc((nval+1)*sizeof(char)); 188 EXIT_IF(lien == NULL,NULL,NULL); 189 190 if ( MEDlienLire(fid, lien, nomlien) < 0 ) { 191 MESSAGE("Erreur a la lecture du lien : "); 192 SSCRUTE(nomlien);SSCRUTE(lien); 193 ret = -1; 194 } else { 195 lien[nval] = '\0'; 196 printf("\t\t|%s|\n\n",lien); 197 } 198 free(lien); 199 } 200 201 /* Interrogation des localisations des points de GAUSS */ 202 nloc = MEDnGauss(fid); 203 204 printf("\nNombre de localisations stockees : %i\n\n",nloc); 205 for (i=1 ; i <= nloc ; i++ ) { 206 if ( MEDgaussInfo(fid, i, locname, &type_geo, &ngauss) < 0) { 207 MESSAGE("Erreur a la demande d'information sur la localisation n° : "); ISCRUTE(i); 208 ret = -1;continue; 209 } 210 printf("\t- Loc. n°%i de nom |%s| et nbr. de pts de GAUSS %i\n",i,locname,ngauss); 211 t1 = (type_geo%100)*(type_geo/100); 212 t2 = ngauss*(type_geo/100); 213 t3 = ngauss; 214 refcoo = (med_float *) malloc(sizeof(med_float)*t1 ); 215 gscoo = (med_float *) malloc(sizeof(med_float)*t2 ); 216 wg = (med_float *) malloc(sizeof(med_float)*t3 ); 217 218 if ( MEDgaussLire(fid, refcoo, gscoo, wg, USER_INTERLACE, locname ) < 0) { 219 MESSAGE("Erreur a la lecture des valeurs de la localisation : "); 220 SSCRUTE(locname); 221 ret = -1; 222 } else { 223 printf("\t Coordonnees de l'element de reference de type %i :\n\t\t",type_geo); 224 for (j=0;j<t1;j++) printf(" %f ",*(refcoo+j)); 225 printf("\n"); 226 printf("\t Localisation des points de GAUSS : \n\t\t"); 227 for (j=0;j<t2;j++) printf(" %f ",*(gscoo+j)); 228 printf("\n"); 229 printf("\t Poids associes aux points de GAUSS :\n\t\t"); 230 for (j=0;j<t3;j++) printf(" %f ",*(wg+j)); 231 printf("\n\n"); 232 } 233 free(refcoo); 234 free(gscoo); 235 free(wg); 236 } 237 238 239 /* fermeture du fichier */ 240 if ( MEDfermer(fid) < 0) return -1; 241 242 return ret; 243 } 244 245 med_err getFieldsOn(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp, 246 med_entite_maillage entite, med_mode_switch stockage) { 247 248 int j,k,l,m,n,nb_geo; 249 med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,ngroup,*vale=NULL,nval; 250 med_int numdt=0,numo=0,lnsize,nbrefmaa; 251 med_float *valr=NULL,dt=0.0; 252 med_err ret=0; 253 med_booleen local; 254 char pflname [MED_TAILLE_NOM+1]=""; 255 char locname [MED_TAILLE_NOM+1]=""; 256 char * lien = NULL; 257 char maa_ass [MED_TAILLE_NOM+1]=""; 258 char dt_unit [MED_TAILLE_PNOM+1]=""; 259 260 261 med_geometrie_element * type_geo; 262 med_geometrie_element typ_noeud[1] = { MED_NONE }; 263 med_geometrie_element typmai[MED_NBR_GEOMETRIE_MAILLE+2] = {MED_POINT1, MED_SEG2, MED_SEG3, MED_TRIA3, 264 MED_QUAD4, MED_TRIA6,MED_QUAD8, MED_TETRA4, 265 MED_PYRA5, MED_PENTA6, MED_HEXA8, MED_TETRA10, 266 MED_PYRA13, MED_PENTA15, MED_HEXA20, 267 MED_POLYGONE, MED_POLYEDRE}; 268 med_geometrie_element typfac[MED_NBR_GEOMETRIE_FACE+1] = {MED_TRIA3,MED_TRIA6, 269 MED_QUAD4,MED_QUAD8, 270 MED_POLYGONE}; 271 med_geometrie_element typare[MED_NBR_GEOMETRIE_ARETE] = {MED_SEG2,MED_SEG3}; 272 273 char ** AFF; 274 275 switch (entite) { 276 case MED_NOEUD : 277 type_geo = typ_noeud; 278 nb_geo = 1; 279 AFF = MED_GEOMETRIE_NOEUD_AFF; 280 break; 281 case MED_MAILLE : 282 case MED_NOEUD_MAILLE : 283 type_geo = typmai; 284 nb_geo = MED_NBR_GEOMETRIE_MAILLE+2; 285 AFF = MED_GEOMETRIE_MAILLE_AFF; 286 break; 287 case MED_FACE : 288 type_geo = typfac; 289 nb_geo = MED_NBR_GEOMETRIE_FACE+1; 290 AFF = MED_GEOMETRIE_FACE_AFF; 291 break; 292 case MED_ARETE : 293 type_geo = typare; 294 nb_geo = MED_NBR_GEOMETRIE_ARETE; 295 AFF = MED_GEOMETRIE_ARETE_AFF; 296 break; 297 } 298 299 300 for (k=0;k<nb_geo;k++) { 301 302 /* Combien de (PDT,NOR) a lire */ 303 nbpdtnor = MEDnPasdetemps(fid,nomcha,entite,type_geo[k]); 304 if (nbpdtnor < 1 ) continue; 305 306 for (j=0;j<nbpdtnor;j++) { 307 308 if ( MEDpasdetempsInfo(fid,nomcha,entite,type_geo[k], 309 j+1, &ngauss, &numdt, &numo, dt_unit, 310 &dt, maa_ass, &local, &nbrefmaa) <0) { 311 MESSAGE("Erreur a la demande d'information sur (pdt,nor) : "); 312 ISCRUTE(numdt); ISCRUTE(numo); 313 ret = -1; continue; 314 }; 315 316 printf("\n +Pas de Temps n.%i (%f) [%s], n. d'ordre %i, avec %i pts de gauss sur le maillage par defaut.\n",numdt,dt,dt_unit,numo,ngauss); 317 318 printf("\tLe maillage par defaut est : |%s|, sur un total de : %i maillages associes\n", 319 maa_ass, nbrefmaa); 320 321 /* Le maillage reference est-il porte par un autre fichier */ 322 if ( !local ) { 323 324 if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 ) { 325 MESSAGE("Erreur a la lecture de la taille du lien : "); 326 SSCRUTE(maa_ass); 327 ret = -1; 328 } else { 329 330 lien = malloc(lnsize*sizeof(char)); 331 EXIT_IF(lien == NULL,NULL,NULL); 332 333 if ( MEDlienLire(fid, lien, maa_ass) < 0 ) { 334 MESSAGE("Erreur a la lecture du lien : "); 335 SSCRUTE(maa_ass);SSCRUTE(lien); 336 ret = -1; 337 } else { 338 printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien); 339 } 340 free(lien); 341 } 342 } 343 344 /* Combien de maillages lies aux (nomcha,ent,geo,numdt,numo) */ 345 /* Notons que cette information est egalement disponible a partir de MEDpasdetempsInfo */ 346 if ( (nbrefmaa = MEDnChampRef(fid,nomcha,entite,type_geo[k],numdt,numo) ) < 0 ) { 347 MESSAGE("Erreur a la demande du nombre de maillages references par le champ : "); 348 SSCRUTE(nomcha); 349 ISCRUTE(entite);ISCRUTE(type_geo[k]); 350 SSCRUTE(MED_ENTITE_MAILLAGE_AFF[(int)entite]);SSCRUTE(AFF[k]); 351 ISCRUTE(numdt); ISCRUTE(numo); 352 ret = -1; continue; 353 }; 354 355 for (l=0;l<nbrefmaa;l++) { 356 357 if ( MEDchampRefInfo(fid,nomcha,entite,type_geo[k], 358 l+1,numdt, numo, maa_ass, &local, &ngauss) <0 ) { 359 MESSAGE("Erreur a la demande d'information sur le maillage utilise par le champ n° : "); 360 SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]); 361 ISCRUTE(l+1);ISCRUTE(numdt); ISCRUTE(numo);SSCRUTE(maa_ass); 362 ret = -1; continue; 363 }; 364 365 /* Prend en compte le nbre de pt de gauss automatiquement */ 366 if ((nval = MEDnVal(fid,nomcha,entite,type_geo[k],numdt,numo,maa_ass,USER_MODE)) <= 0) { 367 MESSAGE("Erreur a la lecture du nombre de valeurs du champ : "); 368 SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]); 369 ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);ISCRUTE(USER_MODE); 370 ret = -1; continue; 371 }; 372 373 374 printf("\t- Il y a %d valeurs en mode %i. Chaque entite %s\ 375 de type geometrique %s associes au maillage |%s| a %i pts de gauss \n", 376 nval,USER_MODE,MED_ENTITE_MAILLAGE_AFF[(int)entite],AFF[k],maa_ass,ngauss); 377 378 /* Le maillage reference est-il porte par un autre fichier */ 379 if ( !local ) { 380 381 if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 ) { 382 MESSAGE("Erreur a la lecture de la taille du lien : "); 383 SSCRUTE(maa_ass); 384 ret = -1; 385 } else { 386 387 lien = malloc(lnsize*sizeof(char) + 1); 388 EXIT_IF(lien == NULL,NULL,NULL); 389 390 if ( MEDlienLire(fid, lien, maa_ass) < 0 ) { 391 MESSAGE("Erreur a la lecture du lien : "); 392 SSCRUTE(maa_ass);SSCRUTE(lien); 393 ret = -1; 394 } else { 395 lien[lnsize] = '\0'; 396 printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien); 397 } 398 free(lien); 399 } 400 } 401 402 /*Lecture des valeurs du champ */ 403 if (typcha == MED_FLOAT64) { 404 405 valr = (med_float*) calloc(ncomp*nval,sizeof(med_float)); 406 EXIT_IF(valr == NULL,NULL,NULL); 407 if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)valr,stockage,MED_ALL,locname, 408 pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) { 409 MESSAGE("Erreur a la lecture du nombre de valeurs du champ : "); 410 SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]); 411 ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass); 412 ret = -1; 413 }; 414 415 } else { 416 417 vale = (med_int*) calloc(ncomp*nval,sizeof(med_int)); 418 EXIT_IF(vale == NULL,NULL,NULL); 419 if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)vale,stockage,MED_ALL,locname, 420 pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) { 421 MESSAGE("Erreur a la lecture des valeurs du champ : "); 422 SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]); 423 ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass); 424 ret = -1; 425 }; 426 427 } 428 429 if ( ngauss > 1 ) 430 printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname); 431 432 if (entite == MED_NOEUD_MAILLE) 433 ngroup = type_geo[k] % 100; 434 else 435 ngroup = ngauss; 436 437 switch (stockage) { 438 439 case MED_FULL_INTERLACE : 440 printf("\t- Valeurs :\n\t"); 441 for (m=0;m<nval/ngroup;m++) { 442 printf("|"); 443 for (n=0;n<ngroup*ncomp;n++) 444 if (typcha == MED_FLOAT64) 445 printf(" %f ",*(valr+(m*ngroup*ncomp)+n)); 446 else 447 printf(" %d ",*(vale+(m*ngroup*ncomp)+n)); 448 449 } 450 break; 451 452 /*Affichage en fonction du profil à traiter*/ 453 case MED_NO_INTERLACE : 454 printf("\t- Valeurs :\n\t"); 455 for (m=0;m<ncomp;m++) { 456 printf("|"); 457 for (n=0;n<nval;n++) 458 if (typcha == MED_FLOAT64) 459 printf(" %f ",*(valr+(m*nval)+n)); 460 else 461 printf(" %d ",*(vale+(m*nval)+n)); 462 } 463 break; 464 } 465 466 printf("|\n"); 467 if (typcha == MED_FLOAT64) { 468 if ( valr ) {free(valr);valr = NULL;}} 469 else 470 if (vale) { free(vale);vale = NULL; } 471 472 /*Lecture du profil associe */ 473 if (strcmp(pflname,MED_NOPFL) == 0 ) 474 printf("\t- Profil : MED_NOPFL\n"); 475 else { 476 477 if ( (pflsize = MEDnValProfil(fid,pflname)) <0 ) { 478 MESSAGE("Erreur a la lecture du nombre de valeurs du profil : "); 479 SSCRUTE(pflname); 480 ret = -1; continue; 481 } 482 483 printf("\t- Profil : |%s| de taille %i\n",pflname,pflsize); 484 485 pflval = (med_int*) malloc(sizeof(med_int)*pflsize); 486 EXIT_IF(pflval == NULL,NULL,NULL); 487 if ( MEDprofilLire(fid,pflval,pflname) <0) { 488 MESSAGE("Erreur a la lecture des valeurs du profil : "); 489 SSCRUTE(pflname); 490 ret = -1; 491 } 492 printf("\t"); 493 for (m=0;m<pflsize;m++) printf(" %i ",*(pflval+m)); 494 printf("\n"); 495 free(pflval); 496 497 } 498 499 } 500 } 501 } /* fin for sur les mailles*/ 502 503 return ret; 504 } 505