SANCinterfaceF.F

      subroutine 
     $  flagset(iqedx,iewx,ibornx,gfschemex,ifggx,ncx,fcx,tlmu2x)
      implicit none!
#include "s2n_declare.h"
      integer iqedx,iewx,ibornx,gfschemex,ifggx
      real*8 ncx,fcx,tlmu2x
*-- Flag settings 
      iqed     = iqedx
      iqcd     = 0      ! may  be something else possible
      iew      = iewx
      iborn    = ibornx
      gfscheme = gfschemex
      ifgg     = ifggx
* Input constants (masses, couplings)
#include "leshw_input.h"
      call s2n_init()
      call ffini()
      nc=ncx
      fc=fcx
      tlmu2=tlmu2x
      end

      subroutine paraget(mtax,conhcx,pix)
      implicit none!
      real*8 mtax,conhcx,pix
#include "s2n_declare.h"
      mtax=mta
      conhcx=conhc
      pix= pi
      end

      subroutine printconsts (pmg)
      implicit none!
      integer pmg,nout
#include "s2n_declare.h"

      if ((pmg.eq.0).or.(pmg.eq.1)) then
         print *
         print *, "flag iqed     = ",iqed
         print *, "flag iew      = ",iew
         print *, "flag qcd      = ",iqcd
         print *, "flag gfscheme = ",gfscheme
         print *, "thmu2(GeV^2)  = ",thmu2
         print *, "tlmu2(GeV^2)  = ",tlmu2
         print *, "mup(GeV)      = ",mup
         print *, "mta(GeV)      = ",mta
       endif
      if ((pmg.eq.0).or.(pmg.eq.2)) then
         print *
         print *, "alphai    = ",alphai
         print *, "alpha     = ",alpha
         print *, "GF(GeV-^2)= ",gf
         print *, "g         = ",g
         print *, "e         = ",e
         print *, "alphas    = ",alpha
         print *, "gs        = ",gs
         print *, "cwm       = ",ctw
         print *, "swm2      = ",stw2
         print *, "mh(GeV)   = ",mh
         print *, "mw(GeV)   = ",mw
         print *, "mz(GeV)   = ",mz
         print *, "wz(GeV)   = ",wz
         print *, "ww(GeV)   = ",ww
         print *, "wh(GeV)   = ",wh
         print *, "wtp(GeV)  = ",wtp
      endif
      
      nout=11
      OPEN(nout   ,file='var.dump')
         write(nout,*) "TAUOLA interace: SANC input parameters:"
      if ((pmg.eq.0).or.(pmg.eq.1)) then
         write(nout,*) " "
         write(nout,*) "flag iqed     = ",iqed
         write(nout,*) "flag iew      = ",iew
         write(nout,*) "flag qcd      = ",iqcd
         write(nout,*) "flag gfscheme = ",gfscheme
         write(nout,*) "thmu2(GeV^2)  = ",thmu2
         write(nout,*) "tlmu2(GeV^2)  = ",tlmu2
         write(nout,*) "mup(GeV)      = ",mup
         write(nout,*) "mta(GeV)      = ",mta
       endif
      if ((pmg.eq.0).or.(pmg.eq.2)) then
         write(nout,*) " "
         write(nout,*) "alphai    = ",alphai
         write(nout,*) "alpha     = ",alpha
         write(nout,*) "GF(GeV-^2)= ",gf
         write(nout,*) "g         = ",g
         write(nout,*) "e         = ",e
         write(nout,*) "alphas    = ",alpha
         write(nout,*) "gs        = ",gs
         write(nout,*) "cwm       = ",ctw
         write(nout,*) "swm2      = ",stw2
         write(nout,*) "mh(GeV)   = ",mh
         write(nout,*) "mw(GeV)   = ",mw
         write(nout,*) "mz(GeV)   = ",mz
         write(nout,*) "wz(GeV)   = ",wz
         write(nout,*) "ww(GeV)   = ",ww
         write(nout,*) "wh(GeV)   = ",wh
         write(nout,*) "wtp(GeV)  = ",wtp
      endif
 
      CLOSE(nout)

      return
      end