1 files changed, 59 insertions, 0 deletions
diff --git a/Boxit.f b/Boxit.f
new file mode 100755
index 0000000..f87ea9a
--- /dev/null
+++ b/Boxit.f
@@ -0,0 +1,59 @@
+ 
+      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 .lt. 2) then
+         do i=1,kount
+            if (xlo .le. xsyn(i)) then
+               lim1obs = i
+               go to 10
+            endif
+         enddo
+         lim1obs = kount
+10       do i=lim1obs,kount
+            if (xhi .lt. 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 .le. xobs(i)) then
+               lim1obs = i
+               go to 20
+            endif
+         enddo
+         lim1obs = lount
+20       do i=lim1obs,lount
+            if (xhi .lt. xobs(i)) then
+               lim2obs = i -1
+               return
+            endif
+         enddo
+         lim2obs = lount
+         return
+      endif
+
+      end
+
+
+
+