src/MEDWrapper/V2_1/MEDgridEcr.cxx

   1 // Copyright (C) 2005  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
   2 // CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS
   3 //
   4 // This library is free software; you can redistribute it and/or
   5 // modify it under the terms of the GNU Lesser General Public
   6 // License as published by the Free Software Foundation; either
   7 // version 2.1 of the License.
   8 //
   9 // This library is distributed in the hope that it will be useful
  10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12 // Lesser General Public License for more details.
  13 //
  14 // You should have received a copy of the GNU Lesser General Public
  15 // License along with this library; if not, write to the Free Software
  16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  17 //
  18 // See http://www.salome-platform.org/
  19 //
  20 /*************************************************************************
  21 * COPYRIGHT (C) 1999 - 2002  EDF R&D
  22 * THIS LIBRARY IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
  23 * IT UNDER THE TERMS OF THE GNU LESSER GENERAL PUBLIC LICENSE
  24 * AS PUBLISHED BY THE FREE SOFTWARE FOUNDATION;
  25 * EITHER VERSION 2.1 OF THE LICENSE, OR (AT YOUR OPTION) ANY LATER VERSION.
  26 *
  27 * THIS LIBRARY IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
  28 * WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
  29 * MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
  30 * LESSER GENERAL PUBLIC LICENSE FOR MORE DETAILS.
  31 *
  32 * YOU SHOULD HAVE RECEIVED A COPY OF THE GNU LESSER GENERAL PUBLIC LICENSE
  33 * ALONG WITH THIS LIBRARY; IF NOT, WRITE TO THE FREE SOFTWARE FOUNDATION,
  34 * INC., 59 TEMPLE PLACE, SUITE 330, BOSTON, MA 02111-1307 USA
  35 *
  36 *************************************************************************/
  37
  38 #include "med_outils.hxx"
  39 #include "med.hxx"
  40
  41 namespace med_2_1{
  42
  43 med_err
  44 MEDgridEcr(med_idt fid, char *maa, med_int mdim, med_float *coo, med_int nb, med_int dim, med_mode_switch mode_coo,
  45            med_repere repere, char *nomcoo, char *unicoo, med_mode_acces mode )
  46 {
  47     /* ecriture des indices */
  48
  49     med_idt  maaid, noeid, ds;
  50     char     chemin[MED_TAILLE_MAA+MED_TAILLE_NOM+1];
  51     med_size dimd[1];
  52     char     *dataset;
  53     med_int  type_rep_int;
  54
  55     /* On inhibe le gestionnaire d'erreur HDF */
  56     _MEDmodeErreurVerrouiller();
  57
  58     /* Si le maillage n'existe pas => erreur */
  59     strcpy(chemin, MED_MAA);
  60     strcat(chemin, maa);
  61     maaid = _MEDdatagroupOuvrir(fid, chemin);
  62     if (maaid < 0) return(-1);
  63
  64     /* Si le Data Group "NOE" n'existe pas on le cree */
  65     if ((noeid = _MEDdatagroupOuvrir(maaid, MED_NOM_NOE)) < 0) {
  66         if ((noeid = _MEDdatagroupCreer(maaid, MED_NOM_NOE)) < 0) {
  67             return(-1);
  68         };
  69     };
  70
  71     switch (dim) {
  72         case 0 : {
  73             dataset = MED_NOM_IN1;
  74             break;
  75         };
  76         case 1 : {
  77             dataset = MED_NOM_IN2;
  78             break;
  79         };
  80         case 2 : {
  81             dataset = MED_NOM_IN3;
  82             break;
  83         };
  84         default : {
  85             return(-1);
  86         };
  87     };
  88
  89     /* Creation du Data Set "IN1" ou "IN2" ou "IN3" */
  90     dimd[0] = nb;
  91     if (_MEDdatasetNumEcrire(noeid, dataset, MED_REEL64, mode_coo, 1, MED_ALL, MED_NOPF, 0, 0, dimd, (unsigned char*)coo, mode) < 0) {
  92         return(-1);
  93     };
  94
  95     /* On re-ouvre le Data Set "IN1" ou "IN2" ou "IN3" pour y placer des attributs */
  96     if ((ds = _MEDdatasetOuvrir(noeid, dataset)) < 0) {
  97         return(-1);
  98     };
  99
 100     /* Attribut NBR (nombre de noeuds) */
 101     if (_MEDattrEntierEcrire(ds, MED_NOM_NBR, &nb, mode) < 0) {
 102         return(-1);
 103     };
 104
 105     /* L'attribut "REP" */
 106     type_rep_int = (med_int)repere;
 107     if (_MEDattrEntierEcrire(ds, MED_NOM_REP, &type_rep_int, mode) < 0) {
 108         return(-1);
 109     };
 110
 111     /* Attribut "NOM" */
 112     if (_MEDattrStringEcrire(ds, MED_NOM_NOM, mdim*MED_TAILLE_PNOM, nomcoo, mode) < 0) {
 113         return(-1);
 114     };
 115
 116     /* Attribut "UNI" */
 117     if (_MEDattrStringEcrire(ds, MED_NOM_UNI, mdim*MED_TAILLE_PNOM, unicoo, mode) < 0) {
 118         return(-1);
 119     };
 120
 121     /* On ferme tout */
 122     if (_MEDdatasetFermer(ds)      < 0) return(-1);
 123     if (_MEDdatagroupFermer(noeid) < 0) return(-1);
 124     if (_MEDdatagroupFermer(maaid) < 0) return(-1);
 125     return(0);
 126 }
 127
 128 }