MED fichier
test31.c
Aller à la documentation de ce fichier.
1 /* This file is part of MED.
2  *
3  * COPYRIGHT (C) 1999 - 2021 EDF R&D, CEA/DEN
4  * MED is free software: you can redistribute it and/or modify
5  * it under the terms of the GNU Lesser General Public License as published by
6  * the Free Software Foundation, either version 3 of the License, or
7  * (at your option) any later version.
8  *
9  * MED is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12  * GNU Lesser General Public License for more details.
13  *
14  * You should have received a copy of the GNU Lesser General Public License
15  * along with MED. If not, see <http://www.gnu.org/licenses/>.
16  */
17 
18 /******************************************************************************
19  * - Nom du fichier : test31.c
20  *
21  * - Description : ecriture d'une numerotation globale dans un maillage MED
22  *
23  *****************************************************************************/
24 
25 #include <med.h>
26 #define MESGERR 1
27 #include <med_utils.h>
28 
29 
30 #ifdef DEF_LECT_ECR
31 #define MODE_ACCES MED_ACC_RDWR
32 #elif DEF_LECT_AJOUT
33 #define MODE_ACCES MED_ACC_RDEXT
34 #else
35 #define MODE_ACCES MED_ACC_CREAT
36 #endif
37 
38 int main (int argc, char **argv)
39 
40 
41 {
42  med_err ret = 0;
43  med_idt fid;
44  med_int mdim,sdim;
45  /* nom du maillage de longueur maxi MED_NAME_SIZE */
46  char maa[MED_NAME_SIZE+1];
47  /* le nombre de noeuds */
48  med_int nnoe = 0;
49  /* table des numeros global */
50  med_int *numglobalnoe=NULL;
51 
52  /* variable de stockage pour reperer le maillage */
53  med_int i;
54  char des[MED_COMMENT_SIZE+1]="";
55  char dtunit[MED_SNAME_SIZE+1]="";
56  char nomcoo[3*MED_SNAME_SIZE+1]="";
57  char unicoo[3*MED_SNAME_SIZE+1]="";
58  med_axis_type rep;
59  med_mesh_type type;
60  med_sorting_type sort;
61  med_int nstep=0;
62  med_bool chgt=MED_FALSE,trsf=MED_FALSE;
63 
64 
65 
66  if (argc != 2) {
67  MESSAGE("Il faut passer un fichier MED en param�tre");
68  return -1;
69  }
70 
71  /* Ouverture du fichier passe en argument */
72  if ((fid = MEDfileOpen(argv[1],MED_ACC_RDWR)) < 0) {
73  MESSAGE("Erreur a l'ouverture du fichier : "); SSCRUTE(argv[1]);
74  return -1;
75  }
76 
77 
78  if ((sdim=MEDmeshnAxis(fid, 1)) <0) {
79  MESSAGE("Erreur � la lecture de la dimension de l'espace du maillage :");
80  SSCRUTE(maa);
81  return -1;
82  }
83 
84  /* Lecture des infos concernant le premier maillage */
85  if ( MEDmeshInfo( fid, 1, maa, &sdim, &mdim, &type, des, dtunit, &sort,
86  &nstep, &rep, nomcoo,unicoo) < 0 ) {
87  MESSAGE("Erreur a la lecture des informations sur le maillage : ");SSCRUTE(maa);
88  return -1;
89  } else {
90  printf("Maillage de nom : |%s| , de dimension : "IFORMAT" , et de type %d\n",maa,mdim,type);
91  printf("\t -Dimension de l'espace : "IFORMAT"\n",sdim);
92  printf("\t -Description du maillage : %s\n",des);
93  printf("\t -Noms des axes : %s\n",nomcoo);
94  printf("\t -Unit�s des axes : %s\n",unicoo);
95  printf("\t -Type de rep�re : %d\n",rep);
96  printf("\t -Nombre d'�tapes de calcul : "IFORMAT"\n",nstep);
97  printf("\t -Unit� des dates : %s\n\n",dtunit);
98  }
99 
100  /* Lecture du nombre de noeuds */
101  if ( (nnoe = MEDmeshnEntity(fid,maa,MED_NO_DT,MED_NO_IT,
103  &chgt,&trsf)) < 0) {
104  MESSAGE("Erreur a la lecture du nombre de noeuds ");
105  return -1;
106  }
107  printf("Nombre de noeuds : "IFORMAT" \n",nnoe);
108 
109  /* Allocations memoires */
110 
111  /* table de la numerotation globale
112  profil : (nombre de noeuds ) */
113  if (nnoe > 0) {
114  numglobalnoe = (med_int*) malloc(sizeof(med_int)*nnoe);
115  for (i=0;i<nnoe;i++) {
116  numglobalnoe[i]=i+100;
117  }
118  }
119 
120  /* ecriture de la numerotation globale */
121 
122  if (MEDmeshGlobalNumberWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_NODE,MED_NONE,nnoe,numglobalnoe)<0) {
123  MESSAGE("Erreur a l''ecriture de la numerotation globale");
124  return -1;
125  }
126 
127  free(numglobalnoe);
128 
129  /* Fermeture du fichier */
130  if (MEDfileClose(fid) < 0) {
131  MESSAGE("Erreur a la fermeture du fichier");
132  return -1;
133  }
134 
135  return 0;
136 }
137 
138 
139 
140 
MED_COMMENT_SIZE
#define MED_COMMENT_SIZE
Definition: med.h:79
MEDmeshnEntity
MEDC_EXPORT med_int MEDmeshnEntity(const med_idt fid, const char *const meshname, const med_int numdt, const med_int numit, const med_entity_type entitype, const med_geometry_type geotype, const med_data_type datatype, const med_connectivity_mode cmode, med_bool *const changement, med_bool *const transformation)
Cette routine permet de lire le nombre d'entités dans un maillage pour une étape de calcul donnée.
Definition: MEDmeshnEntity.c:44
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_err
herr_t med_err
Definition: med.h:332
med_sorting_type
med_sorting_type
Definition: med.h:309
MED_NO_CMODE
Definition: med.h:255
MED_FALSE
Definition: med.h:260
MESSAGE
#define MESSAGE(chaine)
Definition: med_utils.h:324
med_int
int med_int
Definition: med.h:342
MED_ACC_RDWR
Definition: med.h:121
med.h
med_bool
med_bool
Definition: med.h:260
IFORMAT
#define IFORMAT
Definition: med_utils.h:145
MED_NO_DT
#define MED_NO_DT
Definition: med.h:320
MED_NONE
#define MED_NONE
Definition: med.h:231
MED_NODE
Definition: med.h:143
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
MEDmeshGlobalNumberWr
MEDC_EXPORT med_err MEDmeshGlobalNumberWr(const med_idt fid, const char *const meshname, const med_int numdt, const med_int numit, const med_entity_type entitytype, const med_geometry_type geotype, const med_int nentity, const med_int *const number)
Cette routine permet l'écriture d'une numérotation globale sur un maillage pour un type d'entité,...
Definition: MEDmeshGlobalNumberWr.c:41
MED_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:81
MEDfileClose
MEDC_EXPORT med_err MEDfileClose(med_idt fid)
Fermeture d'un fichier MED.
Definition: MEDfileClose.c:30
main
int main(int argc, char **argv)
Definition: test31.c:37
med_axis_type
med_axis_type
Definition: med.h:258
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
MED_NO_IT
#define MED_NO_IT
Definition: med.h:321
MED_COORDINATE
Definition: med.h:149
med_utils.h