1 /*************************************************************************
2 * COPYRIGHT (C) 1999 - 2003 EDF R&D
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 : test6.c
22 *
23 * - Description : ecriture d'elements dans un maillage MED
24 *
25 *****************************************************************************/
26
27 #include <med.h>
28 #include <med_utils.h>
29 #include <stdio.h>
30 #include <string.h>
31 #include <stdlib.h>
32
33 int main (int argc, char **argv)
34 {
35 med_idt fid;
36 med_int nse2 = 5;
37 med_int se2[10] = {1,2,1,3,2,4,3,4,2,3};
38 /* 12345678901234561234567890123456123456789012345612345678901234561234567890123456*/
39 char nomse2[MED_TAILLE_PNOM*5+1]="se1 se2 se3 se4 se5 ";
40 med_int numse2[5] = {1,2,3,4,5};
41 med_int nufase2[5] = {-1,-1,0,-2,-3};
42 med_int ntr3 = 2;
43 med_int tr3[6] = {1,2,-5,-5,3,-4};
44 /* 12345678901234561234567890123456*/
45 char nomtr3[MED_TAILLE_PNOM*2+1] = "tr1 tr2 ";
46 med_int numtr3[2] = {4,5};
47 med_int nufatr3[2] = {0,-1};
48 char maa[MED_TAILLE_NOM+1] = "maa1";
49 med_int mdim = 2;
50
51 /* ouverture du fichier */
52 if ( (fid = MEDouvrir("test6.med",MED_CREATION)) < 0) {
53 MESSAGE("Impossible de creer le fichier test6.med : ");
54 return -1;
55 }
56
57 /* creation du maillage */
58 if ( MEDmaaCr(fid,maa,mdim,MED_NON_STRUCTURE,
59 "un maillage pour test6") < 0 ) {
60 MESSAGE("Impossible de creer le maillage : ");
61 return -1;
62 }
63
64 /* ecriture des connectivites des segments */
65 if ( MEDconnEcr(fid,maa,mdim,se2,MED_FULL_INTERLACE,nse2,MED_ARETE,MED_SEG2,
66 MED_DESC) < 0 ) {
67 MESSAGE("Impossible d'ecrire la connectivité des segments: ");
68 return -1;
69 }
70
71
72 /* ecriture (optionnelle) des noms des segments */
73 if ( MEDnomEcr(fid,maa,nomse2,nse2,MED_ARETE,MED_SEG2) < 0 ) {
74 MESSAGE("Impossible d'ecrire les noms des segments : ");
75 return -1;
76 }
77
78 /* ecriture (optionnelle) des numeros des segments */
79 if ( MEDnumEcr(fid,maa,numse2,nse2,MED_ARETE,MED_SEG2) < 0 ) {
80 MESSAGE("Impossible d'ecrire les numéros des segments : ");
81 return -1;
82 }
83
84 /* ecriture des numeros des familles des segments */
85 if ( MEDfamEcr(fid,maa,nufase2,nse2,MED_ARETE,MED_SEG2) < 0 ) {
86 MESSAGE("Impossible d'ecrire les numéros de famille des segments : ");
87 return -1;
88 }
89
90 /* ecriture des connectivites des triangles */
91 if ( MEDconnEcr(fid,maa,mdim,tr3,MED_FULL_INTERLACE,ntr3,MED_MAILLE,MED_TRIA3,
92 MED_DESC) < 0 ) {
93 MESSAGE("Impossible d'ecrire la connectivité des triangles : ");
94 return -1;
95 }
96
97 /* ecriture (optionnelle) des noms des triangles */
98 if ( MEDnomEcr(fid, maa, nomtr3, ntr3, MED_MAILLE, MED_TRIA3) < 0 ) {
99 MESSAGE("Impossible d'ecrire les noms des triangles : ");
100 return -1;
101 }
102
103 /* ecriture (optionnelle) des numeros des triangles */
104 if ( MEDnumEcr(fid,maa,numtr3,ntr3,MED_MAILLE,MED_TRIA3) < 0 ) {
105 MESSAGE("Impossible d'ecrire les numéros des triangles : ");
106 return -1;
107 }
108
109 /* ecriture des numeros des familles des triangles */
110 if ( MEDfamEcr(fid,maa,nufatr3,ntr3,MED_MAILLE,MED_TRIA3) < 0 ) {
111 MESSAGE("Impossible d'ecrire les numéros de famille des triangles : ");
112 return -1;
113 }
114
115 /* fermeture du fichier */
116 if ( MEDfermer(fid) < 0 ) {
117 MESSAGE("Impossible de fermerle fichier : ");
118 return -1;
119 }
120
121 return 0;
122
123 }
124
125
126