struct apoint *
apoint_scan (FILE *f, struct tm start, struct tm end, char state, char *note)
{
char buf[BUFSIZ], *newline;
time_t tstart, tend, t;
t = time (NULL);
(void)localtime (&t);
/* Read the appointment description */
(void)fgets (buf, sizeof buf, f);
newline = strchr (buf, '\n');
if (newline)
*newline = '\0';
start.tm_sec = end.tm_sec = 0;
start.tm_isdst = end.tm_isdst = -1;
start.tm_year -= 1900;
start.tm_mon--;
end.tm_year -= 1900;
end.tm_mon--;
tstart = mktime (&start);
tend = mktime (&end);
EXIT_IF (tstart == -1 || tend == -1 || tstart > tend,
_("date error in appointment"));
return (apoint_new (buf, note, tstart, tend - tstart, state));
}
/* Free an already created scrollwin. */
void
wins_scrollwin_delete (struct scrollwin *sw)
{
EXIT_IF (sw == NULL, "null pointer");
delwin(sw->win.p);
delwin(sw->pad.p);
}
/*
* Create a new window and its associated pad, which is used to make the
* scrolling faster.
*/
void
wins_scrollwin_init (struct scrollwin *sw)
{
EXIT_IF (sw == NULL, "null pointer");
sw->win.p = newwin (sw->win.h, sw->win.w, sw->win.y, sw->win.x);
sw->pad.p = newpad (sw->pad.h, sw->pad.w);
sw->first_visible_line = 0;
sw->total_lines = 0;
}
static void
apoint_dup (struct apoint *in, struct apoint *bkp)
{
EXIT_IF (!in || !bkp, _("null pointer"));
bkp->start = in->start;
bkp->dur = in->dur;
bkp->state = in->state;
bkp->mesg = mem_strdup (in->mesg);
if (in->note)
bkp->note = mem_strdup (in->note);
}
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;
}
med_err getFieldsOn(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp,
med_entite_maillage entite, med_mode_switch stockage) {
int j,k,l,m,n,nb_geo;
med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,*vale=NULL,nval;
med_int numdt=0,numo=0,lnsize,nbrefmaa;
med_float *valr=NULL,dt=0.0;
med_err ret=0;
med_booleen local;
char pflname [MED_TAILLE_NOM+1]="";
char locname [MED_TAILLE_NOM+1]="";
char * lien = NULL;
char maa_ass [MED_TAILLE_NOM+1]="";
char dt_unit [MED_TAILLE_PNOM+1]="";
med_geometrie_element * type_geo;
med_geometrie_element typ_noeud[1] = { MED_NONE };
med_geometrie_element typmai[MED_NBR_GEOMETRIE_MAILLE+2] = {MED_POINT1, MED_SEG2, MED_SEG3, MED_TRIA3,
MED_QUAD4, MED_TRIA6,MED_QUAD8, MED_TETRA4,
MED_PYRA5, MED_PENTA6, MED_HEXA8, MED_TETRA10,
MED_PYRA13, MED_PENTA15, MED_HEXA20,
MED_POLYGONE, MED_POLYEDRE};
med_geometrie_element typfac[MED_NBR_GEOMETRIE_FACE+1] = {MED_TRIA3,MED_TRIA6,
MED_QUAD4,MED_QUAD8,
MED_POLYGONE};
med_geometrie_element typare[MED_NBR_GEOMETRIE_ARETE] = {MED_SEG2,MED_SEG3};
char ** AFF;
switch (entite) {
case MED_NOEUD :
type_geo = typ_noeud;
nb_geo = 1;
AFF = MED_GEOMETRIE_NOEUD_AFF;
break;
case MED_MAILLE :
type_geo = typmai;
nb_geo = MED_NBR_GEOMETRIE_MAILLE+2;
AFF = MED_GEOMETRIE_MAILLE_AFF;
break;
case MED_FACE :
type_geo = typfac;
nb_geo = MED_NBR_GEOMETRIE_FACE+1;
AFF = MED_GEOMETRIE_FACE_AFF;
break;
case MED_ARETE :
type_geo = typare;
nb_geo = MED_NBR_GEOMETRIE_ARETE;
AFF = MED_GEOMETRIE_ARETE_AFF;
break;
}
for (k=0;k<nb_geo;k++) {
/* Combien de (PDT,NOR) a lire */
nbpdtnor = MEDnPasdetemps(fid,nomcha,entite,type_geo[k]);
if (nbpdtnor < 1 ) continue;
for (j=0;j<nbpdtnor;j++) {
if ( MEDpasdetempsInfo(fid,nomcha,entite,type_geo[k],
j+1, &ngauss, &numdt, &numo, dt_unit,
&dt, maa_ass, &local, &nbrefmaa) <0) {
MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
ISCRUTE(numdt); ISCRUTE(numo);
ret = -1; continue;
};
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);
printf("\tLe maillage par defaut est : |%s|, sur un total de : %i maillages associes\n",
maa_ass, nbrefmaa);
/* Le maillage reference est-il porte par un autre fichier */
if ( !local ) {
if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 ) {
MESSAGE("Erreur a la lecture de la taille du lien : ");
SSCRUTE(maa_ass);
ret = -1;
} else {
lien = malloc(lnsize*sizeof(char));
EXIT_IF(lien == NULL,NULL,NULL);
if ( MEDlienLire(fid, lien, maa_ass) < 0 ) {
MESSAGE("Erreur a la lecture du lien : ");
SSCRUTE(maa_ass);SSCRUTE(lien);
ret = -1;
} else {
printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
}
free(lien);
}
}
/* Combien de maillages lies aux (nomcha,ent,geo,numdt,numo) */
/* Notons que cette information est egalement disponible a partir de MEDpasdetempsInfo */
if ( (nbrefmaa = MEDnChampRef(fid,nomcha,entite,type_geo[k],numdt,numo) ) < 0 ) {
MESSAGE("Erreur a la demande du nombre de maillages references par le champ : ");
SSCRUTE(nomcha); ISCRUTE(numdt); ISCRUTE(numo);
ret = -1; continue;
};
for (l=0;l<nbrefmaa;l++) {
if ( MEDchampRefInfo(fid,nomcha,entite,type_geo[k],
l+1,numdt, numo, maa_ass, &local, &ngauss) <0 ) {
MESSAGE("Erreur a la demande d'information sur le maillage utilise par le champ n° : ");
SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
ISCRUTE(l+1);ISCRUTE(numdt); ISCRUTE(numo);SSCRUTE(maa_ass);
ret = -1; continue;
};
/* Prend en compte le nbre de pt de gauss automatiquement */
if ((nval = MEDnVal(fid,nomcha,entite,type_geo[k],numdt,numo,maa_ass,USER_MODE)) <= 0) {
MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);ISCRUTE(USER_MODE);
ret = -1; continue;
};
printf("\t- Il y a %d valeurs en mode %i. Chaque entite %s\
de type geometrique %s associes au maillage |%s| a %i pts de gauss \n",
nval,USER_MODE,MED_ENTITE_MAILLAGE_AFF[(int)entite],AFF[k],maa_ass,ngauss);
/* Le maillage reference est-il porte par un autre fichier */
if ( !local ) {
if ( (lnsize=MEDnValLien(fid,maa_ass) ) < 0 ) {
MESSAGE("Erreur a la lecture de la taille du lien : ");
SSCRUTE(maa_ass);
ret = -1;
} else {
lien = malloc(lnsize*sizeof(char));
EXIT_IF(lien == NULL,NULL,NULL);
if ( MEDlienLire(fid, lien, maa_ass) < 0 ) {
MESSAGE("Erreur a la lecture du lien : ");
SSCRUTE(maa_ass);SSCRUTE(lien);
ret = -1;
} else {
printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",maa_ass,lien);
}
free(lien);
}
}
/*Lecture des valeurs du champ */
if (typcha == MED_FLOAT64) {
valr = (med_float*) calloc(ncomp*nval,sizeof(med_float));
EXIT_IF(valr == NULL,NULL,NULL);
if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)valr,stockage,MED_ALL,locname,
pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
ret = -1;
};
} else {
vale = (med_int*) calloc(ncomp*nval,sizeof(med_int));
EXIT_IF(vale == NULL,NULL,NULL);
if ( MEDchampLire(fid,maa_ass,nomcha,(unsigned char*)vale,stockage,MED_ALL,locname,
pflname,USER_MODE,entite,type_geo[k],numdt,numo) < 0 ) {
MESSAGE("Erreur a la lecture des valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE(entite);ISCRUTE(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
ret = -1;
};
}
if ( ngauss > 1 )
printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
switch (stockage) {
case MED_FULL_INTERLACE :
printf("\t- Valeurs :\n\t");
for (m=0;m<nval/ngauss;m++) {
printf("|");
for (n=0;n<ngauss*ncomp;n++)
if (typcha == MED_FLOAT64)
printf(" %f ",*(valr+(m*ngauss*ncomp)+n));
else
printf(" %d ",*(vale+(m*ngauss*ncomp)+n));
}
break;
/*Affichage en fonction du profil à traiter*/
case MED_NO_INTERLACE :
printf("\t- Valeurs :\n\t");
for (m=0;m<ncomp;m++) {
printf("|");
for (n=0;n<nval;n++)
if (typcha == MED_FLOAT64)
printf(" %f ",*(valr+(m*nval)+n));
else
printf(" %d ",*(vale+(m*nval)+n));
}
break;
}
printf("|\n");
if (typcha == MED_FLOAT64) {
if ( valr ) {free(valr);valr = NULL;}}
else
if (vale) { free(vale);vale = NULL; }
/*Lecture du profil associe */
if (strcmp(pflname,MED_NOPFL) == 0 )
printf("\t- Profil : MED_NOPFL\n");
else {
if ( (pflsize = MEDnValProfil(fid,pflname)) <0 ) {
MESSAGE("Erreur a la lecture du nombre de valeurs du profil : ");
SSCRUTE(pflname);
ret = -1; continue;
}
printf("\t- Profil : |%s| de taille %i\n",pflname,pflsize);
pflval = (med_int*) malloc(sizeof(med_int)*pflsize);
EXIT_IF(pflval == NULL,NULL,NULL);
if ( MEDprofilLire(fid,pflval,pflname) <0) {
MESSAGE("Erreur a la lecture des valeurs du profil : ");
SSCRUTE(pflname);
ret = -1;
}
printf("\t");
for (m=0;m<pflsize;m++) printf(" %i ",*(pflval+m));
printf("\n");
free(pflval);
}
}
}
} /* fin for sur les mailles*/
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;
}
med_err
MED231champRefInfoEtRenMaa(med_idt fid,char *champ,
med_entite_maillage type_ent, med_geometrie_element type_geo,
int indice, med_int numdt, med_int numo,
char * maa, med_booleen * local, med_int *ngauss)
{
med_err ret=-1;
int num;
med_idt datagroup2=0,datagroup3=0,gid_maa=0,gid_lien=0;
char chemin[(MED_TAILLE_CHA+MED_TAILLE_NOM+1)+(2*MED_TAILLE_NOM_ENTITE+2)+(2*MED_MAX_PARA+1)+(MED_TAILLE_NOM)+1]="";
char chemini[(MED_TAILLE_CHA+MED_TAILLE_NOM+1)+(2*MED_TAILLE_NOM_ENTITE+2)+(2*MED_MAX_PARA+1)+(MED_TAILLE_NOM)+1]="";
char chemin_maa[MED_TAILLE_MAA+MED_TAILLE_NOM+1]="";
char chemin_lien[MED_TAILLE_LIENS+MED_TAILLE_NOM+1]="";
char nomdatagroup1[2*MED_TAILLE_NOM_ENTITE+2]="";
char nomdatagroup2[2*MED_MAX_PARA+1]="";
char tmp1 [MED_TAILLE_NOM_ENTITE+1]="";
char maai [MED_TAILLE_NOM+1];
char maaf [MED_TAILLE_NOM+1]="";
/*
* On inhibe le gestionnaire d'erreur HDF 5
*/
_MEDmodeErreurVerrouiller();
/*
* On construit le nom du datagroup
*/
strcpy(chemin,MED_CHA);
strcat(chemin,champ);
strcat(chemin,"/");
/*
* Si le Data Group de niveau 1 <type_ent>[.<type_geo>] n'existe pas => erreur
*/
/* modif pour la version 2.3.3 */
if (_MEDnomEntite(nomdatagroup1,type_ent) < 0)
goto ERROR;
if ((type_ent != MED_NOEUD)) {
if (_MEDnomGeometrie(tmp1,type_geo) < 0)
goto ERROR;
strcat(nomdatagroup1,".");
strcat(nomdatagroup1,tmp1);
}
strcat(chemin,nomdatagroup1);
strcat(chemin,"/");
/*
* Si le Data Group de niveau 2 <numdtt>.<numoo> n'existe pas => erreur
*/
sprintf(nomdatagroup2,"%*li%*li",MED_MAX_PARA,(long ) numdt,MED_MAX_PARA,(long ) numo);
strcat(chemin,nomdatagroup2);
/*
* Modifie le nom de la première référence à un maillage
* si besoin est
*/
if ( (datagroup2 = _MEDdatagroupOuvrir(fid,chemin)) < 0 ) {
MESSAGE("Erreur à l'ouverture du datagroup : ");
SSCRUTE(chemin); goto ERROR;
}
if ( _MEDattrStringLire(datagroup2,MED_NOM_MAI,MED_TAILLE_NOM,maai) < 0 ) {
MESSAGE("Erreur de lecture de l'attribut MED_NOM_MAI : ");
SSCRUTE(maai); goto ERROR;
}
if ( MAJ_231_232_chaine(maai,maaf) ) {
fprintf(stdout," >>> Normalisation du nom de maillage par défaut [%s] associé au champ [%s] pour (n°dt,n°it) ("IFORMAT","IFORMAT")\n",maai,champ,numdt,numo);
ret = _MEDattrStringEcrire(datagroup2,MED_NOM_MAI,MED_TAILLE_NOM,maaf);
EXIT_IF(ret < 0,"Renommage du maillage en",maaf);
fprintf(stdout," >>> Normalisation du nom du maillage par défaut [%s] ... OK ... \n",maaf);
}
if (_MEDdatagroupFermer(datagroup2) < 0) {
MESSAGE("Impossible de fermer le datagroup : ");
ISCRUTE_int(datagroup2); goto ERROR;
}
strcat(chemin,"/");
/*
* Cherche le datagroup de niveau 3 <maa> correspondant à l'indice <num>
*/
num = indice - 1;
if (_MEDobjetIdentifier(fid,chemin,num,maa) < 0) {
MESSAGE("Impossible de trouver un groupe à l'indice spécifié : ");
SSCRUTE(chemin); ISCRUTE_int(num); goto ERROR;
};
if ( MAJ_231_232_chaine(maa,maaf) ) {
fprintf(stdout," >>> Normalisation du nom de maillage [%s] associé au champ [%s] pour (n°dt,n°it) ("IFORMAT","IFORMAT")\n",maa,champ,numdt,numo);
/* on accede au maillage */
strcpy(chemini,chemin);
strcat(chemini,maa);
strcat(chemin,maaf);
ret = H5Gmove(fid, chemini, chemin );
EXIT_IF(ret < 0,"Renommage du maillage en",maaf);
fprintf(stdout," >>> Normalisation du nom du maillage [%s] ... OK ... \n",maaf);
} else
strcat(chemin,maa);
/*
* Si le Data Group de niveau 3 <maa> n'existe pas => erreur
*/
if ((datagroup3 = _MEDdatagroupOuvrir(fid,chemin)) < 0) {
MESSAGE("Erreur d'ouverture du datagroup lien au maillage : ");
SSCRUTE(chemin); goto ERROR;
};
/* Lire le nbre des points de GAUSS*/
if (_MEDattrEntierLire(datagroup3,MED_NOM_NGA,ngauss) < 0) {
MESSAGE("Erreur à la lecture de l'attribut MED_NOM_NGA : ");
ISCRUTE(*ngauss);goto ERROR;
};
/* Maillage local ou distant */
/* Les noms de maillages n'ayant pas encore été mis à jour
on garde l'ancien nom pour le test local/distant */
strcpy(chemin_maa,MED_MAA);
strcat(chemin_maa,maa);
/* Le maillage est il distant */
if ( (gid_maa = _MEDdatagroupOuvrir(fid,chemin_maa)) < 0) {
/* Verifie que le maillage est bien référencé comme distant */
strcpy(chemin_lien,MED_LIENS);
strcat(chemin_lien,maa);
if ((gid_lien = _MEDdatagroupOuvrir(fid,chemin_lien)) < 0) {
/* MESSAGE("Le maillage n'est ni local, ni distant : "); */
/* SSCRUTE(chemin_maa);SSCRUTE(chemin_lien); goto ERROR; */
*local = MED_FAUX;
}
*local = MED_FAUX;
} else
*local = MED_VRAI;
/*On retourne le nouveau nom de maillage*/
if ( strlen(maaf) ) strcpy(maa,maaf);
/*
* On ferme tout
*/
ret = 0;
ERROR:
if (datagroup3>0) if (_MEDdatagroupFermer(datagroup3) < 0) {
MESSAGE("Impossible de fermer le datagroup : ");
ISCRUTE_int(datagroup3); ret = -1;
}
if (gid_maa>0) if (_MEDdatagroupFermer(gid_maa) < 0) {
MESSAGE("Impossible de fermer le datagroup : ");
ISCRUTE_id(gid_maa); ret = -1;
}
if (gid_lien>0) if (_MEDdatagroupFermer(gid_lien) < 0) {
MESSAGE("Impossible de fermer le datagroup : ");
SSCRUTE(chemin_lien); ret = -1;
}
return ret;
}
void MAJ_231_232_champs(med_idt fid)
{
med_int ncomp,ncha;
med_err lret,rett;
med_type_champ typcha;
char nomcha [MED_TAILLE_NOM+1]="";
char *comp, *unit;
int i;
/* combien de champs dans le fichier */
ncha = MEDnChamp(fid,0);
EXIT_IF(ncha < 0,"lors de la lecture du nombre de champs",NULL);
/****************************************************************************
* LECTURE DES CHAMPS *
****************************************************************************/
/* lecture de tous les champs */
for (i =0;i<ncha;i++) {
lret = 0;
/* Lecture du nombre de composantes */
ncomp = MEDnChamp(fid,i+1);
if (ncomp < 0) {
MESSAGE("Erreur à la lecture du nombre de composantes : "); ISCRUTE(ncomp);
exit(1);
}
/* Lecture du type du champ, des noms des composantes et du nom de l'unité*/
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);
rett = MED231champInfoEtRen(fid,i+1,nomcha,&typcha,comp,unit,ncomp);
if ( rett < 0 ) {
MESSAGE("Erreur à la demande d'information sur les champs ");
exit(1);
}
free(comp);
free(unit);
lret = MED231champNormaliser(fid, nomcha, typcha, ncomp, MED_NOEUD);
if (lret != 0) {
MESSAGE("Erreur à la lecture des champs aux noeuds "); exit(1);
}
lret = MED231champNormaliser(fid, nomcha, typcha, ncomp, MED_MAILLE);
if (lret != 0) {
MESSAGE("Erreur à la lecture des champs aux mailles "); exit(1);
}
lret = MED231champNormaliser(fid, nomcha, typcha, ncomp, MED_FACE);
if (lret != 0) {
MESSAGE("Erreur à la lecture des champs aux faces "); exit(1);
}
lret = MED231champNormaliser(fid, nomcha, typcha, ncomp, MED_ARETE);
if (lret != 0) {
MESSAGE("Erreur à la lecture des champs aux aretes "); exit(1);
}
}
}
void MAJ_236_300_champs(med_idt fid)
{
med_err lret,ret;
/* med_idt _datagroup=0; */
med_field_type typcha;
char nomcha [MED_NAME_SIZE+1]="";
char _meshname [MED_NAME_SIZE+1]="";
char _dtunit [MED_SNAME_SIZE+1]="";
char *comp= NULL, *unit= NULL;
med_int ncomp,ncha;
med_int _ncstp=0;
med_bool _local=MED_FALSE;
htri_t _datasetexist;
char _pathi[MED_TAILLE_CHA+1+MED_NAME_SIZE+1]=MED_CHA;
char _pathf[MED_TAILLE_CHA+2+MED_NAME_SIZE+1]="/CHA_/";
char _pathtmp[MED_TAILLE_CHA+3]="/CHA__/";
int i,j;
char nomlien[MED_NAME_SIZE+1]="";
char * lien = NULL;
med_int nln,nval;
med_int _nloc,_intgeotype,_sdim;
char _pathloc[MED_TAILLE_GAUSS+MED_NAME_SIZE+1]=MED_GAUSS;
med_int _npar,_numdt,_numit;
char _pathpari[MED_TAILLE_NUM_DATA+MED_NAME_SIZE+1+2*MED_MAX_PARA+1+1]=MED_NUM_DATA;
char _pathparf[MED_TAILLE_NUM_DATA+MED_NAME_SIZE+1+2*MED_MAX_PARA+1+1]=MED_NUM_DATA;
int _pathparlen;
med_size _n=0;
char _cpstnamei[2*MED_MAX_PARA+1]="";
char _cpstnamef[2*MED_MAX_PARA+1]="";
char _uniname[MED_SNAME_SIZE+1]="";
/* hid_t _lac_plist_id; */
hid_t _lcp_plist_id;
hid_t _ocp_plist_id;
med_bool _createunt = MED_TRUE;
MAJ_version_num(fid,2,3,6);
/* MAJ des varaibles scalaires */
_npar = MEDnParameter(fid);
if (_npar > 0) {
fprintf(stdout," >>> Normalisation des paramètres scalaires\n");
_pathparf[MED_TAILLE_NUM_DATA-2]='_';
/* _lac_plist_id = H5Pcreate( H5P_LINK_ACCESS ); */
_ocp_plist_id = H5Pcreate( H5P_OBJECT_COPY );
_lcp_plist_id = H5Pcreate( H5P_LINK_CREATE );
H5Pset_create_intermediate_group( _lcp_plist_id, 1 );
H5Pset_copy_object( _ocp_plist_id, H5O_COPY_SHALLOW_HIERARCHY_FLAG);
}
for (i=0 ; i < _npar ; i++ ) {
MED_ERR_EXIT_IF (_MEDobjectGetName(fid, _pathpari ,i, &_pathpari[MED_TAILLE_NUM_DATA]) < 0,
MED_ERR_ACCESS,MED_ERR_DATAGROUP,_pathpari);
strcpy(&_pathparf[MED_TAILLE_NUM_DATA],&_pathpari[MED_TAILLE_NUM_DATA]);
/* SSCRUTE(_pathparf); */
/* SSCRUTE(_pathpari); */
/*Copie le group avec ses attributs et les objets de premier niveau.*/
ret = H5Ocopy(fid,_pathpari,fid,_pathparf,_ocp_plist_id,_lcp_plist_id);
/* ret = H5Lcopy(fid,_pathpari,fid,_pathparf,_lcp_plist_id,_lac_plist_id); */
EXIT_IF(ret < 0,"Copie du datagroup",_pathpari);
_pathparlen=strlen(_pathpari);
_pathpari[_pathparlen]='/';_pathparf[_pathparlen]='/';
++_pathparlen;
_pathpari[_pathparlen]='\0';_pathparf[_pathparlen]='\0';
/* SSCRUTE(_pathparf); */
/* SSCRUTE(_pathpari); */
ret =_MEDnObjects(fid,_pathpari, &_n);
MED_ERR_EXIT_IF( (ret == (MED_ERR_COUNT + MED_ERR_DATAGROUP)), MED_ERR_COUNT,MED_ERR_PARAMETER,_pathpari);
for (j=0 ; j < _n ; ++j ) {
MED_ERR_EXIT_IF (_MEDobjectGetName(fid, _pathpari ,j, &_pathpari[_pathparlen]) < 0,
MED_ERR_ACCESS,MED_ERR_DATAGROUP,_pathpari);
MED_ERR_EXIT_IF (_MEDattributeNumRdByName(fid,_pathpari,MED_NOM_NOR,
MED_INTERNAL_INT,(unsigned char * const ) &_numit) < 0,
MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_NOM_NOR);
MED_ERR_EXIT_IF (_MEDattributeNumRdByName(fid,_pathpari,MED_NOM_NDT,
MED_INTERNAL_INT,(unsigned char * const ) &_numdt) < 0,
MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_NOM_NDT);
MED_ERR_EXIT_IF (_MEDattributeStringRdByName(fid,_pathpari,MED_NOM_UNI,
MED_SNAME_SIZE,_uniname) < 0,
MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_NOM_NDT);
_MEDgetComputationStepName(MED_SORT_DTIT,_numdt,_numit,&_pathparf[_pathparlen]);
/* strcat(_pathpari,"/"); */
/* strcat(_pathparf,"/"); */
/* SSCRUTE(_pathparf); */
/* SSCRUTE(_pathpari); */
/* ret = H5Ocopy(fid,_pathpari,fid,_pathparf,_ocp_plist_id,_lcp_plist_id); */
/* ret = H5Lcopy(fid,_pathpari,fid,_pathparf,_lcp_plist_id,_lac_plist_id); */
/* H5Eprint1(stderr); */
/* EXIT_IF(ret < 0,"Copie d'une étape de calcul du paramètre scalaire ",_pathpari); */
/*On modifie temporairement _pathpari pour pointer dans _pathparf*/
_pathpari[MED_TAILLE_NUM_DATA-2]='_';
/* SSCRUTE(_pathparf); */
/* SSCRUTE(_pathpari); */
ret = H5Gmove(fid, _pathpari, _pathparf );
EXIT_IF(ret < 0,"Renommage de l'étape de calcul",_pathpari);
_pathpari[MED_TAILLE_NUM_DATA-2]='A';
MED_ERR_EXIT_IF(H5Adelete_by_name( fid, _pathparf, MED_NOM_UNI, H5P_DEFAULT ) < 0,
MED_ERR_DELETE,MED_ERR_ATTRIBUTE,_pathparf);
_pathparf[_pathparlen]='\0';
_pathpari[_pathparlen]='\0';
if ( _createunt ) {
MED_ERR_EXIT_IF (_MEDattributeStringWrByName(fid,_pathparf,MED_NOM_UNT, MED_SNAME_SIZE,_uniname) < 0,
MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_NOM_UNT);
_createunt = MED_FALSE;
}
}
_createunt = MED_TRUE;
_pathpari[MED_TAILLE_NUM_DATA]='\0';
}
if ( _npar > 0 ) {
_pathpari[MED_TAILLE_NUM_DATA]='\0';
_pathparf[MED_TAILLE_NUM_DATA]='\0';
MED_ERR_EXIT_IF ( H5Ldelete(fid,_pathpari,H5P_DEFAULT) < 0 ,
MED_ERR_DELETE,MED_ERR_LINK,_pathpari);
ret = H5Gmove(fid, _pathparf, _pathpari );
EXIT_IF(ret < 0,"Renommage du group de paramètres scalaires",_pathparf);
}
/* MAJ des localisations */
_nloc = MEDnLocalization(fid);
/* ISCRUTE(_nloc); */
if (_nloc > 0)
fprintf(stdout," >>> Normalisation des localisations des points d'intégration\n");
for (i=0 ; i < _nloc ; i++ ) {
/* SSCRUTE(_pathloc); */
MED_ERR_EXIT_IF (_MEDobjectGetName(fid, _pathloc ,i, &_pathloc[MED_TAILLE_GAUSS]) < 0,
MED_ERR_ACCESS,MED_ERR_DATAGROUP,_pathloc);
/* SSCRUTE(_pathloc); */
MED_ERR_EXIT_IF (_MEDattributeNumRdByName(fid,_pathloc,MED_NOM_GEO,
MED_INTERNAL_INT,(unsigned char * const ) &_intgeotype) < 0,
MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_NOM_GEO);
_sdim = (_intgeotype/100);
MED_ERR_EXIT_IF (_MEDattributeNumWrByName(fid,_pathloc,MED_NOM_DIM,
MED_INTERNAL_INT,(const unsigned char * const) &_sdim) < 0,
MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_NOM_DIM);
MED_ERR_EXIT_IF ( _MEDattributeStringWrByName(fid,_pathloc,MED_NOM_INM,MED_NAME_SIZE,"") < 0,
MED_ERR_WRITE,MED_ERR_ATTRIBUTE,MED_NOM_INM);
_pathloc[MED_TAILLE_GAUSS]='\0';
}
/* MAJ des liens */
nln = MEDnLink(fid);
if (nln > 0)
fprintf(stdout," >>> Normalisation des liens\n");
for (i=1 ; i <= nln ; i++ ) {
ret = MEDlinkInfo(fid, i, nomlien, &nval);
EXIT_IF(ret,"Erreur a la demande d'information sur le lien",NULL);
/* printf("\t- Lien n°%i de nom |%s| et de taille "IFORMAT"\n",i,nomlien,nval); */
lien = (char *) malloc((nval+1)*sizeof(char));
EXIT_IF(lien == NULL,NULL,NULL);
ret = MEDlinkRd(fid, nomlien, lien );
EXIT_IF(ret,"Erreur a la lecture du lien : ",nomlien);
MAJ_version_num(fid,3,0,8);
ret = MED30linkWr(fid,nomlien,lien);
EXIT_IF(ret,"Erreur a l'écrtiure du lien : ",nomlien);
MAJ_version_num(fid,2,3,6);
lien[nval] = '\0';
/* printf("\t\t|%s|\n\n",lien); */
free(lien);
}
/* combien de champs dans le fichier */
ncha = MEDnField(fid);
EXIT_IF(ncha < 0,"lors de la lecture du nombre de champs",NULL);
/* MAJ des champs */
for (i =0;i<ncha;i++) {
lret = 0;
/* Lecture du nombre de composantes */
ncomp = MEDfieldnComponent(fid,i+1);
if (ncomp < 0) {
MESSAGE("Erreur à la lecture du nombre de composantes : "); ISCRUTE(ncomp);
exit(1);
}
/* Lecture du type du champ, des noms des composantes et du nom de l'unité*/
comp = (char*) malloc(ncomp*MED_SNAME_SIZE+1);
EXIT_IF(comp == NULL,NULL,NULL);
unit = (char*) malloc(ncomp*MED_SNAME_SIZE+1);
EXIT_IF(unit == NULL,NULL,NULL);
/* ret = MED231champInfoEtRen(fid,i+1,nomcha,&typcha,comp,unit,ncomp); */
ret = MEDfieldInfo(fid,i+1,nomcha,_meshname,&_local,&typcha,comp,unit,_dtunit,&_ncstp);
MED_ERR_EXIT_IF(ret,MED_ERR_ACCESS,MED_ERR_FIELD,nomcha);
/* creation du champ destination */
MAJ_version_num(fid,3,0,8);
EXIT_IF( H5Gmove(fid, _pathi, _pathtmp ) < 0,"Switch to ",_pathtmp);
_datasetexist=H5Lexists( fid, _pathf, H5P_DEFAULT );
if (_datasetexist ) { EXIT_IF( (H5Gmove(fid, _pathf, _pathi ) < 0) ,"Switch to ",_pathi); }
MED_ERR_EXIT_IF( MEDfieldCr(fid,nomcha,typcha,ncomp,comp,unit,_dtunit,_meshname ) < 0,
MED_ERR_CREATE,MED_ERR_FIELD,_pathf);
EXIT_IF( H5Gmove(fid, _pathi , _pathf ) < 0,"Switch to ",_pathf);
EXIT_IF( H5Gmove(fid, _pathtmp, _pathi ) < 0,"Switch to ",_pathi);
MAJ_version_num(fid,2,3,6);
free(comp);
free(unit);
lret = MAJ_236_300_fieldOnEntity( fid, nomcha, _meshname, typcha, ncomp, MED_NODE,_ncstp, _pathi, _pathf);
if (lret != 0) {
MESSAGE("Erreur à la lecture des champs aux noeuds "); exit(1);
}
lret = MAJ_236_300_fieldOnEntity( fid, nomcha, _meshname, typcha, ncomp, MED_CELL,_ncstp, _pathi, _pathf);
if (lret != 0) {
MESSAGE("Erreur à la lecture des champs aux mailles "); exit(1);
}
lret = MAJ_236_300_fieldOnEntity( fid, nomcha, _meshname, typcha, ncomp, MED_DESCENDING_FACE,_ncstp, _pathi, _pathf);
if (lret != 0) {
MESSAGE("Erreur à la lecture des champs aux faces "); exit(1);
}
lret = MAJ_236_300_fieldOnEntity( fid, nomcha, _meshname, typcha, ncomp, MED_DESCENDING_EDGE,_ncstp, _pathi, _pathf);
if (lret != 0) {
MESSAGE("Erreur à la lecture des champs aux aretes "); exit(1);
}
lret = MAJ_236_300_fieldOnEntity( fid, nomcha, _meshname, typcha, ncomp, MED_NODE_ELEMENT,_ncstp, _pathi, _pathf);
if (lret != 0) {
MESSAGE("Erreur à la lecture des champs aux aretes "); exit(1);
}
}
_datasetexist=H5Lexists( fid, _pathf, H5P_DEFAULT );
if (_datasetexist ) {
EXIT_IF( (H5Ldelete(fid,_pathi, H5P_DEFAULT) < 0) ,"Delete ",_pathi);
EXIT_IF( (H5Gmove(fid, _pathf, _pathi ) < 0) ,"Switch to ",_pathf);
}
}
int main (int argc, char **argv)
{
med_err ret,lret;
med_idt fid;
char * fichier = NULL;
char maa[MED_NAME_SIZE+1]="";
char desc[MED_COMMENT_SIZE+1]="";
char pflname[MED_NAME_SIZE+1]="",nomlien[MED_NAME_SIZE+1]="";
char _meshname [MED_NAME_SIZE+1]="";
char _dtunit [MED_SNAME_SIZE+1]="";
char locname[MED_NAME_SIZE+1]="";
char * lien = NULL;
char *comp= NULL, *unit= NULL;
char nomcha [MED_NAME_SIZE+1]="";
med_int mdim=0,sdim=0,ncomp,ncha,npro,nln,pflsize,*pflval,nval;
med_int _ncstp=0,ngauss=0,nloc=0,locsdim=0,lnsize=0;
int t1,t2,t3;
med_field_type typcha;
med_geometry_type type_geo;
med_float *refcoo, *gscoo, *wg;
int i,j;
med_bool _local;
char dtunit[MED_SNAME_SIZE+1]="";
char nomcoo[3*MED_SNAME_SIZE+1]="";
char unicoo[3*MED_SNAME_SIZE+1]="";
char geointerpname[MED_SNAME_SIZE+1]="";
char ipointstructmeshname[MED_SNAME_SIZE+1]="";
med_mesh_type type;
med_sorting_type sort;
med_int nstep=0;
med_axis_type rep;
med_int nsectionmeshcell;
med_geometry_type sectiongeotype;
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 = MEDfileOpen(fichier,MED_ACC_RDONLY)) < 0){
MESSAGE("Erreur a l'ouverture du fichier : ");SSCRUTE(fichier);
return -1;
}
ret = 0;
/* Lecture des infos concernant le premier maillage */
if ( MEDmeshInfo( fid, 1, maa, &sdim, &mdim, &type, desc, dtunit, &sort,
&nstep, &rep, nomcoo,unicoo) < 0 ) {
MESSAGE("Erreur a la lecture des informations sur le maillage : ");SSCRUTE(maa);
return -1;
} else {
printf("Maillage de nom : |%s| , de dimension : %d , et de type %d\n",maa,mdim,type);
printf("\t -Dimension de l'espace : %d\n",sdim);
printf("\t -Description du maillage : %s\n",desc);
printf("\t -Noms des axes : |%s|\n",nomcoo);
printf("\t -Unités des axes : |%s|\n",unicoo);
printf("\t -Type de repère : %d\n",rep);
printf("\t -Nombre d'étapes de calcul : %d\n",nstep);
printf("\t -Unité des dates : |%s|\n",dtunit);
}
/* combien de champs dans le fichier */
if ((ncha = MEDnField(fid)) < 0) {
MESSAGE("Impossible de lire le nombre de champs : ");ISCRUTE(ncha);
return ncha;
}
printf("Nombre de champs : "IFORMAT" \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 = MEDfieldnComponent(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_SNAME_SIZE+1);
EXIT_IF(comp == NULL,NULL,NULL);
unit = (char*) malloc(ncomp*MED_SNAME_SIZE+1);
EXIT_IF(unit == NULL,NULL,NULL);
if ( MEDfieldInfo(fid,i+1,nomcha,_meshname,&_local,&typcha,comp,unit,_dtunit,&_ncstp) < 0 ) {
MESSAGE("Erreur a la demande d'information sur les champs : ");
ISCRUTE_int(i+1);SSCRUTE(nomcha);ISCRUTE_int(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);
printf("Unites des dates : |%s| \n",_dtunit);
printf("Le maillage associé est |%s|\n",_meshname);
printf("Nombre de séquences de calcul |%d|\n",_ncstp);
/* Le maillage reference est-il porte par un autre fichier */
if ( !_local ) {
if ( (lnsize=MEDlinkInfoByName(fid,_meshname) ) < 0 ) {
MESSAGE("Erreur a la lecture de la taille du lien : ");
SSCRUTE(_meshname);
ret = -1;
} else {
lien = malloc(lnsize*sizeof(char));
EXIT_IF(lien == NULL,NULL,NULL);
if ( MEDlinkRd(fid,_meshname, lien) < 0 ) {
MESSAGE("Erreur a la lecture du lien : ");
SSCRUTE(_meshname);SSCRUTE(lien);
ret = -1;
} else {
printf("\tLe maillage |%s| est porte par un fichier distant |%s|\n",_meshname,lien);
}
free(lien);
}
}
free(comp);
free(unit);
/*TODO : Créer les API30 spécifiques pour la lecture de champs multi-maillages*/
if (strcmp(nomcha,"champ entier")) {
MESSAGE("There is no API yest for reading field on multiple meshes"); continue;
}
lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NODE, USER_INTERLACE,_ncstp );
if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_CELL, USER_INTERLACE,_ncstp );
else { MESSAGE("Erreur a la lecture des champs aux noeuds "); ret = -1; continue;}
if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_DESCENDING_FACE,USER_INTERLACE,_ncstp);
else { MESSAGE("Erreur a la lecture des champs aux mailles "); ret = -1; continue;}
if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_DESCENDING_EDGE,USER_INTERLACE,_ncstp);
else {MESSAGE("Erreur a la lecture des champs aux faces "); ret = -1; continue;}
if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_NODE_ELEMENT,USER_INTERLACE,_ncstp);
else {MESSAGE("Erreur a la lecture des champs aux aretes"); ret = -1; continue;}
/*TODO */
/* if (lret == 0) lret = getFieldsOn(fid, nomcha, typcha, ncomp, MED_STRUCT_ELEMENT,USER_INTERLACE,_ncstp); */
/* else {MESSAGE("Erreur a la lecture des champs aux éléments de structure"); ret = -1; continue;} */
if (lret != 0) {MESSAGE("Erreur a la lecture des champs aux noeuds des mailles "); ret = -1;};
}
/* Interrogation des profils */
npro = MEDnProfile(fid);
printf("\nNombre de profils stockes : "IFORMAT"\n\n",npro);
for (i=1 ; i <= npro ; i++ ) {
if ( MEDprofileInfo(fid, i, pflname, &nval) < 0) {
MESSAGE("Erreur a la demande d'information sur le profil n° : "); ISCRUTE_int(i);
ret = -1;continue;
}
printf("\t- Profil n°%i de nom |%s| et de taille "IFORMAT"\n",i,pflname,nval);
pflval = (med_int*) malloc(sizeof(med_int)*nval);
if ( MEDprofileRd(fid, pflname, pflval) < 0) {
MESSAGE("Erreur a la lecture des valeurs du profil : ");
SSCRUTE(pflname);
ret = -1;
} else {
printf("\t");
for (j=0;j<nval;j++) printf(" "IFORMAT" ",*(pflval+j));
printf("\n\n");
}
free(pflval);
}
/* Interrogation des liens */
nln = MEDnLink(fid);
printf("\nNombre de liens stockes : "IFORMAT"\n\n",nln);
for (i=1 ; i <= nln ; i++ ) {
if ( MEDlinkInfo(fid, i, nomlien, &nval) < 0) {
MESSAGE("Erreur a la demande d'information sur le lien n° : "); ISCRUTE_int(i);
ret = -1;continue;
}
printf("\t- Lien n°%i de nom |%s| et de taille "IFORMAT"\n",i,nomlien,nval);
lien = malloc((nval+1)*sizeof(char));
EXIT_IF(lien == NULL,NULL,NULL);
if ( MEDlinkRd(fid, nomlien, lien ) < 0 ) {
MESSAGE("Erreur a la lecture du lien : ");
SSCRUTE(nomlien);SSCRUTE(lien);
ret = -1;
} else {
lien[nval] = '\0';
printf("\t\t|%s|\n\n",lien);
}
free(lien);
}
/* Interrogation des localisations des points de GAUSS */
nloc = MEDnLocalization(fid);
printf("\nNombre de localisations stockees : "IFORMAT"\n\n",nloc);
for (i=1 ; i <= nloc ; i++ ) {
if ( MEDlocalizationInfo(fid, i, locname, &type_geo, &locsdim,&ngauss,
geointerpname, ipointstructmeshname,
&nsectionmeshcell,§iongeotype) < 0) {
MESSAGE("Erreur a la demande d'information sur la localisation n° : "); ISCRUTE_int(i);
ret = -1;continue;
}
printf("\t- Loc. n°%i de nom |%s| de dimension %i avec "IFORMAT" pts de GAUSS \n",i,locname,locsdim,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 ( MEDlocalizationRd(fid, locname, USER_INTERLACE, refcoo, gscoo, wg ) < 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 ( MEDfileClose(fid) < 0) return -1;
return ret;
}
med_err getFieldsOn(med_idt fid, char * nomcha, med_field_type typcha, med_int ncomp,
med_entity_type entite, med_switch_mode stockage, med_int ncstp) {
int j,k,l,m,n,nb_geo=0;
med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,ngroup,*vale=NULL,nval;
med_int numdt=0,numo=0,_nprofile;
med_float *valr=NULL,dt=0.0;
med_err ret=0;
char pflname [MED_NAME_SIZE+1]="";
char locname [MED_NAME_SIZE+1]="";
char * lien = NULL;
char dt_unit [MED_SNAME_SIZE+1]="unknown";
med_geometry_type * type_geo;
const char * const * AFF;
const char * const * AFF_ENT=MED_GET_ENTITY_TYPENAME+1;
switch (entite) {
case MED_NODE :
type_geo = MED_GET_NODE_GEOMETRY_TYPE;
nb_geo = MED_N_NODE_FIXED_GEO;
AFF = MED_GET_NODE_GEOMETRY_TYPENAME;
break;
case MED_CELL :
case MED_NODE_ELEMENT :
type_geo = MED_GET_CELL_GEOMETRY_TYPE;
nb_geo = MED_N_CELL_FIXED_GEO;
AFF = MED_GET_CELL_GEOMETRY_TYPENAME;
break;
case MED_DESCENDING_FACE :
type_geo = MED_GET_FACE_GEOMETRY_TYPE;
nb_geo = MED_N_FACE_FIXED_GEO;
AFF = MED_GET_FACE_GEOMETRY_TYPENAME;
break;
case MED_DESCENDING_EDGE :
type_geo = MED_GET_EDGE_GEOMETRY_TYPE;
nb_geo = MED_N_EDGE_FIXED_GEO;
AFF = MED_GET_EDGE_GEOMETRY_TYPENAME;
break;
}
for (k=1;k<=nb_geo;k++) {
/* Combien de (PDT,NOR) a lire */
nbpdtnor = ncstp;
if (nbpdtnor < 1 ) continue;
for (j=0;j<nbpdtnor;j++) {
if ( MEDfieldComputingStepInfo(fid,nomcha,j+1, &numdt, &numo, &dt ) <0) {
MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
ISCRUTE(numdt); ISCRUTE(numo);
ret = -1; continue;
}
if ( (_nprofile = MEDfieldnProfile(fid,nomcha,numdt,numo,entite,type_geo[k],
pflname,locname ) ) < 0 ) {
MESSAGE("Erreur a la demande du nombre de profils referencés par le champ : ");
SSCRUTE(nomcha);
ISCRUTE(numdt); ISCRUTE(numo);
ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
SSCRUTE(AFF_ENT[(int)entite]);SSCRUTE(AFF[k]);
ret = -1; continue;
};
for (l=0;l<_nprofile;l++) {
if ( (nval = MEDfieldnValueWithProfile(fid, nomcha, numdt, numo, entite, type_geo[k],
l+1, USER_MODE, pflname,&pflsize,
locname, &ngauss) ) < 0 ) {
MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE(numdt);ISCRUTE(numo);
ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE_int(USER_MODE);
ret = -1; continue;
};
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);
printf("\t- Il y a "IFORMAT" entités qui portent des valeurs en mode %i. Chaque entite %s\
de type geometrique %s associes au profile |%s| a "IFORMAT" valeurs associées \n",
nval,USER_MODE,AFF_ENT[(int)entite],AFF[k],pflname,ngauss);
/*Lecture des valeurs du champ */
if (typcha == MED_FLOAT64) {
valr = (med_float*) calloc(ncomp*nval*ngauss,sizeof(med_float));
EXIT_IF(valr == NULL,NULL,NULL);
if (MEDfieldValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],
USER_MODE, pflname, stockage,MED_ALL_CONSTITUENT,
(unsigned char*) valr) < 0 ) {
MESSAGE("Erreur a la lecture des valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);
ret = -1;
}
} else {
vale = (med_int*) calloc(ncomp*nval*ngauss,sizeof(med_int));
EXIT_IF(vale == NULL,NULL,NULL);
if (MEDfieldValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],
USER_MODE, pflname, stockage,MED_ALL_CONSTITUENT,
(unsigned char*) vale) < 0 ) {
MESSAGE("Erreur a la lecture des valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);
ret = -1;
};
}
if ( strlen(locname) )
printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname);
if (entite == MED_NODE_ELEMENT)
ngroup = (type_geo[k] % 100);
else
ngroup = ngauss;
switch (stockage) {
case MED_FULL_INTERLACE :
printf("\t- Valeurs :\n\t");
for (m=0;m<(nval*ngauss)/ngroup;m++) {
printf("|");
for (n=0;n<ngroup*ncomp;n++)
if (typcha == MED_FLOAT64)
printf(" %f ",*(valr+(m*ngroup*ncomp)+n));
else
printf(" "IFORMAT" ",*(vale+(m*ngroup*ncomp)+n));
}
break;
/*Affichage en fonction du profil à traiter*/
case MED_NO_INTERLACE :
printf("\t- Valeurs :\n\t");
for (m=0;m<ncomp;m++) {
printf("|");
for (n=0;n<(nval*ngauss);n++)
if (typcha == MED_FLOAT64)
printf(" %f ",*(valr+(m*nval)+n));
else
printf(" "IFORMAT" ",*(vale+(m*nval)+n));
}
break;
}
printf("|\n");
if (typcha == MED_FLOAT64) {
if ( valr ) {free(valr);valr = NULL;}}
else
if (vale) { free(vale);vale = NULL; }
/*Lecture du profil associe */
if (strcmp(pflname,MED_NO_PROFILE) == 0 )
printf("\t- Profil : MED_NO_PROFILE\n");
else {
if ( (pflsize = MEDprofileSizeByName(fid,pflname)) <0 ) {
MESSAGE("Erreur a la lecture du nombre de valeurs du profil : ");
SSCRUTE(pflname);
ret = -1; continue;
}
printf("\t- Profil : |%s| de taille "IFORMAT"\n",pflname,pflsize);
pflval = (med_int*) malloc(sizeof(med_int)*pflsize);
EXIT_IF(pflval == NULL,NULL,NULL);
if ( MEDprofileRd(fid,pflname,pflval) <0) {
MESSAGE("Erreur a la lecture des valeurs du profil : ");
SSCRUTE(pflname);
ret = -1;
}
printf("\t");
for (m=0;m<pflsize;m++) printf(" "IFORMAT" ",*(pflval+m));
printf("\n");
free(pflval);
}
}
}
} /* fin for sur les mailles*/
return ret;
}
med_err MAJ_236_300_fieldOnEntity(med_idt fid, const char * const nomcha, const char * const meshname,
med_field_type typcha, med_int ncomp, med_entity_type entite, med_int ncstp,
char * const _pathi, char * const _pathf) {
med_err ret=-1;
int i,j,k,l,m,n,nb_geo=0;
med_int nbpdtnor=0,pflsize,*pflval,ngauss=0,ngroup,*vale=NULL,nval;
med_int numdt=0,numo=0,_nprofile;
med_int meshnumdt=0,meshnumit=0;
med_float *valr=NULL,dt=0.0;
unsigned char * _val=NULL;
char pflname [MED_NAME_SIZE+1]="";
char locname [MED_NAME_SIZE+1]="";
char _meshname [MED_NAME_SIZE+1]="";
char _fieldname [MED_NAME_SIZE+1]="";
char _pathtmp[MED_FIELD_GRP_SIZE+3]="/CHA__/";
char _pathfb[MED_FIELD_GRP_SIZE+2+MED_NAME_SIZE+1]="/CHA_/";
char * lien = NULL;
char dt_unit [MED_SNAME_SIZE+1]="unknown";
med_bool localmesh;
med_int nmesh=0;
hid_t _ocp_plist_id ;
hid_t _lcp_plist_id ;
int _isavlen=0;
int _fsavlen=0;
int _fieldnamelen=0;
htri_t _groupexist;
char _tmpmeshname[MED_NAME_SIZE+1]="";
med_bool _tmplocal=MED_FALSE;
med_field_type _tmptypcha;
char *_tmpcomp=NULL,*_tmpunit=NULL,_tmpdtunit[MED_SNAME_SIZE+1]="";
med_int _tmpncstp=0;
med_geometry_type * type_geo;
const char * const * AFF;
const char * const * AFF_ENT=MED_GET_ENTITY_TYPENAME+1;
switch (entite) {
case MED_NODE :
type_geo = MED_GET_NODE_GEOMETRY_TYPE;
nb_geo = MED_N_NODE_FIXED_GEO;
AFF = MED_GET_NODE_GEOMETRY_TYPENAME;
break;
case MED_CELL :
case MED_NODE_ELEMENT :
type_geo = MED_GET_CELL_GEOMETRY_TYPE;
nb_geo = MED_N_CELL_FIXED_GEO;
AFF = MED_GET_CELL_GEOMETRY_TYPENAME;
break;
case MED_DESCENDING_FACE :
type_geo = MED_GET_FACE_GEOMETRY_TYPE;
nb_geo = MED_N_FACE_FIXED_GEO;
AFF = MED_GET_FACE_GEOMETRY_TYPENAME;
break;
case MED_DESCENDING_EDGE :
type_geo = MED_GET_EDGE_GEOMETRY_TYPE;
nb_geo = MED_N_EDGE_FIXED_GEO;
AFF = MED_GET_EDGE_GEOMETRY_TYPENAME;
break;
}
strcpy(_fieldname,nomcha);
_isavlen=strlen(_pathi);
_fsavlen=strlen(_pathf);
_fieldnamelen=strlen(_fieldname);
_ocp_plist_id = H5Pcreate( H5P_OBJECT_COPY );
_lcp_plist_id = H5Pcreate( H5P_LINK_CREATE );
/* Reconstruit les objets pointés par un lien symbolique */
H5Pset_copy_object( _ocp_plist_id, H5O_COPY_EXPAND_SOFT_LINK_FLAG );
/*Ne crée pas les liens intermédiaires */
H5Pset_create_intermediate_group(_lcp_plist_id, -1);
/* Ne recopie pas en profondeur*/
H5Pset_copy_object( _ocp_plist_id, H5O_COPY_SHALLOW_HIERARCHY_FLAG);
/* ISCRUTE(nbpdtnor); */
/* SSCRUTE(_fieldname); */
for (k=1;k<=nb_geo;k++) {
/* Combien de (PDT,NOR) a lire */
nbpdtnor = ncstp;
if (nbpdtnor < 1 ) continue;
for (j=0;j<nbpdtnor;j++) {
if ( MEDfield23ComputingStepMeshInfo(fid,nomcha,j+1, &numdt, &numo, &dt,
&nmesh, _meshname,&localmesh, &meshnumdt, &meshnumit ) <0) {
MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
ISCRUTE(numdt); ISCRUTE(numo);ISCRUTE(nmesh);SSCRUTE(_meshname);ISCRUTE_int(localmesh);
ISCRUTE(meshnumdt);ISCRUTE(meshnumit);
ret = -1; continue;
}
/* SSCRUTE(_meshname);SSCRUTE(AFF_ENT[(int)entite]);SSCRUTE(AFF[k]);ISCRUTE(nmesh); */
for (i=0;i< nmesh;++i) {
if ( (_nprofile = MEDfield23nProfile(fid,nomcha,numdt,numo,entite,type_geo[k],i+1,_meshname,
pflname,locname ) ) < 0 ) {
MESSAGE("Erreur a la demande du nombre de profils referencés par le champ : ");
SSCRUTE(nomcha); ISCRUTE(numdt); ISCRUTE(numo);SSCRUTE(_meshname);
ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);SSCRUTE(pflname);SSCRUTE(locname);
SSCRUTE(AFF_ENT[(int)entite]);SSCRUTE(AFF[k]);
ret = -1; continue;
};
/* ISCRUTE(_nprofile); */
/* SSCRUTE(_meshname); */
/* SSCRUTE(meshname); */
/* SSCRUTE(pflname); */
for (l=0;l<_nprofile;l++) {
/* Si le maillage traité dans ce profil n'est pas celui par défaut :
- Vérifie que le nom du champ pour ce maillage secondaire existe.
- Initialise _fieldname au nom du champ qu'il faut traiter par la suite
*/
if (strcmp(_meshname,meshname)) {
/* Lecture du type du champ, des noms des composantes et du nom de l'unité*/
_tmpcomp = (char*) malloc(ncomp*MED_SNAME_SIZE+1);
EXIT_IF(_tmpcomp == NULL,NULL,NULL);
_tmpunit = (char*) malloc(ncomp*MED_SNAME_SIZE+1);
EXIT_IF(_tmpunit == NULL,NULL,NULL);
/* Evite de passer les paramètres _tmpcomp,_tmpunit,_tmpdtunit,&_tmpncstp dans la fct MAJ_236...*/
ret = MEDfieldInfoByName(fid,nomcha,
_tmpmeshname,&_tmplocal,&_tmptypcha,_tmpcomp,_tmpunit,_tmpdtunit,&_tmpncstp);
MED_ERR_EXIT_IF(ret,MED_ERR_ACCESS,MED_ERR_FIELD,nomcha);
/*Ne pose pas de problème de taille de chaîne car MED_NAME_SIZE a doublé en 3.0 */
_fieldname[_fieldnamelen]='_';strcpy(&_fieldname[_fieldnamelen+1],_meshname);
MAJ_version_num(fid,3,0,8);
/*Sauvegarde du champ initial dans _pathtmp*/
EXIT_IF( H5Gmove(fid, _pathi, _pathtmp ) < 0,"Switch to ",_pathtmp);
EXIT_IF( H5Gmove(fid, _pathf, _pathi ) < 0 ,"Switch to ",_pathi);
MED_ERR_EXIT_IF( MEDfieldCr(fid,_fieldname,
_tmptypcha,ncomp,_tmpcomp,_tmpunit,_tmpdtunit,_meshname ) < 0,
MED_ERR_CREATE,MED_ERR_FIELD,_fieldname);
/*Rétablissement des chemins d'accès au champ initial et au champ cible temporaire*/
EXIT_IF( H5Gmove(fid, _pathi , _pathf ) < 0,"Switch to ",_pathf);
EXIT_IF( H5Gmove(fid, _pathtmp, _pathi ) < 0,"Switch to ",_pathi);
MAJ_version_num(fid,2,3,6);
free(_tmpcomp);
free(_tmpunit);
} else {
strcpy(_fieldname,nomcha);
}
/* SSCRUTE(_fieldname); */
if ( (nval = MEDfield23nValueWithProfile(fid, nomcha, numdt, numo, entite, type_geo[k],_meshname,
l+1, MED_COMPACT_PFLMODE, pflname,&pflsize,
locname, &ngauss) ) < 0 ) {
MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(_meshname);
ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ret = -1; continue;
};
/* ISCRUTE(nval); */
/* 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); */
/* printf("\t- Il y a "IFORMAT" entités qui portent des valeurs en mode %i. Chaque entite %s\ */
/* de type geometrique %s associes au profile |%s| a "IFORMAT" valeurs associées \n", */
/* nval,MED_COMPACT_PFLMODE,AFF_ENT[(int)entite],AFF[k],pflname,ngauss); */
/* printf("\t- Le maillage associé est |%s|\n",_meshname); */
/*Lecture des valeurs du champ */
if (typcha == MED_FLOAT64) {
valr = (med_float*) calloc(ncomp*nval*ngauss,sizeof(med_float));
EXIT_IF(valr == NULL,NULL,NULL);
if (MEDfield23ValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],_meshname,
MED_COMPACT_PFLMODE, pflname, MED_NO_INTERLACE,MED_ALL_CONSTITUENT,
(unsigned char*) valr) < 0 ) {
MESSAGE("Erreur a la lecture des valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);
ret = -1;
}
_val = (unsigned char*) valr;
} else {
vale = (med_int*) calloc(ncomp*nval*ngauss,sizeof(med_int));
EXIT_IF(vale == NULL,NULL,NULL);
if (MEDfield23ValueWithProfileRd(fid, nomcha, numdt,numo, entite,type_geo[k],_meshname,
MED_COMPACT_PFLMODE, pflname, MED_NO_INTERLACE,MED_ALL_CONSTITUENT,
(unsigned char*) vale) < 0 ) {
MESSAGE("Erreur a la lecture des valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);
ret = -1;
}
_val = (unsigned char*) vale;
}
/* Ecriture du champ destination */
MAJ_version_num(fid,3,0,8);
/*Sauvegarde du champ initial dans _pathtmp*/
EXIT_IF( H5Gmove(fid, _pathi, _pathtmp ) < 0,"Switch to ",_pathtmp);
_groupexist=H5Lexists( fid, _pathf, H5P_DEFAULT );
EXIT_IF(!_groupexist,"Le champ devrait déjà existé",_pathf);
/*Déplacement du champ cible temporaire au chemin des champs écrits par MED */
if (_groupexist ) { EXIT_IF( (H5Gmove(fid, _pathf, _pathi ) < 0) ,"Switch to ",_pathi); }
/*Ecriture du champ cible au nouveau format*/
/* ISCRUTE(nval); */
MED_ERR_EXIT_IF( MEDfieldValueWithProfileWr(fid, _fieldname, numdt, numo, dt, entite,type_geo[k],
MED_COMPACT_PFLMODE, pflname, locname, MED_NO_INTERLACE, MED_ALL_CONSTITUENT,
nval, _val) < 0,
MED_ERR_WRITE,MED_ERR_FIELD,_fieldname);
free(_val);
/*Rétablissement des chemins d'accès au champ initial et au champ cible temporaire*/
EXIT_IF( H5Gmove(fid, _pathi , _pathf ) < 0,"Switch to ",_pathf);
EXIT_IF( H5Gmove(fid, _pathtmp, _pathi ) < 0,"Switch to ",_pathi);
MAJ_version_num(fid,2,3,6);
/* if ( strlen(locname) ) */
/* printf("\t- Modèle de localisation des points de Gauss de nom |%s|\n",locname); */
/* if (entite == MED_NODE_ELEMENT) */
/* ngroup = (type_geo[k] % 100); */
/* else */
/* ngroup = ngauss; */
/* switch (MED_NO_INTERLACE) { */
/* case MED_FULL_INTERLACE : */
/* printf("\t- Valeurs :\n\t"); */
/* for (m=0;m<(nval*ngauss)/ngroup;m++) { */
/* printf("|"); */
/* for (n=0;n<ngroup*ncomp;n++) */
/* if (typcha == MED_FLOAT64) */
/* printf(" %f ",*(valr+(m*ngroup*ncomp)+n)); */
/* else */
/* printf(" "IFORMAT" ",*(vale+(m*ngroup*ncomp)+n)); */
/* } */
/* break; */
/* case MED_NO_INTERLACE : */
/* printf("\t- Valeurs :\n\t"); */
/* for (m=0;m<ncomp;m++) { */
/* printf("|"); */
/* for (n=0;n<(nval*ngauss);n++) */
/* if (typcha == MED_FLOAT64) */
/* printf(" %f ",*(valr+(m*nval)+n)); */
/* else */
/* printf(" "IFORMAT" ",*(vale+(m*nval)+n)); */
/* } */
/* break; */
/* } */
/* printf("|\n"); */
/* if (typcha == MED_FLOAT64) { */
/* if ( valr ) {free(valr);valr = NULL;}} */
/* else */
/* if (vale) { free(vale);vale = NULL; } */
/* if (strcmp(pflname,MED_NO_PROFILE) == 0 ) */
/* printf("\t- Profil : MED_NO_PROFILE\n"); */
/* else { */
/* if ( (pflsize = MEDprofileSizeByName(fid,pflname)) <0 ) { */
/* MESSAGE("Erreur a la lecture du nombre de valeurs du profil : "); */
/* SSCRUTE(pflname); */
/* ret = -1; continue; */
/* } */
/* printf("\t- Profil : |%s| de taille "IFORMAT"\n",pflname,pflsize); */
/* pflval = (med_int*) malloc(sizeof(med_int)*pflsize); */
/* EXIT_IF(pflval == NULL,NULL,NULL); */
/* if ( MEDprofileRd(fid,pflname,pflval) <0) { */
/* MESSAGE("Erreur a la lecture des valeurs du profil : "); */
/* SSCRUTE(pflname); */
/* ret = -1; */
/* } */
/* printf("\t"); */
/* for (m=0;m<pflsize;m++) printf(" "IFORMAT" ",*(pflval+m)); */
/* printf("\n"); */
/* free(pflval); */
/* } */
}
}
}
} /* fin for sur les mailles*/
ret = 0;
ERROR:
return ret;
}
void MAJ_21_22_maillages(med_idt fid)
{
med_idt gid;
med_err ret;
int n,i;
char nom[MED_TAILLE_NOM+1];
char chemin[MED_TAILLE_MAA+MED_TAILLE_NOM+1];
char description[MED_TAILLE_DESC+1] = "Maillage converti au format MED V2.2";
med_int type = (med_int) MED_NON_STRUCTURE;
med_int dimension;
/* Lecture du nombre de maillages */
n = 0;
_MEDnObjets(fid,(char *) MED_MAA,&n);
EXIT_IF(n < 0,"Erreur a la lecture du nombre de maillage",NULL);
/*
* Mise a jour des maillages :
* - type : MED_NON_STRUCTURE
* - description : "Maillage converti au format V2.2"
*/
for (i=0;i<n;i++) {
/* on recupere le nom du maillage */
ret = _MEDobjetIdentifier(fid,(char *) MED_MAA,i,nom);
EXIT_IF(ret < 0,"Identification d'un maillage",NULL);
fprintf(stdout," >>> Normalisation du maillage [%s] \n",nom);
/* on accede au maillage */
strcpy(chemin,MED_MAA);
strcat(chemin,nom);
gid = _MEDdatagroupOuvrir(fid,chemin);
EXIT_IF(gid < 0,"Accès au maillage",nom);
/* lecture de la dimension du maillage */
ret = _MEDattrEntierLire(gid,(char *)(MED_NOM_DIM),&dimension);
EXIT_IF(ret < 0,"Lecture de la dimension du maillage",nom);
/* Ecriture du type et de la description */
ret = _MEDattrStringEcrire(gid,(char *)(MED_NOM_DES),MED_TAILLE_DESC,description);
EXIT_IF(ret < 0,"Ecriture de la description du maillage ",nom);
ret = _MEDattrEntierEcrire(gid,(char *)(MED_NOM_TYP),&type);
EXIT_IF(ret < 0,"Ecriture de la dimension du maillage ",nom);
/* Mise a jour des noeuds du maillage */
MAJ_21_22_noeuds_maillage(gid,dimension);
fprintf(stdout," ... Normalisation des noeuds effectuée ... \n");
/* Mise a jour des éléments du maillage */
MAJ_21_22_elements_maillage(gid,dimension);
fprintf(stdout," ... Normalisation des éléments effectuée ... \n");
/* Mise a jour des familles du maillage */
MAJ_21_22_familles_maillage(gid);
fprintf(stdout," ... Normalisation des familles effectuée ... \n");
/* On ferme tout */
ret = _MEDdatagroupFermer(gid);
EXIT_IF(ret < 0,"Fermeture de l'accès au maillage",NULL);
fprintf(stdout," >>> Normalisation du maillage [%s] ... OK ... \n",nom);
}
}
void MAJ_21_22_elements_maillage(med_idt mid, med_int dimension)
{
med_idt eid,gid,did,tid;
med_err ret;
int i,j;
med_geometrie_element typmai[MED_NBR_GEOMETRIE_MAILLE] = {MED_POINT1,MED_SEG2,
MED_SEG3,MED_TRIA3,
MED_TRIA6,MED_QUAD4,
MED_QUAD8,MED_TETRA4,
MED_TETRA10,MED_HEXA8,
MED_HEXA20,MED_PENTA6,
MED_PENTA15,MED_PYRA5,
MED_PYRA13};
int taille, edim;
char *nom, *nouvelle_chaine;
char nomgroup[MED_TAILLE_NOM_ENTITE+1];
med_int n;
med_size dimd[1];
med_int *old_conn,*conn;
/* On ne regarde que les mailles et la connectivité nodale */
eid = _MEDdatagroupOuvrir(mid,(char *)(MED_NOM_MAI));
EXIT_IF(eid < 0,"Ouverture du groupe HDF MED_NOM_MAI",NULL);
/* On normalise selon tous les types geometriques */
for (i=0;i<MED_NBR_GEOMETRIE_MAILLE;i++) {
/* On recupere le nom du groupe HDF */
_MEDnomGeometrie(nomgroup,typmai[i]);
/* On accède au type s'il existe dans le fichier */
gid = _MEDdatagroupOuvrir(eid,nomgroup);
if (gid < 0)
continue;
/* Nombre d'element ? */
did = _MEDdatasetOuvrir(gid,(char *)(MED_NOM_NOD));
EXIT_IF(did < 0,"Ouverture du dataset HDF MED_NOM_NOD",NULL);
ret = _MEDattrEntierLire(did,(char *)(MED_NOM_NBR),&n);
EXIT_IF(ret < 0,"Lecture du nombre d'elements",NULL);
ret = _MEDdatasetFermer(did);
EXIT_IF(ret < 0,"Fermeture du dataset HDF MED_NOM_NOD",NULL);
/* on normalise la connectivité si edim < dimension */
edim = typmai[i] / 100;
if (edim < dimension) {
taille = typmai[i]%100 + 1;
old_conn = (med_int *) malloc(sizeof(med_int)*taille*n);
EXIT_IF(old_conn == NULL,NULL,NULL);
#if defined(HAVE_F77INT64)
ret = _MED21datasetNumLire(gid,(char *)(MED_NOM_NOD),MED_INT64,
MED_NO_INTERLACE,(med_size)taille,MED_ALL,
0,NULL,MED_NOPG,
(unsigned char*) old_conn,H5T_NATIVE_INT);
#else
ret = _MED21datasetNumLire(gid,(char *)(MED_NOM_NOD),MED_INT32,
MED_NO_INTERLACE,(med_size) taille,MED_ALL,
0,NULL,MED_NOPG,
(unsigned char*) old_conn,H5T_NATIVE_INT);
#endif
/* On recopie dans le bon tableau */
taille --;
conn = (med_int *) malloc(sizeof(med_int)*taille*n);
EXIT_IF(conn == NULL,NULL,NULL);
for (j=0;j<n*taille;j++)
*(conn+j) = *(old_conn+j);
dimd[0] = n*taille;
#if defined(HAVE_F77INT64)
ret = _MED231datasetNumEcrire(gid,(char *) "TMP",MED_INT64,MED_NO_INTERLACE,
taille,MED_ALL,MED_NOPF,MED_NO_PFLMOD,0,MED_NOPG,dimd,
(unsigned char*) conn);
#else
ret = _MED231datasetNumEcrire(gid,(char *) "TMP",MED_INT32,MED_NO_INTERLACE,
taille,MED_ALL,MED_NOPF,MED_NO_PFLMOD,0,MED_NOPG,dimd,
(unsigned char*) conn);
#endif
EXIT_IF(ret < 0,"Ecriture de la nouvelle connectivité des mailles",NULL);
/* Ecriture du nombre de mailles dans le dataset HDF TMP */
tid = _MEDdatasetOuvrir(gid,"TMP");
EXIT_IF(tid < 0,"Ouverture du dataset HDF TMP",NULL);
ret = _MEDattrEntierEcrire(tid,(char *)(MED_NOM_NBR),&n);
EXIT_IF(ret < 0,"Ecriture du nombre de noeuds dans le dataset HDF TMP",NULL);
ret = _MEDdatasetFermer(tid);
EXIT_IF(ret < 0,"Fermeture du dataset HDF TMP",NULL);
/* Fermeture de l'accès aux connectivites */
ret = H5Gunlink(gid,(char *)(MED_NOM_NOD));
EXIT_IF(ret < 0,"Suppression des anciennes connectivités",NULL);
ret = H5Gmove(gid,"TMP",(char *)(MED_NOM_NOD));
EXIT_IF(ret < 0,"Mise en place des nouvelles connectivités",NULL);
/* on libere la memoire */
free(old_conn);
free(conn);
}
/* Mise a niveau des noms */
nom = (char *) malloc(n*ANCIEN_MED_TAILLE_PNOM+1);
EXIT_IF(nom == NULL,NULL,NULL);
nouvelle_chaine = (char *) malloc(n*MED_TAILLE_PNOM+1);
EXIT_IF(nouvelle_chaine == NULL,NULL,NULL);
ret = _MEDdatasetStringLire(gid,(char *)(MED_NOM_NOM),nom);
if (ret == 0) {
_MED23v30stringConvert(nouvelle_chaine, MED_TAILLE_PNOM,
nom, ANCIEN_MED_TAILLE_PNOM, n );
/* MAJ_21_22_chaine(nom,nouvelle_chaine,n); */
H5Gunlink(gid,(char *)(MED_NOM_NOM));
dimd[0] = n*MED_TAILLE_PNOM+1;
ret = _MEDdatasetStringEcrire(gid,(char *)(MED_NOM_NOM),dimd,nouvelle_chaine);
EXIT_IF(ret < 0,"Ecriture des nouveaux noms des éléments",NULL);
did = _MEDdatasetOuvrir(gid,(char *)(MED_NOM_NOM));
ret = _MEDattrEntierEcrire(did,(char *)(MED_NOM_NBR),&n);
ret = _MEDdatasetFermer(did);
}
free(nom);
free(nouvelle_chaine);
/* on ferme avant de passer au type geometrique suivant */
ret = _MEDdatagroupFermer(gid);
EXIT_IF(ret < 0,"Fermeture de l'accès aux mailles",NULL);
}
/* On ferme tout */
ret = _MEDdatagroupFermer(eid);
EXIT_IF(ret < 0,"Fermeture de l'accès aux mailles",NULL);
}
int main(int argc, char *argv[]) {
med_idt fid =0;
med_int majeur =0, mineur=0, release=0;
med_err ret =-1;
med_bool hdfok =MED_FALSE;
med_bool medok =MED_FALSE;
med_bool fileexist =MED_FALSE;
med_bool accessok =MED_FALSE;
if (argc != 2) {
fprintf(stdout,">> Utilisation : medconforme <nom_de_fichier_med> \n");
return 0;
}
/*
* Quelle version de la bibliotheque MED est utilisee ?
*/
ret=MEDlibraryNumVersion(&majeur, &mineur, &release);
EXIT_IF( ret<0 , "Erreur d'appel de la routine MEDlibraryNumVersion.", NULL);
fprintf(stdout,"- Version de MED-fichier utilisée par medconforme : "IFORMAT"."IFORMAT"."IFORMAT" \n",majeur,mineur,release);
/*
* Le fichier à lire est-il accessible ?
*/
ret = MEDfileExist(argv[1],MED_ACC_RDONLY,&fileexist,&accessok );
MED_ERR_EXIT_IF(ret < 0 , MED_ERR_CALL,MED_ERR_API,"MEDfileExist");
if ( !fileexist ) { fprintf(stdout,"- Le fichier [%s] n'existe pas \n",argv[1]); goto SORTIE; }
if ( !accessok ) { fprintf(stdout,"- Le fichier [%s] n'est pas accessible en lecture \n",argv[1]); goto SORTIE; }
/*
* Le fichier à lire est-il au bon format de fichier HDF ?
*/
ret=MEDfileCompatibility(argv[1],&hdfok,&medok);
MED_ERR_EXIT_IF(ret < 0 , MED_ERR_CALL,MED_ERR_API,"MEDfileCompatibility");
if ( hdfok ) fprintf(stdout,"- Format HDF du fichier MED [%s] conforme au format HDF utilisé par la bibliothèque \n",argv[1]);
else { fprintf(stdout,"- Format HDF du fichier MED [%s] non conforme au format HDF utilisé par la bibliothèque \n",argv[1]); goto SORTIE; }
/*
* Le fichier à lire a-t-il été créé avec une version de la bilbiothèque MED conforme avec celle utilisée ?
* (Numéros majeur identique et mineur de la bibliothèque supérieur à celui du fichier).
*/
if ( medok ) {
fprintf(stdout,"- Version MED du fichier [%s] conforme a la bibliothèque MED utilisée \n",argv[1]);
if ((fid = MEDfileOpen(argv[1],MED_ACC_RDONLY)) < 0) {
MED_ERR_(ret,MED_ERR_OPEN,MED_ERR_FILE,argv[1]);
goto ERROR;
}
/*
* Une fois le fichier ouvert on peut avoir acces au numero de version complet
*/
if ( MEDfileNumVersionRd(fid, &majeur, &mineur, &release) < 0 ) {
MED_ERR_(ret,MED_ERR_CALL,MED_ERR_API,"MEDfileNumVersionRd");
goto ERROR;
}
fprintf(stdout,"- Ce fichier a ete créé avec MED-fichier V"IFORMAT"."IFORMAT"."IFORMAT" \n",majeur,mineur,release);
}
else
fprintf(stdout,"- Version MED du fichier [%s] non conforme avec celle de la bibliothèque utilisée \n",argv[1]);
SORTIE:
ret = 0;
ERROR:
if (fid > 0)
if (MEDfileClose(fid) < 0) {
MED_ERR_(ret,MED_ERR_CLOSE,MED_ERR_FILE,argv[1]);
ret = -1;
}
return ret;
}
med_err
MED231champInfoEtRen(med_idt fid,int indice,char *champ,
med_type_champ *type,char *comp,char *unit,
med_int ncomp)
{
med_err ret=0;
med_idt gid;
char chemin [MED_TAILLE_CHA+MED_TAILLE_NOM+1];
char chemini[MED_TAILLE_CHA+MED_TAILLE_NOM+1];
char champf [MED_TAILLE_NOM+1];
int num;
med_int typechamp;
/*
* On inhibe le gestionnaire d'erreur HDF 5
*/
_MEDmodeErreurVerrouiller();
/*
* On recupere le nom du champ
*/
num = indice - 1;
strcpy(chemin,MED_CHA);
if ((ret = _MEDobjetIdentifier(fid,chemin,num,champ)) < 0)
return -1;
/*
* Si le Data Group cha n'existe pas => erreur
*/
if ( MAJ_231_232_chaine(champ,champf) ) {
fprintf(stdout," >>> Normalisation du nom de champ [%s] \n",champ);
/* on accede au maillage */
strcpy(chemini,chemin);
strcat(chemini,champ);
strcat(chemin,champf);
ret = H5Gmove(fid, chemini, chemin );
EXIT_IF(ret < 0,"Renommage du champ en",champf);
strcpy(champ,champf);
fprintf(stdout," >>> Normalisation du nom du champ [%s] ... OK ... \n",champf);
} else {
strcat(chemin,champ);
}
if ((gid = _MEDdatagroupOuvrir(fid,chemin)) < 0)
return -1;
/*
* La liste des attributs
*/
if ((ret = _MEDattrEntierLire(gid,MED_NOM_TYP,&typechamp)) < 0)
return -1;
*type = (med_type_champ) (typechamp);
if ((ret = _MEDattrStringLire(gid,MED_NOM_NOM,ncomp*MED_TAILLE_PNOM,
comp)) < 0)
return -1;
if ((ret = _MEDattrStringLire(gid,MED_NOM_UNI,ncomp*MED_TAILLE_PNOM,
unit)) < 0)
return -1;
/*
* On ferme tout
*/
if ((ret = _MEDdatagroupFermer(gid)) < 0)
return -1;
return 0;
}
static void
config_file_walk(config_fn_walk_cb_t fn_cb,
config_fn_walk_junk_cb_t fn_junk_cb, void *data)
{
FILE *data_file;
char buf[BUFSIZ], e_conf[BUFSIZ];
char *key, *value;
data_file = fopen(path_conf, "r");
EXIT_IF(data_file == NULL, _("failed to open configuration file"));
pthread_mutex_lock(&nbar.mutex);
for (;;) {
if (fgets(buf, sizeof buf, data_file) == NULL)
break;
io_extract_data(e_conf, buf, sizeof buf);
if (*e_conf == '\0' || *e_conf == '#') {
if (fn_junk_cb)
fn_junk_cb(buf, data);
continue;
}
key = e_conf;
value = strchr(e_conf, '=');
if (value) {
*value = '\0';
value++;
} else {
EXIT(_("invalid configuration directive: \"%s\""),
e_conf);
}
if (strcmp(key, "auto_save") == 0 ||
strcmp(key, "auto_gc") == 0 ||
strcmp(key, "periodic_save") == 0 ||
strcmp(key, "confirm_quit") == 0 ||
strcmp(key, "confirm_delete") == 0 ||
strcmp(key, "skip_system_dialogs") == 0 ||
strcmp(key, "skip_progress_bar") == 0 ||
strcmp(key, "calendar_default_view") == 0 ||
strcmp(key, "week_begins_on_monday") == 0 ||
strcmp(key, "color-theme") == 0 ||
strcmp(key, "layout") == 0 ||
strcmp(key, "side-bar_width") == 0 ||
strcmp(key, "notify-bar_show") == 0 ||
strcmp(key, "notify-bar_date") == 0 ||
strcmp(key, "notify-bar_clock") == 0 ||
strcmp(key, "notify-bar_warning") == 0 ||
strcmp(key, "notify-bar_command") == 0 ||
strcmp(key, "notify-all") == 0 ||
strcmp(key, "output_datefmt") == 0 ||
strcmp(key, "input_datefmt") == 0 ||
strcmp(key, "notify-daemon_enable") == 0 ||
strcmp(key, "notify-daemon_log") == 0) {
WARN_MSG(_("Pre-3.0.0 configuration file format detected, "
"please upgrade running `calcurse-upgrade`."));
}
if (value && (*value == '\0' || *value == '\n')) {
/* Backward compatibility mode. */
if (fgets(buf, sizeof buf, data_file) == NULL)
break;
io_extract_data(e_conf, buf, sizeof buf);
if (*e_conf == '#')
*e_conf = '\0';
value = e_conf;
}
fn_cb(key, value, data);
}
file_close(data_file, __FILE_POS__);
pthread_mutex_unlock(&nbar.mutex);
}
med_err MED231champNormaliser(med_idt fid, char * nomcha, med_type_champ typcha, med_int ncomp,
med_entite_maillage entite) {
int j,k,l,m,n,nb_geo;
med_int nbpdtnor=0,ngauss=0,*vale=NULL,nval;
med_int numdt=0,numo=0,nbrefmaa;
med_float *valr=NULL,dt=0.0;
med_err ret=0, rett;
med_booleen local;
char pflname [MED_TAILLE_NOM+1]="";
char locname [MED_TAILLE_NOM+1]="";
char maa_ass_i [MED_TAILLE_NOM+1]="";
char maa_ass [MED_TAILLE_NOM+1]="";
char dt_unit [MED_TAILLE_PNOM+1]="";
med_geometrie_element * type_geo;
med_geometrie_element typ_noeud[1] = { MED_NONE };
med_geometrie_element typmai[MED_NBR_GEOMETRIE_MAILLE+2] = {MED_POINT1, MED_SEG2, MED_SEG3, MED_TRIA3,
MED_QUAD4, MED_TRIA6,MED_QUAD8, MED_TETRA4,
MED_PYRA5, MED_PENTA6, MED_HEXA8, MED_TETRA10,
MED_PYRA13, MED_PENTA15, MED_HEXA20,
MED_POLYGONE, MED_POLYEDRE};
med_geometrie_element typfac[MED_NBR_GEOMETRIE_FACE+1] = {MED_TRIA3,MED_TRIA6,
MED_QUAD4,MED_QUAD8,
MED_POLYGONE};
med_geometrie_element typare[MED_NBR_GEOMETRIE_ARETE] = {MED_SEG2,MED_SEG3};
char ** AFF;
switch (entite) {
case MED_NOEUD :
type_geo = typ_noeud;
nb_geo = 1;
AFF = MED_GEOMETRIE_NOEUD_AFF;
break;
case MED_MAILLE :
type_geo = typmai;
nb_geo = MED_NBR_GEOMETRIE_MAILLE+2;
AFF = MED_GEOMETRIE_MAILLE_AFF;
break;
case MED_FACE :
type_geo = typfac;
nb_geo = MED_NBR_GEOMETRIE_FACE+1;
AFF = MED_GEOMETRIE_FACE_AFF;
break;
case MED_ARETE :
type_geo = typare;
nb_geo = MED_NBR_GEOMETRIE_ARETE;
AFF = MED_GEOMETRIE_ARETE_AFF;
break;
}
for (k=0;k<nb_geo;k++) {
/* Combien de (PDT,NOR) a lire */
nbpdtnor = MEDnPasdetemps(fid,nomcha,entite,type_geo[k]);
if (nbpdtnor < 1 ) continue;
for (j=0;j<nbpdtnor;j++) {
/*le nom du maillage associé va être modifié dans MED231champRefInfoEtRenMaa */
rett = MEDpasdetempsInfo(fid,nomcha,entite,type_geo[k],
j+1, &ngauss, &numdt, &numo, dt_unit,
&dt, maa_ass_i, &local, &nbrefmaa);
if ( rett <0) {
MESSAGE("Erreur a la demande d'information sur (pdt,nor) : ");
ISCRUTE(numdt); ISCRUTE(numo);
ret = -1; continue;
};
/* Combien de maillages lies aux (nomcha,ent,geo,numdt,numo) */
/* Notons que cette information est egalement disponible a partir de MEDpasdetempsInfo */
nbrefmaa = MEDnChampRef(fid,nomcha,entite,type_geo[k],numdt,numo);
if ( nbrefmaa < 0 ) {
MESSAGE("Erreur a la demande du nombre de maillages references par le champ : ");
SSCRUTE(nomcha); ISCRUTE(numdt); ISCRUTE(numo);
ret = -1; continue;
};
for (l=0;l<nbrefmaa;l++) {
rett = MED231champRefInfoEtRenMaa(fid,nomcha,entite,type_geo[k],
l+1,numdt, numo, maa_ass, &local, &ngauss);
if ( rett < 0 ) {
MESSAGE("Erreur a la demande d'information sur le maillage utilise par le champ n° : ");
ISCRUTE_int(l+1);
ret = -1; continue;
};
/* Prend en compte le nbre de pt de gauss automatiquement */
nval = MEDnVal(fid,nomcha,entite,type_geo[k],numdt,numo,maa_ass,MED_GLOBAL);
if (nval <= 0) {
MESSAGE("Erreur a la lecture du nombre de valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
ret = -1; continue;
};
/*Lecture des valeurs du champ */
if (typcha == MED_FLOAT64) {
valr = (med_float*) calloc(ncomp*nval,sizeof(med_float));
EXIT_IF(valr == NULL,NULL,NULL);
rett = MED231champLireEtUnlink(fid,maa_ass,nomcha,(unsigned char*)valr,MED_NO_INTERLACE,MED_ALL,locname,
pflname,MED_GLOBAL,entite,type_geo[k],numdt,numo);
if ( rett < 0 ) {
MESSAGE("Erreur a la lecture des valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
ret = -1;
};
}
else {
vale = (med_int*) calloc(ncomp*nval,sizeof(med_int));
EXIT_IF(vale == NULL,NULL,NULL);
rett = MED231champLireEtUnlink(fid,maa_ass,nomcha,(unsigned char*)vale,MED_NO_INTERLACE,MED_ALL,locname,
pflname,MED_GLOBAL,entite,type_geo[k],numdt,numo);
if ( rett < 0 ) {
MESSAGE("Erreur a la lecture des valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
ret = -1;
};
}
if (strcmp(pflname,MED_NOPFL)) {
if (typcha == MED_FLOAT64) {
rett = MEDchampEcr(fid,maa_ass,nomcha,(unsigned char*)valr,MED_NO_INTERLACE,nval,locname,
MED_ALL,pflname,MED_GLOBAL,entite,type_geo[k],numdt,dt_unit,dt,numo);
if ( rett < 0 ) {
MESSAGE("Erreur a l'ecriture du nombre de valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
ret = -1;
};
}
else {
rett = MEDchampEcr(fid,maa_ass,nomcha,(unsigned char*)vale,MED_NO_INTERLACE,nval,locname,
MED_ALL,pflname,MED_GLOBAL,entite,type_geo[k],numdt,dt_unit,dt,numo);
if ( rett < 0 ) {
MESSAGE("Erreur a l'ecriture des valeurs du champ : ");
SSCRUTE(nomcha);ISCRUTE_int(entite);ISCRUTE_int(type_geo[k]);
ISCRUTE(numdt);ISCRUTE(numo);SSCRUTE(maa_ass);
ret = -1;
};
}
}
if (typcha == MED_FLOAT64) {
if ( valr ) {free(valr);valr = NULL;}}
else
if (vale) { free(vale);vale = NULL; }
}
}
} /* fin for sur les mailles*/
return ret;
}
