Sources Esope | Cast3M

lispfi
C LISPFI    SOURCE    CHAT      05/01/13    01:22:45     5004
C
C    CALCULE LES FACTEURS DE FORME POUR LE LINESPRING
      SUBROUTINE LISPFI(ASURW,FM,FF)
C=================================================================
C    CALCULE LES FACTEURS DE FORME POUR LE LINESPRING
C
C  ENTREE:
C  ASURW   = RAPPORT PROFONDEUR DE FISSURE SUR EPAISSEUR
C
C  SORTIE:
C  FM  FF = FACTEURS DE FORME
C
C   EBERSOLT MARS 85
C   UHLMANN  OCTOBRE 94
C=================================================================
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
C----------------------------------------------------------
C     INITIALISATION DES COEFFICIENTS CMM CMF CFF
C     REFERENCE TABLEAU 4
C     NOTE TECHNIQUE DRE/STRE/LMA 83/527
C     ACTION P.W.S. 3.79
C     CHEISSOUX  IVARS
C----------------------------------------------------------
C      DIMENSION CM(6),CF(6)
C      DATA CM/1.12,4.5160,-41.62,179.3,-298.4,193.70/
C      DATA CF/1.12,0.5038,-13.06,63.71,-108.6,070.47/
C      FM=0.
C      FF=0.
C      DO 130 IA=1,6
C      IB=7-IA
C      FM=CM(IB)+FM*ASURW
C      FF=CF(IB)+FF*ASURW
C  130 CONTINUE
C-----------------------------------------------------------
      DIMENSION CM(7),CF(7)
      DATA CM/1.12,6.52,-12.39,89.05,-188.61,207.39,-32.05/
      DATA CF/1.12,-1.89,18.01,-87.39,241.90,-319.94,168.01/
      FM=0.
      FF=0.
      DO 130 I1=1,7
      FM = FM + ((ASURW**((I1-1)*2))*CM(I1))
      FF = FF + ((ASURW**((I1-1)*1))*CF(I1))
  130 CONTINUE
C      WRITE (6,*) 'FM = ',FM
C      WRITE (6,*) 'FF = ',FF
C-----------------------------------------------------------
C
      RETURN
      END