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  * - Nom du fichier : test32.c
 21  *
 22  * - Description : lecture nominale d'une numerotation globale dans un maillage MED 
 23  *
 24  *****************************************************************************/
 25 
 26 #include <med.h>
 27 #define MESGERR
 28 #include <med_utils.h>
 29 
 30 #ifdef DEF_LECT_ECR
 31 #define MODE_ACCES MED_LECTURE_ECRITURE
 32 #elif DEF_LECT_AJOUT
 33 #define MODE_ACCES MED_LECTURE_AJOUT
 34 #else
 35 #define MODE_ACCES MED_CREATION
 36 #endif
 37 
 38 int main (int argc, char **argv)
 39 
 40 
 41 {
 42   med_err ret = 0;
 43   med_idt fid;
 44   /* la dimension du maillage */
 45   med_int mdim;
 46   /* nom du maillage de longueur maxi MED_TAILLE_NOM */
 47   char maa[MED_TAILLE_NOM+1];
 48   /* le nombre de noeuds */
 49   med_int nnoe = 0;
 50   /* table des numeros global */
 51   med_int *numglobalnoe;
 52 
 53   /* variable de stockage pour reperer le maillage */
 54   med_int i;
 55   char des[MED_TAILLE_DESC+1];
 56   med_maillage type;
 57 
 58   if (argc != 2) {
 59     MESSAGE("Il faut passer un fichier MED en paramètre");
 60     return -1;
 61   }
 62 
 63   /* Ouverture du fichier passe en argument en lecture seule */
 64   if ((fid = MEDouvrir(argv[1],MED_LECTURE)) < 0) {
 65     MESSAGE("Erreur a l'ouverture du fichier : "); SSCRUTE(argv[1]);
 66     return -1;
 67   }
 68 
 69   /* Lecture des infos sur le premier maillage */
 70   if (MEDmaaInfo(fid,1,maa,&mdim,&type,des) < 0) {
 71     MESSAGE("Erreur a lecture des infos sur le 1er maillage");
 72     return -1;
 73   }
 74   printf("Maillage de nom %s et de dimension %d \n",maa,mdim);
 75 
 76 
 77   /* Lecture du nombre de noeuds */
 78   if ((nnoe = MEDnEntMaa(fid,maa,MED_COOR,MED_NOEUD,0,0)) < 0) {
 79     MESSAGE("Erreur a la lecture du nombre de noeuds ");
 80     return -1;
 81   }
 82   printf("Nombre de noeuds : %d \n",nnoe);
 83 
 84   /* Allocations memoires */
 85 
 86   /* table de la numerotation globale
 87      profil : (nombre de noeuds ) */
 88   if (nnoe > 0) {
 89     numglobalnoe = (med_int*) malloc(sizeof(med_int)*nnoe);
 90   }
 91 
 92   /* lecture de la numerotation globale attachee aux noeuds*/
 93   if (MEDglobalNumLire(fid,maa,numglobalnoe,nnoe,MED_NOEUD,0)<0) {
 94     MESSAGE("Erreur a la lecture de de la numerotation globale");
 95     return -1;
 96   }
 97 
 98   /* ecriture a l'ecran des resultats */
 99   for (i=0;i<nnoe;i++)
100     printf("Numero global du noeud %d : %d \n",i+1,numglobalnoe[i]);
101 
102 
103 
104   free(numglobalnoe);
105 
106   /* Fermeture du fichier */
107   if (MEDfermer(fid) < 0) {
108     MESSAGE("Erreur a la fermeture du fichier ");
109     return -1;
110   }
111 
112   return 0;
113 }
114