# 2.3.6/test28.c

Aller à la documentation de ce fichier.

    00001 /*  This file is part of MED.
    00002  *
    00003  *  COPYRIGHT (C) 1999 - 2011  EDF R&D, CEA/DEN
    00004  *  MED is free software: you can redistribute it and/or modify
    00005  *  it under the terms of the GNU Lesser General Public License as published by
    00006  *  the Free Software Foundation, either version 3 of the License, or
    00007  *  (at your option) any later version.
    00008  *
    00009  *  MED is distributed in the hope that it will be useful,
    00010  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
    00011  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    00012  *  GNU Lesser General Public License for more details.
    00013  *
    00014  *  You should have received a copy of the GNU Lesser General Public License
    00015  *  along with MED.  If not, see <http://www.gnu.org/licenses/>.
    00016  */
    00017
    00018 /******************************************************************************
    00019  * - Nom du fichier : test28.c
    00020  *
    00021  * - Description : lecture d'un maillage structure (grille cartesienne |
    00022  *                 grille polaire) dans le fichier test27.med
    00023  *
    00024  *****************************************************************************/
    00025
    00026 #include <med.h>
    00027 #define MESGERR 1
    00028 #include "med_utils.h"
    00029 #include <string.h>
    00030
    00031 #ifdef DEF_LECT_ECR
    00032 #define MODE_ACCES MED_LECTURE_ECRITURE
    00033 #elif DEF_LECT_AJOUT
    00034 #define MODE_ACCES MED_LECTURE_AJOUT
    00035 #else
    00036 #define MODE_ACCES MED_CREATION
    00037 #endif
    00038
    00039 int main (int argc, char **argv)
    00040
    00041
    00042 {
    00043   int i;
    00044   med_err ret = 0;
    00045   med_idt fid;
    00046   /* la dimension du maillage */
    00047   med_int mdim,axe,j;
    00048   med_float *indices;
    00049   med_int nind,nmaa;
    00050   /* nom du maillage de longueur maxi MED_TAILLE_NOM */
    00051   char maa[MED_TAILLE_NOM+1];
    00052   /* composantes et unites */
    00053   char comp[MED_TAILLE_PNOM+1];
    00054   char unit[MED_TAILLE_PNOM+1];
    00055   char desc[MED_TAILLE_DESC+1];
    00056   med_maillage type_maillage;
    00057   med_type_grille type;
    00058   med_table quoi;
    00059   med_int nnoeuds;
    00060   char *comp2, *unit2;
    00061   med_float *coo;
    00062   med_int *structure_grille;
    00063   med_repere repere;
    00064
    00065   /* Ouverture du fichier test17.med en lecture seule */
    00066   fid = MEDouvrir("test27.med",MED_LECTURE);
    00067   if (fid < 0) {
    00068     MESSAGE("Erreur a l'ouverture du fichier test27.med");
    00069     return -1;
    00070   }
    00071
    00072   /* Lecture du nombre de maillage */
    00073   nmaa = MEDnMaa(fid);
    00074   if (nmaa < 0) {
    00075     MESSAGE("Erreur a la lecture du nombre de maillage");
    00076     return -1;
    00077   }
    00078
    00079   /* On boucle sur les maillages et on ne lit que la grille cartesienne
    00080      et la grille desturcutree */
    00081   for (i=0;i<nmaa;i++) {
    00082
    00083     /* On repere le maillage qui nous interesse */
    00084     if (MEDmaaInfo(fid,i+1,maa,&mdim,&type_maillage,desc) < 0) {
    00085       MESSAGE("Erreur a la lecture des infos sur le maillage");
    00086       return -1;
    00087     }
    00088     printf("Maillage de nom [%s] de dimension "IFORMAT" \n",maa,mdim);
    00089     if (type_maillage == MED_STRUCTURE)
    00090       printf("- Type : Maillage structure \n");
    00091     else
    00092       printf("- Type : Maillage non structure \n");
    00093
    00094     /* On regarde le type de la grille */
    00095     if (type_maillage == MED_STRUCTURE) {
    00096       if (MEDnatureGrilleLire(fid,maa,&type) < 0) {
    00097         MESSAGE("Erreur a la lecture de la nature d'une grille");
    00098         return -1;
    00099       }
    00100       if (type == MED_GRILLE_CARTESIENNE)
    00101         printf("- Grille cartesienne \n");
    00102       if (type == MED_GRILLE_STANDARD)
    00103         printf("- Grille de-structureee \n");
    00104     }
    00105
    00106     /* On regarde les coordonnees de la grille standard */
    00107     if (type_maillage == MED_STRUCTURE && type == MED_GRILLE_STANDARD) {
    00108
    00109       nnoeuds = MEDnEntMaa(fid,maa,MED_COOR,MED_NOEUD,0,0);
    00110       if (nnoeuds < 0) {
    00111         MESSAGE("Erreur a la lecture du nombre de noeuds");
    00112         return -1;
    00113       }
    00114       printf("Nombre de noeuds : "IFORMAT" \n",nnoeuds);
    00115
    00116       structure_grille = (med_int *) malloc(sizeof(med_int)*mdim);
    00117       if (MEDstructureCoordLire(fid,maa,mdim,structure_grille) < 0) {
    00118         MESSAGE("Erreur a la lecture de la structure de la grille");
    00119         ret = -1;
    00120       }
    00121       if (ret == 0) {
    00122         printf("Structure des noeuds de la grille : [ ");
    00123         for (j=0;j<mdim;j++)
    00124           printf(" "IFORMAT" ",*(structure_grille+j));
    00125         printf(" ] \n");
    00126         free(structure_grille);
    00127       }
    00128
    00129       if (ret == 0) {
    00130         coo = (med_float *) malloc(sizeof(med_float)*nnoeuds*mdim);
    00131         comp2 = (char *) malloc(sizeof(char)*(mdim*MED_TAILLE_PNOM+1));
    00132         unit2 = (char *) malloc(sizeof(char)*(mdim*MED_TAILLE_PNOM+1));
    00133
    00134         if ((ret = MEDcoordLire(fid,maa,mdim,coo,MED_FULL_INTERLACE,MED_ALL,NULL,0,
    00135                                 &repere,comp2,unit2)) < 0) {
    00136           MESSAGE("Erreur a la lecture des coordonnees des noeuds");
    00137           ret = -1;
    00138         }
    00139       }
    00140
    00141       if (ret == 0) {
    00142         printf("Coordonnees : [ ");
    00143         for (j=0;j<nnoeuds*mdim;j++)
    00144           printf(" %f ",*(coo+j));
    00145         printf(" ] \n");
    00146         free(coo);
    00147         free(comp2);
    00148         free(unit2);
    00149       }
    00150     }
    00151
    00152     /* On regarde les coordonnees des indices de la grille cartesienne */
    00153     if (type_maillage == MED_STRUCTURE && type == MED_GRILLE_CARTESIENNE)
    00154       for(axe=1;axe<=mdim;axe++) {
    00155         switch(axe) {
    00156         case 1 :
    00157           quoi = MED_COOR_IND1;
    00158           break;
    00159
    00160         case 2 :
    00161           quoi = MED_COOR_IND2;
    00162           break;
    00163
    00164         case 3 :
    00165           quoi = MED_COOR_IND3;
    00166           break;
    00167
    00168         default :
    00169           return -1;
    00170         }
    00171
    00172         nind = MEDnEntMaa(fid,maa,quoi,MED_NOEUD,0,0);
    00173         if (nind < 0) {
    00174           MESSAGE("Erreur a la lecture de la taille de l'indice");
    00175           return -1;
    00176         }
    00177         printf("Lecture de la taille de l'indice : "IFORMAT" \n",nind);
    00178
    00179           /* on lit le tableau des indices */
    00180         indices = (med_float *) malloc(sizeof(med_float)*nind);
    00181         if (MEDindicesCoordLire(fid,maa,mdim,indices,nind,axe,comp,unit) < 0) {
    00182           MESSAGE("Erreur a lecture de indices de coordonnees");
    00183           ret = -1;
    00184         }
    00185         if (ret == 0) {
    00186           printf("Axe %s [%s] : [ ",comp,unit);
    00187           for (j=0;j<nind;j++)
    00188             printf(" %f ",indices[j]);
    00189           printf(" ] \n");
    00190           free(indices);
    00191         }
    00192       }
    00193   }
    00194
    00195   /* On ferme le fichier */
    00196   if (MEDfermer(fid) < 0) {
    00197     MESSAGE("Erreur a la fermeture du fichier");
    00198     return -1;
    00199   }
    00200   printf("Fermeture du fichier \n");
    00201
    00202   return ret;
    00203 }