subroutine lmstssr (kode,is,ie)
c
c a lmstst routine, that does testing to the sample resolution
c of the given instrument data.   compares the start and stop
c times of the data presently in the data common area with
c the start and stop times currently defined in the /comlim/
c common area.   if any of the samples are with the current processing
c limits kode is return=1.   summarising kode:
c kode = -1 not reached required data yet
c kode =  0 at least one of data samples in inclusive start, stop times
c kode =  1 this data past current defined limits.
c
c v1.0 jfe johnson       31 oct 81
c v1.1 jfej  3 nov 81       uses current time in /blktim/
c v1.2 jfej 13 nov 81       split into separate call to lmstst
c
c
c see lmstst routine for following comments
c len is length of data records in millisecs
c nsamp is the number of data samples in that record
c f is computed from f=nsamp/len (samples/milliseconds)
c recipf is ifix(1/f)
c
      real*8    strip_times(4),rec_yd,rec_t,idms
c
      integer*4 iyd,it,iflag,me,ms
c
      common    /blktim/iyd,it
      common    /comlim/strip_times
c
      data f,nsamp/0.064,512/
c
      call lmststr (kode)
      if (kode .ne. 0) return
c
c here there must be more checking for time differences greater
c than one day, otherwise diftms could return silly value
c now added (jfej), iflag =0 idms ok, iflag=1 idms not defined=
c gap greater than 1 day
c
      is=1
      rec_yd=iyd
      rec_t=float(it)*1000.0d0
      call diftmsr (strip_times(1),strip_times(2),rec_yd,rec_t,
     *                    iflag,idms)
      if (iflag .ne. 0) go to 20
      ms=1.0d0+(idms/1000.0d0)*f
      if (ms .gt. 0) is=ms
      if (ms .lt. 1) is=1
   20 continue
      ie=nsamp
      call diftmsr (strip_times(3),strip_times(4),rec_yd,rec_t,iflag,
     *                    idms)
      if (iflag .ne. 0) return
      me=1.0d0+(idms/1000.0d0)*f
      if (me .le. nsamp) ie=me
      if (me. gt. nsamp) ie=nsamp
      return
      end