Advertisement
Guest User

Untitled

a guest
Mar 10th, 2019
78
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1.       program test
  2.           implicit none
  3.           REAL:: age, y, z, feh
  4.           REAL, dimension(10):: age1
  5.           real, dimension(3):: y1
  6.           real, dimension(9):: z1
  7.           real, dimension(6):: feh1
  8.           integer i
  9.           OPEN(1,file='input.txt')
  10.           read(1,*) age, y, z
  11.           close(1)
  12.           feh = z/0.017
  13.           OPEN(2,file='result.txt',status='replace')
  14.           WRITE(2,*) "age=",age, "y=",y, "z=",z, "Fe/H=",feh
  15.           close(2)
  16.           DATA age1/9.7,9.75,9.8,9.85,9.9,9.95,10,10.05,10.1,10.15/
  17.           DATA y1/0.23,0.26,0.3/
  18.           DATA z1/0.0001,0.0004,0.002,0.002,0.004,0.008,0.017,0.04,0.07/
  19.           DATA feh1/0.00316,0.01,0.0316,0.1,0.316,1/
  20.           open(3,file='border.txt')
  21.           do i = 1, 10
  22.           if(age1(i).lt.age.and.age1(i+1).gt.age) then
  23.           write(3,*) "age1=", age1(i), "age2=", age1(i+1)
  24.           endif
  25.           enddo
  26.           do i=1,3
  27.           if(y1(i).lt.y.and.y1(i+1).gt.y) then
  28.           write(3,*) "y1=", y1(i), "y2=", y1(i+1)
  29.           end if
  30.           end do
  31.           do i=1,9
  32.           if(z1(i).lt.z.and.z1(i+1).gt.z) then
  33.           write(3,*) "z1=", z1(i), "z2=", z1(i+1)
  34.           end if
  35.           end do
  36.           do i=1,9
  37.           if(feh1(i).lt.feh.and.feh1(i+1).gt.feh) then
  38.           write(3,*) "feh1=", feh1(i), "feh2=", feh1(i+1)
  39.           end if
  40.           end do
  41.           close(3)
  42.  
  43.       end program test
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement