# MEDmeshGridStructWr.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
    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 }