MED fichier

c/test3.c
/* This file is part of MED.
*
* COPYRIGHT (C) 1999 - 2021 EDF R&D, CEA/DEN
* MED is free software: you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* MED is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with MED. If not, see <http://www.gnu.org/licenses/>.
*/
/******************************************************************************
* - Nom du fichier : test3.c
*
* - Description : lecture des informations sur les maillages d'un fichier MED.
*
*****************************************************************************/
#include <med.h>
#define MESGERR 1
#include "med_utils.h"
#include <string.h>
#ifdef DEF_LECT_ECR
#define MODE_ACCES MED_ACC_RDWR
#elif DEF_LECT_AJOUT
#define MODE_ACCES MED_ACC_RDEXT
#else
#define MODE_ACCES MED_ACC_CREAT
#endif
#define OBJEXIST(oname)\
if ( MEDfileObjectExist(fid,MED_MESH, #oname ,& oname##exist) < 0) { \
MESSAGE("Erreur de test de présence du maillage "#oname); \
return -1; \
} \
if ( oname##exist) { MESSAGE("Le maillage "#oname" existe."); } else \
{ MESSAGE("Le maillage "#oname" n'existe pas.");}
int main (int argc, char **argv)
{
med_err ret = 0;
med_idt fid = 0;
med_int nmaa=0,mdim=0,sdim=0,nstep=0;
int i=0;
/*Pour les outils de type memchecker */
/* char maa[MED_NAME_SIZE+1]=""; */
/* char nomu[MED_LNAME_SIZE+1]=""; */
/* char des [MED_COMMENT_SIZE+1]=""; */
/* char dtunit[MED_SNAME_SIZE+1]=""; */
char *maa ;
char *nomu ;
char *des ;
char *dtunit;
char *axisname=0,*axisunit=0;
med_mesh_type meshtype;
med_err inomu;
med_sorting_type sortingtype;
med_axis_type axistype;
med_bool maa1exist=MED_FALSE,maa2exist=MED_FALSE;
med_bool maa3exist=MED_FALSE,maa4exist=MED_FALSE;
/*Pour les outils de type memchecker */
maa = (char *) malloc(sizeof(char)*(MED_NAME_SIZE+1 ));
nomu = (char *) malloc(sizeof(char)*(MED_LNAME_SIZE+1 ));
des = (char *) malloc(sizeof(char)*(MED_COMMENT_SIZE+1));
dtunit = (char *) malloc(sizeof(char)*(MED_SNAME_SIZE+1 ));
/* Ouverture du fichier "test2.med" en lecture seule */
fid = MEDfileOpen("test2.med",MED_ACC_RDONLY);
if (fid < 0) {
MESSAGE("Erreur a l'ouverture du fichier test2.med");
return -1;
}
/* Teste la présence des maillages */
OBJEXIST(maa1); OBJEXIST(maa2); OBJEXIST(maa3); OBJEXIST(maa4);
/* Lecture du nombre de maillage dans le fichier */
nmaa = MEDnMesh(fid);
if (nmaa < 0) {
MESSAGE("Erreur a la lecture du nombre de maillage");
return -1;
}
/* Boucle sur tous les maillages, pour chaque maillage, on lit :
- Le nom.
- Le type
- La dimension
- La description
- La dimension de l'espace si elle existe
- Le nom universel s'il existe
*/
printf("- Nombre de maillage dans test2.med = "IFORMAT"\n",nmaa);
for (i=0;i< nmaa;i++) {
/* lecture des informations */
if ((sdim=MEDmeshnAxis(fid, i+1)) <0) {
MESSAGE("Erreur a la lecture de la dimension de l'espace du maillage :");
ret = -1;
}
axisname = (char*) malloc(MED_SNAME_SIZE*sdim+1);
axisunit = (char*) malloc(MED_SNAME_SIZE*sdim+1);
if (MEDmeshInfo(fid,i+1, maa, &sdim, &mdim, &meshtype, des, dtunit, &sortingtype, &nstep,
&axistype, axisname, axisunit) < 0) {
MESSAGE("Erreur a la lecture des informations du maillage :"); SSCRUTE(maa);
ret = -1;
}
/* lecture du nom universel */
inomu = MEDmeshUniversalNameRd(fid,maa,nomu);
/* affichage des donnees lues */
if (inomu < 0)
printf("maillage %d de nom %s, de dimension "IFORMAT" \n",i+1,maa,mdim);
else
printf("maillage %d de nom %s, de dimension "IFORMAT" et de nom univ. %s\n",i+1,maa,mdim,nomu);
printf("La dimension de l'espace est "IFORMAT" \n",sdim);
if (meshtype == MED_STRUCTURED_MESH)
printf("Il s'agit d'un maillage structure \n");
else
printf("Il s'agit d'un maillage non structure \n");
printf("Description associee au maillage : %s \n\n",des);
printf("\t -Noms des axes : %s\n", axisname);
printf("\t -Unités des axes : %s\n",axisunit);
printf("\t -Type de repère : %d\n", axistype);
printf("\t -Nombre d'étape de calcul : "IFORMAT"\n",nstep);
printf("\t -Unité des dates : %s\n",dtunit);
free(axisname);
free(axisunit);
}
free( maa );
free( nomu );
free( des );
free( dtunit );
/* Fermeture du fichier */
if ( MEDfileClose(fid) < 0) {
MESSAGE("Erreur a la fermeture du fichier test2.med");
return -1;
}
return ret;
}
MED_COMMENT_SIZE
#define MED_COMMENT_SIZE
Definition: med.h:79
MEDfileOpen
MEDC_EXPORT med_idt MEDfileOpen(const char *const filename, const med_access_mode accessmode)
Ouverture d'un fichier MED.
Definition: MEDfileOpen.c:42
MED_SNAME_SIZE
#define MED_SNAME_SIZE
Definition: med.h:82
med_idt
hid_t med_idt
Definition: med.h:331
MED_ACC_RDONLY
Definition: med.h:120
med_err
herr_t med_err
Definition: med.h:332
med_sorting_type
med_sorting_type
Definition: med.h:309
MED_FALSE
Definition: med.h:260
MESSAGE
#define MESSAGE(chaine)
Definition: med_utils.h:324
MEDmeshUniversalNameRd
MEDC_EXPORT med_err MEDmeshUniversalNameRd(const med_idt fid, const char *const meshname, char *const univname)
Cette routine permet la lecture du nom universel d'un maillage.
Definition: MEDmeshUniversalNameRd.c:38
med_int
int med_int
Definition: med.h:342
med.h
med_bool
med_bool
Definition: med.h:260
IFORMAT
#define IFORMAT
Definition: med_utils.h:145
MEDmeshnAxis
MEDC_EXPORT med_int MEDmeshnAxis(const med_idt fid, const int meshit)
Cette routine permet de lire dans un maillage le nombre d'axes du repère des coordonnées des noeuds.
Definition: MEDmeshnAxis.c:35
med_mesh_type
med_mesh_type
Definition: med.h:131
SSCRUTE
#define SSCRUTE(chaine)
Definition: med_utils.h:323
MED_LNAME_SIZE
#define MED_LNAME_SIZE
Definition: med.h:83
MED_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:81
MEDnMesh
MEDC_EXPORT med_int MEDnMesh(const med_idt fid)
Cette routine permet de lire le nombre de maillages dans un fichier.
Definition: MEDnMesh.c:34
MEDfileClose
MEDC_EXPORT med_err MEDfileClose(med_idt fid)
Fermeture d'un fichier MED.
Definition: MEDfileClose.c:30
OBJEXIST
#define OBJEXIST(oname)
Definition: test3.c:39
med_axis_type
med_axis_type
Definition: med.h:258
MED_STRUCTURED_MESH
Definition: med.h:131
MEDmeshInfo
MEDC_EXPORT med_err MEDmeshInfo(const med_idt fid, const int meshit, char *const meshname, med_int *const spacedim, med_int *const meshdim, med_mesh_type *const meshtype, char *const description, char *const dtunit, med_sorting_type *const sortingtype, med_int *const nstep, med_axis_type *const axistype, char *const axisname, char *const axisunit)
Cette routine permet de lire les informations relatives à un maillage dans un fichier.
Definition: MEDmeshInfo.c:43
main
int main(int argc, char **argv)
Definition: 3.0.8/test10.c:50
med_utils.h