# _MEDmeshGridIndexCoordinateRd30.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 _MEDmeshGridIndexCoordinateRd30(int dummy, ...) {
    00026
    00027
    00028   med_idt         _meshid=0;
    00029   med_err         _ret         = -1;
    00030   med_data_type   _datatype    = MED_UNDEF_DATATYPE;
    00031   med_grid_type   _gridtype    = MED_UNDEF_GRID_TYPE;
    00032   med_int         _intgridtype = 0;
    00033   med_int         _meshdim     = 0;
    00034   char            _meshpath[MED_TAILLE_MAA+MED_NAME_SIZE+1]=MED_MAA;
    00035   med_int  _intaxistype=0;
    00036
    00037
    00038   MED_VARGS_DECL(const, med_idt               , , fid       );
    00039   MED_VARGS_DECL(const, char*  , const          , meshname  );
    00040   MED_VARGS_DECL(const, med_int               , , numdt     );
    00041   MED_VARGS_DECL(const, med_int               , , numit     );
    00042   MED_VARGS_DECL(const, med_int               , , axis      );
    00043   MED_VARGS_DECL(, med_float *, const           , gridindex );
    00044   MED_VARGS_DECL(, med_err *                   ,, fret      );
    00045
    00046   va_list params;
    00047   va_start(params,dummy);
    00048
    00049   MED_VARGS_DEF(const, med_idt               , , fid       );
    00050   MED_VARGS_DEF(const, char*  , const          , meshname  );
    00051   MED_VARGS_DEF(const, med_int               , , numdt     );
    00052   MED_VARGS_DEF(const, med_int               , , numit     );
    00053   MED_VARGS_DEF(const, med_int               , , axis      );
    00054   MED_VARGS_DEF(, med_float *, const           , gridindex );
    00055   MED_VARGS_DEF(, med_err *                   ,, fret      );
    00056
    00057   /*
    00058    * On inhibe le gestionnaire d'erreur HDF 5
    00059    */
    00060   _MEDmodeErreurVerrouiller();
    00061
    00062   strcat(_meshpath,meshname);
    00063   if ((_meshid = _MEDdatagroupOuvrir(fid,_meshpath)) < 0) {
    00064     MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,MED_ERR_MESH_MSG);
    00065     SSCRUTE(meshname);SSCRUTE(_meshpath); goto ERROR;
    00066   }
    00067
    00068   /*
    00069    * Lecture du type de grille (attribut MED_NOM_GTY)
    00070    */
    00071   if (_MEDattrEntierLire(_meshid,MED_NOM_GTY,&_intgridtype) < 0) {
    00072     MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
    00073     SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE(_intgridtype);goto ERROR;
    00074   }
    00075   _gridtype=(med_grid_type) _intgridtype;
    00076
    00077   if ((_gridtype != MED_CARTESIAN_GRID) && (_gridtype != MED_POLAR_GRID)) {
    00078     MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
    00079     SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
    00080   }
    00081
    00082
    00083  /* Lecture de l'attribut MED_NOM_REP */
    00084   if ( _MEDattrEntierLire(_meshid,MED_NOM_REP,&_intaxistype) < 0) {
    00085     MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
    00086     SSCRUTE(meshname);SSCRUTE(MED_NOM_REP);
    00087     ISCRUTE(_intaxistype);goto ERROR;
    00088   }
    00089
    00090   if ((med_mesh_type)( _intaxistype) != MED_CARTESIAN ) {
    00091     MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_AXISTYPE,MED_ERR_MESH_MSG);
    00092     SSCRUTE(meshname);ISCRUTE(_intaxistype);goto ERROR;
    00093   }
    00094
    00095   /* Lecture de l'attribut MED_NOM_DIM  */
    00096   if (_MEDattrEntierLire(_meshid,MED_NOM_DIM,&_meshdim) < 0) {
    00097     MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
    00098     SSCRUTE(meshname);SSCRUTE(MED_NOM_DIM);ISCRUTE(_meshdim);goto ERROR;
    00099   }
    00100
    00101   if (axis > _meshdim ) {
    00102     MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
    00103     SSCRUTE(meshname);ISCRUTE(_meshdim);ISCRUTE(axis);goto ERROR;
    00104   }
    00105
    00106   /*
    00107    * On ecrit le tableau d'indice dans un dataset HDF
    00108    */
    00109   switch(axis)
    00110     {
    00111     case 1 :
    00112       _datatype = MED_COORDINATE_AXIS1;
    00113       break;
    00114     case 2 :
    00115       _datatype = MED_COORDINATE_AXIS2;
    00116       break;
    00117     case 3 :
    00118       _datatype = MED_COORDINATE_AXIS3;
    00119       break;
    00120     default :
    00121       MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
    00122       SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
    00123     }
    00124
    00125   if (_MEDmeshAdvancedRd(fid,
    00126                          meshname,
    00127                          _datatype,
    00128                          MED_NO_NAME,
    00129                          MED_INTERNAL_UNDEF,
    00130                          numdt,
    00131                          numit,
    00132                          MED_NODE,
    00133                          MED_NONE,
    00134                          MED_NO_CMODE,
    00135                          MED_UNDEF_PFLMODE,
    00136                          MED_NO_PROFILE,
    00137                          MED_FULL_INTERLACE,
    00138                          MED_ALL_CONSTITUENT,
    00139                          NULL,
    00140                          (unsigned char *) gridindex) < 0 ) {
    00141     MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAdvancedRd");
    00142     goto ERROR;
    00143   }
    00144
    00145   _ret = 0;
    00146
    00147  ERROR:
    00148
    00149   if (_meshid>0)            if (_MEDdatagroupFermer(_meshid) < 0) {
    00150     MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_meshpath);
    00151     ISCRUTE_id(_meshid);
    00152   }
    00153
    00154
    00155   va_end(params);
    00156   *fret = _ret;
    00157   return;
    00158 }