Télécharger chmest.eso

Retour à la liste

Numérotation des lignes :

chmest
  1. C CHMEST SOURCE CHAT 05/01/12 21:58:59 5004
  2. SUBROUTINE CHMEST(IDX,IONZ,TOT,NXDIM,XMUNEW)
  3. C======================================================================
  4. C ISSU DE TRIOEF(TRIEST)
  5. C OBJET: EVALUATION DE LA FORCE IONIQUE D'APRES LES CONCENTRATIONS
  6. C TOTALES EN COMPOSANTS
  7. C
  8. C APPELLEE PAR CHIMI2
  9. C
  10. C METHODE: LA FORCE IONIQUE EST FORCEE A 0 SI L'ELECTRONEUTRALITE N'EST
  11. C PAS RESPECTEE A PLUS DE 5 %
  12. C
  13. C FORTRAN 77
  14. C
  15. C======================================================================
  16. IMPLICIT INTEGER(I-N)
  17. IMPLICIT REAL*8(A-H,O-Z)
  18. DIMENSION IDX(NXDIM),IONZ(NXDIM),TOT(NXDIM)
  19. XMUNEW=0.D0
  20. VI=0.D0
  21. VJ=0.D0
  22. DO 45 J=1,NXDIM
  23. IDXJ=IDX(J)
  24.  
  25. IF (IDXJ.EQ.99) GOTO 40
  26. IF (IDXJ.EQ.80) GOTO 40
  27. IF (IDXJ.EQ.81) GOTO 40
  28. IF (IDXJ.EQ.90) GOTO 40
  29. IF (IDXJ.EQ.91) GOTO 40
  30. IF (IDXJ.EQ.92) GOTO 40
  31. IF (IDXJ.EQ.93) GOTO 40
  32. IF (IDXJ.EQ.94) GOTO 40
  33. IF (IDXJ.EQ.95) GOTO 40
  34. IF (IDXJ.EQ.96) GOTO 40
  35. C GF=FLOAT(IONZ( J))
  36. GF=IONZ( J)
  37. C
  38. F1=ABS(TOT(J))*GF
  39. IF (F1.LT.0.D0)THEN
  40. VI=VI+F1
  41. XMUNEW=XMUNEW+0.5D0*F1*GF
  42. ELSEIF(F1.GT.0.D0)THEN
  43. VJ=VJ+F1
  44. XMUNEW=XMUNEW+0.5D0*F1*GF
  45. ENDIF
  46. 40 CONTINUE
  47. 45 CONTINUE
  48.  
  49. IF (XMUNEW.NE.0.D0) THEN
  50. F1=(VI+VJ)*2.D0/(ABS(VI)+VJ)
  51. IF (F1.GT.0.05D0)THEN
  52. XMUNEW=0.D0
  53. ENDIF
  54. ENDIF
  55. RETURN
  56. END
  57. C
  58.  
  59.  
  60.  
  61.  

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