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