subroutine solread(wave, flux, crs, nsol)
        implicit double precision (a-h,o-z)
        include "parameters.dat"

        dimension wave(nd5),flux(nd5),crs(nd5,0:3)

 1      continue
        write(unit=*,fmt='(26h Enter 0=solmax,1=solmin >,$)')
        read (unit=*,fmt=*) icycle
        if(icycle .eq. 0) then
           open(unit=60,file='soluvmax.dat',status='old')
        else if(icycle .eq. 1) then
           open(unit=60,file='soluvmin.dat',status='old')
        else
           goto 1
        endif

        read(unit=60,fmt=*) nsol
        do n = 1, nsol
           read(unit=60,fmt=*) wave(n),flux(n),(crs(n,j),j=0,3)
        end do

 10     close(unit=60)

        return
        end