1 : C CHAREM SOURCE CHAT 05/01/12 21:55:52 5004 2 : SUBROUTINE CHAREM(IDIAM,R0,RM,B,R,EPSTAR,EPST,DEPS,SIGEL, 3 : .X1,X2,S,IBOU,ICENT2,JA,ALPHA1,ALPHA2,PHI,PSI,OME) 4 : IMPLICIT INTEGER(I-N) 5 : IMPLICIT REAL*8(A-H,O-Z) 6 : DIMENSION SIGEL(*),S(*),X1(*),X2(*),ALPHA1(*),ALPHA2(*) 7 : EPST=EPSTAR 8 : DEPS=0.D00 9 : IF(IDIAM.NE.0) R=RM-(RM-R0)*EXP(-B*EPSTAR) 10 : PHI=1.D00 11 : IF(PSI.NE.1.D00) PHI=1.D00+(PSI-1.D00)*EXP(-OME*EPSTAR) 12 : JX=JA 13 : DO 1 I=1,IBOU 14 : S(I)=SIGEL(I) 15 : X1(I)=ALPHA1(JX) 16 : IF(ICENT2.EQ.0) GO TO 1 17 : X1(I)=X1(I)-ALPHA2(JX) 18 : X2(I)=ALPHA2(JX) 19 : 1 JX=JX+1 20 : RETURN 21 : END 22 :
© Cast3M 2003 - All rights reserved.
Disclaimer