1 : C FZ4 SOURCE CHAT 05/01/13 00:14:20 5004 2 : FUNCTION FZ4(T) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8 (A-H,O-Z) 5 : COMMON/CMCHAM/ D1,D2,R,A,U,DD1,DD2,LSOMD,LD1,LD2 6 : COMMON/CMARGT/SINT,COST,RSINT,RCOST,RSIN2,XA,XU,XXA,XXU,SQA2,SQA1 7 : +,SQU2,SQU1 8 : LOGICAL LSOMD,LD1,LD2 9 : FON(DMI,D,ABSD,SQA,SQU)=SIGN(DMI,D)*LOG(((SQA-ABSD)*(SQU+ABSD))/ 10 : +((SQU-ABSD)*(SQA+ABSD))) 10 : CALL ARGUME(T) 11 : IF(ABS(SINT).LT.0.00001D0) GOTO 11 12 : DMI=0.5D0 13 : ABSD1=ABS(D1) 14 : ABSD2=ABS(D2) 15 : IF(LSOMD) GOTO 3 16 : 4 IF(LD1) GOTO 6 17 : 7 IF(LD2) GOTO 8 18 : GOTO 9 19 : 3 FZ= FON(DMI,D2,ABSD2,SQA2,SQU2)*2.D0 20 : GOTO 10 21 : 6 FZ= FON(DMI,D2,ABSD2,SQA2,SQU2) 22 : GOTO 10 23 : 8 FZ=-FON(DMI,D1,ABSD1,SQA1,SQU1) 24 : GOTO 10 25 : 9 FZ= FON(DMI,D2,ABSD2,SQA2,SQU2)-FON(DMI,D1,ABSD1,SQA1,SQU1) 26 : 10 FZ4=RSINT*FZ 27 : RETURN 28 : 11 FZ4=0.D0 29 : RETURN 30 : END 31 :
© Cast3M 2003 - All rights reserved.
Disclaimer