Sources Esope | Cast3M

chmest
C CHMEST    SOURCE    CHAT      05/01/12    21:58:59     5004
      SUBROUTINE CHMEST(IDX,IONZ,TOT,NXDIM,XMUNEW)
C======================================================================
C       ISSU DE TRIOEF(TRIEST)
C OBJET: EVALUATION DE LA FORCE IONIQUE D'APRES LES CONCENTRATIONS
C TOTALES EN COMPOSANTS
C
C APPELLEE PAR CHIMI2
C
C METHODE: LA FORCE IONIQUE EST FORCEE A 0 SI L'ELECTRONEUTRALITE N'EST
C PAS RESPECTEE A PLUS DE 5 %
C
C   FORTRAN 77
C
C======================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION IDX(NXDIM),IONZ(NXDIM),TOT(NXDIM)
      XMUNEW=0.D0
      VI=0.D0
      VJ=0.D0
      DO 45 J=1,NXDIM
      IDXJ=IDX(J)
 
      IF (IDXJ.EQ.99) GOTO 40
      IF (IDXJ.EQ.80) GOTO 40
      IF (IDXJ.EQ.81) GOTO 40
      IF (IDXJ.EQ.90) GOTO 40
      IF (IDXJ.EQ.91) GOTO 40
      IF (IDXJ.EQ.92) GOTO 40
      IF (IDXJ.EQ.93) GOTO 40
      IF (IDXJ.EQ.94) GOTO 40
      IF (IDXJ.EQ.95) GOTO 40
      IF (IDXJ.EQ.96) GOTO 40
C     GF=FLOAT(IONZ( J))
      GF=IONZ( J)
C
      F1=ABS(TOT(J))*GF
      IF (F1.LT.0.D0)THEN
            VI=VI+F1
            XMUNEW=XMUNEW+0.5D0*F1*GF
      ELSEIF(F1.GT.0.D0)THEN
            VJ=VJ+F1
            XMUNEW=XMUNEW+0.5D0*F1*GF
      ENDIF
   40 CONTINUE
   45 CONTINUE
 
      IF (XMUNEW.NE.0.D0) THEN
            F1=(VI+VJ)*2.D0/(ABS(VI)+VJ)
            IF (F1.GT.0.05D0)THEN
                  XMUNEW=0.D0
            ENDIF
      ENDIF
      RETURN
      END
C