subroutine fortinit(fname_in, fname_out)

      implicit none

#include "skhead.h"
#include "skparm.h"
#include "sktq.h"
#include "skbadc.h"
#include "geopmt.h"

      INTEGER fortread

      integer iret, lenchr

      character*(*) fname_in, fname_out

c*** paw and zebra variables
      INTEGER NH
      PARAMETER (NH=1000000)
      REAL     paw(NH),wallsk
      COMMON /PAWC/PAW

      integer idim
      parameter (idim=4000000)
      real zbs(idim)
      common/kzbs/zbs

c*** init.
      call kzinit
      call hlimit(-nh)

c      call skbadopt(23)
c      call SKOPTN("31,30,25")
      CALL SKOPTN('31,30,27,26,25,23,16')

      if (lenchr(fname_out).ge.1) then
c*** open output file
         call set_rflist(20,fname_out,'LOCAL',' ','WRT',' ',' ',
     &        'recl=5670 status=new',' ',' ')
         print *, 'open as LOCAL: ', fname_out(1:lenchr(fname_out))
         call skopenf(20,1,'Z',iret)
         print *,'out iret = ',iret
         if (iret.ne.0) goto 9999
      endif

c*** open input rfm file
      if (fname_in(1:5) .eq. '/disk') then
         call set_rflist(10,fname_in,'DISK',' ','RED',' ',' ',
     &        'recl=5670 status=old',' ',' ')
         print *, 'open as DISK: ', fname_in(1:lenchr(fname_in))
      else
         call set_rflist(10,fname_in,'LOCAL',' ','RED',' ',' ',
     &        'recl=5670 status=old',' ',' ')
         print *, 'open as LOCAL: ', fname_in(1:lenchr(fname_in))
      endif

      call skopenf(10,1,'Z',iret) ! open input file
      print *,'in iret = ',iret
      if (iret.ne.0) goto 9999

cc*** read an event
      iret=fortread()
      if (iret.ne.0) goto 9999

      return

 9999 continue
      stop
      end