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 : test4.c
22 *
23 * - Description : ecriture des noeuds d'un maillage MED.
24 *
25 *****************************************************************************/
26
27 #include <med.h>
28 #define MESGERR
29 #include <med_utils.h>
30
31 #ifdef DEF_LECT_ECR
32 #define MODE_ACCES MED_LECTURE_ECRITURE
33 #elif DEF_LECT_AJOUT
34 #define MODE_ACCES MED_LECTURE_AJOUT
35 #else
36 #define MODE_ACCES MED_CREATION
37 #endif
38
39 int main (int argc, char **argv)
40
41
42 {
43 med_err ret = 0;
44 med_idt fid;
45 /* la dimension du maillage */
46 med_int mdim = 2;
47 /* nom du maillage de longueur maxi MED_TAILLE_NOM */
48 char maa[MED_TAILLE_NOM+1] = "maa1";
49 /* le nombre de noeuds */
50 med_int nnoe = 4;
51 /* table des coordonnees
52 profil : (dimension * nombre de noeuds) */
53 med_float coo[8] = {0.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0};
54 /* tables des noms et des unites des coordonnees
55 profil : (dimension*MED_TAILLE_PNOM+1) */
56 /* 12345678901234561234567890123456*/
57 char nomcoo[2*MED_TAILLE_PNOM+1] = "x y ";
58 char unicoo[2*MED_TAILLE_PNOM+1] = "cm cm ";
59 /* tables des noms, numeros, numeros de familles des noeuds
60 autant d'elements que de noeuds - les noms ont pout longueur
61 MED_TAILLE_PNOM */
62 /* 1234567890123456123456789012345612345678901234561234567890123456*/
63 char nomnoe[4*MED_TAILLE_PNOM+1] = "nom1 nom2 nom3 nom4 ";
64 med_int numnoe[4] = {1,2,3,4};
65 med_int nufano[4] = {0,1,2,2};
66
67 /* Creation du fichier "test4.med" */
68 fid = MEDouvrir("test4.med",MODE_ACCES);
69 if (fid < 0) {
70 MESSAGE("Erreur a la creation du fichier test4.med");
71 return -1;
72 }
73
74 /* Creation du maillage "maa" de type MED_NON_STRUCURE
75 et de dimension 2 */
76 if (MEDmaaCr(fid,maa,mdim,MED_NON_STRUCTURE,
77 "un maillage pour test4") < 0) {
78 MESSAGE("Erreur a la creation du maillage : "); SSCRUTE(maa);
79 ret = -1;
80 }
81
82 /* Ecriture des coordonnees des noeuds en mode MED_FULL_INTERLACE :
83 (X1,Y1, X2,Y2, X3,Y3, ...) dans un repere cartesien */
84 if (MEDcoordEcr(fid,maa,mdim,coo,MED_FULL_INTERLACE,nnoe,
85 MED_CART,nomcoo,unicoo) < 0) {
86 MESSAGE("Erreur a l'ecriture des coordonnees des noeuds");
87 ret = -1;
88 }
89
90 /* Ecriture des noms des noeuds (optionnel dans un maillage MED) */
91 if (MEDnomEcr(fid,maa,nomnoe,nnoe,MED_NOEUD,0) < 0) {
92 MESSAGE("Erreur a l'ecriture des noms des noeuds");
93 ret = -1;
94 }
95
96 /* Ecriture des numeros des noeuds (optionnel dans un maillage MED) */
97 if (MEDnumEcr(fid,maa,numnoe,nnoe,MED_NOEUD,0) < 0) {
98 MESSAGE("Erreur a l'ecriture des numeros des noeuds");
99 ret = -1;
100 }
101
102 /* Ecriture des numeros de famille des noeuds */
103 if (MEDfamEcr(fid,maa,nufano,nnoe,MED_NOEUD,0) < 0) {
104 MESSAGE("Erreur a l'ecriture des numeros de familles des noeuds");
105 ret = -1;
106 }
107
108 /* Fermeture du fichier */
109 if (MEDfermer(fid) < 0) {
110 MESSAGE("Erreur a la fermeture du fichier test4.med");
111 return -1;
112 }
113
114 return ret;
115 }
116
117
118
119