*
*...Simple main-program example for using the GRV routines :
*
      PROGRAM MAIN
      IMPLICIT DOUBLE PRECISION (A-H, O-Z)
*
*...TABLES OF X AND Q2 VALUES :
      DIMENSION QS(16), XB(26)
      DATA QS / 1.E0, 2.E0, 5.E0, 1.E1, 2.E1, 5.E1, 1.E2, 2.E2, 
     1          5.E2, 1.E3, 3.E3, 1.E4, 3.E4, 1.E5, 3.E5, 1.E6 /
      DATA XB / 1.E-8, 3.E-8, 1.E-7, 3.E-7, 1.E-6, 3.E-6,
     1          1.E-5, 2.E-5, 5.E-5, 1.E-4, 2.E-4, 5.E-4, 
     1          1.E-3, 2.E-3, 5.E-3, 1.E-2, 2.E-2, 5.E-2,
     2          0.1, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5, 0.7 /
      DATA PI / 3.1415 92654
*
*...INTERPOLATION PROGRAM INITIALIZATION PARAMETER :
      COMMON / INTINIP / IINIP
      COMMON / INTINIF / IINIF
      IINIP = 0
      IINIF = 0
*
*...SELECTION OF THE SET OF DISTRIBUTIONS :
      ISET = 2
*
*...X AND Q2 LOOPS :
      DO 1 K1 = 1, 12
        Q2  = QS(K1)
        IALP = 1
        IF (ISET .GT. 1) IALP = 2 
        ALS = ALPHAS (Q2, IALP) * 4.* PI
        WRITE (6,10) Q2, ALS
  10    FORMAT (/,8X,' Q2(GEV2) =',1PE9.1,4X,' ALPHA_S = '0PF6.4,/)
      DO 2 K2 = 1, 26
        X = XB(K2)
*
*...CALL OF THE INTERPOLATION PROGRAM :
c      CALL GRV98PA (ISET, X, Q2, UV, DV, US, DS, SS, GL)
       CALL GRV98F2 (ISET, X, Q2, F2L, F2C, F2B, F2)
*
*...OUTPUT :
c      WRITE (6,11) X, UV, DV, US, DS, SS, GL
       WRITE (6,11) X, F2L, F2C, F2B, F2
  11   FORMAT (1X,1PE11.3,1X,6(1X,1PE10.3))
   2  CONTINUE
   1  CONTINUE
      STOP
      END
*
      INCLUDE 'grv98.f'