* WTO.20 version of July 24 1997
*
* Needs the following files
*
* WTO_20_IO_T
* WTO_20_MAIN_T
* WTO_20_CCA_T
* WTO_20_NCA_T
* WTO_20_HCA_T
* WTO_20_LIB_T
*
* 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)
*
*------------------------------------------------------------------------
*
      oflag3= 'ofl'
      oval= 'c'
      call wtocflag(rs,oflag3,oval,-1,-1,-1.d0,-1.d0)
      oflag4= 'ospe'
      oval= 'y'
      call wtocflag(rs,oflag4,oval,-1,-1,-1.d0,-1.d0)
*
*------------------------------------------------------------------------
*
      oflag5= 'osinw'
      oval= 'f'
      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)
      separa= 0.1d0
      asccsw= -1.d0
      eusw= 0.d0
      edsw= 0.d0
      alsw= 1.d0
      sthqsw= -1.d0
      call wtoswcutset(separa,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