1 : C INTFIL SOURCE PV 05/04/25 21:15:30 5077 2 : SUBROUTINE INTFIL(GP1,GP2,DKSI,TETA1,TETA2,DK1,DK2) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : C 6 : C======================================================================= 7 : C SOUS-PROGRAMME FORTRAN QUI CALCULE LES TERMES DES MATRICES ANBN 8 : C 9 : C TETA2 DK2 10 : C 1 / / SIN (K * (TETA-DKSI)) 11 : C GP1(T) = ___ / / _____________________ DK DTETA 12 : C PI / / K 13 : C TETA1 DK1 14 : C TETA2 DK2 15 : C 1 / / 16 : C GP2(T) = ___ / / COS(K*(TETA-DKSI)) DK DTETA 17 : C PI / / 18 : C TETA1 DK1 19 : C 20 : C 21 : C----------------------------------------------------------------------- 22 : C CREATION: J-M GAUTIER AVRIL 87 23 : C======================================================================= 24 : C 25 : -INC CCREEL 26 : -INC CCOPTIO 27 : C 28 : IF (DK1.EQ.DK2) RETURN 29 : IF ((DK1.LT.0.D0).OR.(DK2.LT.0.D0)) RETURN 30 : C----------------------------------------------------------------------- 31 : C INITIALISATION 32 : C----------------------------------------------------------------------- 33 : SINT1=0 34 : SINT2=0 35 : SINT3=0 36 : SINT4=0 37 : C 38 : GP1=0 39 : GP2=0 40 : C 41 : DELTET=TETA2-TETA1 42 : C 43 : DO 10 IETA=-1,1,2 44 : TETKS1=TETA1+IETA*DKSI 45 : TETKS2=TETA2+IETA*DKSI 46 : C----------------------------------------------------------------------- 47 : X1=TETKS2*DK1 48 : X2=TETKS2*DK2 49 : X3=TETKS1*DK2 50 : X4=TETKS1*DK1 51 : C----------------------------------------------------------------------- 51 : CALL SINI(X1,SINT1) 51 : CALL SINI(X2,SINT2) 51 : CALL SINI(X3,SINT3) 51 : CALL SINI(X4,SINT4) 52 : C----------------------------------------------------------------------- 53 : C CALCUL DE GP1 ET GP2 54 : C 55 : GP1=GP1+SINT2-SINT1 56 : GP2=GP2+SINT2-SINT3-SINT1+SINT4 57 : 10 CONTINUE 58 : C 59 : XPI2=2.D0*XPI 60 : GP1=DELTET*GP1/XPI2 61 : GP2=GP2/XPI2 62 : RETURN 63 : END 64 : 65 : 66 : 67 :
© Cast3M 2003 - All rights reserved.
Disclaimer