c
c   sdcdata.fh
c
c  Data for controlling SD correction optimisation
c
c  bin_ibin(max_ibin,2)   intensity bin limits, for fulls & partials
c        
c  ksdtype               switch for type of mapping of refined parameters
c                         permanent version from input
c  lsdtype               switch for type of mapping of refined parameters
c                         to sdfacc parameters
c                         temporary version for recycling
c        lsdtype = -1    no refinement of parameters
c                =  0    same parameters for all runs, fulls & partials
c                = +1    independent parameters for each run,
c                         and for fulls & partials if present
c                = +2    same SdB, SdAdd for all, different Sdfac
c
c  nsdparam              number of SD parameters
c  nibin(2)              number of intensity bins
c  kfullpart(maxrun)      for each run, 
c                         =  0  fulls only
c                         = +1  partials only
c                         = +2  partials + "significant" number of fulls
c  fixSdB                .true. to fix SdB parameters
c  resid                 residual
c  residrun(irun, mfp)   residual by run & full/part
c  sddelta(irun,ibin,mfp) SD of delta by run, intensity, full|part
c  nsddelta               corresponding number of contributions
c 

c  Maximum number of intensity bins
c (allow more than in the analysis but not too many)
      integer max_ibin
      parameter (max_ibin = 21)

      common /sdcdata/
     $   resid, residrun(maxrun,2), swrun(maxrun,2),
c
     $   sddelta(maxrun,max_ibin,2),
     $   bin_ibin(max_ibin,2),
c
     $   ksdtype, lsdtype, nsdparam, nibin(2),
     $   kfullpart(maxrun), nsddelta(maxrun,max_ibin,2),
     $   fixSdB
      
      double precision resid, residrun, swrun
      real bin_ibin, sddelta
      integer  ksdtype, lsdtype, nsdparam, nibin, kfullpart, nsddelta
      logical fixSdB

      save /sdcdata/