# MEDindicesCoordEcr.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
    00023 #include <stdlib.h>
    00024 #include <string.h>
    00025
    00026 med_err
    00027 MEDindicesCoordEcr(med_idt fid,char *maillage,med_int mdim,med_float *indices,
    00028                    med_int n,med_int axe,char *comp,char *unit)
    00029 {
    00030   med_idt maaid, noeid, dataset;
    00031   med_err ret;
    00032   med_size dimd[1];
    00033   char chemin[MED_TAILLE_MAA+MED_TAILLE_NOM+1];
    00034   char nom_dataset[MED_TAILLE_NOM_ENTITE+1];
    00035   med_int att;
    00036   med_maillage maillage_type;
    00037   med_type_grille type;
    00038
    00039   /*
    00040    * Si axe > mdim => erreur
    00041    */
    00042   if (axe > mdim)
    00043     return -1;
    00044
    00045   /*
    00046    * On inhibe le gestionnaire d'erreur HDF
    00047    */
    00048   _MEDmodeErreurVerrouiller();
    00049 if (MEDcheckVersion(fid) < 0) return -1;
    00050
    00051
    00052   /*
    00053    * Si le maillage n'existe pas => erreur
    00054    */
    00055   strcpy(chemin,MED_MAA);
    00056   strcat(chemin,maillage);
    00057   if ((maaid = _MEDdatagroupOuvrir(fid,chemin)) < 0)
    00058       return -1;
    00059
    00060   /*
    00061    * Si le maillage est de type MED_NON_STRUCTURE => erreur
    00062    */
    00063   if ((ret = _MEDattrEntierLire(maaid,MED_NOM_TYP,&att)) < 0)
    00064     return -1;
    00065   maillage_type = (med_maillage) att;
    00066   if (maillage_type == MED_NON_STRUCTURE)
    00067     return -1;
    00068
    00069   /*
    00070    * Si la grille n'est pas de type MED_GRILLE_CARTESIENNE ou
    00071    * MED_GRILLE_POLAIRE => erreur
    00072    */
    00073   if ((ret = _MEDattrEntierLire(maaid,MED_NOM_GTY,&att)) < 0)
    00074     return -1;
    00075   /* sizeof(enum) tjrs = sizeof(int) en C, or
    00076      sur machines 64 bits par défaut med_int==long (normalement 64bits),
    00077      du coup sur  machines 64 bits _MEDattrEntierLire utilise
    00078      le type hdf NATIVE_LONG, ce genere parfois un warning
    00079   */
    00080   type = (med_type_grille) att;
    00081   if ((type != MED_GRILLE_CARTESIENNE) && (type != MED_GRILLE_POLAIRE))
    00082     return -1;
    00083
    00084   /*
    00085    * Si le groupe HDF "NOE" n'existe pas => erreur
    00086    */
    00087   if ((noeid = _MEDdatagroupOuvrir(maaid,MED_NOM_NOE)) < 0)
    00088     if ((noeid = _MEDdatagroupCreer(maaid,MED_NOM_NOE)) < 0)
    00089       return -1;
    00090
    00091   /*
    00092    * On ecrit le tableau d'indice dans un dataset HDF
    00093    */
    00094   switch(axe) {
    00095
    00096   case 1 :
    00097     strcpy(nom_dataset,MED_NOM_IN1);
    00098     break;
    00099
    00100   case 2 :
    00101     strcpy(nom_dataset,MED_NOM_IN2);
    00102     break;
    00103
    00104   case 3 :
    00105     strcpy(nom_dataset,MED_NOM_IN3);
    00106     break;
    00107
    00108   default :
    00109     return -1;
    00110
    00111   }
    00112
    00113   dimd[0] = n;
    00114   if ((ret = _MEDdatasetNumEcrire(noeid,nom_dataset,MED_FLOAT64,MED_FULL_INTERLACE,1,MED_ALL,MED_NOPF,MED_NO_PFLMOD,0,0,
    00115                                   MED_NOPG,dimd,(unsigned char*) indices)) < 0)
    00116     return -1;
    00117
    00118   /*
    00119    * On re-ouvre le dataset HDF pour y placer des attributs
    00120    */
    00121   if ((dataset = _MEDdatasetOuvrir(noeid,nom_dataset)) < 0)
    00122     return -1;
    00123
    00124   /*
    00125    * Attribut NBR (taille du tableau d'indices)
    00126    */
    00127   if ((ret = _MEDattrEntierEcrire(dataset,MED_NOM_NBR,&n)) < 0)
    00128     return -1;
    00129
    00130   /*
    00131    * Attribut "NOM" (nom de la composante)
    00132    */
    00133   if ((ret = _MEDattrStringEcrire(dataset,MED_NOM_NOM,MED_TAILLE_PNOM,comp)) < 0)
    00134     return -1;
    00135
    00136   /*
    00137    * Attribut "UNI" (unite de la composante)
    00138    */
    00139   if ((ret = _MEDattrStringEcrire(dataset,MED_NOM_UNI,MED_TAILLE_PNOM,unit)) < 0)
    00140     return -1;
    00141
    00142   /*
    00143    * On ferme tout
    00144    */
    00145   if ((ret = _MEDdatasetFermer(dataset)) < 0)
    00146     return -1;
    00147   if ((ret = _MEDdatagroupFermer(noeid)) < 0)
    00148     return -1;
    00149   if ((ret = _MEDdatagroupFermer(maaid)) < 0)
    00150     return -1;
    00151
    00152   return 0;
    00153 }