1 /************************************************************************* 2 * COPYRIGHT (C) 1999 - 2007 EDF R&D, CEA/DEN 3 * THIS LIBRARY IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY 4 * IT UNDER THE TERMS OF THE GNU LESSER GENERAL PUBLIC LICENSE 5 * AS PUBLISHED BY THE FREE SOFTWARE FOUNDATION; 6 * EITHER VERSION 2.1 OF THE LICENSE, OR (AT YOUR OPTION) ANY LATER VERSION. 7 * 8 * THIS LIBRARY IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT 9 * WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF 10 * MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU 11 * LESSER GENERAL PUBLIC LICENSE FOR MORE DETAILS. 12 * 13 * YOU SHOULD HAVE RECEIVED A COPY OF THE GNU LESSER GENERAL PUBLIC LICENSE 14 * ALONG WITH THIS LIBRARY; IF NOT, WRITE TO THE FREE SOFTWARE FOUNDATION, 15 * INC., 59 TEMPLE PLACE, SUITE 330, BOSTON, MA 02111-1307 USA 16 * 17 *************************************************************************/ 18 19 20 /****************************************************************************** 21 * - Nom du fichier : test27.c 22 * 23 * - Description : creation de maillages structures (grille cartesienne | 24 * grille standard ) dans le fichier test27.med 25 * 26 *****************************************************************************/ 27 28 #include <med.h> 29 #define MESGERR 30 #include <med_utils.h> 31 32 #ifdef DEF_LECT_ECR 33 #define MODE_ACCES MED_LECTURE_ECRITURE 34 #elif DEF_LECT_AJOUT 35 #define MODE_ACCES MED_LECTURE_AJOUT 36 #else 37 #define MODE_ACCES MED_CREATION 38 #endif 39 40 int main (int argc, char **argv) 41 42 43 { 44 med_idt fid; 45 /* la dimension du maillage */ 46 med_int mdim,axe; 47 med_float indice[4]; 48 med_float coo[8] = {0.0,0.0,1.0,0.0,0.0,1.0,1.0,1.0}; 49 med_int nnoeuds = 4; 50 med_int structure_grille[2] = {2,2}; 51 med_int nind; 52 /* nom du maillage de longueur maxi MED_TAILLE_NOM */ 53 char maa[MED_TAILLE_NOM+1]; 54 /* composantes et unites */ 55 char comp[MED_TAILLE_PNOM+1]; 56 char unit[MED_TAILLE_PNOM+1]; 57 /* 12345678901234561234567890123456 */ 58 char comp2[2*MED_TAILLE_PNOM+1] = "X Y "; 59 char unit2[2*MED_TAILLE_PNOM+1] = "cm cm "; 60 61 /* Creation du fichier test27.med */ 62 fid = MEDouvrir("test27.med",MODE_ACCES); 63 if (fid < 0) { 64 MESSAGE("Erreur a la creation du fichier test27.med"); 65 return -1; 66 } 67 printf("Creation du fichier test27.med \n"); 68 69 /* Creation d'un maillage MED_NON_STRUCTURE */ 70 if (MEDmaaCr(fid,"maillage vide",3,MED_NON_STRUCTURE,"un maillage vide") < 0) { 71 MESSAGE("Erreur a la creation du maillage MED_NON_STRUCTURE"); 72 return -1; 73 } 74 75 /* creation d'une grille cartesienne de dimension 2 */ 76 strcpy(maa,"grille_cartesienne"); 77 mdim = 2; 78 /* on commence par definir un maillage MED_STRUCTURE 79 de dimension 2 */ 80 if (MEDmaaCr(fid,maa,mdim,MED_STRUCTURE, 81 "un exemple de grille cartesienne") < 0) { 82 MESSAGE("Erreur a la creation de la grille"); 83 return -1; 84 } 85 printf("Creation d'un maillage structure MED_STRUCTURE \n"); 86 87 /* On specifie la nature du maillage structure : MED_GRILLE_CARTESIENNE */ 88 if (MEDnatureGrilleEcr(fid,maa,MED_GRILLE_CARTESIENNE) < 0) { 89 MESSAGE("Erreur a l'ecriture de la nature de la grille"); 90 return -1; 91 } 92 printf("On definit la nature du maillage structure : MED_GRILLE_CARTESIENNE \n"); 93 94 /* on definit les indices des coordonnees de la grille selon chaque dimension */ 95 /* axe des "X" */ 96 indice[0] = 1.1; 97 indice[1] = 1.2; 98 indice[2] = 1.3; 99 indice[3] = 1.4; 100 nind = 4; 101 axe = 1; 102 strcpy(comp,"X"); 103 strcpy(unit,"cm"); 104 if (MEDindicesCoordEcr(fid,maa,mdim,indice,nind,axe,comp,unit) < 0) { 105 MESSAGE("Erreur a l'ecriture de l'axe X"); 106 return -1; 107 } 108 printf("Ecriture des indices des coordonnees selon l'axe des X \n"); 109 110 /* axe des "Y" */ 111 indice[0] = 2.1; 112 indice[1] = 2.2; 113 indice[2] = 2.3; 114 indice[3] = 2.4; 115 nind = 4; 116 strcpy(comp,"Y"); 117 strcpy(unit,"cm"); 118 axe = 2; 119 if (MEDindicesCoordEcr(fid,maa,mdim,indice,nind,axe,comp,unit) < 0) { 120 MESSAGE("Erreur a l'ecriture de l'axe Y"); 121 return -1; 122 } 123 printf("Ecriture des indices des coordonnees selon l'axe des Y \n"); 124 125 /* Creation d'une grille MED_STANDARD de dimension 2 */ 126 strcpy(maa,"grille_standard"); 127 mdim = 2; 128 /* on commence par definir un maillage MED_STRUCTURE 129 de dimension 2 */ 130 if (MEDmaaCr(fid,maa,mdim,MED_STRUCTURE, 131 "un exemple de grille standard") < 0) { 132 MESSAGE("Erreur a la creation de la 2e grille"); 133 return -1; 134 } 135 printf("Creation d'un maillage structure MED_STRUCTURE \n"); 136 137 /* On specifie la nature du maillage structure : MED_GRILLE_STANDARD */ 138 if (MEDnatureGrilleEcr(fid,maa,MED_GRILLE_STANDARD) < 0) { 139 MESSAGE("Erreur a l'ecriture du type de la grille"); 140 return -1; 141 } 142 printf("On definit la nature du maillage structure : MED_GRILLE_STANDARD \n"); 143 144 /* On ecrit les coordonnees de la grille */ 145 if (MEDcoordEcr(fid,maa,mdim,coo,MED_FULL_INTERLACE,nnoeuds,MED_CART,comp2,unit2) < 0) { 146 MESSAGE("Erreur a l'ecriture des noeuds de la grille MED_GRILLE_STANDARD"); 147 return -1; 148 } 149 printf("Ecriture des coordonnees des noeuds \n"); 150 151 /* On definit la structure de la grille */ 152 if (MEDstructureCoordEcr(fid,maa,mdim,structure_grille) < 0) { 153 MESSAGE("Erreur a l'ecriture de la structure de la grille"); 154 return -1; 155 } 156 printf("Ecriture de la structure de la grille : / 2,2 / \n"); 157 158 /* On ferme le fichier */ 159 if (MEDfermer(fid) < 0) { 160 MESSAGE("Erreur a la fermeture du fichier"); 161 return -1; 162 } 163 printf("Fermeture du fichier \n"); 164 165 return 0; 166 }