c      integer b,e
      character*2 e
      character*80 line
      character*54 a
      character*10 l
      real g,h,i,j,k,occup,bfac
c      character*4 c,d,f,g,h,i,j,k
c
c look for command line arguments
        n_arg = iargc()
        if (n_arg .EQ. 0)  STOP 'no value given for occupancy!'
c           recover month
            call getarg(1, line)
            read(line, *, end=666, err=666) temp
            occup = temp

c   5  read(5,7,end=200)a,b,c,d,e,f,g,h,i,j,k
c   7  format(A4,I7,A5,A4,A2,I4,4X,3F8.3,2F6.2)
   5  read(5,7,end=200)a,j,k,l
   7  format(A54,2F6.2,A10)
       if (a(1:4) .eq. "END ")    goto 300
       write(6,10)a,occup,k,l
  10  format(A54,2F6.2,A10)


      goto 5
200   continue
300   continue
c 300  write(6,310)'END'
c      if(k.eq.0.0) write(stderr.*)'no atoms read or last B=0!'
310   format(A3)
      stop
666    STOP 'error reading command line argument'
      end