**********************************************************
      program example
**********************************************************
*
*    calculation of LSS'05 NLO(LO) PPDs and g1(p,n,d)
*
*     example program
*
**********************************************************
*
c      implicit none
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
*
      integer nq2,nxbj
      parameter (nq2=1 , nxbj=25)
*
      integer iini,iset
      integer iq2,ixbj
*
      double precision q2ref(nq2)
      double precision xbj(nxbj)
      double precision x,q2
      double precision uv,dv,gl,sea,qb,g1pLT,g1p,g1nLT,g1n,g1d
      double precision duv,ddv,dgl,dsea,dqb,dg1pLT,dg1p,dg1nLT,dg1n,dg1d
*    
      common /intini/ iini

*
*     set data:
*
        data q2ref  / 3.0d0 /


       DATA xbj /
     1           
     3           1.D-2, 1.5D-2, 2.2D-2, 3.2D-2, 5.0D-2, 7.5D-2,
     4           0.1, 0.125, 0.15, 0.175, 0.2, 0.225, 0.25, 0.275,
     5           0.3, 0.325, 0.35, 0.375, 0.4, 0.45,  0.5, 0.55,
     6           0.6, 0.65,0.7 /

      OPEN(6,FILE='test_LSS2005_q3.dat'
     1 ,STATUS='NEW',FORM='FORMATTED')

*
*     start reading data
*
      iini = 0                  ! reading data first time
*
*     choose order of parton sets and fit mode
*            ISET = 1   Set 1 NEXT-TO-LEADING ORDER (MS-bar)        *
*                      (DATA FILE 'Set1_NLO_MS.grid' UNIT=11,       *
*                       TO BE DEFINED BY THE USER )                 *
*            ISET = 2   Set 1 NEXT-TO-LEADING ORDER (JET)           *
*                      (DATA FILE 'Set1_NLO_JET.grid' UNIT=22,      *
*                       TO BE DEFINED BY THE USER )                 *
*            ISET = 3   Set 2 LEADING ORDER                         *
*                      (DATA FILE 'Set2_LO.grid' UNIT=33,           *
*                       TO BE DEFINED BY THE USER )                 *
*            ISET = 4   Set 2 NEXT-TO-LEADING ORDER (MS-bar)        *
*                      (DATA FILE 'Set2_NLO_MS.grid' UNIT=44,       *
*                       TO BE DEFINED BY THE USER )                 *
*            ISET = 5   Set 2 NEXT-TO-LEADING ORDER (JET)           *
*                      (DATA FILE 'Set2_NLO_JET.grid' UNIT=55,      *
*                       TO BE DEFINED BY THE USER )             


*
*     g1 and pdf vs x at Q2ref
*


*         MS-bar
           iset=4
  
      PRINT '(''iset = '' ,i3)'
     1 ,iset	  
         write(6,2000)	 
	   
	   q2=q2ref(1)  
         do ixbj = 1,25
*
            x = xbj(ixbj)

c     To obtain the values of PPDs and g1 outside the 
c     resonanse region use the cut W**2>4 GeV**2  
    
      call LSS2005 (ISET,X,Q2,UUB,DDB,SSB,GL,UV,DV,UB,DB,ST,
     1g1pLT,g1p,g1nLT,g1n)

          g1d = 0.5d0 * (g1p+g1n) * (1.0d0-1.5d0*0.05d0)

        write(6,2001) q2,x,UUB,DDB,SSB,GL,UV,DV,UB,DB,ST
	1,g1pLT,g1p,g1nLT,g1n
*

         enddo


 2000 format('* Q2 x UUB DDB SSB GL UV DV UB DB ST g1pLT g1p g1nLT g1n')
 2001   FORMAT (2e9.3,13(1pe12.4))
*
      end

*


      include 'LSS2005pdf_g1.f'        

*=======================================================================