Sources Esope | Cast3M

chmmod
C CHMMOD    SOURCE    CHAT      05/01/12    21:59:32     5004
      SUBROUTINE CHMMOD(IDSCHI,XMU,XMUNEW,TMPNEW,GNEW)
C=======================================================================
C            ISSU DE TRIOEF(TRIMOD)
C
C - OBJET: MODIFICATION DES LOGK POUR UNE FORCE IONIQUE DONNEE (DE XMU A
C XMUNEW) A TEMPERATURE FIXEE TMP.
C
C
C - APPELLEE PAR CHMKMD()
C
C - METHODES: LE COEFFICIENT GFO POUR UN ION UNIVALENT EST CELUI DE
C L'APPROXIMATION DE DAVIES PRENANT EN COMPTE LA TEMPERATURE; ON UTILISE
C UNE DEPENDANCE LINEAIRE EN TEMPERATURE POUR LA CONSTANTE DIELECTRIQUE
C DE L'EAU
C
C LES ESPECES DE TYPE III SONT PRISES EN COMPTE DE TELLE SORTE QUE LEUR
C ACTIVITE INITIALE SOIT MAINTENUE.
C
C======================================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
      SEGMENT IDSCHI
           REAL*8 GK(NYDIM),AA(NYDIM,NXDIM),FF(NZDIM,NPDIM)
           INTEGER IDX(NXDIM),IDY(NYDIM),IDZ(NZDIM),IDP(NPDIM),NN(6)
           INTEGER IDECY(NYDIM),IONZ(NXDIM)
           CHARACTER*32 NAME(NXDIM),NAMESP(NYDIM)
      ENDSEGMENT
 
 
C     write(6,*)' XMUNEW    XMU ',XMUNEW,XMU
      IF (XMUNEW.EQ.XMU) RETURN
C
      NC=NN(1)+NN(2)+NN(3)+NN(4)+NN(5)+NN(6)
      I1=NN(1)+NN(2)+1
      I2=NN(1)+NN(2)+NN(3)
      NX=IDX(/1)
C
C
        T1=TMPNEW+273.15D0
      IF(TMPNEW.LE.100.D0) THEN
        TE=TMPNEW-25.D0
        ET=78.54D0*(1.D0-4.579D-3*TE+1.19D-5*TE*TE
     *  -2.8D-8*TE*TE*TE)
      ELSEIF((TMPNEW.GT.100.D0).AND.(TMPNEW.LE.370.D0)) THEN
        ET=5321.D0/T1+233.76D0-0.9297D0*T1+1.417D-3*T1*T1
     *     -8.292D-7*T1*T1*T1
      ENDIF
C
C
C
      AT=1.82D6/SQRT((ET*T1)**3)
      SI=SQRT(XMU)
      GOLD=-AT*(SI/(1.D0+SI)-0.2D0*XMU)
      SI=SQRT(XMUNEW)
      GNEW=-AT*(SI/(1.D0+SI)-0.2D0*XMUNEW)
      IF (ABS(GNEW).EQ.0.D0) GNEW=0.D0
      DGF=GNEW-GOLD
C      WRITE(6,*)'  '
C      WRITE(6,*)' CHMMOD ', ' NC ',NC,' NX ',NX
C      WRITE(6,*)' ET ',ET,' T1 ',T1,' AT ',AT,' SI ',SI
C      WRITE(6,*)' DGF ',DGF,' GNEW ',GNEW
      DO 20 I=1,NC
         VJ=0.D0
         VI=0.D0
         DO 30 J=1,NX
            VJ=VJ+AA(I,J)*IONZ( J)*IONZ( J)
            VI=VI+AA(I,J)*IONZ( J)
30       CONTINUE
         GKT=GK(I)+DGF*(VJ-VI*VI)
C        write(6,*)'idy',idy(i),'gkt',gkt,gk(i),vi,vj,i,idecy(i)
         IF (I.GE.I1.AND.I.LE.I2) GKT=GK(I)+DGF*VJ
         IF ( IDECY(I)   .EQ. 1 ) GKT=GK(I)+DGF*VJ
         GK(I)=GKT
20    CONTINUE
      XMU=XMUNEW
      RETURN
      END
C
* LE COEFFICIENT D'ACTIVITE D'UN COMPLEXE DE SURFACE EST CONSTANT
* ET EGAL A 1/CONCENTRATION DE SITES DE SURFACE
C
C
C
C