subroutine readin
c
c     author:  c. seftor
c     date:    28-jan-1991
c     purpose: reads in all necessary input for the mateer code
c
      implicit integer*4(i-n), real*8(a-h,o-z)
      character file_name*200
c
      include "input.inc"
      include "inchr.inc"
      include "log.inc"
      include "contrl.inc"
      include "in.inc"
      include "prints.inc"
c
c     new statement
      include "depolt.inc"
c     end of statement
c
      if (iargc() .eq. 0) then
        file_name = 'PROF'
      else
        call getarg(1, file_name)
      end if
      open (unit=22, file=file_name, status='old')
      read (22,'(a)') prfnam
      read (22,*) pres
      read (22,*) nthet
      read (22,*) (theta(num),num=1,nthet)
      read (22,*) nscan
      read (22,*) (scan(num), num=1, nscan)
      read (22,*) naz
      read (22,*) (azmth(num), num=1, naz)
      read (22,*) nalb
      read (22,*) (alb(num), num=1, nalb)

      read (22, *) start_wave, end_wave
      read (22,*) (xprf(num),num=11,1,-1)
      read (22,*) (tmpprf(num),num=11,1,-1)
      read (22,*) (jprint(i),i=1,10)
      read (22,*) num_iter
      if (num_iter .gt. max_num_iter) then
        print*,'Number of iteration ranges requested was ', num_iter
        print*,'Maximum allowed is ', max_num_iter,' Program stopped.'
        stop 'Too many iteration ranges requested.'
      end if
      read (22,*) (wave_iter(i),i=1,num_iter)
      read (22,*) (iter(i),i=1,num_iter)
      read (22,*) ipol
      close (22)
      imuz = nthet
      imu = nscan
      iazmth =naz
      do i=1, nthet
        thnot(i) = theta(i)
      end do

      do i = 1, nscan
        thta(i) = scan(i)
      end do
      if (jprint(9).eq.1) then
         write (*, '(''PRFNAM = '',a)') prfnam
         write (*, *) 'PRES = ', pres
         write (*,*) 'NTHET = ', nthet
         write (*,*) 'THETA = ', (theta(num),num=1,nthet)
         write (*,*) 'NSCAN = ',nscan
         write (*,*) 'SCAN = ', (scan(num),num=1,nscan)
         write (*,*) 'NAZ = ',naz
         write (*,*) 'azmth = ', (azmth(num),num=1, naz)
         write (*,*) 'NALB = ',nalb
         write (*,*) 'ALB = ', (alb(num),num=1,nalb)
         write (*, *) 'START, stop = ', start_wave, end_wave
         write (*,*) 'Ozone profile =',xprf
         write (*,*) 'Temp  profile =',tmpprf
         write (*,*) 'jprint =', (jprint(i),i=1,8)
         write (*,*) 'NUM_ITER = ', num_iter
         write (*,*) (wave_iter(i), i=1,num_iter)
         write (*,*) (iter(i),i=1,num_iter)
         write (*,*) 'IPOL = ', ipol
      endif
      return
      end