CTITLE SARINIT -- COEFFICIENT INITIALIZATION FOR STATIC APPLICATIONS
C***********************************************************************
C                 COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991            *
C***********************************************************************
CA
CA    DESIGNER               J. V. S. HARVEY
CA    LANGUAGE               VS FORTRAN (77)
CA    SYSTEM                 IBM
CA    REWRITTEN              29 MAR 1988
C     REVISED                03-03-92 JJC CHANGED COFGEN TO MCOFGN
CA
CA
CA            CALL SARINIT( IPR, KPBUGF, COEF )
CA
CA
CA    IN/OUT  ARGUMENT  TYPE  DESCRIPTION
CA
CA    IN      IPR         I4  DEBUG OUTPUT UNIT NUMBER
CA    IN      KPBUGF      I4  DEBUG FLAG
CA
CA    OUT     COEF        R4  COEFFICIENTS FOR TRACE INTERPOLATION
CA               ( 2-D ARRAY DIMENSIONED:  12 BY 0:100 )
CA
C***********************************************************************
C
C     SUBROUTINES CALLED:  MCOFGN -- TRACE INTERP. COEF. DETERMINATION
C
C***********************************************************************
C
      SUBROUTINE SARINIT( IPR, KPBUGF, COEF )
C
      REAL COEF(12,0:100)
C
      REAL RW(12), G(14), FF(14), A(14)
C***********************************************************************
C
C                       COEFFICIENT CALCALULATION
C
      CALL MCOFGN( 12, 0.8, RW, G, FF, A, COEF(1,1) )
C
      DO 275 I = 1, 12
         COEF(I,  0) = 0.0
         COEF(I,100) = 0.0
  275    CONTINUE
C
      COEF(6,  0) = 1.0
      COEF(7,100) = 1.0
C======================================================================
C
C                            DEBUG PRINT
C
      IF( KPBUGF .GE. 2 ) THEN
         WRITE( IPR, 8000 )
         DO 900 IFRAC = 0, 50
            WRITE( IPR, 8010 ) IFRAC, ( COEF(I,IFRAC), I=1,12 )
  900       CONTINUE
C
         WRITE( IPR, 8000 )
         DO 950 IFRAC = 51, 100
            WRITE( IPR, 8010 ) IFRAC, ( COEF(I,IFRAC), I=1,12 )
  950       CONTINUE
         ENDIF
      RETURN
C======================================================================
C
C                         FORMAT STATEMENTS
C
 8000 FORMAT('1',45X,'COEFFICIENTS FOR STATICS APPLICATION'//
     *       ' FRAC.'//)
 8010 FORMAT(I5,5X,12(F10.5))
      END