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   * - Nom du fichier : test28.c
 21   *
 22   * - Description : lecture d'un maillage structure (grille cartesienne |
 23   *                 grille polaire) dans le fichier test27.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    int i;
 34    med_err ret = 0;
 35    med_idt fid;
 36    /* la dimension du maillage */
 37    med_int mdim,axe,j;
 38    med_float *indices;
 39    med_int nind,nmaa;
 40    /* nom du maillage de longueur maxi MED_TAILLE_NOM */
 41    char maa[MED_TAILLE_NOM+1];
 42    /* composantes et unites */
 43    char comp[MED_TAILLE_PNOM+1];
 44    char unit[MED_TAILLE_PNOM+1];
 45    char desc[MED_TAILLE_DESC+1];
 46    med_maillage type_maillage;
 47    med_type_grille type;
 48    med_table quoi;
 49    med_int nnoeuds;
 50    char *comp2, *unit2;
 51    med_float *coo;
 52    med_int *structure_grille;
 53    med_repere repere;
 54
 55    /* Ouverture du fichier test17.med en lecture seule */
 56    fid = MEDouvrir("test27.med",MED_LECTURE);
 57    if (fid < 0) {
 58      MESSAGE("Erreur a l'ouverture du fichier test27.med");
 59      return -1;
 60    }
 61
 62    /* Lecture du nombre de maillage */
 63    nmaa = MEDnMaa(fid);
 64    if (nmaa < 0) {
 65      MESSAGE("Erreur a la lecture du nombre de maillage");
 66      return -1;
 67    }
 68
 69    /* On boucle sur les maillages et on ne lit que la grille cartesienne 
 70       et la grille desturcutree */
 71    for (i=0;i<nmaa;i++) {
 72
 73      /* On repere le maillage qui nous interesse */
 74      if (MEDmaaInfo(fid,i+1,maa,&mdim,&type_maillage,desc) < 0) {
 75        MESSAGE("Erreur a la lecture des infos sur le maillage");
 76        return -1;
 77      }
 78      printf("Maillage de nom [%s] de dimension %d \n",maa,mdim);
 79      if (type_maillage == MED_STRUCTURE)
 80        printf("- Type : Maillage structure \n");
 81      else
 82        printf("- Type : Maillage non structure \n");
 83
 84      /* On regarde le type de la grille */
 85      if (type_maillage == MED_STRUCTURE) {
 86        if (MEDnatureGrilleLire(fid,maa,&type) < 0) {
 87      MESSAGE("Erreur a la lecture de la nature d'une grille");
 88      return -1;
 89        }
 90        if (type == MED_GRILLE_CARTESIENNE)
 91      printf("- Grille cartesienne \n");
 92        if (type == MED_GRILLE_STANDARD)
 93      printf("- Grille de-structureee \n");
 94      }
 95
 96      /* On regarde les coordonnees de la grille standard */
 97      if (type_maillage == MED_STRUCTURE && type == MED_GRILLE_STANDARD) {
 98
 99        nnoeuds = MEDnEntMaa(fid,maa,MED_COOR,MED_NOEUD,0,0);
100        if (nnoeuds < 0) {
101      MESSAGE("Erreur a la lecture du nombre de noeuds");
102      return -1;
103        }
104        printf("Nombre de noeuds : %d \n",nnoeuds);
105
106        structure_grille = (med_int *) malloc(sizeof(med_int)*mdim);
107        if (MEDstructureCoordLire(fid,maa,mdim,structure_grille) < 0) {
108      MESSAGE("Erreur a la lecture de la structure de la grille");
109      ret = -1;
110        }
111        if (ret == 0) {
112      printf("Structure des noeuds de la grille : [ ");
113      for (j=0;j<mdim;j++)
114        printf(" %d ",*(structure_grille+j));
115      printf(" ] \n");
116      free(structure_grille);
117        }
118
119        if (ret == 0) {
120      coo = (med_float *) malloc(sizeof(med_float)*nnoeuds*mdim);
121      comp2 = (char *) malloc(sizeof(char)*(mdim*MED_TAILLE_PNOM+1));
122      unit2 = (char *) malloc(sizeof(char)*(mdim*MED_TAILLE_PNOM+1));
123
124      if (ret = MEDcoordLire(fid,maa,mdim,coo,MED_FULL_INTERLACE,MED_ALL,NULL,0,
125                     &repere,comp2,unit2) < 0) {
126        MESSAGE("Erreur a la lecture des coordonnees des noeuds");
127        ret = -1;
128      }
129        }
130
131        if (ret == 0) {
132      printf("Coordonnees : [ ");
133      for (j=0;j<nnoeuds*mdim;j++)
134        printf(" %f ",*(coo+j));
135      printf(" ] \n");
136      free(coo);
137      free(comp2);
138      free(unit2);
139        }
140      }
141
142      /* On regarde les coordonnees des indices de la grille cartesienne */
143      if (type_maillage == MED_STRUCTURE && type == MED_GRILLE_CARTESIENNE)
144        for(axe=1;axe<=mdim;axe++) {
145      switch(axe) {
146      case 1 :
147        quoi = MED_COOR_IND1;
148        break;
149
150      case 2 :
151        quoi = MED_COOR_IND2;
152        break;
153
154      case 3 :
155        quoi = MED_COOR_IND3;
156        break;
157
158      default :
159        return -1;
160      }
161
162      nind = MEDnEntMaa(fid,maa,quoi,MED_NOEUD,0,0);
163      if (nind < 0) {
164        MESSAGE("Erreur a la lecture de la taille de l'indice");
165        return -1;
166      }
167      printf("Lecture de la taille de l'indice : %d \n",nind);
168
169        /* on lit le tableau des indices */
170      indices = (med_float *) malloc(sizeof(med_float)*nind);
171      if (MEDindicesCoordLire(fid,maa,mdim,indices,nind,axe,comp,unit) < 0) {
172        MESSAGE("Erreur a lecture de indices de coordonnees");
173        ret = -1;
174      }
175      if (ret == 0) {
176        printf("Axe %s [%s] : [ ",comp,unit);
177        for (j=0;j<nind;j++)
178          printf(" %f ",indices[j]);
179        printf(" ] \n");
180        free(indices);
181      }
182        }
183    }
184
185    /* On ferme le fichier */
186    if (MEDfermer(fid) < 0) {
187      MESSAGE("Erreur a la fermeture du fichier");
188      return -1;
189    }
190    printf("Fermeture du fichier \n");
191
192    return ret;
193  }