.v1和.v2 求第三边长

zlzshubao

有两个数据文件，分别是以.v1和.v2作为文件扩展名，里面是很复杂的数据，如下图所示：


现在想要求出两份类似文件的每个对应数据的平方和（相当于每两个对应数据是直角三角形两条直角边，求第三条边长），之后将结果存储在一个新的文件中。

Fortran小白求指导代码应该怎样编写~

谢谢~~

fcode

你的需求描述得不清楚，没法解答。

文件里各列是什么含义？哪些是你所需的？读出来之后，如何判断“对应”？

zlzshubao

fcode 发表于 2015-6-2 09:35
你的需求描述得不清楚，没法解答。

文件里各列是什么含义？哪些是你所需的？读出来之后，如何判断“对应” ...

谢谢您的回复。抱歉问题问的不清楚，我补充一下：
我有两个文件，打开均如上图所示，一堆数据，每一个数据都是一个时刻的某个方向的流速（两个文件分别是两个方向的流速），现在想要求两个文件中对应时刻的流速的合速度，也就是两个速度的平方和然后开根号，得到合速度。“对应”指的就是两个文件中的第一个数据对应第一个数据，第二个数据对应第二个数据……
之后将求得的合速度的数据文件单独存起来，成为一个新的dat或txt文件。
不知道有没有描述清楚，麻烦你啦~

fcode

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

Program www_fcode_cn
  Implicit None
  Integer , parameter :: N = 20 !// 每行20个数
  Real :: a( N ) , b( N )
  Integer :: err
  Open( 12 , File = "a.v1" )
  Open( 13 , File = "b.v2" )
  Open( 14 , File = "c.v3" )
  Do
    Read( 12 , * , iostat = err ) a
    if ( err /= 0 ) exit
    Read( 13 , * , iostat = err ) b
    if ( err /= 0 ) exit
    Write( 14 , * ) sqrt( a*a + b*b )
  End Do
  Close( 12 )
  Close( 13 )
End Program www_fcode_cn

zlzshubao

fcode 发表于 2015-6-2 22:50
[mw_shl_code=fortran,true]Program www_fcode_cn
  Implicit None
  Integer , parameter :: N = 20 !//  ...

不知道为什么，还是有这样的错误

fcode

windows 的绝对路径大概是这个格式：
D:\\Adp_M2_1215to1218.v1

注意是冒号，不是分号。两个斜线 \\
最后，看错误提示要在黑色窗口看。你帖的最后一张图看似是错误，但其实不是（无视它）