FHOutput - write FeynHiggs inputs and outputs to a file

       integer error
       character*(*) filename
       integer key
       double precision sqrts

       subroutine FHOutput(error, filename, key, sqrts)

       FHOutput[filename, key, sqrts]

       FHOutput  runs  the  FeynHiggs  functions encoded in key and writes the
       output to the given file.

       error (output)
              zero if successful, otherwise the line number in  Output.F  from
              which the error message was emitted.

       filename (input)
              the  name  of  the output file, "stdout" and "stderr" are recog-
              nized as special filenames.

       key (input)
              a bit-wise encoding of which output to give:
              bit 0 (1) = input parameters (FHRetrievePara),
              bit 1 (2) = derived parameters (FHGetPara),
              bit 2 (4) = Higgs masses and mixings (FHHiggsCorr),
              bit 3 (8) = uncertainties (FHUncertainties),
              bit 4 (16) = decay widths and branching ratios (FHCouplings),
              bit 5 (32) = the 'fast' flag of FHCouplings,
              bit 6 (64) = electroweak constraints (FHConstraints),
              bit 7 (128) = flavour observables (FHFlavour).
              In case of doubt put key = -1.

       sqrts (input)
              the value of sqrts for running FHHiggsProd.  FHHiggsProd is  run
              only if a positive value of sqrts is given.


                                  14-May-2011                      FHOUTPUT(1)