1 : C DRUCO2 SOURCE CHAT 05/01/12 22:55:17 5004 2 : C DRUCO2 SOURCE JFG 02/07/04 3 : SUBROUTINE DRUCO2(S1,SEQ,DPHI,DLAM,VEC,BETINSA) 4 : C 5 : IMPLICIT INTEGER(I-N) 6 : IMPLICIT REAL*8(A-H,O-Z) 7 : DIMENSION S1(6),VEC(6) 8 : C 9 : SEGMENT BETINSA 10 : REAL*8 RT,RC,YOUN,XNU,GFT,GFC,CAR 11 : REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC 12 : INTEGER IFIS,IPLA,IBB,IGAU1 13 : ENDSEGMENT 14 : C 15 : BETA = 1.2D0 16 : ALC=(BETA-1.)/(2.*BETA-1.) 17 : BEC=BETA/(2.*BETA-1.) 18 : C 19 : AC = YOUN / (1.- 2.*XNU) 20 : AD = (3 * DLAM * ALC * AC ) 21 : C 22 : DPHI = BEC * SEQ - ALC*(S1(1)+S1(2)+S1(3)+AD) 23 : C 23 : CALL ZERO(VEC,6,1) 24 : VEC(1) = ALC * AC 25 : VEC(2) = ALC * AC 26 : VEC(3) = ALC * AC 27 : C 28 : RETURN 29 : END 30 : 31 : 32 : 33 :
© Cast3M 2003 - All rights reserved.
Disclaimer