利用fortran程序，打开并读取200707_Jul.dat中的数据

zznwode

yinlingr 发表于 2012-9-27 20:44

                               
登录/注册后可看大图

简单处理了一下，看看是不是比较容易fortran读取呢？

好处理多了！谢谢奥！

yinlingr

zznwode 发表于 2012-9-27 20:50

                               
登录/注册后可看大图

好处理多了！谢谢奥！

客气啦！也就是提供了一个思路而已。

zznwode

yinlingr 发表于 2012-9-27 20:53

                               
登录/注册后可看大图

客气啦！也就是提供了一个思路而已。

你看看这个行不？
Program lv
Real jp(1:1488)
Integer i,j,m,n,k
Real tt(1:12,1:1488),Sv,Tv,Bc
Real(KIND=8) S,T,SUM
Open(1,File='200707_Jul.dat')
Read (1,*)((tt(m,n),m=1,12),n=1,1488)

    S=0.0
  Do i=1,1488,1
    S=S+tt(1,i)
  EndDo
  Sv=S/1488

SUM=0.0
Do k=1,1488
   jp(k)=tt(1,k)-Sv
   jp(k)=jp(k)*jp(k)/1488
   SUM=SUM+jp(K)
   Bc=sqrt(SUM)

Enddo
   
  Do j=1,1488,1
    T=T+tt(12,j)
  EndDo
  Tv=T/1488
Write(*,*) Bc,Tv
Close(1)
End
有个学长说read处有错，我改了，还是不通

zznwode

做个霸气的木头 发表于 2012-9-27 20:39

                               
登录/注册后可看大图

别人帮你写程序不可能的，通过自己努力，遇到问题了在来吧~

您给看看
Program lv
Real jp(1:1488)
Integer i,j,m,n,k
Real tt(1:12,1:1488),Sv,Tv,Bc
Real(KIND=8) S,T,SUM
Open(1,File='200707_Jul.dat')
Read (1,*)((tt(m,n),m=1,12),n=1,1488)

    S=0.0
  Do i=1,1488,1
    S=S+tt(1,i)
  EndDo
  Sv=S/1488

SUM=0.0
Do k=1,1488
   jp(k)=tt(1,k)-Sv
   jp(k)=jp(k)*jp(k)/1488
   SUM=SUM+jp(K)
   Bc=sqrt(SUM)

Enddo
   
  Do j=1,1488,1
    T=T+tt(12,j)
  EndDo
  Tv=T/1488
Write(*,*) Bc,Tv
Close(1)
End
那里错了。

做个霸气的木头

zznwode 发表于 2012-9-27 22:14

                               
登录/注册后可看大图

您给看看
Program lv
Real jp(1:1488)

你bc是求的标准差吧~，标准差的公式应该是ip全部加在一起后再开方吧，你开方放在里面不就是加一个开方一次吗

zznwode

做个霸气的木头 发表于 2012-9-28 07:48

                               
登录/注册后可看大图

你bc是求的标准差吧~，标准差的公式应该是ip全部加在一起后再开方吧，你开方放在里面不就是加一个开方一次 ...

jp(k)=jp(k)*jp(k)/1488
   SUM=SUM+jp(K)
   Bc=sqrt(SUM)
bc中的sum不是在上一步有限制么？

做个霸气的木头

zznwode 发表于 2012-9-28 17:54

                               
登录/注册后可看大图

jp(k)=jp(k)*jp(k)/1488
   SUM=SUM+jp(K)
   Bc=sqrt(SUM)

不是应该求完和再平方的吗，我是说这个开方应该放在循环后面的~

zznwode

做个霸气的木头 发表于 2012-9-28 18:01

                               
登录/注册后可看大图

不是应该求完和再平方的吗，我是说这个开方应该放在循环后面的~

您给改下！我试试。

做个霸气的木头

zznwode 发表于 2012-9-28 18:03

                               
登录/注册后可看大图

您给改下！我试试。

这还不明白？把那个bc=那一行放在这个循环后面~

zznwode

做个霸气的木头 发表于 2012-9-28 18:08

                               
登录/注册后可看大图

这还不明白？把那个bc=那一行放在这个循环后面~

嗯，改了。多谢！