* fichier : solsoltest.dgibi ************************************************************************ ************************************************************************ SAUT PAGE ; * ************************************************************************** * UTILISATION DES OPERATEURS CHI1 ET CHI2 * EN PRESENCE DE SOLUTIONS SOLIDES ************************************************************************** * * repertoire des fichiers "divers" * option echo 0 ; a= 0. 0. ; b= 0. 5. ; option elem qua4 ; ab= a droit 1 b ; vec1= 100. 0. ; mt = ab trans 1 vec1 ; * ************************************************************************* ************************************************************************* * 2 SOLUTIONS SOLIDES : * * LA SOLUTION SOLIDE TERNAIRE CALCITE-MAGNESITE-SIDERITE , ET * LA SOLUTION SOLIDE D'ORDRE 27 DE FORMULE GENERALE : * (Si Al )(Al Fe Mg ) O (OH) (K Na Ca Mg ) * 4-x x 2-y-2/3z y z 10 2 x-a-2b-2g a b g * * 8 MINERAUX : * * PHLOGOPITE-H, GIBBSITE, MICROCLINE, SAPONITE, QUARTZ, * ALBITE, ILLITE, DOLOMITE ************************************************************************* ************************************************************************* * ************************************************************************* * CAS 1 : LES 2 SOLUTIONS SOLIDES AVEC PHLOGOPITE-H, GIBBSITE, DOLOMITE, * MICROCLINE, QUARTZ, SAPONITE, ALBITE ET ON IMPOSE LE PH ************************************************************************* * * ------------------------- * TABLE DES DONNEES DE CHI1 * ------------------------- * tabdon = TABLE ; 2148 2149 2150 2168 2169 2313 2314 2315 2316 2317 2369 2379 2380 2381 2171 2280 2278 2279 2284 2285 2282 2283 2366 2368 2370 2371 2372 2376 2377 2378 ; tabdon.NVSOSO = TABLE ; tabdon.NVSOSO.1 = TABLE ; tabdon.NVSOSO.1 . IDEN = 30000 ; tabdon.NVSOSO.1 . ITYP = 5 ; 2369 2379 2380 2381 2171 2280 2278 2279 2284 2285 2282 2283 2366 2368 2370 2371 2372 2376 2377 2378 ; tabdon.NVSOSO.2 = TABLE ; tabdon.NVSOSO.2 . IDEN = 31000 ; tabdon.NVSOSO.2 . ITYP = 5 ; tabdon.CLIM = TABLE ; tabdon.BDD = 'STRASBG' ; tb1 = CHI1 tabdon *IMPCHI1 tb1 'TB1' ; * *----------------------------- * TABLE DES PARAMETRES DE CHI2 *----------------------------- * tbparm= TABLE ; tbparm.ITMAX= 200 ; tbparm.ITERSOLI = 20 ; tbparm.IAFFICHE = 2 ; tbparm.NFI= 8 ; tbparm.EPS= 1.D-4 ; tbpar3= TABLE ; tbpar3.'SOUSTYPE'= 'DONNEES_CHIMIQUES' ; X061 -3.50 X050 -3.27 X001 -2.77 X005 -1.22 X101 -5.19 X007 -8.54 ; totco= totco1 + totco2 ; totca= totca1 + totca2 ; totna= totna1 + totna2 ; totk= totk1 + totk2 ; total= total1 + total2 ; totmg= totmg1 + totmg2 ; totfe2i= totfe2i1 + totfe2i2 ; totsi= totsi1 + totsi2 ; tbpar3.TOT= totca + totna + totk + total + totmg + totsi +totfe2i + toth + totco ; climph = climph1 + climph2 ; tbpar3.CLIM= climph ; *option echo 1 ; *IMPCHI2 tb3 'TB3' ; option echo 0 ; * *----------------------- * CONTROLE DES RESULTATS *----------------------- * 'NATURE' DISCRET ; FIONTES= FIONTE1 + FIONTE2 ; *LIST VERR1 ; W027 0. W028 1.31176E-05 W029 1.54413E-04 W030 8.44764E-04 W031 0. W032 4.81838E-06 W033 1.43656E-03 'NATURE' DISCRET ; W026 (0. - 0.) W027 (0. - 0.) W028 (1.31176E-05 - 1.31176E-05) W029 (1.54414E-04 - 1.54413E-04) W030 (7.34384E-04 - 8.44764E-04) W031 (3.77770E-05 - 0.) W032 (3.83440E-06 - 4.81838E-06) W033 (1.43653E-03 - 1.43656E-03) 'NATURE' DISCRET ; PRECTES= PRECTE1+PRECTE2 ; *LIST PRECTES ; *LIST TB3.PREC ; PRECD=PRECTES / 50. ; *LIST VERR2 ; *LIST TB3.TYP5 ; W027 2.73971E-03 W028 0.00000E+00 W029 0.00000E+00 W030 0.00000E+00 W031 9.24659E-01 W032 0.00000E+00 W033 0.00000E+00 'NATURE' DISCRET ; W026 (7.05903E-02 - 7.05934E-02) W027 (2.73959E-03 - 2.73971E-03) W028 (0.00000E+00 - 0.00000E+00) W029 (0.00000E+00 - 0.00000E+00) W030 (0.00000E+00 - 0.00000E+00) W031 (0.00000E+00 - 9.24659E-01) W032 (0.00000E+00 - 0.00000E+00) W033 (0.00000E+00 - 0.00000E+00) 'NATURE' DISCRET ; TY5TES= TY5TE1+TY5TE2 ; *LIST TY5TES ; TY5D= TY5TES/ 50. ; VERR3= ( ABS ( TY5TES - TB3.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ; *LIST VERR3 ; W036 6.85256E-15 W037 1.71692E-11 W038 7.52657E-10 W039 9.00205E-10 W040 2.99667E-09 W041 5.91304E-14 W042 7.07221E-14 W043 1.34885E-15 W044 2.35425E-13 W045 3.60628E-03 W046 3.66581E-03 W047 6.39823E-01 'NATURE' DISCRET; W050 4.71032E-16 W051 4.13296E-19 W052 4.30296E-09 W053 8.24023E-14 W054 2.46968E-14 W055 2.06489E-14 W056 7.21357E-17 W057 2.16697E-17 W058 1.81179E-17 W059 7.51028E-07 W060 2.25610E-07 W061 1.88631E-07 W063 6.09031E-12 W064 5.98672E-01 'NATURE' DISCRET; W035 (4.42490E-04 - 4.42526E-04) W036 (6.85067E-15 - 6.85256E-15) W037 (1.85587E-11 - 1.71692E-11) W038 (7.52490E-10 - 7.52657E-10) W039 (9.00004E-10 - 9.00205E-10) W040 (2.99604E-09 - 2.99667E-09) W041 (5.91135E-14 - 5.91304E-14) 'NATURE' DISCRET; W043 (1.45792E-15 - 1.34885E-15) W044 (2.35360E-13 - 2.35425E-13) W045 (3.60522E-03 - 3.60628E-03) W046 (3.96216E-03 - 3.66581E-03) W047 (6.39634E-01 - 6.39823E-01) W048 (1.60652E-01 - 1.60701E-01) 'NATURE' DISCRET; W050 (5.09151E-16 - 4.71032E-16) W051 (4.46720E-19 - 4.13296E-19) W052 (4.65087E-09 - 4.30296E-09) W053 (8.23848E-14 - 8.24023E-14) W054 (2.46913E-14 - 2.46968E-14) W055 (2.06443E-14 - 2.06489E-14) W056 (7.21166E-17 - 7.21357E-17) 'NATURE' DISCRET; W058 (1.81129E-17 - 1.81179E-17) W059 (7.50817E-07 - 7.51028E-07) W060 (2.25544E-07 - 2.25610E-07) W061 (1.88577E-07 - 1.88631E-07) W063 (6.08901E-12 - 6.09031E-12) W064 (5.98672E-01 - 5.98672E-01) 'NATURE' DISCRET; FRACTES= FRACTE1+FRACTE2+FRACTE3+FRACTE4+FRACTE5+FRACTE6 ; *LIST FRACTES ; FRACD= FRACTES/50 ; VERR4= (ABS ( FRACTES - TB3.POLE)) MASQUE SUPERIEUR FRACD SOMME ; *LIST VERR4 ; * *************************************************************************** * CAS 2 : SYSTEME SANS CONTRAINTE AYANT LE PH PRECEDENT ET DE L' ILLITE *************************************************************************** * *--------------------------- * TABLE DES DONNEES DES CHI1 *--------------------------- * tabdon2 = TABLE ; tabdon2.IDEN = tabdon.IDEN ; 2148 2149 2150 2168 2169 2313 2314 2315 2316 2317 2369 2379 2380 2381 2171 2280 2278 2279 2284 2285 2282 2283 2366 2368 2370 2371 2372 2376 2377 2378 ; tabdon2.NVSOSO = tabdon.NVSOSO ; tabdon2.BDD = 'STRASBG' ; tb2 = CHI1 tabdon2 *IMPCHI1 tb1 'TB1' ; * *----------------------------- * TABLE DES PARAMETRES DE CHI2 *----------------------------- * tbpar4= TABLE ; tbpar4.'SOUSTYPE'= 'DONNEES_CHIMIQUES' ; X061 -4.00 X050 -2.45 X001 -2.74 X005 -1.22 X101 -5.23 X007 -8.18 ; toto = (tb3.TYP3) * (-1.) ; tbpar4.TOT = totca + totna + totk + total + totmg + totfe2i + totsi + toth + totco ; *LIST tbpar4.TOT ; tbparm2 = TABLE ; tbparm2.ITMAX= 200 ; tbparm2.ITERSOLI = 20 ; tbparm2.IAFFICHE = 2 ; tbparm2.NFI= 8 ; tbparm2.EPS= 1.D-4 ; option echo 1 ; *IMPCHI2 tb5 'TB5' ; option echo 0 ; * *----------------------- * CONTROLE DES RESULTATS *----------------------- * 'NATURE' DISCRET ; FIONTES= FIONTE1 + FIONTE2 ; *LIST VERR1 ; W027 0. W028 1.31140E-05 W029 1.36945E-04 W030 8.02850E-04 W031 0. W032 0.00000E+00 W033 1.61111E-05 W034 1.45930E-03 'NATURE' DISCRET ; W026 (0.00000E+00 - 0.) W027 (0.00000E+00 - 0.) W028 (1.31034E-05 - 1.31140E-05) W029 (8.35208E-05 - 1.36945E-04) W030 (5.48759E-04 - 8.02850E-04) W031 (0.00000E+00 - 0.) W032 (0.00000E+00 - 0.00000E+00) W033 (8.87255E-05 - 1.61111E-05) W034 (1.52577E-03 - 1.45930E-03) 'NATURE' DISCRET ; PRECTES= PRECTE1+PRECTE2 ; *LIST PRECTES ; *LIST TB3.PREC ; PRECD=PRECTES / 50. ; *LIST VERR2 ; *LIST TB3.TYP5 ; W027 6.75986E-02 W028 0.00000E+00 W029 0.00000E+00 W030 0.00000E+00 W031 8.58262E-01 W032 9.92157E-01 W033 0.00000E+00 W034 0. 'NATURE' DISCRET ; W026 (2.25337E-03 - 2.61667E-03) W027 (5.86999E-02 - 6.75986E-02) W028 (0.00000E+00 - 0.00000E+00) W029 (0.00000E+00 - 0.00000E+00) W030 (0.00000E+00 - 0.00000E+00) W031 (8.75253E-01 - 8.58262E-01) W032 (9.67437E-01 - 9.92157E-01) W033 (0.00000E+00 - 0.00000E+00) W034 (0.00000E+00 - 0.00000E+00) 'NATURE' DISCRET ; TY5TES= TY5TE1+TY5TE2 ; *LIST TY5TES ; TY5D= TY5TES/ 50. ; VERR7= ( ABS ( TY5TES - TB5.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ; *LIST VERR3 ; W037 6.83202E-15 W038 1.61230E-11 W039 7.51750E-10 W040 8.98560E-10 W041 2.99190E-09 W042 5.89761E-14 W043 7.04936E-14 W044 1.26488E-15 W045 2.34720E-13 W046 3.60700E-03 W047 3.44863E-03 W048 6.39951E-01 'NATURE' DISCRET; W051 4.42330E-16 W052 3.87658E-19 W053 4.04681E-09 W054 8.22710E-14 W055 2.46517E-14 W056 2.06240E-14 W057 7.19364E-17 W058 2.16048E-17 W059 1.80749E-17 W060 7.50954E-07 W061 2.25535E-07 W062 1.88686E-07 W064 6.08060E-12 W065 5.98974E-01 'NATURE' DISCRET; W036 (4.43949E-04 - 4.42837E-04) W037 (6.73199E-15 - 6.83202E-15) W038 (1.62507E-11 - 1.61230E-11) W039 (7.55828E-10 - 7.51750E-10) W040 (9.00193E-10 - 8.98560E-10) W041 (3.00728E-09 - 2.99190E-09) W042 (5.81290E-14 - 5.89761E-14) 'NATURE' DISCRET; W044 (1.24981E-15 - 1.26488E-15) W045 (2.31283E-13 - 2.34720E-13) W046 (3.60910E-03 - 3.60700E-03) W047 (3.46017E-03 - 3.44863E-03) W048 (6.40323E-01 - 6.39951E-01) W049 (1.60934E-01 - 1.60795E-01) 'NATURE' DISCRET; W051 (4.45834E-16 - 4.42330E-16) W052 (3.84309E-19 - 3.87658E-19) W053 (4.06344E-09 - 4.04681E-09) W054 (8.26940E-14 - 8.22710E-14) W055 (2.46965E-14 - 2.46517E-14) W056 (2.07359E-14 - 2.06240E-14) W057 (7.11184E-17 - 7.19364E-17) 'NATURE' DISCRET; W059 (1.78744E-17 - 1.80749E-17) W060 (7.51960E-07 - 7.50954E-07) W061 (2.25090E-07 - 2.25535E-07) W062 (1.88992E-07 - 1.88686E-07) W064 (6.11187E-12 - 6.08060E-12) W065 (6.00700E-01 - 5.98974E-01) 'NATURE' DISCRET; FRACTES= FRACTE1+FRACTE2+FRACTE3+FRACTE4+FRACTE5+FRACTE6 ; *LIST FRACTES ; FRACD= FRACTES/50 ; VERR8= (ABS ( FRACTES - TB5.POLE)) MASQUE SUPERIEUR FRACD SOMME ; *LIST VERR4 ; * * VERR= VERR1+VERR2+VERR3+VERR4+VERR5+VERR6+VERR7+VERR8 ; VERR= VERR7+VERR8 ; SI (VERR EGA 0 ) ; SINO ; FINSI ; FIN ;
© Cast3M 2003 - Tous droits réservés.
Mentions légales