ADAS Subroutine hapixv

       subroutine hapixv( ndwvl  , ndpix  , fcrit  ,
     &                    nwvl   , npix   , wvmin  , wvmax  ,
     &                    iwvrg  , cpixmx ,
     &                    wvl    , tev    , amssno , pec    ,
     &                    cpixa  , ind1   , ind2
     &                  )

c-----------------------------------------------------------------------
c
c  ****************** fortran77 subroutine: hapixv *********************
c
c  purpose:  distribute Doppler broadened line emission into pixel range
c
c  calling program: hapecf
c
c
c  subroutine:
c
c  input : (i*4)  ndwvl   = maximum number of wavelength intervals
c  input : (i*4)  ndpix   = maximum number of pixels per wvln. interval
c  input : (r*8)  fcrit   = pixel counts for the selected line below
c                           this fraction of the largest pixel count are
c                           discounted.
c
c  input : (i*4)  nwvl    = number of wavelength intervals
c  input : (i*4)  npix()  = number of pixels assigned to wavelength interval
c  input : (r*8)  wvmin() = lower limit of wavelength interval (ang)
c  input : (r*8)  wvmax() = upper limit of wavelength interval (ang)
c
c  input : (i*4)  iwvrg   = index of wavelength range in which line lies
c  input : (r*8)  cpixmx  = largest pixel count currently found
c                           for the wavelength range
c
c  input : (r*8)  wvl     = input line wavelength for test(ang)
c  input : (r*8)  tev     = electron temperature (eV)
c  input : (r*8)  amssno  = atomic mass number
c  input : (r*8)  pec     = photon emissivity coefficient for line
c
c  output: (r*8)  cpixa() = counts in each pixel for the line
c  output: (r*8)  ind1    = first pixel with non-negligible count
c  output: (r*8)  ind2    = last pixel with non-negligible count
c
c routines:
c          routine    source    brief description
c          -------------------------------------------------------------
c          i4unit     adas      fetch unit number for output of messages
c          r8erfc     adas      returns erfc(x) function value
c
c author:  Hugh Summers, University of Strathclyde
c          JA7.08
c          tel. 0141-548-4196
c
c date:    15/01/02
c
c version  : 1.1
c date     : 15-01-2002
c modified : H P Summers
c             - first version.
c
c version  : 1.2
c date     : 18-06-2007
c modified : H P Summers
c             - corrected error in ind1 & ind2 return.
c
c-----------------------------------------------------------------------
c-----------------------------------------------------------------------
      INTEGER             IND1,        IND2,        IWVRG,       NDPIX
      INTEGER             NDWVL,       NPIX(NDWVL), NWVL
      REAL*8              AMSSNO,      CPIXA(NDPIX),             CPIXMX
      REAL*8              FCRIT,       PEC,         TEV,         WVL
      REAL*8              WVMAX(NDWVL),             WVMIN(NDWVL)