Esope source files

   1 : C LISPFI    SOURCE    CHAT      05/01/13    01:22:45     5004           
   2 : C
   3 : C    CALCULE LES FACTEURS DE FORME POUR LE LINESPRING                           
   4 :       SUBROUTINE LISPFI(ASURW,FM,FF)                                            
   5 : C=================================================================              
   6 : C    CALCULE LES FACTEURS DE FORME POUR LE LINESPRING                           
   7 : C                                                                               
   8 : C  ENTREE:
   9 : C  ASURW   = RAPPORT PROFONDEUR DE FISSURE SUR EPAISSEUR  
  10 : C                      
  11 : C  SORTIE:                                                                     
  12 : C  FM  FF = FACTEURS DE FORME                                                  
  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     INITIALISATION DES COEFFICIENTS CMM CMF CFF                 
  23 : C     REFERENCE TABLEAU 4                        
  24 : C     NOTE TECHNIQUE DRE/STRE/LMA 83/527                            
  25 : C     ACTION P.W.S. 3.79                                              
  26 : C     CHEISSOUX  IVARS                                                 
  27 : C----------------------------------------------------------
  28 : C      DIMENSION CM(6),CF(6)  
  29 : C      DATA CM/1.12,4.5160,-41.62,179.3,-298.4,193.70/              
  30 : C      DATA CF/1.12,0.5038,-13.06,63.71,-108.6,070.47/       
  31 : C      FM=0.                                                        
  32 : C      FF=0.                                                         
  33 : C      DO 130 IA=1,6                                 
  34 : C      IB=7-IA                                                 
  35 : C      FM=CM(IB)+FM*ASURW                                           
  36 : C      FF=CF(IB)+FF*ASURW                                               
  37 : C  130 CONTINUE                                                     
  38 : C-----------------------------------------------------------
  39 :       DIMENSION CM(7),CF(7)                                             
  40 :       DATA CM/1.12,6.52,-12.39,89.05,-188.61,207.39,-32.05/        
  41 :       DATA CF/1.12,-1.89,18.01,-87.39,241.90,-319.94,168.01/
  42 :       FM=0.                                                           
  43 :       FF=0.                                                               
  44 :       DO 130 I1=1,7                 
  45 :       FM = FM + ((ASURW**((I1-1)*2))*CM(I1))
  46 :       FF = FF + ((ASURW**((I1-1)*1))*CF(I1))
  47 :   130 CONTINUE                                                      
  48 : C      WRITE (6,*) 'FM = ',FM
  49 : C      WRITE (6,*) 'FF = ',FF
  50 : C-----------------------------------------------------------
  51 : C
  52 :       RETURN                                                        
  53 :       END                                                          
  54 :