fsila.for

Для каждой характеристики генерируется структурной похожий кусок кода. Пример кода для ГПР (гидропневматической рессоры):

...      
      REAL*8 FUNCTION FLX020(X,Y,XO,YO)
      IMPLICIT REAL*8 (A-H,O-Z)
      COMMON /PAR/ gv(90000)
C    ГПР1'                                   
      EQUIVALENCE (PAR01,GV(  1822))
      EQUIVALENCE (PAR02,GV(  1823))
      EQUIVALENCE (PAR03,GV(  1824))
      EQUIVALENCE (PAR04,GV(  1825))
      EQUIVALENCE (PAR05,GV(  1826))
      EQUIVALENCE (PAR06,GV(  1827))
      EQUIVALENCE (PAR07,GV(  1828))
      EQUIVALENCE (PAR08,GV(  1829))
      EQUIVALENCE (PAR09,GV(  1830))
      EQUIVALENCE (PAR10,GV(  1831))
      EQUIVALENCE (PAR11,GV(  1832))
      EQUIVALENCE (PAR12,GV(  1833))
      IF(Y.GE.0.0) GOTO 100
      IF(Y.GE.PAR05) GOTO 101
      YO=Y*PAR01+PAR07
      GOTO 150
 101  CONTINUE
      YO=Y*PAR02
      GOTO 150
 100  CONTINUE
      IF(Y.GE.PAR06) GOTO 102
      YO=Y*PAR03
      GOTO 150
 102  CONTINUE
      YO=Y*PAR04+PAR08
 150  CONTINUE
      DLT=PAR09-X*PAR11
      IF(DLT.LE.0.0)CALL OTKAZ(13,1,1)
      XO=PAR10*PAR11*(PAR09/(PAR09-X*PAR11))**PAR12
      FLX020=XO+YO
      IF(X.LE.0.0)XO=0.0
      FLX020=XO+YO
      RETURN
      END
...

В массиве GV сначала описываются параметры диссипативной функции:

а затем параметры упругой функции

Адреса для каждой модели свои!