      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  }