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