program lire_evol
!     Version: 19/08/94
!     Lecture de fichiers d'evolution et reecriture sous une autre forme
!     Auteur: T. Guillot

      implicit none
      integer i,j,nb
      real m,r,logp,logt,logro,logl,s
      character*60 nom_fich1,nom_fich2,caractere

      write(*,*)'Entrer les noms des deux fichiers a lire'
      write(*,*)'Attention: ces fichiers doivent etre complets'
      read(*,'(a)')nom_fich1
      read(*,'(a)')nom_fich2

      open(unit=22,file='lire_evol.des',status='unknown')
      open(unit=11,file=nom_fich1,status='old')
      read(11,*)
      do j=1,9
       read(11,'(a)')caractere
       read(11,*)nb
!       goto 902
! 901   write(*,*)'Lecture d''un caractere'
!       goto 900

       if ((j.ge.6).and.(j.le.9)) then
        write(22,*)caractere
        write(22,*)nb
       endif
       do i=1,nb
        read(11,*)m,r,logp,logt,logro,logl,s
        if ((j.ge.6).and.(j.le.9)) then
!        On passe en Mbar
         logp=logp-12
         write(22,'(3(1x,1pd12.5))')m,10**logp,10**logt
        endif
       enddo

      enddo
      close(11)
      
      write(22,*)'C *************** 2EME FICHIER ***********'
      write(*,*)'C *************** 2EME FICHIER ***********'
      open(unit=11,file=nom_fich2,status='old')
      read(11,*)
      do j=1,9
       read(11,'(a)')caractere
       read(11,*)nb
!       goto 902
! 901   write(*,*)'Lecture d''un caractere'
!       goto 900

       if ((j.ge.6).and.(j.le.9)) then
        write(22,*)caractere
        write(22,*)nb
       endif
       do i=1,nb
        read(11,*)m,r,logp,logt,logro,logl,s
        if ((j.ge.6).and.(j.le.9)) then
         logp=logp-12
         write(22,'(3(1x,1pd12.5))')m,10**logp,10**logt
        endif
       enddo

      enddo
      close(11)

      close(22)
      write(*,*)'J''ai cree le fichier "lire_evol.des"'

      end