1 /*************************************************************************
2 * COPYRIGHT (C) 1999 - 2007 EDF R&D, CEA/DEN
3 * THIS LIBRARY IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
4 * IT UNDER THE TERMS OF THE GNU LESSER GENERAL PUBLIC LICENSE
5 * AS PUBLISHED BY THE FREE SOFTWARE FOUNDATION;
6 * EITHER VERSION 2.1 OF THE LICENSE, OR (AT YOUR OPTION) ANY LATER VERSION.
7 *
8 * THIS LIBRARY IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
9 * WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
10 * MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
11 * LESSER GENERAL PUBLIC LICENSE FOR MORE DETAILS.
12 *
13 * YOU SHOULD HAVE RECEIVED A COPY OF THE GNU LESSER GENERAL PUBLIC LICENSE
14 * ALONG WITH THIS LIBRARY; IF NOT, WRITE TO THE FREE SOFTWARE FOUNDATION,
15 * INC., 59 TEMPLE PLACE, SUITE 330, BOSTON, MA 02111-1307 USA
16 *
17 *************************************************************************/
18
19
20 /******************************************************************************
21 * - Nom du fichier : test27.c
22 *
23 * - Description : creation de maillages structures (grille cartesienne |
24 * grille standard ) dans le fichier test27.med
25 *
26 *****************************************************************************/
27
28 #include <med.h>
29 #define MESGERR
30 #include <med_utils.h>
31
32 #ifdef DEF_LECT_ECR
33 #define MODE_ACCES MED_LECTURE_ECRITURE
34 #elif DEF_LECT_AJOUT
35 #define MODE_ACCES MED_LECTURE_AJOUT
36 #else
37 #define MODE_ACCES MED_CREATION
38 #endif
39
40 int main (int argc, char **argv)
41
42
43 {
44 med_idt fid;
45 /* la dimension du maillage */
46 med_int mdim,axe;
47 med_float indice[4];
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 med_int nind;
52 /* nom du maillage de longueur maxi MED_TAILLE_NOM */
53 char maa[MED_TAILLE_NOM+1];
54 /* composantes et unites */
55 char comp[MED_TAILLE_PNOM+1];
56 char unit[MED_TAILLE_PNOM+1];
57 /* 12345678901234561234567890123456 */
58 char comp2[2*MED_TAILLE_PNOM+1] = "X Y ";
59 char unit2[2*MED_TAILLE_PNOM+1] = "cm cm ";
60
61 /* Creation du fichier test27.med */
62 fid = MEDouvrir("test27.med",MODE_ACCES);
63 if (fid < 0) {
64 MESSAGE("Erreur a la creation du fichier test27.med");
65 return -1;
66 }
67 printf("Creation du fichier test27.med \n");
68
69 /* Creation d'un maillage MED_NON_STRUCTURE */
70 if (MEDmaaCr(fid,"maillage vide",3,MED_NON_STRUCTURE,"un maillage vide") < 0) {
71 MESSAGE("Erreur a la creation du maillage MED_NON_STRUCTURE");
72 return -1;
73 }
74
75 /* creation d'une grille cartesienne de dimension 2 */
76 strcpy(maa,"grille_cartesienne");
77 mdim = 2;
78 /* on commence par definir un maillage MED_STRUCTURE
79 de dimension 2 */
80 if (MEDmaaCr(fid,maa,mdim,MED_STRUCTURE,
81 "un exemple de grille cartesienne") < 0) {
82 MESSAGE("Erreur a la creation de la grille");
83 return -1;
84 }
85 printf("Creation d'un maillage structure MED_STRUCTURE \n");
86
87 /* On specifie la nature du maillage structure : MED_GRILLE_CARTESIENNE */
88 if (MEDnatureGrilleEcr(fid,maa,MED_GRILLE_CARTESIENNE) < 0) {
89 MESSAGE("Erreur a l'ecriture de la nature de la grille");
90 return -1;
91 }
92 printf("On definit la nature du maillage structure : MED_GRILLE_CARTESIENNE \n");
93
94 /* on definit les indices des coordonnees de la grille selon chaque dimension */
95 /* axe des "X" */
96 indice[0] = 1.1;
97 indice[1] = 1.2;
98 indice[2] = 1.3;
99 indice[3] = 1.4;
100 nind = 4;
101 axe = 1;
102 strcpy(comp,"X");
103 strcpy(unit,"cm");
104 if (MEDindicesCoordEcr(fid,maa,mdim,indice,nind,axe,comp,unit) < 0) {
105 MESSAGE("Erreur a l'ecriture de l'axe X");
106 return -1;
107 }
108 printf("Ecriture des indices des coordonnees selon l'axe des X \n");
109
110 /* axe des "Y" */
111 indice[0] = 2.1;
112 indice[1] = 2.2;
113 indice[2] = 2.3;
114 indice[3] = 2.4;
115 nind = 4;
116 strcpy(comp,"Y");
117 strcpy(unit,"cm");
118 axe = 2;
119 if (MEDindicesCoordEcr(fid,maa,mdim,indice,nind,axe,comp,unit) < 0) {
120 MESSAGE("Erreur a l'ecriture de l'axe Y");
121 return -1;
122 }
123 printf("Ecriture des indices des coordonnees selon l'axe des Y \n");
124
125 /* Creation d'une grille MED_STANDARD de dimension 2 */
126 strcpy(maa,"grille_standard");
127 mdim = 2;
128 /* on commence par definir un maillage MED_STRUCTURE
129 de dimension 2 */
130 if (MEDmaaCr(fid,maa,mdim,MED_STRUCTURE,
131 "un exemple de grille standard") < 0) {
132 MESSAGE("Erreur a la creation de la 2e grille");
133 return -1;
134 }
135 printf("Creation d'un maillage structure MED_STRUCTURE \n");
136
137 /* On specifie la nature du maillage structure : MED_GRILLE_STANDARD */
138 if (MEDnatureGrilleEcr(fid,maa,MED_GRILLE_STANDARD) < 0) {
139 MESSAGE("Erreur a l'ecriture du type de la grille");
140 return -1;
141 }
142 printf("On definit la nature du maillage structure : MED_GRILLE_STANDARD \n");
143
144 /* On ecrit les coordonnees de la grille */
145 if (MEDcoordEcr(fid,maa,mdim,coo,MED_FULL_INTERLACE,nnoeuds,MED_CART,comp2,unit2) < 0) {
146 MESSAGE("Erreur a l'ecriture des noeuds de la grille MED_GRILLE_STANDARD");
147 return -1;
148 }
149 printf("Ecriture des coordonnees des noeuds \n");
150
151 /* On definit la structure de la grille */
152 if (MEDstructureCoordEcr(fid,maa,mdim,structure_grille) < 0) {
153 MESSAGE("Erreur a l'ecriture de la structure de la grille");
154 return -1;
155 }
156 printf("Ecriture de la structure de la grille : / 2,2 / \n");
157
158 /* On ferme le fichier */
159 if (MEDfermer(fid) < 0) {
160 MESSAGE("Erreur a la fermeture du fichier");
161 return -1;
162 }
163 printf("Fermeture du fichier \n");
164
165 return 0;
166 }