program readbin
!     Version: 07/01/97
!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
!     Reads files of the kind "evol.bin" from CEPAM                    c
!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      implicit none
      logical ppt
      integer n,i
      real*8 a(11),b(7)
      character*40 fic

      real*8 rgaz,rjup
      data rgaz/8.31451d7/,rjup/7d9/

!      write(*,*)'File name (ex: evol.bin)'
!      read(*,'(a)')fic
!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
!     Reads the input file "readbin.don"                               c
!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      open(unit=11,file='readbin.don',status='old')
      read(11,'(a)')fic
      close(11)
      write(*,*)'Input file: ',fic
      open(unit=35,file=fic,status='old',form='unformatted')
      read(35)b,ppt
      read(35)n
      write(*,101)b,ppt,n
 101  format('Mass=',1pd10.3,'; Mcore=',d10.3,'; Pice=',d10.3,          &
     &     '; Omega=',d10.3,/,'Alpha=',d10.3,'; X0=',0p,f8.5,'; Teq=',  &
     &     f8.2,'; PPT=',l1,/,'Number of steps (0=unknown):',i5,/)
      write(*,102)
 102  format(t1,'Age',t13,'R/Rjup',t23,'Teff',t33,'Tsurf',t43,'Tint',   &
     &     t53,'Ssurf/R',t63,'Sint/R')
 103  format(f12.5,f10.6,3f10.3,2f10.4)
      do i=0,999
       read(35,err=333)a
       a(2)=a(2)/rjup
       a(5)=a(5)/rgaz
       a(11)=a(11)/rgaz
       write(*,103)a(1),a(2),a(8),a(4),a(10),a(5),a(11)
      enddo
      write(*,*)'Problem: end of file'
      pause
 333  write(*,*)'Number of steps:',i
      end