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, dimension(2):: age2(2), y2(2), z2(2), feh2(2)
          real, dimension(416):: age3(416), y3(416), z3(416), feh3(416)
      real, dimension(1000000)::flx2(2,2,2,2,40000),flx(40000,416),
     * wavn(40000,1)

          integer i, i1, j1, j, k, l, nl
          DATA age1/9.7,9.75,9.8,9.85,9.9,9.95,10.0,10.05,10.1,10.15/
          DATA y1/0.23,0.26,0.30/
          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.0/
          OPEN(1,file='input.txt')
          read(1,*) age, y, z
          close(1)
          feh = z/0.017
          OPEN(2,file='feh.txt',status='replace')
          WRITE(2,6) "age=",age, "y=",y, "z=",z, "Fe/H=",feh
          close(2)
          open(3,file='border.txt')
          do i = 1, 10
          if(age1(i).lt.age.and.age1(i+1).gt.age) then
          write(3,6) "age1=", age1(i), "age2=", age1(i+1)
          age2(1)=age1(i)
          age2(2)=age1(i+1)
          endif
          enddo
          do i=1,3
          if(y1(i).lt.y.and.y1(i+1).gt.y) then
          write(3,6) "y1=", y1(i), "y2=", y1(i+1)
          y2(1)=y1(i)
          y2(2)=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,6) "z1=", z1(i), "z2=", z1(i+1)
          z2(1)=z1(i)
          z2(2)=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,6) "Fe/H1=", feh1(i),"Fe/H2=", feh1(i+1)
          feh2(1)=feh1(i)
          feh2(2)=feh1(i+1)
          end if
          end do
          close(3)

          open(4,file='FLUXTREC.dat')
          open(5,file='wtf.txt')
          open(6,file='result.txt')
          nl=40000
          read(4,5) (wavn(i,1),i=1,nl)
          write(5,*) (wavn(i,1),i=1,nl)
          do j1=1,415
          read(4,7,end=100) age3(j1),y3(j1),z3(j1),feh3(j1),(flx(i1,j1),
     *i1=1,nl)
          do i=1,2
          do j=1,2
          do k=1,2
          do l=1,2
          if(age3(j1).eq.age2(i).and.y3(j1).eq.y2(j).and.z3(j1).eq.z2(k)
     *.and.feh3(j1).eq.feh2(l)) then
          write(6,8) age3(j1),y3(j1),z3(j1),feh3(j1), (flx(i1,j1),
     *i1=1,nl)
  !        flx2(i,j,k,l,i1)=flx(i1)
          endif
          end do
          end do
          end do
          end do
          end do

 100     continue
          close(4)
          close(5)
          close(6)
7         format(2x,F8.5,3x,F7.5,3x,F7.5,3x,F7.5,40000E11.4)
5         format(40x,414F11.4)
6         format(2X,A6,F9.6,4X,A6,F9.6,A6,F9.6)
8         format(f6.3,3x,f5.3,3x,f7.5,3x,f8.6,3x,40000E11.4)
       end program test