MED fichier: Fichier source de test27.c

MED fichier

- Documentation MED
- Guides d'utilisation
- Guides de référence

test27.c
Aller à la documentation de ce fichier.
1 /* This file is part of MED.
2 *
3 * COPYRIGHT (C) 1999 - 2019 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 /******************************************************************************
20 * - Nom du fichier : test27.c
21 *
22 * - Description : creation de maillages structures (grille cartesienne |
23 * grille standard ) dans le fichier test27.med
24 *
25 *****************************************************************************/
26
27 #include <med.h>
28 #define MESGERR 1
29 #include "med_utils.h"
30 #include <string.h>
31
32 #ifdef DEF_LECT_ECR
33 #define MODE_ACCES MED_ACC_RDWR
34 #elif DEF_LECT_AJOUT
35 #define MODE_ACCES MED_ACC_RDEXT
36 #else
37 #define MODE_ACCES MED_ACC_CREAT
38 #endif
39
40 int main (int argc, char **argv)
41
42 /*TODO : Tester l'écriture des attributs famille,numéros optionnels, noms optionnels*/
43 {
44 med_idt fid=0;
45 med_int mdim=2,axe=0,nind=0;
46 med_float indiceX[4] = {1.0,1.1,1.2,1.3};
47 med_float indiceY[4] = {2.0,2.1,2.2,2.3};
48 med_float coo[8] = {0.0,0.0,1.0,0.0,0.0,1.0,1.0,1.0};
49 med_int nnoeuds = 4;
50 med_int structure_grille[2] = {2,2};
51 char maa [MED_NAME_SIZE+1]= "grille_cartesian";
52 char maa2[MED_NAME_SIZE+1]= "grille_curvilinear";
53 /* composantes et unites */
54 /* 12345678901234561234567890123456 */
55 char comp[2*MED_SNAME_SIZE+1] = "X Y ";
56 char unit[2*MED_SNAME_SIZE+1] = "cm cm ";
57
58 /* Creation du fichier test27.med */
59 fid = MEDfileOpen("test27.med",MODE_ACCES);
60 if (fid < 0) {
61 MESSAGE("Erreur a la creation du fichier test27.med");
62 return -1;
63 }
64 printf("Creation du fichier test27.med \n");
65
66 /* Creation du maillage "maa" de type MED_NON_STRUCURE et de dimension 2 */
67 if (MEDmeshCr( fid, "maillage vide",2, 2, MED_UNSTRUCTURED_MESH,
68 "un maillage vide","s", MED_SORT_DTIT,
69 MED_CARTESIAN, comp, unit) < 0) {
70 MESSAGE("Erreur a la creation du maillage MED_UNSTRUCTURED_MESH : "); SSCRUTE(maa);
71 return -1;
72 }
73 /* creation d'une grille cartesienne de dimension 2 */
74 /* on commence par definir un maillage MED_STRUCTURED_MESH
75 de dimension 2 */
76 if (MEDmeshCr( fid, maa,mdim, mdim, MED_STRUCTURED_MESH,
77 "un exemple de grille cartesienne","s", MED_SORT_DTIT,
78 MED_CARTESIAN, comp, unit) < 0) {
79 MESSAGE("Erreur a la creation de la grille");
80 return -1;
81 }
82 printf("Creation d'un maillage structure MED_STRUCTURED_MESH \n");
83
84 /* On specifie la nature du maillage structure : MED_GRILLE_CARTESIENNE */
85 if (MEDmeshGridTypeWr(fid,maa, MED_CARTESIAN_GRID) < 0) {
86 MESSAGE("Erreur a l'ecriture de la nature de la grille");
87 return -1;
88 }
89
90 printf("On definit la nature du maillage structure : MED_GRILLE_CARTESIENNE \n");
91
92 /* on definit les indices des coordonnees de la grille selon chaque dimension */
93 /* axe des "X" */
94 nind = 4;
95 axe = 1;
96 if (MEDmeshGridIndexCoordinateWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,
97 axe,nind,indiceX) < 0) {
98 MESSAGE("Erreur a l'ecriture de l'axe X");
99 return -1;
100 }
101 printf("Ecriture des indices des coordonnees selon l'axe des X \n");
102
103 /* axe des "Y" */
104 nind = 4;
105 axe = 2;
106 if (MEDmeshGridIndexCoordinateWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,
107 axe,nind,indiceY) < 0) {
108 MESSAGE("Erreur a l'ecriture de l'axe Y");
109 return -1;
110 }
111 printf("Ecriture des indices des coordonnees selon l'axe des Y \n");
112
113 /* Creation d'une grille MED_CURVILINEAR_GRID de dimension 2 */
114 /* on commence par definir un maillage MED_STRUCTURED_MESH
115 de dimension 2 */
116 if (MEDmeshCr( fid, maa2,mdim, mdim, MED_STRUCTURED_MESH,
117 "un exemple de grille standard","s", MED_SORT_DTIT,
118 MED_CARTESIAN, comp, unit) < 0) {
119 MESSAGE("Erreur a la creation de la 2e grille");
120 return -1;
121 }
122 printf("Creation d'un maillage structure MED_STRUCTURED_MESH \n");
123
124 /* On specifie la nature du maillage structure : MED_CURVILINEAR_GRID */
125 if (MEDmeshGridTypeWr(fid,maa2, MED_CURVILINEAR_GRID) < 0) {
126 MESSAGE("Erreur a l'ecriture de la nature de la grille");
127 return -1;
128 }
129 printf("On definit la nature du maillage structure : MED_CURVILINEAR_GRID \n");
130
131
132 if (MEDmeshNodeCoordinateWr(fid,maa2,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,
133 MED_FULL_INTERLACE,nnoeuds, coo) < 0) {
134 MESSAGE("Erreur a l'ecriture des noeuds de la grille MED_CURVILINEAR_GRID");
135 return -1;
136 }
137 printf("Ecriture des coordonnees des noeuds \n");
138
139 /* On definit la structure de la grille */
140 if ( MEDmeshGridStructWr(fid,maa2,MED_NO_DT,MED_NO_IT, MED_UNDEF_DT, structure_grille ) < 0) {
141 MESSAGE("Erreur a l'ecriture de la structure de la grille");
142 return -1;
143 }
144 printf("Ecriture de la structure de la grille : / 2,2 / \n");
145
146 /* On ferme le fichier */
147 if (MEDfileClose(fid) < 0) {
148 MESSAGE("Erreur a la fermeture du fichier");
149 return -1;
150 }
151 printf("Fermeture du fichier \n");
152
153 return 0;
154 }
MED_NO_IT
#define MED_NO_IT
Definition: med.h:316
med_utils.h
MED_CURVILINEAR_GRID
Definition: med.h:137
MEDfileOpen
MEDC_EXPORT med_idt MEDfileOpen(const char *const filename, const med_access_mode accessmode)
Ouverture d'un fichier MED.
Definition: MEDfileOpen.c:41
MED_FULL_INTERLACE
Definition: med.h:94
MED_SORT_DTIT
Definition: med.h:306
med_float
double med_float
Definition: med.h:331
MEDmeshGridTypeWr
MEDC_EXPORT med_err MEDmeshGridTypeWr(const med_idt fid, const char *const meshname, const med_grid_type gridtype)
Cette routine permet de définir le type d'un maillage structuré (MED_STRUCTURED_MESH).
Definition: MEDmeshGridTypeWr.c:37
MED_UNDEF_DT
#define MED_UNDEF_DT
Definition: med.h:317
MESSAGE
#define MESSAGE(chaine)
Definition: med_utils.h:316
MEDfileClose
MEDC_EXPORT med_err MEDfileClose(med_idt fid)
Fermeture d'un fichier MED.
Definition: MEDfileClose.c:30
MEDmeshGridIndexCoordinateWr
MEDC_EXPORT med_err MEDmeshGridIndexCoordinateWr(const med_idt fid, const char *const meshname, const med_int numdt, const med_int numit, const med_float dt, const med_int axis, const med_int indexsize, const med_float *const gridindex)
Cette routine permet l'écriture des coordonnées des noeuds d'un maillage structuré selon un axe du re...
Definition: MEDmeshGridIndexCoordinateWr.c:43
MED_UNSTRUCTURED_MESH
Definition: med.h:129
MODE_ACCES
#define MODE_ACCES
Definition: test27.c:37
med_int
int med_int
Definition: med.h:337
MEDmeshNodeCoordinateWr
MEDC_EXPORT med_err MEDmeshNodeCoordinateWr(const med_idt fid, const char *const meshname, const med_int numdt, const med_int numit, const med_float dt, const med_switch_mode switchmode, const med_int nentity, const med_float *const coordinates)
Cette routine permet d'écrire dans un maillage le tableau des coordonnées des noeuds, selon une séquence de calcul donnée.
Definition: MEDmeshNodeCoordinateWr.c:45
MEDmeshCr
MEDC_EXPORT med_err MEDmeshCr(const med_idt fid, const char *const meshname, const med_int spacedim, const med_int meshdim, const med_mesh_type meshtype, const char *const description, const char *const dtunit, const med_sorting_type sortingtype, const med_axis_type axistype, const char *const axisname, const char *const axisunit)
Cette routine permet de créer un maillage dans un fichier.
Definition: MEDmeshCr.c:45
MED_CARTESIAN
Definition: med.h:255
MED_CARTESIAN_GRID
Definition: med.h:135
med_idt
hid_t med_idt
Definition: med.h:326
MED_NO_DT
#define MED_NO_DT
Definition: med.h:315
MEDmeshGridStructWr
MEDC_EXPORT med_err MEDmeshGridStructWr(const med_idt fid, const char *const meshname, const med_int numdt, const med_int numit, const med_float dt, const med_int *const gridstruct)
Cette routine définit la structure (nombre de points sur chaque axe du repère) d'un maillage structur...
Definition: MEDmeshGridStructWr.c:41
MED_NAME_SIZE
#define MED_NAME_SIZE
Definition: med.h:79
MED_SNAME_SIZE
#define MED_SNAME_SIZE
Definition: med.h:80
med.h
main
int main(int argc, char **argv)
Definition: test27.c:40
SSCRUTE
#define SSCRUTE(chaine)
Definition: med_utils.h:315
MED_STRUCTURED_MESH
Definition: med.h:129