00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027 #include <med.h>
00028 #define MESGERR 1
00029 #include "med_utils.h"
00030 #include <string.h>
00031
00032 #ifdef DEF_LECT_ECR
00033 #define MODE_ACCES MED_ACC_RDWR
00034 #elif DEF_LECT_AJOUT
00035 #define MODE_ACCES MED_ACC_RDEXT
00036 #else
00037 #define MODE_ACCES MED_ACC_CREAT
00038 #endif
00039
00040 int main (int argc, char **argv)
00041
00042
00043 {
00044 med_idt fid=0;
00045 med_int mdim=2,axe=0,nind=0;
00046 med_float indiceX[4] = {1.0,1.1,1.2,1.3};
00047 med_float indiceY[4] = {2.0,2.1,2.2,2.3};
00048 med_float coo[8] = {0.0,0.0,1.0,0.0,0.0,1.0,1.0,1.0};
00049 med_int nnoeuds = 4;
00050 med_int structure_grille[2] = {2,2};
00051 char maa [MED_NAME_SIZE+1]= "grille_cartesian";
00052 char maa2[MED_NAME_SIZE+1]= "grille_curvilinear";
00053
00054
00055 char comp[2*MED_SNAME_SIZE+1] = "X Y ";
00056 char unit[2*MED_SNAME_SIZE+1] = "cm cm ";
00057
00058
00059 fid = MEDfileOpen("test27.med",MODE_ACCES);
00060 if (fid < 0) {
00061 MESSAGE("Erreur a la creation du fichier test27.med");
00062 return -1;
00063 }
00064 printf("Creation du fichier test27.med \n");
00065
00066
00067 if (MEDmeshCr( fid, "maillage vide",2, 2, MED_UNSTRUCTURED_MESH,
00068 "un maillage vide","s", MED_SORT_DTIT,
00069 MED_CARTESIAN, comp, unit) < 0) {
00070 MESSAGE("Erreur a la creation du maillage MED_UNSTRUCTURED_MESH : "); SSCRUTE(maa);
00071 return -1;
00072 }
00073
00074
00075
00076 if (MEDmeshCr( fid, maa,mdim, mdim, MED_STRUCTURED_MESH,
00077 "un exemple de grille cartesienne","s", MED_SORT_DTIT,
00078 MED_CARTESIAN, comp, unit) < 0) {
00079 MESSAGE("Erreur a la creation de la grille");
00080 return -1;
00081 }
00082 printf("Creation d'un maillage structure MED_STRUCTURED_MESH \n");
00083
00084
00085 if (MEDmeshGridTypeWr(fid,maa, MED_CARTESIAN_GRID) < 0) {
00086 MESSAGE("Erreur a l'ecriture de la nature de la grille");
00087 return -1;
00088 }
00089
00090 printf("On definit la nature du maillage structure : MED_GRILLE_CARTESIENNE \n");
00091
00092
00093
00094 nind = 4;
00095 axe = 1;
00096 if (MEDmeshGridIndexCoordinateWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,
00097 axe,nind,indiceX) < 0) {
00098 MESSAGE("Erreur a l'ecriture de l'axe X");
00099 return -1;
00100 }
00101 printf("Ecriture des indices des coordonnees selon l'axe des X \n");
00102
00103
00104 nind = 4;
00105 axe = 2;
00106 if (MEDmeshGridIndexCoordinateWr(fid,maa,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,
00107 axe,nind,indiceY) < 0) {
00108 MESSAGE("Erreur a l'ecriture de l'axe Y");
00109 return -1;
00110 }
00111 printf("Ecriture des indices des coordonnees selon l'axe des Y \n");
00112
00113
00114
00115
00116 if (MEDmeshCr( fid, maa2,mdim, mdim, MED_STRUCTURED_MESH,
00117 "un exemple de grille standard","s", MED_SORT_DTIT,
00118 MED_CARTESIAN, comp, unit) < 0) {
00119 MESSAGE("Erreur a la creation de la 2e grille");
00120 return -1;
00121 }
00122 printf("Creation d'un maillage structure MED_STRUCTURED_MESH \n");
00123
00124
00125 if (MEDmeshGridTypeWr(fid,maa2, MED_CURVILINEAR_GRID) < 0) {
00126 MESSAGE("Erreur a l'ecriture de la nature de la grille");
00127 return -1;
00128 }
00129 printf("On definit la nature du maillage structure : MED_CURVILINEAR_GRID \n");
00130
00131
00132 if (MEDmeshNodeCoordinateWr(fid,maa2,MED_NO_DT,MED_NO_IT,MED_UNDEF_DT,
00133 MED_FULL_INTERLACE,nnoeuds, coo) < 0) {
00134 MESSAGE("Erreur a l'ecriture des noeuds de la grille MED_CURVILINEAR_GRID");
00135 return -1;
00136 }
00137 printf("Ecriture des coordonnees des noeuds \n");
00138
00139
00140 if ( MEDmeshGridStructWr(fid,maa2,MED_NO_DT,MED_NO_IT, MED_UNDEF_DT, structure_grille ) < 0) {
00141 MESSAGE("Erreur a l'ecriture de la structure de la grille");
00142 return -1;
00143 }
00144 printf("Ecriture de la structure de la grille : / 2,2 / \n");
00145
00146
00147 if (MEDfileClose(fid) < 0) {
00148 MESSAGE("Erreur a la fermeture du fichier");
00149 return -1;
00150 }
00151 printf("Fermeture du fichier \n");
00152
00153 return 0;
00154 }