src/MEDWrapper/V2_1/MEDgridEcr.cxx

   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.
   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 #include "med_outils.hxx"
  20 #include "med.hxx"
  21
  22 namespace med_2_1{
  23
  24 med_err
  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 )
  27 {
  28     /* ecriture des indices */
  29
  30     med_idt  maaid, noeid, ds;
  31     char     chemin[MED_TAILLE_MAA+MED_TAILLE_NOM+1];
  32     med_size dimd[1];
  33     char     *dataset;
  34     med_int  type_rep_int;
  35
  36     /* On inhibe le gestionnaire d'erreur HDF */
  37     _MEDmodeErreurVerrouiller();
  38
  39     /* Si le maillage n'existe pas => erreur */
  40     strcpy(chemin, MED_MAA);
  41     strcat(chemin, maa);
  42     maaid = _MEDdatagroupOuvrir(fid, chemin);
  43     if (maaid < 0) return(-1);
  44
  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) {
  48             return(-1);
  49         };
  50     };
  51
  52     switch (dim) {
  53         case 0 : {
  54             dataset = MED_NOM_IN1;
  55             break;
  56         };
  57         case 1 : {
  58             dataset = MED_NOM_IN2;
  59             break;
  60         };
  61         case 2 : {
  62             dataset = MED_NOM_IN3;
  63             break;
  64         };
  65         default : {
  66             return(-1);
  67         };
  68     };
  69
  70     /* Creation du Data Set "IN1" ou "IN2" ou "IN3" */
  71     dimd[0] = nb;
  72     if (_MEDdatasetNumEcrire(noeid, dataset, MED_REEL64, mode_coo, 1, MED_ALL, MED_NOPF, 0, 0, dimd, (unsigned char*)coo, mode) < 0) {
  73         return(-1);
  74     };
  75
  76     /* On re-ouvre le Data Set "IN1" ou "IN2" ou "IN3" pour y placer des attributs */
  77     if ((ds = _MEDdatasetOuvrir(noeid, dataset)) < 0) {
  78         return(-1);
  79     };
  80
  81     /* Attribut NBR (nombre de noeuds) */
  82     if (_MEDattrEntierEcrire(ds, MED_NOM_NBR, &nb, mode) < 0) {
  83         return(-1);
  84     };
  85
  86     /* L'attribut "REP" */
  87     type_rep_int = (med_int)repere;
  88     if (_MEDattrEntierEcrire(ds, MED_NOM_REP, &type_rep_int, mode) < 0) {
  89         return(-1);
  90     };
  91
  92     /* Attribut "NOM" */
  93     if (_MEDattrStringEcrire(ds, MED_NOM_NOM, mdim*MED_TAILLE_PNOM, nomcoo, mode) < 0) {
  94         return(-1);
  95     };
  96
  97     /* Attribut "UNI" */
  98     if (_MEDattrStringEcrire(ds, MED_NOM_UNI, mdim*MED_TAILLE_PNOM, unicoo, mode) < 0) {
  99         return(-1);
 100     };
 101
 102     /* On ferme tout */
 103     if (_MEDdatasetFermer(ds)      < 0) return(-1);
 104     if (_MEDdatagroupFermer(noeid) < 0) return(-1);
 105     if (_MEDdatagroupFermer(maaid) < 0) return(-1);
 106     return(0);
 107 }
 108
 109 }