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 : test27.c
 22   *
 23   * - Description : creation de maillages structures (grille cartesienne |
 24   *                 grille standard ) dans le fichier test27.med
 25   *
 26   *****************************************************************************/
 27
 28  #include <med.h>
 29  #include <med_utils.h>
 30  #include <string.h>
 31
 32  int main (int argc, char **argv)
 33  {
 34    med_idt fid;
 35    /* la dimension du maillage */
 36    med_int mdim,axe;
 37    med_float indice[4];
 38    med_float coo[8] = {0.0,0.0,1.0,0.0,0.0,1.0,1.0,1.0};
 39    med_int nnoeuds = 4;
 40    med_int structure_grille[2] = {2,2};
 41    med_int nind;
 42    /* nom du maillage de longueur maxi MED_TAILLE_NOM */
 43    char maa[MED_TAILLE_NOM+1];
 44    /* composantes et unites */
 45    char comp[MED_TAILLE_PNOM+1];
 46    char unit[MED_TAILLE_PNOM+1];
 47    /*                               12345678901234561234567890123456 */
 48    char comp2[2*MED_TAILLE_PNOM+1] = "X               Y               ";
 49    char unit2[2*MED_TAILLE_PNOM+1] = "cm              cm              ";
 50
 51    /* Creation du fichier test27.med */
 52    fid = MEDouvrir("test27.med",MED_CREATION);
 53    if (fid < 0) {
 54      MESSAGE("Erreur a la creation du fichier test27.med");
 55      return -1;
 56    }
 57    printf("Creation du fichier test27.med \n");
 58
 59    /* Creation d'un maillage MED_NON_STRUCTURE */
 60    if (MEDmaaCr(fid,"maillage vide",3,MED_NON_STRUCTURE,"un maillage vide") < 0) {
 61      MESSAGE("Erreur a la creation du maillage MED_NON_STRUCTURE");
 62      return -1;
 63    }
 64
 65    /* creation d'une grille cartesienne de dimension 2 */
 66    strcpy(maa,"grille_cartesienne");
 67    mdim = 2;
 68    /* on commence par definir un maillage MED_STRUCTURE 
 69       de dimension 2 */
 70    if (MEDmaaCr(fid,maa,mdim,MED_STRUCTURE,
 71             "un exemple de grille cartesienne") < 0) {
 72      MESSAGE("Erreur a la creation de la grille");
 73      return -1;
 74    }
 75    printf("Creation d'un maillage structure MED_STRUCTURE \n");
 76
 77    /* On specifie la nature du maillage structure : MED_GRILLE_CARTESIENNE */
 78    if (MEDnatureGrilleEcr(fid,maa,MED_GRILLE_CARTESIENNE) < 0) {
 79      MESSAGE("Erreur a l'ecriture de la nature de la grille");
 80      return -1;
 81    }
 82    printf("On definit la nature du maillage structure : MED_GRILLE_CARTESIENNE \n");
 83
 84    /* on definit les indices des coordonnees de la grille selon chaque dimension  */
 85    /* axe des "X" */
 86    indice[0] = 1.1;
 87    indice[1] = 1.2;
 88    indice[2] = 1.3;
 89    indice[3] = 1.4;
 90    nind = 4;
 91    axe = 1;
 92    strcpy(comp,"X");
 93    strcpy(unit,"cm");
 94    if (MEDindicesCoordEcr(fid,maa,mdim,indice,nind,axe,comp,unit) < 0) {
 95      MESSAGE("Erreur a l'ecriture de l'axe X");
 96      return -1;
 97    }
 98    printf("Ecriture des indices des coordonnees selon l'axe des X \n");
 99
100    /* axe des "Y" */
101    indice[0] = 2.1;
102    indice[1] = 2.2;
103    indice[2] = 2.3;
104    indice[3] = 2.4;
105    nind = 4;
106    strcpy(comp,"Y");
107    strcpy(unit,"cm");
108    axe = 2;
109    if (MEDindicesCoordEcr(fid,maa,mdim,indice,nind,axe,comp,unit) < 0) {
110      MESSAGE("Erreur a l'ecriture de l'axe Y");
111      return -1;
112    }
113    printf("Ecriture des indices des coordonnees selon l'axe des Y \n");
114
115    /* Creation d'une grille MED_STANDARD de dimension 2 */
116    strcpy(maa,"grille_standard");
117    mdim = 2;
118    /* on commence par definir un maillage MED_STRUCTURE 
119       de dimension 2 */
120    if (MEDmaaCr(fid,maa,mdim,MED_STRUCTURE,
121             "un exemple de grille standard") < 0) {
122      MESSAGE("Erreur a la creation de la 2e grille");
123      return -1;
124    }
125    printf("Creation d'un maillage structure MED_STRUCTURE \n");
126
127    /* On specifie la nature du maillage structure : MED_GRILLE_STANDARD */
128    if (MEDnatureGrilleEcr(fid,maa,MED_GRILLE_STANDARD) < 0) {
129      MESSAGE("Erreur a l'ecriture du type de la grille");
130      return -1;
131    }
132    printf("On definit la nature du maillage structure : MED_GRILLE_STANDARD \n");
133
134    /* On ecrit les coordonnees de la grille */
135    if (MEDcoordEcr(fid,maa,mdim,coo,MED_FULL_INTERLACE,nnoeuds,MED_CART,comp2,unit2) < 0) {
136      MESSAGE("Erreur a l'ecriture des noeuds de la grille MED_GRILLE_STANDARD");
137      return -1;
138    }
139    printf("Ecriture des coordonnees des noeuds \n");
140
141    /* On definit la structure de la grille */
142    if (MEDstructureCoordEcr(fid,maa,mdim,structure_grille) < 0) {
143      MESSAGE("Erreur a l'ecriture de la structure de la grille");
144      return -1;
145    }
146    printf("Ecriture de la structure de la grille : / 2,2 / \n");
147
148    /* On ferme le fichier */
149    if (MEDfermer(fid) < 0) {
150      MESSAGE("Erreur a la fermeture du fichier");
151      return -1;
152    }
153    printf("Fermeture du fichier \n");
154
155    return 0;
156  }