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