* WTO.21 version of July 24 1997
*
* Needs the following files
*
* WTO_21_IO
* WTO_21_MAIN
* WTO_21_CCA
* WTO_21_NCA
* WTO_21_HCA
* WTO_21_LIB
*
* 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= 14
      oxc= 'n'
      rs= 190.d0
      wmi= 80.23d0
      zmi= 91.1888d0
      zgi= 2.4974d0
      tqmi= 173.8d0
      hmi= 100.d0
      alszi= 0.12d0
*
*-----------------------------------------------------------------------
* Initialization
*
      call wtoinit(ipro,oxc,rs,wmi,zmi)
*-----------------------------------------------------------------------
* Changing flags
*
      oflag5= 'ocoul'
      oval= 'y'
      call wtocflag(rs,oflag5,oval,-1,-1,-1.d0,-1.d0)
      oflag3= 'ipm'
      call wtocopt(oflag3,8)
      oflag3= 'irm'
      call wtocopt(oflag3,8)
      oflag4= 'oqcd'
      oval= 'y'
      call wtocflag(rs,oflag4,oval,-1,-1,-1.d0,-1.d0)
*
*-----------------------------------------------------------------------
* Calling the code and getting some output
*
      call wto(rs,wmi,zmi,zgi,tqmi,hmi,alszi,sig,esig)
      call wtooutput(ipro)
*
      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