program test
          implicit none
          REAL:: age, y, z, feh
          REAL, dimension(10):: age1
          real, dimension(3):: y1
          real, dimension(9):: z1
          real, dimension(6):: feh1
          real, allocatable, dimension(2:,:)::array
          integer i
          OPEN(1,file='input.txt')
          read(1,*) age, y, z
          close(1)
          feh = z/0.017
          OPEN(2,file='result.txt',status='replace')
          WRITE(2,*) "age=",age, "y=",y, "z=",z, "Fe/H=",feh
          close(2)
          DATA age1/9.7,9.75,9.8,9.85,9.9,9.95,10,10.05,10.1,10.15/
          DATA y1/0.23,0.26,0.3/
          DATA z1/0.0001,0.0004,0.002,0.002,0.004,0.008,0.017,0.04,0.07/
          DATA feh1/0.00316,0.01,0.0316,0.1,0.316,1/
          open(3,file='border.txt')
          do i = 1, 10
          if(age1(i).lt.age.and.age1(i+1).gt.age) then
          write(3,*) "age1=", age1(i), "age2=", age1(i+1)
          endif
          enddo
          do i=1,3
          if(y1(i).lt.y.and.y1(i+1).gt.y) then
          write(3,*) "y1=", y1(i), "y2=", y1(i+1)
          end if
          end do
          do i=1,9
          if(z1(i).lt.z.and.z1(i+1).gt.z) then
          write(3,*) "z1=", z1(i), "z2=", z1(i+1)
          end if
          end do
          do i=1,9
          if(feh1(i).lt.feh.and.feh1(i+1).gt.feh) then
          write(3,*) "feh1=", feh1(i), "feh2=", feh1(i+1)
          end if
          end do
          close(3)
          open(4, File='FLUXTREC.data')

      end program test