Fortran Coder

查看: 1084|回复: 4
打印 上一主题 下一主题

寻求帮助

[复制链接]

12

帖子

4

主题

0

精华

入门

F 币
52 元
贡献
26 点
跳转到指定楼层
楼主
发表于 2023-9-12 11:19:47 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
本帖最后由 IDCY 于 2023-9-12 15:47 编辑

假设两个txt文件
txt1
t=    x=   y=
0     1      2
0     2      3
0      3      4
txt2
t=    x=   y=
0     1      9
0     2      4
0      3      1

计算txt2中y的值减去txt1中y的值的绝对误差,即abs(y1-y2)
并将结果保存到一个新的txt文件

代码在下面
数据在附上


结果运行出不来


是这两个dat文件中的数据存在问题嘛
如果要修改的话怎么修改


下载.png (17.04 KB, 下载次数: 108)

下载.png

Riemann.dat

1.68 KB, 下载次数: 1

result.dat

1.27 KB, 下载次数: 0

分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

12

帖子

4

主题

0

精华

入门

F 币
52 元
贡献
26 点
沙发
 楼主| 发表于 2023-9-12 15:43:05 | 只看该作者
program main
implicit none
integer,parameter:: nx=21,ny=4
real::er(nx,ny),ef(nx,ny),e(nx,ny)
integer i,j
call read_data("C:\Users\20455\Desktop\error\result.dat",er)
call read_data("C:\Users\20455\Desktop\error\Riemann.dat",ef)
do i=1,nx
do j=2,ny
e(i,j)= abs(er(i,j)-ef(i,j))
end do
end do
e(:,1)=er(:,1)

call write_data("C:\Users\20455\Desktop\error\error.dat",e)
end program main

subroutine read_data(filename,data)
implicit none
integer i,j
integer,parameter:: nx=21,ny=4
character*(*),intent(in)::filename
real,intent(out)::data(nx,ny)
open(10,file=filename)
do i=1,nx
do j=2,ny
read(10,*)(data(i,j))
end do
end do
close(10)
end subroutine read_data

subroutine write_data(filename,data)
implicit none
integer i,j
integer,parameter:: nx=21,ny=4
character*(*),intent(in)::filename
real,intent(in)::data(nx,ny)
open(10,file=filename)
do i=1,nx
do j=1,ny
write(10,*)(data(i,j))
end do
end do
close(10)
end subroutine write_data

12

帖子

4

主题

0

精华

入门

F 币
52 元
贡献
26 点
板凳
 楼主| 发表于 2023-9-12 16:05:22 | 只看该作者
好了 ,我已经改好了
是其中的j=2,ny的问题 应该是j=1,ny

12

帖子

4

主题

0

精华

入门

F 币
52 元
贡献
26 点
地板
 楼主| 发表于 2023-9-12 16:06:06 | 只看该作者
IDCY 发表于 2023-9-12 15:43
program main
implicit none
integer,parameter:: nx=21,ny=4

subroutine read_data(filename,data)
implicit none
integer i,j
integer,parameter:: nx=21,ny=4
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)
end do

close(10)
end subroutine read_data

subroutine write_data(filename,data)
implicit none
integer i,j
integer,parameter:: nx=21,ny=4
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)
end do

close(10)
end subroutine write_data

12

帖子

4

主题

0

精华

入门

F 币
52 元
贡献
26 点
5#
 楼主| 发表于 2023-9-12 16:07:50 | 只看该作者
已解决,修改后的代码附上
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-1 08:03

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表