* WTO.22 version of October 18 1999
*
* Needs the following files
*
* WTO_22_IO
* WTO_22_MAIN
* WTO_22_CCA
* WTO_22_NCA
* WTO_22_HCA
* WTO_22_SWA.FOR
* WTO_22_LIB
* WTO_22_JETSET.FOR
*
* To be linked with NAGLIB
*
*-----------------------------------------------------------------------
* For help use call wtohelp() and look for wto_help.dat
*-----------------------------------------------------------------------
*
      implicit real*8(a-h,o-z)
      character*1 oxc,oval
      character*2 oflag2
      character*3 oflag3
      character*4 oflag4
      character*5 oflag5
      character*6 oflag6
*
      parameter(nin=21)
*
      dimension afsm(4)
*
*-----------------------------------------------------------------------
* Input parameters
*-----------------------------------------------------------------------
*
      ipro= 5
      oxc= 'c'
      rs= 183.d0
      wmi= 80.350d0
      zmi= 91.1867d0
      zgi= 2.4974d0
      tqmi= 173.8d0
      hmi= 100.d0
      alszi= 0.12d0
*
*----------------------------------------------------------------------
*
      call wtoinit(ipro,oxc,rs,wmi,zmi)
      oflag3= 'oww'
      oval= 'f'
      call wtocflag(rs,oflag3,oval,-1,-1,-1.d0,-1.d0)
      oflag3= 'ozz'
      oval= 'f'
      call wtocflag(rs,oflag3,oval,-1,-1,-1.d0,-1.d0)
      oflag5= 'osinw'
      oval= 'c'
      call wtocflag(rs,oflag5,oval,-1,-1,-1.d0,-1.d0)
      oflag3= 'ipm'
      call wtocopt(oflag3,6)
      oflag3= 'irm'
      call wtocopt(oflag3,6)
      oflag4= 'iosf'
      call wtocopt(oflag4,0)
      asccsw= 0.997d0
      eusw= 15.d0
      edsw= 15.d0
      alsw= 1.d0
      sthqsw= 0.d0
      call wtoswcutset(asccsw,sthqsw,eusw,edsw,alsw)
      call wto(rs,wmi,zmi,zgi,tqmi,hmi,alszi,sig,esig)
*
      print 1,sig,esig
      re= 1.d2*esig/sig
      print 2,re
    1 format(' sigma (nb)    =',e20.5,3x,'+-',3x,e20.5)
    2 format(/' Rel. error of ',f10.3,' %',///)
*
*-----------------------------------------------------------------------
*
      stop
      end