Télécharger la plaquette Cast3M


fr en

Télécharger chmion.eso

Retour à la liste

Numérotation des lignes : 

chmion

  1. C CHMION    SOURCE    CHAT      05/01/12    21:59:14     5004
  2.       SUBROUTINE CHMION(IDSCHI,SP2,XMUNEW)
  3. C=======================================================================
  4. C            ISSU DE TRIO-EF (TRICAL)
  5. C
  6. C - OBJET: CALCUL DE FORCE IONIQUE
  7. C
  8. C - ARGUMENTS: POINTEURS DES SEGMENTS IDSCHI SP2
  9. C              XMUNEW = NOUVELLE FORCE IONIQUE
  10. C
  11. C - METHODES: CALCULE LA FORCE IONIQUE D'APRES LES CHARGES DES COMPOSANTS
  12. C ET LES VALEURS DE CONCENTRATIONS "C" DES ESPECES DISSOUTES DE TYPE I ET
  13. C II
  14. C
  15. C PEUT ETRE UTILISE POUR AFFINER LA CORRECTION DES LOGK AU SEIN D'UNE
  16. C BOUCLE, EN FONCTION D'UN CRITERE:
  17. C
  18. C               10 CALL CHMKMD()
  19. C                  ...
  20. C                  CALL CHMION()
  21. C                  ...CRITERE DE SORTIE DE BOUCLE
  22. C                  GOTO 10
  23. C
  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. C
  39.       I1=NN(1)+1
  40.       I2=NN(1)+NN(2)
  41.       NXDIM=IDX(/1)
  42.       XMUNEW=0.D0
  43.       DO 50 I=I1,I2
  44.       VI=0.D0
  45.       DO 51 J=1,NXDIM
  46.       VI=VI+AA(I,J)*IONZ(J)
  47. 51    CONTINUE
  48.       XMUNEW=XMUNEW+VI*VI*CC(I)
  49. 50    CONTINUE
  50.  
  51. *
  52. * SEULES LES ESPECES DE TYPE 1 PARTICIPENT A LA FORCE IONIQUE
  53. * L'ESPECE SIMPLE D'UN COMPOSANT MIS EN TYPE 6 N'INTERVIENT PAS
  54. * (CF ELECTRON ,SITE DE SURFACE
  55.       DO 54 J=1,NN(1)
  56.       IDYJ=IDY(J)
  57.       CALL CHIADY(IDX,NXDIM,IDYJ,IYJ)
  58.       VI=IONZ( IYJ)*IONZ( IYJ)
  59.       XMUNEW=XMUNEW+VI*ABS(CC(J))
  60.  54   CONTINUE
  61.       XMUNEW=XMUNEW/2.D0
  62.       RETURN
  63.       END
  64. C
  65. C
  66. C
  67. C
  68.  
  69.  
  70.  
  71.  
  72.  

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