# _MEDmeshGridIndexCoordinateRd236.c

Aller à la documentation de ce fichier.

    00001 /*  This file is part of MED.
    00002  *
    00003  *  COPYRIGHT (C) 1999 - 2011  EDF R&D, CEA/DEN
    00004  *  MED is free software: you can redistribute it and/or modify
    00005  *  it under the terms of the GNU Lesser General Public License as published by
    00006  *  the Free Software Foundation, either version 3 of the License, or
    00007  *  (at your option) any later version.
    00008  *
    00009  *  MED is distributed in the hope that it will be useful,
    00010  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
    00011  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    00012  *  GNU Lesser General Public License for more details.
    00013  *
    00014  *  You should have received a copy of the GNU Lesser General Public License
    00015  *  along with MED.  If not, see <http://www.gnu.org/licenses/>.
    00016  */
    00017
    00018
    00019 #include <med.h>
    00020 #include <med_config.h>
    00021 #include <med_outils.h>
    00022 #include <string.h>
    00023 #include <stdlib.h>
    00024
    00025 void _MEDmeshGridIndexCoordinateRd236(int dummy, ...) {
    00026
    00027
    00028   med_access_mode _MED_ACCESS_MODE;
    00029   med_idt         _meshid=0;
    00030   med_err         _ret         = -1;
    00031   med_data_type   _datatype    = MED_UNDEF_DATATYPE;
    00032   med_grid_type   _gridtype    = MED_UNDEF_GRID_TYPE;
    00033   med_int         _intgridtype = 0;
    00034   med_int         _meshdim     = 0;
    00035   char            _meshpath[MED_TAILLE_MAA+MED_NAME_SIZE+1]=MED_MAA;
    00036   med_int  _intaxistype=0;
    00037
    00038
    00039   MED_VARGS_DECL(const, med_idt               , , fid       );
    00040   MED_VARGS_DECL(const, char*  , const          , meshname  );
    00041   MED_VARGS_DECL(const, med_int               , , numdt     );
    00042   MED_VARGS_DECL(const, med_int               , , numit     );
    00043   MED_VARGS_DECL(const, med_int               , , axis      );
    00044   MED_VARGS_DECL(, med_float *, const           , gridindex );
    00045   MED_VARGS_DECL(, med_err *                   ,, fret      );
    00046
    00047   va_list params;
    00048   va_start(params,dummy);
    00049
    00050   MED_VARGS_DEF(const, med_idt               , , fid       );
    00051   MED_VARGS_DEF(const, char*  , const          , meshname  );
    00052   MED_VARGS_DEF(const, med_int               , , numdt     );
    00053   MED_VARGS_DEF(const, med_int               , , numit     );
    00054   MED_VARGS_DEF(const, med_int               , , axis      );
    00055   MED_VARGS_DEF(, med_float *, const           , gridindex );
    00056   MED_VARGS_DEF(, med_err *                   ,, fret      );
    00057
    00058   /*
    00059    * On inhibe le gestionnaire d'erreur HDF 5
    00060    */
    00061   _MEDmodeErreurVerrouiller();
    00062
    00063   strcat(_meshpath,meshname);
    00064   if ((_meshid = _MEDdatagroupOuvrir(fid,_meshpath)) < 0) {
    00065     MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,MED_ERR_MESH_MSG);
    00066     SSCRUTE(meshname);SSCRUTE(_meshpath); goto ERROR;
    00067   }
    00068
    00069   /*
    00070    * Lecture du type de grille (attribut MED_NOM_GTY)
    00071    */
    00072   if (_MEDattrEntierLire(_meshid,MED_NOM_GTY,&_intgridtype) < 0) {
    00073     MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
    00074     SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE(_intgridtype);goto ERROR;
    00075   }
    00076   _gridtype=(med_grid_type) _intgridtype;
    00077
    00078   if ((_gridtype != MED_CARTESIAN_GRID) && (_gridtype != MED_POLAR_GRID)) {
    00079     MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
    00080     SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
    00081   }
    00082
    00083
    00084  /* Lecture de l'attribut MED_NOM_REP */
    00085   /* En 2.3.6, il n'y en a pas pour les grilles !=  MED_GRILLE_STANDARD*/
    00086 /*   if ( _MEDattrEntierLire(_meshid,MED_NOM_REP,&_intaxistype) < 0) { */
    00087 /*     MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG); */
    00088 /*     SSCRUTE(meshname);SSCRUTE(MED_NOM_REP); */
    00089 /*     SSCRUTE(_intaxistype);goto ERROR; */
    00090 /*   } */
    00091
    00092 /*   if ((med_mesh_type)( _intaxistype) != MED_CARTESIAN ) { */
    00093 /*     MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_AXISTYPE,MED_ERR_MESH_MSG); */
    00094 /*     SSCRUTE(meshname);ISCRUTE_int(_intaxistype);goto ERROR; */
    00095 /*   } */
    00096
    00097   /* Lecture de l'attribut MED_NOM_DIM  */
    00098   if (_MEDattrEntierLire(_meshid,MED_NOM_DIM,&_meshdim) < 0) {
    00099     MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
    00100     SSCRUTE(meshname);SSCRUTE(MED_NOM_DIM);ISCRUTE(_meshdim);goto ERROR;
    00101   }
    00102
    00103   if (axis > _meshdim ) {
    00104     MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
    00105     SSCRUTE(meshname);ISCRUTE(_meshdim);ISCRUTE(axis);goto ERROR;
    00106   }
    00107
    00108   /*
    00109    * On ecrit le tableau d'indice dans un dataset HDF
    00110    */
    00111   switch(axis)
    00112     {
    00113     case 1 :
    00114       _datatype = MED_COORDINATE_AXIS1;
    00115       break;
    00116     case 2 :
    00117       _datatype = MED_COORDINATE_AXIS2;
    00118       break;
    00119     case 3 :
    00120       _datatype = MED_COORDINATE_AXIS3;
    00121       break;
    00122     default :
    00123       MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
    00124       SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
    00125     }
    00126
    00127   if (_MEDmeshAdvancedRd(fid,
    00128                          meshname,
    00129                          _datatype,
    00130                          MED_NO_NAME,
    00131                          MED_INTERNAL_UNDEF,
    00132                          numdt,
    00133                          numit,
    00134                          MED_NODE,
    00135                          MED_NONE,
    00136                          MED_NO_CMODE,
    00137                          MED_UNDEF_PFLMODE,
    00138                          MED_NO_PROFILE,
    00139                          MED_FULL_INTERLACE,
    00140                          MED_ALL_CONSTITUENT,
    00141                          NULL,
    00142                          (unsigned char * const) gridindex) < 0 ) {
    00143     MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAdvancedRd");
    00144     goto ERROR;
    00145   }
    00146
    00147   _ret = 0;
    00148
    00149  ERROR:
    00150
    00151   if (_meshid>0)            if (_MEDdatagroupFermer(_meshid) < 0) {
    00152     MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_meshpath);
    00153     ISCRUTE_id(_meshid);
    00154   }
    00155
    00156
    00157   va_end(params);
    00158   *fret = _ret;
    00159   return;
    00160 }