int main (int argc, char **argv)
{
med_err ret,lret;
med_idt fid;
char * fichier = NULL;
char maa[MED_TAILLE_NOM+1]="";
char desc[MED_TAILLE_DESC+1]="";
char pflname[MED_TAILLE_NOM+1]="",nomlien[MED_TAILLE_NOM+1]="";
char locname[MED_TAILLE_NOM+1]="";
char * lien = NULL;
char *comp= NULL, *unit= NULL;
char nomcha [MED_TAILLE_NOM+1]="";
med_int mdim,ncomp,ncha,npro,nln,pflsize,*pflval,nval;
med_int ngauss,nloc;
int t1,t2,t3;
med_type_champ typcha;
med_maillage type;
med_geometrie_element type_geo;
med_float *refcoo, *gscoo, *wg;
int i,j;
if (argc != 2) {
MESSAGE("Aucun nom de fichier precise, fichier test10.med utilise ");
fichier = "test10.med";
} else {
fichier = argv[1];
};
/* Ouverture du fichier med */
if ((fid = MEDouvrir(fichier,MED_LECTURE)) < 0){
MESSAGE("Erreur a l'ouverture du fichier : ");SSCRUTE(fichier);
return -1;
}
ret = 0;
/* infos sur le premier maillage */
if ( MEDmaaInfo(fid,1,maa,&mdim,&type,desc) < 0 ) {
MESSAGE("Erreur a la lecture des informations sur le maillage : ");
SSCRUTE(maa);ISCRUTE(mdim);ISCRUTE(type);SSCRUTE(desc);
return -1;
}
printf("Maillage de nom |%s| et de dimension %d \n",maa,mdim);
/* combien de champs dans le fichier */
if ((ncha = MEDnChamp(fid,0)) < 0) {
MESSAGE("Impossible de lire le nombre de champs : ");ISCRUTE(ncha);
return ncha;
}
printf("Nombre de champs : %d \n",ncha);
/* lecture de tous les champs */
for (i =0;i<ncha;i++) {
lret = 0;
printf("\nChamp numero : %d \n",i+1);
/* Lecture du nombre de composantes */
if ((ncomp = MEDnChamp(fid,i+1)) < 0) {
MESSAGE("Erreur a la lecture du nombre de composantes : "); ISCRUTE(ncomp);
ret = -1; continue;
}
/* Lecture du type du champ, des noms des composantes et du nom de l'unite*/
comp = (char*) malloc(ncomp*MED_TAILLE_PNOM+1);
EXIT_IF(comp == NULL,NULL,NULL);
unit = (char*) malloc(ncomp*MED_TAILLE_PNOM+1);
EXIT_IF(unit == NULL,NULL,NULL);
if ( MEDchampInfo(fid,i+1,nomcha,&typcha,comp,unit,ncomp) < 0 ) {
MESSAGE("Erreur a la demande d'information sur les champs : ");
ISCRUTE(i+1);SSCRUTE(nomcha);ISCRUTE(typcha);SSCRUTE(comp);SSCRUTE(unit);
ISCRUTE(ncomp);
ret = -1; continue;
}
printf("Nom du champ : |%s| de type %d\n",nomcha,typcha);
printf("Nom des composantes : |%s|\n",comp);
printf("Unites des composantes : |%s| \n",unit);
free(comp);
free(unit);
lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NOEUD, USER_INTERLACE );
if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_MAILLE, USER_INTERLACE );
else { MESSAGE("Erreur a la lecture des champs aux noeuds "); ret = -1; continue;}
if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_FACE,USER_INTERLACE);
else { MESSAGE("Erreur a la lecture des champs aux mailles "); ret = -1; continue;}
if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_ARETE,USER_INTERLACE);
else {MESSAGE("Erreur a la lecture des champs aux faces "); ret = -1; continue;}
if (lret != 0) {MESSAGE("Erreur a la lecture des champs aux aretes "); ret = -1;};
}
/* Interrogation des profils */
npro = MEDnProfil(fid);
printf("\nNombre de profils stockes : %i\n\n",npro);
for (i=1 ; i <= npro ; i++ ) {
if ( MEDprofilInfo(fid, i, pflname, &nval) < 0) {
MESSAGE("Erreur a la demande d'information sur le profil n° : "); ISCRUTE(i);
ret = -1;continue;
}
printf("\t- Profil n°%i de nom |%s| et de taille %i\n",i,pflname,nval);
pflval = (med_int*) malloc(sizeof(med_int)*nval);
if ( MEDprofilLire(fid, pflval, pflname) < 0) {
MESSAGE("Erreur a la lecture des valeurs du profil : ");
SSCRUTE(pflname);
ret = -1;
} else {
printf("\t");
for (j=0;j<nval;j++) printf(" %i ",*(pflval+j));
printf("\n\n");
}
free(pflval);
}
/* Interrogation des liens */
nln = MEDnLien(fid);
printf("\nNombre de liens stockes : %i\n\n",nln);
for (i=1 ; i <= nln ; i++ ) {
if ( MEDlienInfo(fid, i, nomlien, &nval) < 0) {
MESSAGE("Erreur a la demande d'information sur le lien n° : "); ISCRUTE(i);
ret = -1;continue;
}
printf("\t- Lien n°%i de nom |%s| et de taille %i\n",i,nomlien,nval);
lien = malloc((nval+1)*sizeof(char));
EXIT_IF(lien == NULL,NULL,NULL);
if ( MEDlienLire(fid, lien, nomlien) < 0 ) {
MESSAGE("Erreur a la lecture du lien : ");
SSCRUTE(nomlien);SSCRUTE(lien);
ret = -1;
} else
printf("\t\t|%s|\n\n",lien);
free(lien);
}
/* Interrogation des localisations des points de GAUSS */
nloc = MEDnGauss(fid);
printf("\nNombre de localisations stockees : %i\n\n",nloc);
for (i=1 ; i <= nloc ; i++ ) {
if ( MEDgaussInfo(fid, i, locname, &type_geo, &ngauss) < 0) {
MESSAGE("Erreur a la demande d'information sur la localisation n° : "); ISCRUTE(i);
ret = -1;continue;
}
printf("\t- Loc. n°%i de nom |%s| et nbr. de pts de GAUSS %i\n",i,locname,ngauss);
t1 = (type_geo%100)*(type_geo/100);
t2 = ngauss*(type_geo/100);
t3 = ngauss;
refcoo = (med_float *) malloc(sizeof(med_float)*t1 );
gscoo = (med_float *) malloc(sizeof(med_float)*t2 );
wg = (med_float *) malloc(sizeof(med_float)*t3 );
if ( MEDgaussLire(fid, refcoo, gscoo, wg, USER_INTERLACE, locname ) < 0) {
MESSAGE("Erreur a la lecture des valeurs de la localisation : ");
SSCRUTE(locname);
ret = -1;
} else {
printf("\t Coordonnees de l'element de reference de type %i :\n\t\t",type_geo);
for (j=0;j<t1;j++) printf(" %f ",*(refcoo+j));
printf("\n");
printf("\t Localisation des points de GAUSS : \n\t\t");
for (j=0;j<t2;j++) printf(" %f ",*(gscoo+j));
printf("\n");
printf("\t Poids associes aux points de GAUSS :\n\t\t");
for (j=0;j<t3;j++) printf(" %f ",*(wg+j));
printf("\n\n");
}
free(refcoo);
free(gscoo);
free(wg);
}
/* fermeture du fichier */
if ( MEDfermer(fid) < 0) return -1;
return ret;
}
int MEDimport(char * filein, char * fileout) {
med_idt fid, gid;
med_err ret;
med_int majeur, mineur, release;
med_bool hdfok=MED_FALSE;
med_bool medok=MED_FALSE;
char *_fileout,*tmp=NULL;
int _fileoutsize;
bool hasfileout=false;
char *commande;
med_int nprofil;
char chemin_profils[MED_TAILLE_PROFILS+1];
char chemin_liens[MED_TAILLE_LIENS+1];
char version[9];
int MAJ_21_22 = 0, MAJ_231_232 = 0, MAJ_236_300 = 0, MAJ_300_310 = 0, MAJ_310_320 = 0 ;
#ifdef PPRO_NT
char *drive, *dir, *ext;
#endif
unsigned char reponse='o';
med_bool _noversion=MED_FALSE;
EXIT_IF(filein == NULL,"Le nom du fichier d'entrée est vide : ", filein);
hasfileout = strcmp(fileout,"");
if ( hasfileout ) {
_fileoutsize = strlen(fileout);
_fileout = fileout;
} else {
_fileoutsize = strlen(filein)+strlen(PACKAGE_VERSION);
tmp = (char *) malloc(sizeof(char)*(_fileoutsize+1));
strcpy(tmp,filein);
strcat(tmp,PACKAGE_VERSION);
#ifdef PPRO_NT
_splitpath( tmp, drive, dir, _fileout, ext );
#else
_fileout = basename(tmp);
#endif
_fileoutsize = strlen(_fileout);
}
/* Test du format du fichier */
ret = MEDfileCompatibility(filein,&hdfok,&medok);
if (ret < 0 ) {
fprintf(stdout,">>> Attention le fichier %s ne contient pas de numéro de version. \n",filein);
fprintf(stdout,">>> Le fichier %s est supposé être en version 2.1.1. \n",filein);
/* PAs d'interactif dans une bibliothèque !*/
/* fprintf(stdout,">>> Voulez-vous essayer une conversion d'un fichier < 2.2 (o/n) ? "); */
/* scanf("%c",&reponse); */
if ( (reponse != 'o') && (reponse != 'O') && (reponse != 'y') && (reponse != 'Y') ) {
EXIT_IF(MEDfileCompatibility(filein,&hdfok,&medok) < 0,
"Erreur d'appel de MEDfileCompatibility : ", filein);
}
_noversion = MED_TRUE;
}
EXIT_IF( !hdfok ,
"Le fichier d'entrée n'est pas dans un format HDF compatible : ", filein);
/* EXIT_IF( !medok , */
/* "MEDimport ne gère pas le format MED de ce fichier : ", filein); */
/* Creation et ouverture du fichier que l'on va convertir au format MED actuel */
commande = (char *) malloc(sizeof(char)*(strlen("cp ")+strlen(filein)+
strlen(" ")+_fileoutsize + 4 +1 ) );
EXIT_IF(commande == NULL,NULL,NULL);
strcpy(commande,"cp \"");
strcat(commande,filein);
strcat(commande,"\" \"");
strcat(commande,_fileout);
strcat(commande,"\"");
fprintf(stdout,">>> Creation du fichier %s : %s \n",_fileout,commande);
system(commande);
free(commande);
commande = (char *) malloc(sizeof(char)*(strlen("chmod u+w \"") + _fileoutsize +1 +1 ) );
EXIT_IF(commande == NULL,NULL,NULL);
strcpy(commande,"chmod u+w \"");
strcat(commande,_fileout);
strcat(commande,"\"");
fprintf(stdout,">>> Chmod +w du fichier %s : %s \n",_fileout,commande);
system(commande);
free(commande);
fid = MEDfileOpen(_fileout,MED_ACC_RDWR);
EXIT_IF(fid < 0,"Ouverture du fichier : ", _fileout);
/* Verification du numero de version */
if (! _noversion)
ret = MEDfileNumVersionRd(fid,&majeur,&mineur,&release);
else {
ret=0;
majeur=2;
mineur=1;
release=1;
}
sprintf(version, IFORMAT"_"IFORMAT"_"IFORMAT, majeur, mineur, release);
EXIT_IF(ret < 0,"Lecture du numero de version de MED-fichier",NULL);
if (strcmp(version, "2_2_0") < 0)
MAJ_21_22 = 1;
if (strcmp(version, "2_3_2") < 0)
MAJ_231_232 = 1;
if (strcmp(version, "3_0_0") < 0)
MAJ_236_300 = 1;
if (strcmp(version, "3_1_0") < 0)
MAJ_300_310 = 1;
if (strcmp(version, "3_2_0") < 0)
MAJ_310_320 = 1;
/* Ne pas oublier que la version cible du fichier à convertir est celui de la bibliothèque. */
if (MAJ_310_320 == 0) {
fprintf(stdout,"Le fichier %s est déjà au bon format !!! \n",_fileout);
ret = MEDfileClose(fid);
EXIT_IF(ret < 0,"Fermeture du fichier",filein);
return 0;
}
/* On avertit qu'on commence la conversion */
fprintf(stdout,">>> Lancement de la normalisation du fichier selon le format " PACKAGE_VERSION " ...\n");
/* On inhibe le gestionnaire d'erreur HDF5 */
_MEDmodeErreurVerrouiller();
/* Mise a jour du numero de version */
fprintf(stdout,"- Lancement de la mise à jour du numéro de version ... \n");
/* La mise à jour MAJ_version(fid); doit être différée pour que les majs des fichiers anciens
fonctionnent correctement*/
/* MAJ_version(fid); */
MAJ_write_version_num(fid,2,3,6);
fprintf(stdout," Numéro de version : ... OK ... \n");
if (MAJ_21_22) {
/* Mise a jour des maillages : type = MED_NON_STRUCTURE, description, ... */
fprintf(stdout,"- Lancement de la mise à jour des maillages (21_22)... \n");
MAJ_21_22_maillages(fid);
fprintf(stdout," Maillage(s) : ... OK ...\n");
/* Mise a jour des champs */
fprintf(stdout,"- Lancement de la mise à jour des champs de résultats (21_22)... \n");
MAJ_21_22_champs(fid);
fprintf(stdout," Champs(s) : ... OK ...\n");
/* Mise a jour des profils eventuels */
nprofil = MEDnProfil(fid);
if (nprofil > 0) {
fprintf(stdout,"- Lancement de la mise à jour des profils (21_22)... \n");
MAJ_21_22_profils(fid,nprofil);
fprintf(stdout," Profils(s) : ... OK ...\n");
} else {
strncpy(chemin_profils,MED_PROFILS,MED_TAILLE_PROFILS-1);
chemin_profils[MED_TAILLE_PROFILS-1] = '\0';
gid = _MEDdatagroupCreer(fid,chemin_profils);
EXIT_IF(gid < 0,"Creation du groupe HDF sur les profils",chemin_profils);
ret = _MEDdatagroupFermer(gid);
EXIT_IF(ret < 0,"Fermeture du groupe HDF sur les profils",chemin_profils);
}
/* On cree le groupe HDF pour les liens */
strncpy(chemin_liens,MED_LIENS,MED_TAILLE_LIENS-1);
chemin_liens[MED_TAILLE_LIENS-1] = '\0';
gid = _MEDdatagroupCreer(fid,chemin_liens);
EXIT_IF(gid < 0,"Creation du groupe HDF sur les liens",chemin_liens);
ret = _MEDdatagroupFermer(gid);
EXIT_IF(ret < 0,"Fermeture du groupe HDF sur les liens",chemin_liens);
}
if (MAJ_231_232) {
/* Mise a jour des champs */
fprintf(stdout,"- Lancement de la mise à jour des champs de résultats (231_232)... \n");
MAJ_231_232_champs(fid);
fprintf(stdout," Champs(s) : ... OK ...\n");
fprintf(stdout,"- Lancement de la mise à jour des noms de maillages (231_232)... \n");
MAJ_231_232_maillages(fid);
fprintf(stdout," Noms(s) de maillage(s): ... OK ...\n");
}
if (MAJ_236_300) {
/* Le système de cache de version a été developpé à partir de la 3.0*/
/* Initialise le cache sur une 2.3.6 (cas d'absence d'INFO)*/
_MEDfileVersion(fid);
/* Mise a jour des champs */
fprintf(stdout,"- Lancement de la mise à jour des champs de résultats (236_300)... \n");
MAJ_236_300_champs(fid);
fprintf(stdout," Champs(s) : ... OK ...\n");
/* MAJ_version(fid); */
fprintf(stdout,"- Lancement de la mise à jour des maillages (236_300)... \n");
MAJ_236_300_maillages(fid);
fprintf(stdout," Maillage(s): ... OK ...\n");
/* MAJ_version(fid); */
}
if (MAJ_300_310) {
/* Le système de cache de version a été developpé à partir de la 3.0*/
/* Initialise le cache sur une 3.0.8 (cas d'absence d'INFO)*/
/* s'il n'a pas déjà été instanciée ds les MAJ précédentes */
MAJ_write_version_num(fid,3,0,8);
_MEDfileVersion(fid);
/* Si le cache était dèjà instancié, met à jour le cache */
MAJ_version_num(fid,3,0,8);
/* Mise a jour des champs */
fprintf(stdout,"- Lancement de la mise à jour des champs de résultats (300_310) ... \n");
MAJ_300_310_champs(fid);
fprintf(stdout," Champs(s) : ... OK ...\n");
}
if (MAJ_310_320) {
/* Le système de cache de version a été developpé à partir de la 3.0*/
/* Initialise le cache sur une 3.0.8 (cas d'absence d'INFO)*/
/* s'il n'a pas été déjà été instanciée ds les MAJ_ précédentes */
MAJ_write_version_num(fid,3,1,0);
_MEDfileVersion(fid);
/* Si le cache était dèjà instancié, met à jour le cache */
MAJ_version_num(fid,3,1,0);
/* Mise a jour des familles/groupes */
fprintf(stdout,"- Lancement de la mise à jour des familles/groupes (310_320) ... \n");
MAJ_310_320_familles(fid);
fprintf(stdout," Famille(s)/Groupe(s) : ... OK ...\n");
}
/* A l'écriture d'une nouvelle version de MAJ ex 310_320,
il est necessaire de revisiter les appels à MAJ_version(fid) pour
les remplacer par les appels MAJ_version(fid,3,1,Lastest31z) */
MAJ_version(fid);
MAJ_write_version_num(fid,MED_NUM_MAJEUR,MED_NUM_MINEUR,MED_NUM_RELEASE);
/* Fermeture du fichier */
ret = MEDfileClose(fid);
EXIT_IF(ret < 0,"Fermeture du fichier",_fileout);
/* On avertit que c'est fini */
fprintf(stdout,">>> Conversion du fichier %s au format MED V" PACKAGE_VERSION " terminée\n",
_fileout);
/* On libere la memoire */
if (!hasfileout) free(tmp);
return 0;
}
