subroutine lmssetr (kode)
c
c read start and stop limits and subperiod definitions (ms)
c return kode is as follows kode=-1 end of file in input
c                           kode=0  successful input
c                           kode=1  syntactically incorrect input
c
c v1.1      jfej      15 nov 81      now uses thms for input
c v1.2      jfej       6 dec 81      debug removed, stop time=0 allowed
c v1.3      jfej       7 dec 81      deluns common added
c v1.4      jfej      23 jan 82      correction
c v1.5  jfej    26 feb 82      allows negative ilims(1,2) to specify
c                        start/stop record nos.
c v1.6  jfej    23 jul 82       subperiod definition now subset/fraset
c v2.0  rlw     03 apr 86       changed from i*4 to real*8
c
      integer*4 kode,ih,im,is
c
      integer*4 hhmmss
c
      real*8      run_times(4),strip_times(4),ilims(8),inc
c
      common    /comtms/run_times,inc,dur
     *                /comlim/strip_times
     *                /deluns/lundat,lunin,lunout
c
      kode=0
   5  write (*,500)
  500 format (' start yr,dn,hhmmss,msec stop yr,dn,hhmmss,msec')
      read (*,*) ilims
c
c            ! see if end year and day same as start year and day
c
      if (ilims(5) .eq. 0) ilims(5)=ilims(1)
      if (ilims(6) .eq. 0) ilims(6)=ilims(2)
c
c            ! convert start yy,ddd to yddd
c
c      ilims(1)=mod(ilims(1),10.0d0)
      ilims(1)=ilims(1) - 80.0d0
      ilims(1)=ilims(1)*1000.0d0+ilims(2)
      run_times(1)=80000.0d0 + ilims(1)
c
c            ! convert start hhmmss to sec * 1000000.
c
      hhmmss=ilims(3)
      call thms(hhmmss,ih,im,is)
      ilims(2)=ih
      ilims(2)=1000000.0d0*(ilims(2)*3600.0d0+im*60.0d0+is)
      run_times(2)=ilims(2)+ilims(4)*1000.0d0
c
c            ! convert end yy,ddd to yddd
c
c      ilims(3)=mod(ilims(5),10.0d0)
      ilims(3)=ilims(5) - 80.0d0
      ilims(3)=ilims(3)*1000.0d0+ilims(6)
      run_times(3)=80000.0d0 + ilims(3)
c
c            ! convert end hhmmss to sec * 1000000.
c
      hhmmss=ilims(7)
      call thms(hhmmss,ih,im,is)
      ilims(4)=ih
      ilims(4)=1000000.0d0*(ilims(4)*3600.0d0+im*60.0d0+is)
      run_times(4)=ilims(4)+ilims(8)*1000.0d0
c
c update current time limits
c
      do i=1,4
         strip_times(i)=run_times(i)
      end do
      inc=0.0d0
      dur=0.0d0
      return
c
   99 kode=-1
      return
      end