Télécharger lispal.eso

Retour à la liste

Numérotation des lignes :

lispal
  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.  
  71.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales