subroutine boxit c****************************************************************************** c This subroutine figures out the point numbers of the wavelength c boundaries for synthetic spectrum plots c****************************************************************************** implicit real*8 (a-h,o-z) include 'Atmos.com' include 'Pstuff.com' include 'Linex.com' include 'Equivs.com' c*****if only a synthetic spectrum is being plotted, use beginning and c end points of the synthetic spectrum instead of those of the observed c spectrum if (plotopt < 2) then do i=1,kount if (xlo <= xsyn(i)) then lim1obs = i go to 10 endif enddo lim1obs = kount 10 do i=lim1obs,kount if (xhi < xsyn(i)) then lim2obs = i -1 return endif enddo lim2obs = kount return c*****and here is the same logic when an observed spectrum exists else do i=1,lount if (xlo <= xobs(i)) then lim1obs = i go to 20 endif enddo lim1obs = lount 20 do i=lim1obs,lount if (xhi < xobs(i)) then lim2obs = i -1 return endif enddo lim2obs = lount return endif end