00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019 #include <med.h>
00020 #include <med_config.h>
00021 #include <med_outils.h>
00022
00023 #include <string.h>
00024 #include <stdlib.h>
00025
00026 med_err
00027 MEDconnEcr(med_idt fid,char *maa, med_int mdim, med_int *connectivite,med_mode_switch mode_switch,
00028 med_int nbre, med_entite_maillage type_ent,
00029 med_geometrie_element type_geo,med_connectivite type_conn)
00030 {
00031 med_idt maaid, entid, geoid, dataset;
00032 med_err ret;
00033 med_size dimd[1];
00034 char chemin[MED_TAILLE_MAA+MED_TAILLE_NOM+1];
00035 char nom_ent[MED_TAILLE_NOM_ENTITE+1];
00036 char nom_geo[MED_TAILLE_NOM_ENTITE+1];
00037 char nom_dataset[MED_TAILLE_NOM_ENTITE+1];
00038 int dim, nnoe, ndes;
00039 int taille;
00040 med_entite_maillage _type_ent=type_ent;
00041
00042 if ( type_ent == MED_NOEUD_MAILLE ) _type_ent=MED_NOEUD ;
00043
00044
00045
00046
00047 _MEDmodeErreurVerrouiller();
00048 if (MEDcheckVersion(fid) < 0) return -1;
00049
00050
00051
00052
00053
00054 strcpy(chemin,MED_MAA);
00055 strcat(chemin,maa);
00056 if ((maaid = _MEDdatagroupOuvrir(fid,chemin)) < 0)
00057 return -1;
00058
00059
00060
00061
00062
00063 if ((ret = _MEDnomEntite(nom_ent,_type_ent)) < 0)
00064 return -1;
00065
00066
00067
00068 if ((entid = _MEDdatagroupOuvrir(maaid,nom_ent)) < 0)
00069 if ((entid = _MEDdatagroupCreer(maaid,nom_ent)) < 0)
00070 return -1;
00071
00072
00073
00074
00075 if ((ret = _MEDnomGeometrie(nom_geo,type_geo)) < 0)
00076 return -1;
00077 if ((geoid = _MEDdatagroupOuvrir(entid,nom_geo)) < 0)
00078 if ((geoid = _MEDdatagroupCreer(entid,nom_geo)) < 0)
00079 return -1;
00080
00081
00082
00083
00084 if ((ret=_MEDparametresGeometrie(_type_ent,type_geo,&dim,&nnoe,&ndes))<0)
00085 return -1;
00086
00087 switch(type_conn)
00088 {
00089 case MED_NOD :
00090 strcpy(nom_dataset,MED_NOM_NOD);
00091 taille = nnoe;
00092 break;
00093
00094 case MED_DESC :
00095 strcpy(nom_dataset,MED_NOM_DES);
00096 taille = ndes;
00097 break;
00098
00099 default :
00100 return -1;
00101 }
00102 dimd[0] = nbre*taille;
00103 #if defined(HAVE_F77INT64)
00104 ret = _MEDdatasetNumEcrire(geoid,nom_dataset,MED_INT64,mode_switch,(med_size)taille,MED_ALL,
00105 MED_NOPF,MED_NO_PFLMOD,0,0,MED_NOPG,dimd,
00106 (unsigned char*) connectivite);
00107 #else
00108 ret = _MEDdatasetNumEcrire(geoid,nom_dataset,MED_INT32,mode_switch,(med_size)taille,MED_ALL,
00109 MED_NOPF,MED_NO_PFLMOD,0,0,MED_NOPG,dimd,
00110 (unsigned char*) connectivite);
00111 #endif
00112 if (ret < 0 ) {
00113 MESSAGE("Impossible d'ecrire le dataset coon de taille : ");
00114 ISCRUTE_int(taille); return -1;
00115 }
00116
00117
00118
00119
00120 if ((dataset = _MEDdatasetOuvrir(geoid,nom_dataset)) < 0)
00121 return -1;
00122 if ((ret = _MEDattrEntierEcrire(dataset,MED_NOM_NBR,&nbre)) < 0)
00123 return -1;
00124
00125
00126
00127
00128 if ((ret = _MEDdatasetFermer(dataset)) < 0)
00129 return -1;
00130 if ((ret = _MEDdatagroupFermer(geoid)) < 0)
00131 return -1;
00132 if ((ret = _MEDdatagroupFermer(entid)) < 0)
00133 return -1;
00134 if ((ret = _MEDdatagroupFermer(maaid)) < 0)
00135 return -1;
00136
00137 return 0;
00138 }
00139
00140