Текст программы к 1 лабораторной (допуск к зачету)

parameter (N=49,M=98,MO3=77,MO2=49)
!real F(N,M), SO2(M), SO3(M), O2(N), O3(N)
real, dimension (N,M)::F
real, dimension (M)::SO2, SO3
real AL(M),SO2R(MO2),SO3R(MO3),ALO2(MO2),ALO3(MO3)
real, dimension (N)::O2, O3, z
!F-солнечная рад.
!SO2-сечение поглощения O2
!O2-концентрация
character AAA*10
open(unit=11,file='D:\Fortran\lis_shm_year.mon')
!read(11,*)F
do I=1,17
read (11,*) AAA
print *,AAA
end do
do J=1,M
read(11,*)AL1,AL2,F(1,J)
AL(J)=(AL1+AL2)/2
print *,AL(J),AL1,AL2,F(1,J)
end do
close(11)
open(unit=12,file='D:\Fortran\x_sec_00.dat')
do J=1,20
read (12,*) AAA
print *,AAA
end do
do J=1,17
read(12,*)nom,ALO2(J),SO2R(J)
print *,nom,ALO2(J),SO2R(J)
end do
J=1
read(12,*) AAA
print *,AAA
do J=18,49
read(12,*)nom,ALO2(J),SO2R(J)
print *, nom,ALO2(J),SO2R(J)
end do
J=1
read(12,*) AAA
print *,AAA
do J=1,77
read(12,*)AL1,AL2,SO3R(J)
ALO3(J)=(AL1+AL2)/2
print *,ALO3(J),AL1,AL2,SO3R(J)
end do
close(12)
open(unit=13,file='D:\Fortran\trop_centers.txt')
do J=1,2
read (13,*) AAA
print *,AAA
end do
do J=1,49
read(13,100)z(J),O3(J),O2(J)
100 format(e11.4,38x,e10.4,e12.4)
print 50,z(J),O3(J),O2(J)
50 format(e11.4,e12.4,e12.4)
end do
do J=1,M
SO2(J)=0
do I=1,MO2
if(AL(J)<=ALO2(I))goto 1
end do
! или.le.<= .lt.< .ge.>= .gt.> .eq.==
goto 2
1 JR=I
JL=I-1
SO2(J)=SO2R(JL)+(SO2R(JR)-SO2R(JL))/(ALO2(JR)-ALO2(JL))*(AL(J)-ALO2(JL))
2 continue
print *,J,AL(J),SO2(J),JR,JL
end do
end