src/MEDWrapper/V2_1/MEDnVal.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.hxx"
  20 #include <string.h>
  21 #include <stdlib.h>
  22 #include "med_outils.hxx"
  23
  24 namespace med_2_1{
  25
  26 med_int
  27 MEDnVal(med_idt fid, char *champ, med_entite_maillage type_ent,
  28         med_geometrie_element type_geo,med_int numdt, med_int numo)
  29 {
  30   med_int n;
  31   med_idt datagroup;
  32   med_err ret;
  33   char nomdatagroup1[2*MED_TAILLE_NOM_ENTITE+2],nomdatagroup2[2*MED_MAX_PARA+1];
  34   char tmp1   [MED_TAILLE_NOM_ENTITE+1];
  35   char chemin [MED_TAILLE_CHA+(MED_TAILLE_NOM+1)+(2*MED_TAILLE_NOM_ENTITE+2)+(2*MED_MAX_PARA)+1+100];
  36
  37   /*
  38    * On inhibe le gestionnaire d'erreur HDF
  39    */
  40   _MEDmodeErreurVerrouiller();
  41
  42   /*
  43    * On cree le chemin d'accès
  44    */
  45   strcpy(chemin,MED_CHA);
  46   strcat(chemin,champ);
  47   strcat(chemin,"/");
  48
  49   /* On cree le nom du datagroup de niveau 1 */
  50   if ((ret = _MEDnomEntite(nomdatagroup1,type_ent)) < 0)
  51     return -1;
  52   if ((type_ent != MED_NOEUD))
  53     {
  54       if ((ret = _MEDnomGeometrie(tmp1,type_geo)) < 0)
  55         return -1;
  56       strcat(nomdatagroup1,".");
  57       strcat(nomdatagroup1,tmp1);
  58     }
  59   strcat(chemin,nomdatagroup1);
  60   strcat(chemin,"/");
  61
  62   /* Creation du datagroup de niveau 2 <numdt>.<numoo> */
  63   sprintf(nomdatagroup2,"%*li%*li",MED_MAX_PARA,(long ) numdt,MED_MAX_PARA,(long ) numo);
  64   strcat(chemin,nomdatagroup2);
  65
  66  /*
  67    * Acces au champ
  68   */
  69   if ((datagroup = _MEDdatagroupOuvrir(fid,chemin)) < 0)
  70     return 0;
  71   if ((ret = _MEDattrEntierLire(datagroup,MED_NOM_NBR,&n)) < 0)
  72     return -1;
  73
  74   /*
  75    * fermetures
  76    */
  77   if ((ret = _MEDdatagroupFermer(datagroup)) < 0)
  78     return -1;
  79
  80   return n;
  81 }
  82
  83 }