Télécharger solsoltest.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : solsoltest.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. SAUT PAGE ;
  5. *
  6. **************************************************************************
  7. * UTILISATION DES OPERATEURS CHI1 ET CHI2
  8. * EN PRESENCE DE SOLUTIONS SOLIDES
  9. **************************************************************************
  10. *
  11. option echo 0 ;
  12. option dime 2 ;
  13. a= 0. 0. ;
  14. b= 0. 5. ;
  15. option elem qua4 ;
  16. ab= a droit 1 b ;
  17. vec1= 100. 0. ;
  18. mt = ab trans 1 vec1 ;
  19. *
  20. *************************************************************************
  21. *************************************************************************
  22. * 2 SOLUTIONS SOLIDES :
  23. *
  24. * LA SOLUTION SOLIDE TERNAIRE CALCITE-MAGNESITE-SIDERITE , ET
  25. * LA SOLUTION SOLIDE D'ORDRE 27 DE FORMULE GENERALE :
  26. * (Si Al )(Al Fe Mg ) O (OH) (K Na Ca Mg )
  27. * 4-x x 2-y-2/3z y z 10 2 x-a-2b-2g a b g
  28. *
  29. * 8 MINERAUX :
  30. *
  31. * PHLOGOPITE-H, GIBBSITE, MICROCLINE, SAPONITE, QUARTZ,
  32. * ALBITE, ILLITE, DOLOMITE
  33. *************************************************************************
  34. *************************************************************************
  35. *
  36. *************************************************************************
  37. * CAS 1 : LES 2 SOLUTIONS SOLIDES AVEC PHLOGOPITE-H, GIBBSITE, DOLOMITE,
  38. * MICROCLINE, QUARTZ, SAPONITE, ALBITE ET ON IMPOSE LE PH
  39. *************************************************************************
  40. *
  41. * -------------------------
  42. * TABLE DES DONNEES DE CHI1
  43. * -------------------------
  44. *
  45. tabdon = TABLE ;
  46. tabdon.IDEN = LECT 1 2 4 5 7 50 61 101 165 ;
  47. tabdon.CHXMX = LECT 2286 2144 2158 2184 2157 2152 2159
  48. 2148 2149 2150
  49. 2168 2169 2313 2314 2315 2316 2317
  50. 2369 2379 2380 2381
  51. 2171 2280 2278 2279
  52. 2284 2285 2282 2283
  53. 2366 2368 2370 2371 2372
  54. 2376 2377 2378 ;
  55. tabdon.NVSOSO = TABLE ;
  56. tabdon.NVSOSO.1 = TABLE ;
  57. tabdon.NVSOSO.1 . IDEN = 30000 ;
  58. tabdon.NVSOSO.1 . ITYP = 5 ;
  59. tabdon.NVSOSO.1 . SOLID = LECT 2168 2169 2313 2314 2315 2316 2317
  60. 2369 2379 2380 2381
  61. 2171 2280 2278 2279
  62. 2284 2285 2282 2283
  63. 2366 2368 2370 2371 2372
  64. 2376 2377 2378 ;
  65. tabdon.NVSOSO.2 = TABLE ;
  66. tabdon.NVSOSO.2 . IDEN = 31000 ;
  67. tabdon.NVSOSO.2 . ITYP = 5 ;
  68. tabdon.NVSOSO.2 . SOLID = LECT 2148 2149 2150 ;
  69. tabdon.CLIM = TABLE ;
  70. tabdon.CLIM.TYP3 = LECT 50 ;
  71. tabdon.BDD = 'STRASBG' ;
  72.  
  73. tb1 = CHI1 tabdon
  74. COMP '/u/castem/divers/COMPSM' ;
  75.  
  76. *IMPCHI1 tb1 'TB1' ;
  77. *
  78. *-----------------------------
  79. * TABLE DES PARAMETRES DE CHI2
  80. *-----------------------------
  81. *
  82. tbparm= TABLE ;
  83. tbparm.ITMAX= 200 ;
  84. tbparm.ITERSOLI = 20 ;
  85. tbparm.IAFFICHE = 2 ;
  86. tbparm.NFI= 8 ;
  87. tbparm.EPS= 1.D-4 ;
  88. tbparm.IMPRIM= LECT 1 4 ;
  89. tbparm.SORTIE= MOTS 'PREC' 'SOLU' 'FION' 'TYP3' 'TYP5' 'POLE' 'LOGK' ;
  90. tbpar3= TABLE ;
  91. tbpar3.'SOUSTYPE'= 'DONNEES_CHIMIQUES' ;
  92. tbpar3.LOGC= MANU CHPO mt 9 X165 -7.68 X002 -2.59 X004 -4.86
  93. X061 -3.50 X050 -3.27 X001 -2.77
  94. X005 -1.22 X101 -5.19 X007 -8.54 ;
  95. totco1= MANU CHPO mt 1 X101 3.6D-3 ;
  96. totco2= MANU CHPO ab 1 X101 (3.6D-3 - 3.6D-3) ;
  97. totco= totco1 + totco2 ;
  98. totca1= MANU CHPO mt 1 X001 3.D-3 ;
  99. totca2= MANU CHPO ab 1 X001 (3.D-3 - 3.D-3) ;
  100. totca= totca1 + totca2 ;
  101. totna1= MANU CHPO mt 1 X005 1.676D-1 ;
  102. totna2= MANU CHPO ab 1 X005 (1.676D-1 - 1.676D-1) ;
  103. totna= totna1 + totna2 ;
  104. totk1= MANU CHPO mt 1 X004 1.5D-4 ;
  105. totk2= MANU CHPO ab 1 X004 (2.0D-4 - 1.5D-4) ;
  106. totk= totk1 + totk2 ;
  107. total1= MANU CHPO mt 1 X165 8.7D-3 ;
  108. total2= MANU CHPO ab 1 X165 (4.0D-2 - 8.7D-3) ;
  109. total= total1 + total2 ;
  110. totmg1= MANU CHPO mt 1 X002 3.9D-3 ;
  111. totmg2= MANU CHPO ab 1 X002 (3.9D-3 - 3.9D-3) ;
  112. totmg= totmg1 + totmg2 ;
  113. totfe2i1= MANU CHPO mt 1 X007 4.0D-5 ;
  114. totfe2i2= MANU CHPO ab 1 X007 (4.0D-5 - 4.0D-5) ;
  115. totfe2i= totfe2i1 + totfe2i2 ;
  116. totsi1= MANU CHPO mt 1 X061 1.0D-3 ;
  117. totsi2= MANU CHPO ab 1 X061 (1.0D-3 - 1.0D-3) ;
  118. totsi= totsi1 + totsi2 ;
  119. toth = MANU CHPO mt 1 X050 0.D0 ;
  120. tbpar3.TOT= totca + totna + totk + total + totmg + totsi +totfe2i +
  121. toth + totco ;
  122. LIST tbpar3.TOT ;
  123. climph1= MANU CHPO mt 1 W024 7.75 ;
  124. climph2= MANU CHPO ab 1 W024 ( 7.75 - 7.75 ) ;
  125. climph = climph1 + climph2 ;
  126. tbpar3.CLIM= climph ;
  127.  
  128. tb3= CHI2 tb1 tbparm tbpar3 ;
  129.  
  130. *option echo 1 ;
  131.  
  132. *IMPCHI2 tb3 'TB3' ;
  133.  
  134. option echo 0 ;
  135. *
  136. *-----------------------
  137. * CONTROLE DES RESULTATS
  138. *-----------------------
  139. *
  140. FIONTE1=MANU CHPO mt 1 SCAL 9.48864D-2 'NATURE' DISCRET ;
  141. FIONTE2=MANU CHPO ab 1 SCAL (9.48931D-2 - 9.48864D-2 )
  142. 'NATURE' DISCRET ;
  143. FIONTES= FIONTE1 + FIONTE2 ;
  144. VERR1= ( ABS ( FIONTES - TB3.FION )) MASQUE SUPERIEUR 5.D-8 SOMME ;
  145. *LIST VERR1 ;
  146. PRECTE2= MANU CHPO mt 9 W025 8.68203E-03 W026 0.
  147. W027 0. W028 1.31176E-05
  148. W029 1.54413E-04 W030 8.44764E-04
  149. W031 0. W032 4.81838E-06
  150. W033 1.43656E-03 'NATURE' DISCRET ;
  151.  
  152. PRECTE1= MANU CHPO ab 9 W025 (3.99452E-02 - 8.68203E-03)
  153. W026 (0. - 0.)
  154. W027 (0. - 0.)
  155. W028 (1.31176E-05 - 1.31176E-05)
  156. W029 (1.54414E-04 - 1.54413E-04)
  157. W030 (7.34384E-04 - 8.44764E-04)
  158. W031 (3.77770E-05 - 0.)
  159. W032 (3.83440E-06 - 4.81838E-06)
  160. W033 (1.43653E-03 - 1.43656E-03)
  161. 'NATURE' DISCRET ;
  162. PRECTES= PRECTE1+PRECTE2 ;
  163. *LIST PRECTES ;
  164. *LIST TB3.PREC ;
  165. PRECD=PRECTES / 50. ;
  166. VERR2= ( ABS ( PRECTES - TB3.PREC )) MASQUE SUPERIEUR PRECD SOMME ;
  167. *LIST VERR2 ;
  168. *LIST TB3.TYP5 ;
  169. TY5TE1= MANU CHPO mt 9 W025 0.00000E+00 W026 7.05934E-02
  170. W027 2.73971E-03 W028 0.00000E+00
  171. W029 0.00000E+00 W030 0.00000E+00
  172. W031 9.24659E-01 W032 0.00000E+00
  173. W033 0.00000E+00 'NATURE' DISCRET ;
  174.  
  175. TY5TE2= MANU CHPO ab 9 W025 (0.00000E+00 - 0.00000E+00)
  176. W026 (7.05903E-02 - 7.05934E-02)
  177. W027 (2.73959E-03 - 2.73971E-03)
  178. W028 (0.00000E+00 - 0.00000E+00)
  179. W029 (0.00000E+00 - 0.00000E+00)
  180. W030 (0.00000E+00 - 0.00000E+00)
  181. W031 (0.00000E+00 - 9.24659E-01)
  182. W032 (0.00000E+00 - 0.00000E+00)
  183. W033 (0.00000E+00 - 0.00000E+00)
  184. 'NATURE' DISCRET ;
  185. TY5TES= TY5TE1+TY5TE2 ;
  186. *LIST TY5TES ;
  187. TY5D= TY5TES/ 50. ;
  188. VERR3= ( ABS ( TY5TES - TB3.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  189. *LIST VERR3 ;
  190.  
  191. FRACTE1= MANU CHPO mt 14 W034 4.00886E-01 W035 4.42526E-04
  192. W036 6.85256E-15 W037 1.71692E-11 W038 7.52657E-10
  193. W039 9.00205E-10 W040 2.99667E-09 W041 5.91304E-14
  194. W042 7.07221E-14 W043 1.34885E-15 W044 2.35425E-13
  195. W045 3.60628E-03 W046 3.66581E-03 W047 6.39823E-01
  196. 'NATURE' DISCRET;
  197.  
  198. FRACTE2= MANU CHPO mt 16 W048 1.60701E-01 W049 1.92204E-01
  199. W050 4.71032E-16 W051 4.13296E-19 W052 4.30296E-09
  200. W053 8.24023E-14 W054 2.46968E-14 W055 2.06489E-14
  201. W056 7.21357E-17 W057 2.16697E-17 W058 1.81179E-17
  202. W059 7.51028E-07 W060 2.25610E-07 W061 1.88631E-07
  203. W063 6.09031E-12 W064 5.98672E-01 'NATURE' DISCRET;
  204.  
  205. FRACTE3= MANU CHPO ab 8 W034 (4.00885E-01 - 4.00886E-01)
  206. W035 (4.42490E-04 - 4.42526E-04)
  207. W036 (6.85067E-15 - 6.85256E-15)
  208. W037 (1.85587E-11 - 1.71692E-11)
  209. W038 (7.52490E-10 - 7.52657E-10)
  210. W039 (9.00004E-10 - 9.00205E-10)
  211. W040 (2.99604E-09 - 2.99667E-09)
  212. W041 (5.91135E-14 - 5.91304E-14)
  213. 'NATURE' DISCRET;
  214.  
  215. FRACTE4= MANU CHPO ab 7 W042 (7.07018E-14 - 7.07221E-14)
  216. W043 (1.45792E-15 - 1.34885E-15)
  217. W044 (2.35360E-13 - 2.35425E-13)
  218. W045 (3.60522E-03 - 3.60628E-03)
  219. W046 (3.96216E-03 - 3.66581E-03)
  220. W047 (6.39634E-01 - 6.39823E-01)
  221. W048 (1.60652E-01 - 1.60701E-01)
  222. 'NATURE' DISCRET;
  223.  
  224. FRACTE5= MANU CHPO ab 8 W049 (1.92145E-01 - 1.92204E-01)
  225. W050 (5.09151E-16 - 4.71032E-16)
  226. W051 (4.46720E-19 - 4.13296E-19)
  227. W052 (4.65087E-09 - 4.30296E-09)
  228. W053 (8.23848E-14 - 8.24023E-14)
  229. W054 (2.46913E-14 - 2.46968E-14)
  230. W055 (2.06443E-14 - 2.06489E-14)
  231. W056 (7.21166E-17 - 7.21357E-17)
  232. 'NATURE' DISCRET;
  233.  
  234. FRACTE6= MANU CHPO ab 7 W057 (2.16637E-17 - 2.16697E-17)
  235. W058 (1.81129E-17 - 1.81179E-17)
  236. W059 (7.50817E-07 - 7.51028E-07)
  237. W060 (2.25544E-07 - 2.25610E-07)
  238. W061 (1.88577E-07 - 1.88631E-07)
  239. W063 (6.08901E-12 - 6.09031E-12)
  240. W064 (5.98672E-01 - 5.98672E-01)
  241. 'NATURE' DISCRET;
  242.  
  243. FRACTES= FRACTE1+FRACTE2+FRACTE3+FRACTE4+FRACTE5+FRACTE6 ;
  244. *LIST FRACTES ;
  245. FRACD= FRACTES/50 ;
  246. VERR4= (ABS ( FRACTES - TB3.POLE)) MASQUE SUPERIEUR FRACD SOMME ;
  247. *LIST VERR4 ;
  248. *
  249. ***************************************************************************
  250. * CAS 2 : SYSTEME SANS CONTRAINTE AYANT LE PH PRECEDENT ET DE L' ILLITE
  251. ***************************************************************************
  252. *
  253. *---------------------------
  254. * TABLE DES DONNEES DES CHI1
  255. *---------------------------
  256. *
  257. tabdon2 = TABLE ;
  258. tabdon2.IDEN = tabdon.IDEN ;
  259. tabdon2.CHXMX = LECT 2286 2144 2158 2184 2157 2152 2159 2172
  260. 2148 2149 2150
  261. 2168 2169 2313 2314 2315 2316 2317
  262. 2369 2379 2380 2381
  263. 2171 2280 2278 2279
  264. 2284 2285 2282 2283
  265. 2366 2368 2370 2371 2372
  266. 2376 2377 2378 ;
  267. tabdon2.NVSOSO = tabdon.NVSOSO ;
  268. tabdon2.BDD = 'STRASBG' ;
  269.  
  270. tb2 = CHI1 tabdon2
  271. COMP '/u/castem/divers/COMPSM' ;
  272.  
  273. *IMPCHI1 tb1 'TB1' ;
  274.  
  275. *
  276. *-----------------------------
  277. * TABLE DES PARAMETRES DE CHI2
  278. *-----------------------------
  279. *
  280. tbpar4= TABLE ;
  281. tbpar4.'SOUSTYPE'= 'DONNEES_CHIMIQUES' ;
  282. tbpar4.LOGC= MANU CHPO mt 9 X165 -7.21 X002 -2.53 X004 -3.82
  283. X061 -4.00 X050 -2.45 X001 -2.74
  284. X005 -1.22 X101 -5.23 X007 -8.18 ;
  285. toto = (tb3.TYP3) * (-1.) ;
  286. toth = EXCO W024 toto X050 ;
  287. tbpar4.TOT = totca + totna + totk + total + totmg + totfe2i + totsi +
  288. toth + totco ;
  289. *LIST tbpar4.TOT ;
  290. tbparm2 = TABLE ;
  291. tbparm2.ITMAX= 200 ;
  292. tbparm2.ITERSOLI = 20 ;
  293. tbparm2.IAFFICHE = 2 ;
  294. tbparm2.NFI= 8 ;
  295. tbparm2.EPS= 1.D-4 ;
  296. tbparm2.IMPRIM= LECT 1 4 ;
  297. tbparm2.SORTIE= MOTS 'PREC' 'SOLU' 'FION' 'TYP5' 'POLE' 'LOGK' ;
  298.  
  299. tb5= CHI2 tb2 tbparm2 tbpar4 ;
  300.  
  301. option echo 1 ;
  302.  
  303. *IMPCHI2 tb5 'TB5' ;
  304.  
  305. option echo 0 ;
  306. *
  307. *-----------------------
  308. * CONTROLE DES RESULTATS
  309. *-----------------------
  310. *
  311. FIONTE1=MANU CHPO mt 1 SCAL 9.49043E-02 'NATURE' DISCRET ;
  312. FIONTE2=MANU CHPO ab 1 SCAL (9.49673E-02 - 9.49043E-02 )
  313. 'NATURE' DISCRET ;
  314. FIONTES= FIONTE1 + FIONTE2 ;
  315. VERR5= ( ABS ( FIONTES - TB5.FION )) MASQUE SUPERIEUR 5.D-8 SOMME ;
  316. *LIST VERR1 ;
  317. PRECTE2= MANU CHPO mt 10 W025 8.64980E-03 W026 0.
  318. W027 0. W028 1.31140E-05
  319. W029 1.36945E-04 W030 8.02850E-04
  320. W031 0. W032 0.00000E+00
  321. W033 1.61111E-05 W034 1.45930E-03
  322. 'NATURE' DISCRET ;
  323.  
  324. PRECTE1= MANU CHPO ab 10 W025 (3.97828E-02 - 8.64980E-03)
  325. W026 (0.00000E+00 - 0.)
  326. W027 (0.00000E+00 - 0.)
  327. W028 (1.31034E-05 - 1.31140E-05)
  328. W029 (8.35208E-05 - 1.36945E-04)
  329. W030 (5.48759E-04 - 8.02850E-04)
  330. W031 (0.00000E+00 - 0.)
  331. W032 (0.00000E+00 - 0.00000E+00)
  332. W033 (8.87255E-05 - 1.61111E-05)
  333. W034 (1.52577E-03 - 1.45930E-03)
  334. 'NATURE' DISCRET ;
  335. PRECTES= PRECTE1+PRECTE2 ;
  336. *LIST PRECTES ;
  337. *LIST TB3.PREC ;
  338. PRECD=PRECTES / 50. ;
  339. VERR6= ( ABS ( PRECTES - TB5.PREC )) MASQUE SUPERIEUR PRECD SOMME ;
  340. *LIST VERR2 ;
  341. *LIST TB3.TYP5 ;
  342. TY5TE1= MANU CHPO mt 10 W025 0.00000E+00 W026 2.61667E-03
  343. W027 6.75986E-02 W028 0.00000E+00
  344. W029 0.00000E+00 W030 0.00000E+00
  345. W031 8.58262E-01 W032 9.92157E-01
  346. W033 0.00000E+00 W034 0.
  347. 'NATURE' DISCRET ;
  348.  
  349. TY5TE2= MANU CHPO ab 10 W025 (0.00000E+00 - 0.00000E+00)
  350. W026 (2.25337E-03 - 2.61667E-03)
  351. W027 (5.86999E-02 - 6.75986E-02)
  352. W028 (0.00000E+00 - 0.00000E+00)
  353. W029 (0.00000E+00 - 0.00000E+00)
  354. W030 (0.00000E+00 - 0.00000E+00)
  355. W031 (8.75253E-01 - 8.58262E-01)
  356. W032 (9.67437E-01 - 9.92157E-01)
  357. W033 (0.00000E+00 - 0.00000E+00)
  358. W034 (0.00000E+00 - 0.00000E+00)
  359. 'NATURE' DISCRET ;
  360. TY5TES= TY5TE1+TY5TE2 ;
  361. *LIST TY5TES ;
  362. TY5D= TY5TES/ 50. ;
  363. VERR7= ( ABS ( TY5TES - TB5.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  364. *LIST VERR3 ;
  365.  
  366. FRACTE1= MANU CHPO mt 14 W035 4.00583E-01 W036 4.42837E-04
  367. W037 6.83202E-15 W038 1.61230E-11 W039 7.51750E-10
  368. W040 8.98560E-10 W041 2.99190E-09 W042 5.89761E-14
  369. W043 7.04936E-14 W044 1.26488E-15 W045 2.34720E-13
  370. W046 3.60700E-03 W047 3.44863E-03 W048 6.39951E-01
  371. 'NATURE' DISCRET;
  372.  
  373. FRACTE2= MANU CHPO mt 16 W049 1.60795E-01 W050 1.92197E-01
  374. W051 4.42330E-16 W052 3.87658E-19 W053 4.04681E-09
  375. W054 8.22710E-14 W055 2.46517E-14 W056 2.06240E-14
  376. W057 7.19364E-17 W058 2.16048E-17 W059 1.80749E-17
  377. W060 7.50954E-07 W061 2.25535E-07 W062 1.88686E-07
  378. W064 6.08060E-12 W065 5.98974E-01 'NATURE' DISCRET;
  379.  
  380. FRACTE3= MANU CHPO ab 8 W035 (3.98856E-01 - 4.00583E-01)
  381. W036 (4.43949E-04 - 4.42837E-04)
  382. W037 (6.73199E-15 - 6.83202E-15)
  383. W038 (1.62507E-11 - 1.61230E-11)
  384. W039 (7.55828E-10 - 7.51750E-10)
  385. W040 (9.00193E-10 - 8.98560E-10)
  386. W041 (3.00728E-09 - 2.99190E-09)
  387. W042 (5.81290E-14 - 5.89761E-14)
  388. 'NATURE' DISCRET;
  389.  
  390. FRACTE4= MANU CHPO ab 7 W043 (6.92318E-14 - 7.04936E-14)
  391. W044 (1.24981E-15 - 1.26488E-15)
  392. W045 (2.31283E-13 - 2.34720E-13)
  393. W046 (3.60910E-03 - 3.60700E-03)
  394. W047 (3.46017E-03 - 3.44863E-03)
  395. W048 (6.40323E-01 - 6.39951E-01)
  396. W049 (1.60934E-01 - 1.60795E-01)
  397. 'NATURE' DISCRET;
  398.  
  399. FRACTE5= MANU CHPO ab 8 W050 (1.91673E-01 - 1.92197E-01)
  400. W051 (4.45834E-16 - 4.42330E-16)
  401. W052 (3.84309E-19 - 3.87658E-19)
  402. W053 (4.06344E-09 - 4.04681E-09)
  403. W054 (8.26940E-14 - 8.22710E-14)
  404. W055 (2.46965E-14 - 2.46517E-14)
  405. W056 (2.07359E-14 - 2.06240E-14)
  406. W057 (7.11184E-17 - 7.19364E-17)
  407. 'NATURE' DISCRET;
  408.  
  409. FRACTE6= MANU CHPO ab 7 W058 (2.12884E-17 - 2.16048E-17)
  410. W059 (1.78744E-17 - 1.80749E-17)
  411. W060 (7.51960E-07 - 7.50954E-07)
  412. W061 (2.25090E-07 - 2.25535E-07)
  413. W062 (1.88992E-07 - 1.88686E-07)
  414. W064 (6.11187E-12 - 6.08060E-12)
  415. W065 (6.00700E-01 - 5.98974E-01)
  416. 'NATURE' DISCRET;
  417.  
  418. FRACTES= FRACTE1+FRACTE2+FRACTE3+FRACTE4+FRACTE5+FRACTE6 ;
  419. *LIST FRACTES ;
  420. FRACD= FRACTES/50 ;
  421. VERR8= (ABS ( FRACTES - TB5.POLE)) MASQUE SUPERIEUR FRACD SOMME ;
  422. *LIST VERR4 ;
  423. *
  424. * VERR= VERR1+VERR2+VERR3+VERR4+VERR5+VERR6+VERR7+VERR8 ;
  425. VERR= VERR7+VERR8 ;
  426. SI (VERR EGA 0 ) ;
  427. ERRE 0 ;
  428. SINO ;
  429. ERRE 5 ;
  430. FINSI ;
  431. FIN ;
  432.  
  433.  
  434.  
  435.  
  436.  
  437.  
  438.  
  439.  
  440.  
  441.  
  442.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales