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