Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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
- 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)
- end program test
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement