1 /*************************************************************************
2 * COPYRIGHT (C) 1999 - 2002 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.
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.
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
17 *************************************************************************/
19 #include "med_outils.hxx"
25 MEDgridEcr(med_idt fid, char *maa, med_int mdim, med_float *coo, med_int nb, med_int dim, med_mode_switch mode_coo,
26 med_repere repere, char *nomcoo, char *unicoo, med_mode_acces mode )
28 /* ecriture des indices */
30 med_idt maaid, noeid, ds;
31 char chemin[MED_TAILLE_MAA+MED_TAILLE_NOM+1];
36 /* On inhibe le gestionnaire d'erreur HDF */
37 _MEDmodeErreurVerrouiller();
39 /* Si le maillage n'existe pas => erreur */
40 strcpy(chemin, MED_MAA);
42 maaid = _MEDdatagroupOuvrir(fid, chemin);
43 if (maaid < 0) return(-1);
45 /* Si le Data Group "NOE" n'existe pas on le cree */
46 if ((noeid = _MEDdatagroupOuvrir(maaid, MED_NOM_NOE)) < 0) {
47 if ((noeid = _MEDdatagroupCreer(maaid, MED_NOM_NOE)) < 0) {
54 dataset = MED_NOM_IN1;
58 dataset = MED_NOM_IN2;
62 dataset = MED_NOM_IN3;
70 /* Creation du Data Set "IN1" ou "IN2" ou "IN3" */
72 if (_MEDdatasetNumEcrire(noeid, dataset, MED_REEL64, mode_coo, 1, MED_ALL, MED_NOPF, 0, 0, dimd, (unsigned char*)coo, mode) < 0) {
76 /* On re-ouvre le Data Set "IN1" ou "IN2" ou "IN3" pour y placer des attributs */
77 if ((ds = _MEDdatasetOuvrir(noeid, dataset)) < 0) {
81 /* Attribut NBR (nombre de noeuds) */
82 if (_MEDattrEntierEcrire(ds, MED_NOM_NBR, &nb, mode) < 0) {
86 /* L'attribut "REP" */
87 type_rep_int = (med_int)repere;
88 if (_MEDattrEntierEcrire(ds, MED_NOM_REP, &type_rep_int, mode) < 0) {
93 if (_MEDattrStringEcrire(ds, MED_NOM_NOM, mdim*MED_TAILLE_PNOM, nomcoo, mode) < 0) {
98 if (_MEDattrStringEcrire(ds, MED_NOM_UNI, mdim*MED_TAILLE_PNOM, unicoo, mode) < 0) {
103 if (_MEDdatasetFermer(ds) < 0) return(-1);
104 if (_MEDdatagroupFermer(noeid) < 0) return(-1);
105 if (_MEDdatagroupFermer(maaid) < 0) return(-1);