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

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