Télécharger chmion.eso

Retour à la liste

Numérotation des lignes :

  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