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