Télécharger la plaquette Cast3M


fr en

Télécharger uferdx.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  uferdx.dgibi
  2.  SAUT PAGE ;
  3. *********************************************************************
  4. *    Utilisation des modules CHI1 et CHI2
  5. *                    avec presence de redox
  6. *********************************************************************
  7. * repertoire des fichiers "divers"
  8. DIVERS = VENV 'CASTEM_DIVERS';
  9. *
  10. OPTION DIME 2 ;                                                         
  11. *
  12. *
  13. *                 DEFINITION DU MAILLAGE
  14. *
  15. n1 = 1 ;
  16. n2 = 1   ;
  17. a = 0.0 0.0 ;
  18. b = 0.0 0.1 ;
  19. c = 2. 0.1 ;
  20. d = 2. 0. ;
  21. *
  22. *
  23. option elem qua4 ;
  24. *
  25. ab = a droit n1 b ; 
  26. bc = b droit n2 c ; 
  27. cd = c droit n1 d ; 
  28. da = d droit n2 a ; 
  29. *
  30. *
  31. GP = AB BC CD DA DALL 'PLAN' ;
  32. ELIM 0.001 GP ;
  33. AB= CHANG POI1 AB ;
  34. G1= CHANG POI1 GP ;
  35. G0= DIFF G1 AB ;
  36. *
  37. *                        TABLE ENTREE DE CHI1
  38. TABDON=TABLE ;                                                          
  39. TABDON.IDEN= LECT 5 4 103 99 6 7 32 39 40 41 50 ;
  40.  TABDON.CLIM= TABLE ;
  41.  TABDON.CLIM . TYP6 = LECT 99 ;
  42.  TABDON.CHXMX = LECT 40007 20310 40501 32530 32520 21460
  43.                 32009 32011 ;
  44. *                     on rajoute 3 espèces ne figurant pas dans
  45. *                    la base de données
  46.   TABDON.NVESP= TABLE   ;
  47.   TABDON.NVESP.1 =TABLE ;
  48.   TABDON.NVESP.1 .IDEN= 40501 ;
  49.   TABDON.NVESP.1 .LOGK= -92.81 ;
  50.   TABDON.NVESP.1 .ITYP= 5 ;
  51.   TABDON.NVESP.1 .COMP= LECT 40 99 ;
  52.   TABDON.NVESP.1 .STOECH= PROG 1. 4. ;
  53.   TABDON.NVESP.2 =TABLE ;
  54.   TABDON.NVESP.2 .IDEN= 32530 ;
  55.   TABDON.NVESP.2 .LOGK= 13.03 ;
  56.   TABDON.NVESP.2 .ITYP= 5 ;
  57.   TABDON.NVESP.2 .COMP= LECT 32 40 50 ;
  58.   TABDON.NVESP.2 .STOECH= PROG 1. 3. -14. ;
  59.   TABDON.NVESP.3 =TABLE ;
  60.   TABDON.NVESP.3 .IDEN= 32520 ;
  61.   TABDON.NVESP.3 .LOGK= -2.46 ;
  62.   TABDON.NVESP.3 .ITYP= 5 ;
  63.   TABDON.NVESP.3 .COMP= LECT 32 40 50 ;
  64.   TABDON.NVESP.3 .STOECH= PROG 2. 1. -8. ;
  65. *OPTION DONN 5 ;
  66. *
  67. TB1=CHI1 TABDON
  68.     COMP ('CHAINE' DIVERS '/COMPOM')  ;
  69. *
  70. *                          TABLE DES PARAMETRES DE CHI2
  71. *OPTION DONN 5 ;
  72.  TBPAR2= TABLE ;
  73.  TBPAR2.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  74.  TBPAR2.FIONI= MANU CHPO G1 1 SCAL 0.001 ;
  75.  TBPAR2.LOGC= MANU CHPO G1 11 X005 -3. X004 -11. X103 -3. X099 -5.
  76.          X006 -10. X007 -10. X032 -10. X039 -10. X040 -10. X041 -10.
  77.          X050 -7.  ;
  78. *
  79. *                    on forme le CHPOIN des concentrations totales
  80. *
  81.   TOTNA1= MANU CHPO G1 1 X005 1.0008D-3 ;
  82.   TOTNA2= MANU CHPO AB 1 X005 (1.D-2 - 1.0008D-3) ;
  83.    TOTNA= TOTNA1 + TOTNA2 ;
  84.   TOTK1=MANU CHPO G1 1 X004 1.D-11 ;
  85.   TOTK2=MANU CHPO AB 1 X004 (1.D-5 - 1.D-11) ;
  86.    TOTK= TOTK1 + TOTK2 ;
  87.   TOTCL1=MANU CHPO G1 1 X103 1.0915D-3 ;
  88.   TOTCL2=MANU CHPO AB 1 X103 (1.0022057D-2 - 1.0915D-3) ;
  89.     TOTCL=TOTCL1 + TOTCL2 ;
  90.   TOTE1= MANU CHPO G1 1 X099 4.5619624684D-5 ;
  91.   TOTE2= MANU CHPO AB 1 X099 (1.49708535D-12 - 4.5619624684D-5) ;
  92.     TOTE= TOTE1 + TOTE2 ;
  93.   TOTFE3= MANU CHPO G1 1 X006 4.559E-5 ;
  94.   TOTFE5= MANU CHPO AB 1 X006 (1.D-8 - 4.559E-5) ;
  95.    TOTFE3= TOTFE3 + TOTFE5 ;
  96.    TOTFE2= MANU CHPO G1 1 X007 0.D0 ;
  97.   TOTU6 = MANU CHPO G1 1 X032 2.D-8 ;
  98.   TOTU  = MANU CHPO AB 1 X032 ( 1.D-6 - 2.D-8) ;
  99.   TOTU6 = TOTU6 + TOTU ;
  100.   TOTU3 = MANU CHPO G1 1 X039 0.D0 ;
  101.   TOTU4 = MANU CHPO G1 1 X040 0.D0 ;
  102.   TOTU5 = MANU CHPO G1 1 X041 0.D0 ;
  103.   TOTH1= MANU CHPO G1 1 X050 -4.8578064D-7 ;
  104.   TOTH2= MANU CHPO AB 1 X050 (1.0028126D-5 + 4.8578064D-7) ;
  105.    TOTH = TOTH1 + TOTH2 ;
  106.  TBPAR2.TOT= TOTNA + TOTK + TOTCL + TOTE + TOTFE3 + TOTFE2 + TOTU6 +
  107.                TOTU5 + TOTU4 + TOTU3 + TOTH ;
  108.  TBPARM= TABLE ;
  109.  TBPARM.ITMAX = 80;
  110.  TBPARM.EPS= 1.D-6 ;
  111.  TBPARM.NFI= 4 ;
  112. *TBPARM.ITERSOLI=15 ;
  113. *
  114. *                     on précise les elements de la table CHIMI2
  115. *              que l'on souhaite conserver
  116. *
  117.  TBPARM.SORTIE= MOTS 'TYP5' 'PREC' 'FION' ;
  118. *
  119.  TB3= CHI2 TB1 TBPARM TBPAR2 ;
  120. *                             CONTROLE DES RESULTATS
  121.  
  122.  FIONTE1=MANU CHPO G1 1 SCAL 1.13636E-03 'NATURE' DISCRET ;
  123.  FIONTE2=MANU CHPO AB 1 SCAL (1.00229E-02 - 1.13636E-03 )
  124.  'NATURE' DISCRET ;
  125.  FIONTES= FIONTE1 + FIONTE2 ;
  126.  VERR1= ( ABS ( FIONTES - TB3.FION )) MASQUE SUPERIEUR 5.D-8 SOMME ;
  127.  PRECTE1= MANU CHPO G1 8 W050 0. W051 0. W052 5.18766E-09 W053 0.
  128.   W054 0. W055 1.99605E-08 W056 0. W057 0.  'NATURE' DISCRET ;
  129.  PRECTE2= MANU CHPO AB 8 W050 0. W051 0. W052 -5.18766E-09 W053 0.
  130.   W054 0. W055 -1.99605E-08 W056 0. W057 0.  'NATURE' DISCRET ;
  131.  PRECTES= PRECTE1+PRECTE2 ;
  132.  PRECD=PRECTES /  50. ;
  133.  VERR2= ( ABS ( PRECTES - TB3.PREC )) MASQUE SUPERIEUR PRECD SOMME ;
  134.  TY5TE1=  MANU CHPO AB 8 W050 8.28223E-16 W051 4.85086E-02
  135.  W052 2.07294E-03 W053 1.31814E-04 W054 1.65945E-02 W055 2.17882E-16
  136.  W056 2.31600E-180    W057 1.16045E-44 'NATURE' DISCRET ;
  137.  TY5TE2=  MANU CHPO G0 8 W050 1.15633E-36 W051 8.45130E-07 W052 0.
  138.  W053 7.27012E-23 W054 9.15253E-21 W055 0. W056 7.36067E-97
  139.  W057 6.18787E-16 'NATURE' DISCRET ;
  140.  TY5TES= TY5TE1+TY5TE2 ;
  141.  TY5D= TY5TES/ 50. ;
  142. *          DANS LE CAS DES CONCENTRATIONS TRES FAIBLES
  143. *    ON ADMETTRA UNE ERREUR PLUS IMPORTANTE SUR LE RESULTAT
  144.  TT1= TY5TES MASQUE DIFFERENT 0. ;
  145.  TTT= TT1 MASQUE INFERIEUR 1.E-40 ;
  146. TY5D= TY5D + ( TTT / 2. ) ;
  147.  VERR3= ( ABS ( TY5TES - TB3.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  148. *OPTION DONN 5 ;
  149. *                                TABLE ENTREE DE CHI1
  150.   TABLIM2= TABLE ;
  151.   TABLIM2. TYP3 = LECT 50 99  ;
  152.   TABDON.CLIM = TABLIM2  ;
  153. *OPTION DONN 5 ;
  154. TB4=CHI1 TABDON
  155.     COMP ('CHAINE' DIVERS '/COMPOM')
  156.     LOGK ('CHAINE' DIVERS '/COMPOM') ;
  157. *
  158. *                                 TABLE ENTREE DE CHI2
  159. *
  160.  TBPAR3= TABLE ;
  161.  TBPAR3.'SOUSTYPE'='DONNEES_CHIMIQUES' ;
  162.  TBPAR2.FIONI= MANU CHPO G1 1 SCAL 0.001 ;
  163.  TBPAR3.LOGC= MANU CHPO G1 11 X005 -3. X004 -11. X103 -3. X099 -5.
  164.          X006 -10. X007 -10. X032 -10. X039 -10. X040 -10. X041 -10.
  165.          X050 -7.  ;
  166.  TBPAR3.FIONI= TB3.FION ;
  167.    PH1=MANU CHPO G1 1 W046 7.D0 ;
  168.    PH2=MANU CHPO AB 1 W046 (5.D0 - 7.D0) ;
  169.     PH=PH1+PH2 ;
  170.    PE1=MANU CHPO G1 1 W045 -7.5D0 ;
  171.    PE2=MANU CHPO AB 1 W045 (12.D0 + 7.5D0) ;
  172.         PE=PE1+PE2 ;
  173.  TBPAR3.CLIM= PH + PE  ;
  174.   TOTE=MANU CHPO G1 1 X099 0.D0 ;
  175.   TOTH=MANU CHPO G1 1 X050 0.D0 ;
  176.  TBPAR3.TOT= TOTNA + TOTK + TOTCL + TOTE + TOTFE3 + TOTFE2 + TOTU6 +
  177.                TOTU5 + TOTU4 + TOTU3 + TOTH ;
  178.  TB5= CHI2 TB4 TBPARM TBPAR3 ;
  179. *                             CONTROLE DES RESULTATS
  180.  
  181.  FIONTE1=MANU CHPO G1 1 SCAL 1.13652E-03 'NATURE' DISCRET ;
  182.  FIONTE2=MANU CHPO AB 1 SCAL (1.00178E-02 - 1.13652E-03 )
  183.  'NATURE' DISCRET ;
  184.  FIONTES= FIONTE1 + FIONTE2 ;
  185.  FIOND= FIONTES/ 2.D5 ;
  186.  VERR4= ( ABS ( FIONTES - TB5.FION )) MASQUE SUPERIEUR FIOND SOMME ;
  187.  PRECTE1= MANU CHPO G1 8  W051 0. W052 0. W053 0. W054 0.
  188.  W055 1.99815E-08  W056 0. W057 0. W058 0.   'NATURE' DISCRET ;
  189.  PRECTE2= MANU CHPO AB 8  W051 0. W052 0. W053 0. W054 0.
  190.  W055 -1.99815E-08  W056 0. W057 0. W058 0.   'NATURE' DISCRET ;
  191.  PRECTES= PRECTE1+PRECTE2 ;
  192.  PRECD=PRECTES /  50. ;
  193.  VERR5= ( ABS ( PRECTES - TB5.PREC )) MASQUE SUPERIEUR PRECD SOMME ;
  194.  TY5TE1=  MANU CHPO AB 8 W051 4.45529E-02  W052 7.02472E-04
  195.  W053  1.13068E-04  W054 1.42344E-02  W055 3.57552E-17
  196.  W056 1.39104E-182  W057 4.39903E-47 W058 1.00004E-16
  197.   'NATURE' DISCRET ;
  198.  TY5TE2=  MANU CHPO G0 8 W051 1.22641E-07 W052 4.63350E-03
  199.  W053 3.16228E-23 W054 3.98107E-21   W055 0. W056  3.89045E-96
  200.  W057  2.69153E-16   W058 2.18776E-37  'NATURE' DISCRET ;
  201.  TY5TES= TY5TE1+TY5TE2 ;
  202.  TY5D= TY5TES/ 50. ;
  203. *          DANS LE CAS DES CONCENTRATIONS TRES FAIBLES
  204. *    ON ADMETTRA UNE ERREUR PLUS IMPORTANTE SUR LE RESULTAT
  205.  TT1= TY5TES MASQUE DIFFERENT 0. ;
  206.  TTT= TT1 MASQUE INFERIEUR 1.E-40 ;
  207. TY5D= TY5D + ( TTT / 2. ) ;
  208.  VERR6= ( ABS ( TY5TES - TB5.TYP5 )) MASQUE SUPERIEUR TY5D SOMME ;
  209. *
  210.  VERR= VERR1+VERR2+VERR3+VERR4+VERR5+VERR6 ;
  211.  SI (VERR EGA 0 ) ;
  212.   ERRE 0 ;
  213.  SINO  ;
  214.   ERRE 5 ;
  215.  FINSI ;
  216. FIN ;
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  

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