MEDmeshGridStructWr.c

Aller à la documentation de ce fichier.
00001 /*  This file is part of MED.
00002  *
00003  *  COPYRIGHT (C) 1999 - 2015  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 
00037 med_err
00038 MEDmeshGridStructWr(const med_idt               fid,
00039                     const char*  const          meshname,
00040                     const med_int               numdt,
00041                     const med_int               numit,
00042                     const med_float             dt,
00043                     const med_int * const       gridstruct)
00044 {
00045   med_access_mode _MED_ACCESS_MODE;
00046   med_idt         _meshid=0;
00047   med_err         _ret         = -1;
00048   med_data_type   _datatype    = MED_UNDEF_DATATYPE;
00049   med_grid_type   _gridtype    = MED_UNDEF_GRID_TYPE;
00050   med_int         _intgridtype = 0;
00051   med_int         _meshdim     = 0;
00052   char            _meshpath[MED_TAILLE_MAA+MED_NAME_SIZE+1]=MED_MAA;
00053   int             _i=0;
00054 
00055   /*
00056    * On inhibe le gestionnaire d'erreur HDF 5
00057    */
00058   _MEDmodeErreurVerrouiller();
00059  if (_MEDcheckVersion30(fid) < 0) goto ERROR;
00060 
00061   if ( (_MED_ACCESS_MODE = _MEDmodeAcces(fid) ) == MED_ACC_UNDEF ) {
00062     MED_ERR_(_ret,MED_ERR_UNRECOGNIZED,MED_ERR_ACCESSMODE,MED_ERR_FILE_MSG);
00063     goto ERROR;
00064   }
00065 
00066   if ( _MED_ACCESS_MODE == MED_ACC_RDONLY) {
00067     MED_ERR_(_ret,MED_ERR_INVALID,MED_ERR_ACCESSMODE,MED_ERR_FILE_MSG);
00068     ISCRUTE_int(_MED_ACCESS_MODE);
00069     goto ERROR;
00070   }
00071 
00072   strcat(_meshpath,meshname);
00073   if ((_meshid = _MEDdatagroupOuvrir(fid,_meshpath)) < 0) {
00074     MED_ERR_(_ret,MED_ERR_OPEN,MED_ERR_DATAGROUP,MED_ERR_MESH_MSG);
00075     SSCRUTE(meshname);SSCRUTE(_meshpath); goto ERROR;
00076   }
00077 
00078   /* Lecture de l'attribut MED_NOM_DIM  */
00079   if (_MEDattrEntierLire(_meshid,MED_NOM_DIM,&_meshdim) < 0) {
00080     MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
00081     SSCRUTE(meshname);SSCRUTE(MED_NOM_DIM);ISCRUTE(_meshdim);goto ERROR;
00082   }
00083 
00084   /*
00085    * Lecture du type de grille (attribut MED_NOM_GTY)
00086    */
00087   if (_MEDattrEntierLire(_meshid,MED_NOM_GTY,&_intgridtype) < 0) {
00088     MED_ERR_(_ret,MED_ERR_READ,MED_ERR_ATTRIBUTE,MED_ERR_MESH_MSG);
00089     SSCRUTE(meshname);SSCRUTE(MED_NOM_GTY);ISCRUTE_int(_gridtype);goto ERROR;
00090   }
00091   _gridtype=(med_grid_type) _intgridtype;
00092 
00093   if ((_gridtype != MED_CURVILINEAR_GRID)) {
00094     MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
00095     SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
00096   }
00097 
00098   /*
00099    * Attribut NBR (taille du tableau d'indices)
00100    */
00101   for (_i=0;_i<_meshdim;++_i) {
00102 
00103     switch(_i) {
00104      case 0 :
00105       _datatype = MED_COORDINATE_AXIS1;
00106       break;
00107     case 1 :
00108       _datatype = MED_COORDINATE_AXIS2;
00109       break;
00110     case 2 :
00111       _datatype = MED_COORDINATE_AXIS3;
00112       break;
00113     default :
00114       MED_ERR_(_ret,MED_ERR_RANGE,MED_ERR_GRIDTYPE,MED_ERR_MESH_MSG);
00115       SSCRUTE(meshname);ISCRUTE_int(_gridtype);goto ERROR;
00116     }
00117 
00118     if (_MEDmeshAdvancedWr(fid,
00119                            meshname,
00120                            _datatype,
00121                            MED_NO_NAME,
00122                            MED_INTERNAL_INT,
00123                            numdt,
00124                            numit,
00125                            dt,
00126                            MED_NODE,
00127                            MED_NONE,
00128                            MED_NO_CMODE,
00129                            MED_UNDEF_PFLMODE,
00130                            MED_NO_PROFILE,
00131                            MED_FULL_INTERLACE,
00132                            MED_ALL_CONSTITUENT,
00133                            NULL,
00134                            1,
00135                            (gridstruct+_i)) < 0 ) {
00136       MED_ERR_(_ret,MED_ERR_CALL,MED_ERR_API,"_MEDmeshAdvancedWr");
00137       goto ERROR;
00138     }
00139 
00140   }
00141 
00142  _ret = 0;
00143 
00144  ERROR:
00145 
00146   if (_meshid>0)            if (_MEDdatagroupFermer(_meshid) < 0) {
00147     MED_ERR_(_ret,MED_ERR_CLOSE,MED_ERR_DATAGROUP,_meshpath);
00148     ISCRUTE_id(_meshid);
00149   }
00150 
00151   return _ret;
00152 
00153 }

Généré le Thu Oct 8 14:26:16 2015 pour MED fichier par  doxygen 1.6.1