Fortran Coder

查看: 9133|回复: 2
打印 上一主题 下一主题

[求助] fortran程序能运行,生成结果是空的?

[复制链接]

3

帖子

2

主题

0

精华

新人

F 币
15 元
贡献
10 点
跳转到指定楼层
楼主
发表于 2020-4-1 14:21:19 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
本帖最后由 是来来啊 于 2020-4-1 15:20 编辑

附上程序
[Fortran] 纯文本查看 复制代码
program main
  integer ist,dd 
  parameter(ist=50,dd=12) 
        dimension cor(ist)
        character*8 istid(ist)
  real istr(ist),istr19(ist,dd)
!    ---------------读实况值 istr(i)
    open(20,file='F:\shuju\sk_06R.txt')
       do  i=1,ist
       read(20,*)            
       enddo
    close(20)
!    ---------------读预报值 istr19(i,j)
    open(20,file='F:\shuju\ms_06R_26.txt')
         do  i=1,50
         read(20,*)(istr19(i,j),j=1,dd)
         enddo
    close(20)
open(2,file='F:\shuju\result\ACC_06R.txt')  !!!!输出ACC文件
     do j=1,ist
        write(2,'(I5,2x,<6>(f10.3,2X))')cor(j)
     enddo
   close(2)
   end

    Subroutine cor(x,y,n,coefficient)
        dimension x(n),y(n)
        real coefficient,avex,avey,sum1,sum2
        sum1=0.0;sum2=0.0
        do i=1,n
        sum1=x(i)+sum1
        sum2=y(i)+sum2
        enddo
        avex=sum1/n;avey=sum2/n
        sum1=0.0;sum2=0.0;sum3=0.0
        do j=1,n
        sum1=sum1+(x(j)-avex)*(y(j)-avey)
        sum2=sum2+(x(j)-avex)**2
        sum3=sum3+(y(j)-avey)**2
        enddo
        coefficient=sum1/sqrt(sum2*sum3)
        return
        end

sk_06R.txt

3.06 KB, 下载次数: 1

ms_06R_26.txt

3.47 KB, 下载次数: 1

acc_20200401151749.pdf

95.02 KB, 下载次数: 0

要求两组数据的距平先关系数

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

2033

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1641 元
贡献
709 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

沙发
发表于 2020-4-1 14:35:50 | 只看该作者
这代码逻辑太差了,完全看不懂你要做什么。

主程序里打开了 sk_06R.txt,然后啥也没读,溜了一圈就走了。

然后打开 ms_06R_26.txt 文件,读入到 istr19,然后就啥也没了。

然后打开 ACC_06R.txt,输出了 cor数组,这个cor数组没有值。

cor函数需要 x y n,计算 coefficient
但是和上面的 istr19,cor数组没有任何关联。

3

帖子

2

主题

0

精华

新人

F 币
15 元
贡献
10 点
板凳
 楼主| 发表于 2020-4-1 14:54:43 | 只看该作者
fcode 发表于 2020-4-1 14:35
这代码逻辑太差了,完全看不懂你要做什么。

主程序里打开了 sk_06R.txt,然后啥也没读,溜了一圈就走了。

捂脸,是要求两个文本数据的距平相关系数,刚开始接触fortran,不是很懂
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2024-12-23 22:18

Powered by Tencent X3.4

© 2013-2024 Tencent

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