Télécharger la plaquette Cast3M


fr en

Télécharger chmele.eso

Retour à la liste

Numérotation des lignes : 

chmele

  1. C CHMELE    SOURCE    CHAT      05/01/12    21:58:55     5004
  2.       SUBROUTINE CHMELE(IDSCHI,LGKMOD,LGKTMP,SP2,IOXYDO,Q99,PEMIN,
  3.      * PEMAX,DE,PE)
  4. C
  5. C- VERIFICATION DE LA PRESENCE DE L'ELECTRON PARMI LES COMPOSANTS
  6. C Modif PhM : Possibilite d'utiliser d'autre composants que l'electron
  7. C pour les reactions redox.
  8. C L'identifiant du composant redox est forcement 99.
  9. C On fait les distinction suivante en affectant IOXYDO
  10. C                IOXYDO = 0 -> Le composant redox (99) n'est pas present
  11. C                              donc pas de raction redoxPas de reaction redox
  12. C                IOXYDO = 1 -> Presence du composant redox en type
  13. C                              3. Reaction redox avec activite fixe
  14. C                              pour le composant redox (Pe fixe par ex )
  15. C                IOXYDO = 2 -> Presence du composant redox en type >3,.
  16. C                              On considere que c'est l'electron mis en
  17. C                              type 6. On le place en type 3 et
  18. c                              on utilisera la dichotomie.
  19. C                IOXYDO = 3 -> Presence du composant redox mais en type
  20. C                              1. Il y a reaction redox avec un autre
  21. C                              composantque l'electron (pas de
  22. C                              dichotomie,presence de ce composant
  23. C                               dans les bilans de masse)
  24. C
  25.       IMPLICIT INTEGER(I-N)
  26.       IMPLICIT REAL*8(A-H,O-Z)
  27.       SEGMENT IDSCHI
  28.            REAL*8 GK(NYDIM),AA(NYDIM,NXDIM),FF(NZDIM,NPDIM)
  29.            INTEGER IDX(NXDIM),IDY(NYDIM),IDZ(NZDIM),IDP(NPDIM),NN(6)
  30.            INTEGER IDECY(NYDIM),IONZ(NXDIM)
  31.            CHARACTER*32 NAME(NXDIM),NAMESP(NYDIM)
  32.       ENDSEGMENT
  33.       SEGMENT SP2
  34.            REAL*8 GX(NXDIM),XX(NXDIM),GS(NZDIM),SS(NZDIM)
  35.            REAL*8 TOT(NXDIM),TOTAQ(NXDIM),TOTFIX(NXDIM),GKS(NZDIM)
  36.            REAL*8 YY(NXDIM),ZZ(NXDIM,NXDIM),CC(NYDIM),GC(NYDIM)
  37.       ENDSEGMENT
  38.       IOXYDO = 0
  39.       NXDIM=IDX(/1)
  40.       NYDIM=IDY(/1)
  41.       NZDIM=IDZ(/1)
  42.       NPDIM=IDP(/1)
  43.       CALL CHIADY(IDX,NXDIM,99,IPOX99)
  44.       IF(IPOX99.EQ.0)RETURN
  45.       Q99 = TOT(IPOX99)
  46. C
  47. C- MISE EN TYPE 3 (IE PE IMPOSÉ) DE L'ELECTRON SI IL EST EN TYPE 6
  48. C- (IE Q IMPOSE) AFIN DE FAIRE UNE MÉTHODE DE DICHOTOMIE.
  49. C- IL FAUT SORTIR L'ELECTRON DU SYSTEME AFIN D'AMELIORER LA STABILITE
  50. C- NUMERIQUE DU SCHEMA. C'EST CE QUI SE PASSE SI L'E- EST EN TYPE 3,
  51. C- D'OU L'IDEE D'UN DICHOTOMIE SUR LE PE POUR RETROUVER LE Q IMPOSE.
  52. C
  53.       CALL CHIADY(IDY,NYDIM,99,IPOS99)
  54.       N3FIN = NN(1) + NN(2) + NN(3)
  55.       IF (IPOS99.GT.N3FIN) THEN
  56.          IOXYDO = 2
  57.          Q99    = TOT(IPOX99)
  58.          ID     = 99
  59.          LIN    = 6
  60.          LEN    = 3
  61.          CALL CHMREX(IDSCHI,LGKMOD,LGKTMP,ID,LIN,LEN)
  62.          N3DEB  = NN(1) + NN(2) + 1
  63.          N3FIN  = NN(1) + NN(2) + NN(3)
  64.          CALL CHIADY(IDY,NYDIM,99,IPOS99)
  65.  
  66.          PEMIN  = -GX(IPOX99) -0.5D0*DE
  67.          PEMAX  = PEMIN + DE
  68.  
  69.          PE     = PEMAX
  70.          CALL CHIADY(IDX,NXDIM,99,IDXT)
  71.          GK(IPOS99)  = PE
  72.          XX(IDXT)   = 10.D0 ** GX(IDXT)
  73.          TOT(IDXT) = 0.D0
  74. C
  75.       ELSEIF(IPOS99.GT.(NN(1)+NN(2)))THEN
  76.          IOXYDO=1
  77. C MODIF PhM
  78.       ELSEIF (IPOS99.LE.NN(1)) THEN
  79.          IOXYDO=3
  80. C MODIF PhM
  81.       ENDIF
  82.       RETURN
  83.       END
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  

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