两组数据分别在两个文件中。格式一摸一样。第一行是半径坐标,第一列是角度坐标,剩下的是对应的数值,共721行,161列
要算的是把两组数据的平方和开根号生成在一个新的文件中,坐标轴不变,写的代码问题很大,求大佬指导·。感觉还是没学明白循环和定义,有没有大佬指条明路啊
[Fortran] 纯文本查看 复制代码 program main
implicit none
real,dimension(721,161)::er
real,dimension(721,161)::ef
real,dimension(721,161)::e
integer i,j
open(unit=10,file="C:\Users\zklz\Desktop\er.csv")
open(unit=20,file="C:\Users\zklz\Desktop\ef.csv")
open(unit=30,file="C:\Users\zklz\Desktop\e.txt")
do i=1,721
read(10,*)(er(i,j),j=1,161)
read(20,*)(ef(i,j),j=1,161)
end do
close(10)
close(20)
do j=1,161
do i=1,721
if(i==1.or.j==1)then
e(i,j)=er(i,j)
else if(i/=1.and.j/=1)then
e(i,j)=sqrt(er(i,j)**2+ef(i,j)**2)
end if
end do
write(30,*)(e(i,j),i=1,721)
end do
close(30)
end
|