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);
}
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);
}
}
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;
}
/* Rename a group. */
int
NC4_rename_grp(int grpid, const char *name)
{
NC_GRP_INFO_T *grp, *g;
NC_HDF5_FILE_INFO_T *h5;
char norm_name[NC_MAX_NAME + 1];
int retval;
LOG((2, "nc_rename_grp: grpid 0x%x name %s", grpid, name));
/* Find info for this file and group, and set pointer to each. */
if ((retval = nc4_find_grp_h5(grpid, &grp, &h5)))
return retval;
if (!h5)
return NC_ENOTNC4;
if (h5->no_write)
return NC_EPERM; /* attempt to write to a read-only file */
/* Do not allow renaming the root group */
if(grp->parent == NULL)
return NC_EBADGRPID;
/* Check and normalize the name. */
if ((retval = nc4_check_name(name, norm_name)))
return retval;
/* Check that this name is not in use as a var, grp, or type. */
if ((retval = nc4_check_dup_name(grp, norm_name)))
return retval;
/* If it's not in define mode, switch to define mode. */
if (!(h5->flags & NC_INDEF))
if ((retval = NC4_redef(grpid)))
return retval;
/* Rename the group, if it exists in the file */
if (grp->hdf_grpid)
{
/* Close the group */
if (H5Gclose(grp->hdf_grpid) < 0)
return NC_EHDFERR;
grp->hdf_grpid = 0;
/* Attempt to rename & re-open the group, if the parent group is open */
if (grp->parent->hdf_grpid)
{
/* Rename the group */
if (H5Gmove(grp->parent->hdf_grpid, grp->name, name) < 0)
return NC_EHDFERR;
/* Reopen the group, with the new name */
if ((grp->hdf_grpid = H5Gopen2(grp->parent->hdf_grpid, name, H5P_DEFAULT)) < 0)
return NC_EHDFERR;
}
}
/* Give the group its new name in metadata. UTF8 normalization
* has been done. */
free(grp->name);
if (!(grp->name = malloc((strlen(norm_name) + 1) * sizeof(char))))
return NC_ENOMEM;
strcpy(grp->name, norm_name);
return NC_NOERR;
}
int
main()
{
printf("\n*** Checking HDF5 group functions.\n");
printf("*** Checking out root group...");
{
hid_t fileid, grpid, access_plistid;
/* Open the root group of a new file. */
if ((access_plistid = H5Pcreate(H5P_FILE_ACCESS)) < 0) ERR;
if (H5Pset_fclose_degree(access_plistid, H5F_CLOSE_SEMI)) ERR;
if ((fileid = H5Fcreate(FILE_NAME, H5F_ACC_TRUNC, H5P_DEFAULT,
access_plistid)) < 0) ERR;
if ((grpid = H5Gopen(fileid, "/")) < 0) ERR;
if (H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0) ERR;
/* Reopen file and root group. */
if ((fileid = H5Fopen(FILE_NAME, H5F_ACC_RDWR,
access_plistid)) < 0) ERR;
if ((grpid = H5Gopen(fileid, "/")) < 0) ERR;
if (H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0)
ERR;
}
SUMMARIZE_ERR;
printf("*** Checking out H5Gmove...");
{
hid_t fileid, grpid;
hid_t datasetid, spaceid;
/* Create file with one dataset. */
if ((fileid = H5Fcreate(FILE_NAME, H5F_ACC_TRUNC, H5P_DEFAULT,
H5P_DEFAULT)) < 0) ERR;
if ((grpid = H5Gopen(fileid, "/")) < 0) ERR;
if ((spaceid = H5Screate(H5S_SCALAR)) < 0) ERR;
if ((datasetid = H5Dcreate(grpid, DATASET_NAME, H5T_NATIVE_INT,
spaceid, H5P_DEFAULT)) < 0) ERR;
if (H5Dclose(datasetid) < 0 ||
H5Sclose(spaceid) < 0 ||
H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0) ERR;
/* Reopen file and check, then rename dataset. */
if ((fileid = H5Fopen(FILE_NAME, H5F_ACC_RDWR,
H5P_DEFAULT)) < 0) ERR;
if ((grpid = H5Gopen(fileid, "/")) < 0) ERR;
if ((datasetid = H5Dopen1(grpid, DATASET_NAME)) < 0) ERR;
if (H5Dclose(datasetid) < 0) ERR;
if (H5Gmove(grpid, DATASET_NAME, NEW_NAME) < 0) ERR;
if ((datasetid = H5Dopen1(grpid, NEW_NAME)) < 0) ERR;
if (H5Dclose(datasetid) < 0 ||
H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0)
ERR;
}
SUMMARIZE_ERR;
printf("*** Checking out sub-groups...");
{
hid_t fileid, grpid, subgrpid;
/* Create file with some nested groups. */
if ((fileid = H5Fcreate(FILE_NAME, H5F_ACC_TRUNC, H5P_DEFAULT,
H5P_DEFAULT)) < 0) ERR;
if ((grpid = H5Gcreate(fileid, GRP_NAME, 0)) < 0) ERR;
if ((subgrpid = H5Gcreate(grpid, SUB_GRP_NAME, 0)) < 0) ERR;
if (H5Gclose(subgrpid) < 0 ||
H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0) ERR;
/* Reopen file and discover groups. */
if ((fileid = H5Fopen(FILE_NAME, H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) ERR;
if ((grpid = H5Gopen(fileid, GRP_NAME)) < 0) ERR;
if ((subgrpid = H5Gopen(grpid, SUB_GRP_NAME)) < 0) ERR;
if (H5Gclose(subgrpid) < 0 ||
H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0) ERR;
}
SUMMARIZE_ERR;
printf("*** Checking out UTF8 named sub-group...");
{
hid_t fileid, grpid, subgrpid;
/* Create file with nested group. */
if ((fileid = H5Fcreate(FILE_NAME, H5F_ACC_TRUNC, H5P_DEFAULT,
H5P_DEFAULT)) < 0) ERR;
if ((grpid = H5Gcreate(fileid, (char *)norm_utf8, 0)) < 0) ERR;
if ((subgrpid = H5Gcreate(grpid, SUB_GRP_NAME, 0)) < 0) ERR;
if (H5Gclose(subgrpid) < 0 ||
H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0) ERR;
/* Reopen file and discover groups. */
if ((fileid = H5Fopen(FILE_NAME, H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) ERR;
if ((grpid = H5Gopen(fileid, (char *)norm_utf8)) < 0) ERR;
if ((subgrpid = H5Gopen(grpid, SUB_GRP_NAME)) < 0) ERR;
if (H5Gclose(subgrpid) < 0 ||
H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0) ERR;
}
SUMMARIZE_ERR;
printf("*** Checking out UTF8 named sub-group with group creation ordering...");
{
hid_t fileid, grpid, subgrpid;
hid_t fapl_id, fcpl_id, gcpl_id;
/* Create file with nested group. */
if ((fapl_id = H5Pcreate(H5P_FILE_ACCESS)) < 0) ERR;
if (H5Pset_libver_bounds(fapl_id, H5F_LIBVER_LATEST, H5F_LIBVER_LATEST) < 0) ERR;
if ((fcpl_id = H5Pcreate(H5P_FILE_CREATE)) < 0) ERR;
if (H5Pset_link_creation_order(fcpl_id, H5P_CRT_ORDER_TRACKED|H5P_CRT_ORDER_INDEXED) < 0) ERR;
if ((fileid = H5Fcreate(FILE_NAME, H5F_ACC_TRUNC, fcpl_id, fapl_id)) < 0) ERR;
if ((gcpl_id = H5Pcreate(H5P_GROUP_CREATE)) < 0) ERR;
if (H5Pset_link_creation_order(gcpl_id, H5P_CRT_ORDER_TRACKED|H5P_CRT_ORDER_INDEXED) < 0) ERR;
if ((grpid = H5Gcreate_anon(fileid, gcpl_id, H5P_DEFAULT)) < 0) ERR;
if ((H5Olink(grpid, fileid, (char *)norm_utf8, H5P_DEFAULT, H5P_DEFAULT)) < 0) ERR;
if ((subgrpid = H5Gcreate(grpid, SUB_GRP_NAME, 0)) < 0) ERR;
if (H5Gclose(subgrpid) < 0 ||
H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0) ERR;
/* Reopen file and discover groups. */
if ((fileid = H5Fopen(FILE_NAME, H5F_ACC_RDONLY, H5P_DEFAULT)) < 0) ERR;
if ((grpid = H5Gopen(fileid, (char *)norm_utf8)) < 0) ERR;
if ((subgrpid = H5Gopen(grpid, SUB_GRP_NAME)) < 0) ERR;
if (H5Gclose(subgrpid) < 0 ||
H5Gclose(grpid) < 0 ||
H5Fclose(fileid) < 0) ERR;
}
SUMMARIZE_ERR;
FINAL_RESULTS;
}
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;
}
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;
}
