1 : C LISPAL SOURCE CHAT 05/01/13 01:22:29 5004 2 : C CALCULE LES COEFFS ALPHA DU LINESPRING 3 : SUBROUTINE LISPAL(ASURW,ALMM,ALMF,ALFF,DELTA) 4 : C 5 : C================================================================ 6 : C 7 : C ENTREES: 8 : C ASURW = RAPPORT PROFONDEUR DE LA FISSURE SUR EPAISSEUR 9 : C 10 : C SORTIES: 11 : C ALMM ALMF ALFF COEFF ALPHA 12 : C DELTA = ALMM * ALFF - ALMF * ALMF 13 : C 14 : C EBERSOLT MARS 85 15 : C UHLMANN OCTOBRE 94 16 : C================================================================ 17 : C 18 : IMPLICIT INTEGER(I-N) 19 : IMPLICIT REAL*8(A-H,O-Z) 20 : C 21 : C------------------------------------------------------------- 22 : C DIMENSION CMM(5),CMF(5),CFF(5) 23 : C DATA CMM/24.00,-232.0,896.1,-1470.0,922.6/ 24 : C DATA CMF/10.02,-86.71,336.0,-549.60,346.8/ 25 : C DATA CFF/4.817,-33.39,127.2,-206.10,130.3/ 26 : C DATA XZER/0.D0/ 27 : C ALMM=XZER 28 : C ALMF=XZER 29 : C ALFF=XZER 30 : C DO 120 IA=1,5 31 : C IB=6-IA 32 : C ALMM=ALMM*ASURW+CMM(IB) 33 : C ALMF=ALMF*ASURW+CMF(IB) 34 : C ALFF=ALFF*ASURW+CFF(IB) 35 : C 120 CONTINUE 36 : C------------------------------------------------------------- 37 : DIMENSION CMM(13),CMF(19),CFF(13) 38 : DATA CMM/1.9761,11.4870,7.7086,15.0143,280.1207,-1099.72, 39 : 1 3418.9795,-7686.9237,12794.1279,-13185.0403,7868.2682, 40 : 1 -1740.2463,124.1360/ 41 : DATA CMF/1.9735,-2.2166,21.6051,-69.3133,196.3,-406.2608, 42 : 1 644.9350,-408.9569,-159.6927,-988.9879,4266.5487, 43 : 1 -2997.1408,-6050.7849,8855.3615,3515.4345,-11744.1116, 44 : 1 4727.9784,1695.6087,-845.8958/ 45 : DATA CFF/1.9710,-4.4277,34.4952,-165.7321,626.3926, 46 : 1 -2144.4651,7043.4169,-19003.2199,37853.3028,-52595.4681, 47 : 1 48079.2948,-25980.1559,6334.2425/ 48 : DATA XZER/0.D0/ 49 : ALMM=XZER 50 : ALMF=XZER 51 : ALFF=XZER 52 : DO 120 I1=1,13 53 : ALMM = ALMM + ((ASURW**((I1-1)*2))*CMM(I1)) 54 : 120 CONTINUE 55 : DO 121 I2=1,19 56 : ALMF = ALMF + ((ASURW**(I2-1))*CMF(I2)) 57 : 121 CONTINUE 58 : DO 122 I3=1,13 59 : ALFF = ALFF + ((ASURW**(I3-1))*CFF(I3)) 60 : 122 CONTINUE 61 : C------------------------------------------------------------- 62 : C 63 : ALMM=ALMM*ASURW*ASURW 64 : ALMF=ALMF*ASURW*ASURW 65 : ALFF=ALFF*ASURW*ASURW 66 : DELTA=ALMM*ALFF-ALMF*ALMF 67 : C 68 : RETURN 69 : END 70 :
© Cast3M 2003 - All rights reserved.
Disclaimer