两个数组的计算问题

necrohan

[Fortran] 纯文本查看 复制代码

program main
implicit none
integer,parameter::nx=721,ny=604
real::er(nx,ny),ef(nx,ny),e(nx,ny)
integer i,j

! 读入数据
call read_data("C:\Users\zklz\Desktop\er.csv",er)
call read_data("C:\Users\zklz\Desktop\ef.csv",ef)

! 计算
do i=2,721
  do j=2,161
    e(i,j)=sqrt(er(i,j)**2+ef(i,j)**2)
  enddo
enddo
e(1,:)=er(1,:); e(:,1)=er(:,1) ! 保存第一行、第一列

! 输出
call write_data("C:\Users\zklz\Desktop\e.txt",e)
end program main
!-----------------------------
subroutine read_data(filename,data)
implicit none
integer,parameter::nx=721,ny=604
character*(*),intent(in)::filename
real,intent(out)::data(nx,ny)

open(10,file=filename)
do i=1,nx
  read(10,*)(data(i,j),j=1,ny)
enddo
close(10)
end subroutine read_data
!-----------------------------
subroutine write_data(filename,data)
implicit none
integer,parameter::nx=721,ny=604
character*(*),intent(in)::filename
real,intent(in)::data(nx,ny)

open(10,file=filename)
do i=1,nx
  write(10,*)(data(i,j),j=1,ny)
enddo
close(10)
end subroutine write_data
!-----------------------------

necrohan

循环不明白，你连教材都没看吧，静下心来看半小时就足够足够的了

necrohan

wawewen 发表于 2019-12-2 16:37
感觉连基本的认识都没有啊，看教材感觉有时候看不太明白

如果不能看进去，以后问题只会更多

necrohan

wawewen 发表于 2019-12-3 15:29
又除了彭国伦的书，大佬有没有推荐的教材

随便找一本
http://pan.fcode.cn/