      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 }