00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018 #include <med.h>
00019 #include <med_config.h>
00020 #include <med_outils.h>
00021
00022 #include <string.h>
00023 #include <stdlib.h>
00024
00025 med_err
00026 MEDgaussLire(med_idt fid, med_float *refcoo, med_float *gscoo, med_float * wg, med_mode_switch mode_coo, char *locname)
00027 {
00028 med_err ret = -1;
00029 med_idt locid=0;
00030 char chemin[MED_TAILLE_GAUSS+MED_TAILLE_NOM+1]="";
00031 med_int type_geo;
00032
00033
00034
00035
00036 _MEDmodeErreurVerrouiller();
00037 if (MEDcheckVersion(fid) < 0) return -1;
00038
00039
00040
00041
00042
00043 strcpy(chemin,MED_GAUSS);
00044 strcat(chemin,locname);
00045 if ((locid = _MEDdatagroupOuvrir(fid,chemin)) < 0) {
00046 MESSAGE("Impossible d'ouvrir le datagroup : ");
00047 SSCRUTE(chemin); goto ERROR;
00048 }
00049
00050
00051
00052 if (_MEDattrEntierLire(locid,MED_NOM_GEO,&type_geo) < 0) {
00053 MESSAGE("Erreur à la lecture de l'attribut MED_NOM_GEO : ");
00054 ISCRUTE(type_geo);goto ERROR;
00055 };
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065 if ( _MEDdatasetNumLire(locid,MED_NOM_COO,MED_FLOAT64,mode_coo,(type_geo/100),MED_ALL,MED_NOPF,MED_NO_PFLMOD,MED_PFL_NON_COMPACT,0,MED_NOPG,0,
00066 (unsigned char*) refcoo) < 0 ) {
00067 MESSAGE("Impossible de lire le dataset : ");SSCRUTE(MED_NOM_COO);
00068 goto ERROR;
00069 }
00070
00071
00072
00073
00074
00075 if ( _MEDdatasetNumLire(locid,MED_NOM_GAU,MED_FLOAT64,mode_coo,(type_geo/100),MED_ALL,MED_NOPF,MED_NO_PFLMOD,MED_PFL_NON_COMPACT,0,MED_NOPG,0,
00076 (unsigned char*) gscoo) < 0 ) {
00077 MESSAGE("Impossible de lire le dataset : ");SSCRUTE(MED_NOM_GAU);
00078 goto ERROR;
00079 }
00080
00081
00082
00083
00084
00085 if ( _MEDdatasetNumLire(locid,MED_NOM_VAL,MED_FLOAT64,mode_coo,1,MED_ALL,MED_NOPF,MED_NO_PFLMOD,MED_PFL_NON_COMPACT,0,MED_NOPG,0,
00086 (unsigned char*) wg) < 0 ) {
00087 MESSAGE("Impossible de lire le dataset : ");SSCRUTE(MED_NOM_VAL);
00088 goto ERROR;
00089 }
00090
00091
00092
00093
00094
00095
00096 ret = 0;
00097 ERROR:
00098
00099 if ( locid > 0 ) if ( _MEDdatagroupFermer(locid) < 0) {
00100 MESSAGE("Impossible de fermer le datagroup : ");
00101 ISCRUTE_id(locid); ret = -1;
00102 }
00103
00104 return ret;
00105 }
00106