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