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