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