Fortran Coder

标题: 如何按行进行数据操作 [打印本页]
作者: wxy    时间: 2017-1-29 21:32
标题: 如何按行进行数据操作
我在做一个计算积温的小程序，根据彭国伦计算学生成绩的例子做的，已经实现一整年积温的计算 ，现在我想根据第一列的天数 1.实现输入起止日期即起止行数，灵活的计算两个指定行之间的数据而不是简单的一整列的计算，不知道该怎末给指定行赋值 2.就是temp我现在是手动赋值 ，能不能根据文件行数自动赋值给它，我试着把获取文件行数的代码放进去，一放进去就提示mod不能用，大过年的打扰大家了
[Fortran] 纯文本查看 复制代码
module typedef

    type tdd

        integer :: stat=0

        real ::  tmax , tmin ,  total , average

    end type

end module



PROGRAM MAIN

    USE typedef

    implicit none

    real::dy,rain

    INTEGER,PARAMETER::temp=365

    type(tdd)::t(temp)! 保存气温值

    type(tdd):: total !日均温累加

    integer i



    open(7,file="GDNH94")

    read(7,*)



    do i=1,temp

        read(7,*)dy , t(i)%tmax , t(i)%tmin , rain



        t(i)%average=(t(i)%tmax+t(i)%tmin)/2    !日均温

        total%average=total%average+t(i)%average    !日均温累加=积温

    enddo

    print *,total%average    !

    STOP

end PROGRAM

QQ截图20170129210133.png (5.26 KB, 下载次数: 431)

QQ截图20170129210133.png
作者: vvt    时间: 2017-1-29 23:16
本帖最后由 vvt 于 2017-1-29 23:26 编辑

问题1. 在你的例子里，已经把所有文件中的365行，都读入到内存里（t数组）了，所以你可以自由的操作内存（t数组），而不再需要管文件的第几行了。
   此外，行是文件的内容，如果你不需要修改文件内容。就不存在“怎么给行赋值”的问题。（行也不能，没法赋值）
   注意赋值是针对内存中的变量。文件内容只能（从文件）读取（到内存）或（从内存）写入（到文件）

问题2. 请给出代码和错误提示。

作者: wxy    时间: 2017-1-31 10:20
本帖最后由 wxy 于 2017-1-31 11:32 编辑
vvt 发表于 2017-1-29 23:16
问题1. 在你的例子里，已经把所有文件中的365行，都读入到内存里（t数组）了，所以你可以自由的操作内存（t ...

谢谢 ，我还有一个问题，帮忙看一下第23行怎么改能够把start，end输入的值赋给i,没有错误提示。但是运行时还没赋值就输出了错误的结果
[Fortran] 纯文本查看 复制代码
module typedef

    type tdd

        integer :: stat=0

        real ::  tmax , tmin ,  total , average

    end type

end module



PROGRAM MAIN

    USE typedef

    implicit none

    real::rain

    INTEGER,PARAMETER::temp=365

    type(tdd)::t(temp)! 保存气温值

    type(tdd):: total !日均温累加

    integer i,dy,start,end



    open(7,file="1988.txt")



    write(*,*)"start:"

    read(7,*) start

    write(*,*)"end:"

    read(7,*) end



    do i=start,end



        read(7,*)dy , t(i)%tmax , t(i)%tmin , rain



        t(i)%average=(t(i)%tmax+t(i)%tmin)/2    !日均温

        total%average=total%average+t(i)%average    !日均温累加=积温

    print *,dy , t(i)%tmax , t(i)%tmin , rain , total%average    !



    enddo

    STOP

end PROGRAM



QQ截图20170131113201.png (19.33 KB, 下载次数: 366)

QQ截图20170131113201.png
作者: vvt    时间: 2017-1-31 15:42
你可以：
1. 只读取 start 到 end 之间的数据，然后计算。
2. 读取所有数据，然后只计算 start 到 end 之间的数据。

请选择其中一种方法。

如果你选用第一种方法，请先跳过第 1 到 第 start-1 的数据。
如果你选用第二种方法，请先完成所有数据的读取，然后再计算。（而不是一边读取一边计算）
作者: wxy    时间: 2017-1-31 18:07
本帖最后由 wxy 于 2017-1-31 18:40 编辑
vvt 发表于 2017-1-31 15:42
你可以：
1. 只读取 start 到 end 之间的数据，然后计算。
2. 读取所有数据，然后只计算 start 到 end 之间 ...

成功了 谢谢，但是还有两个问题 1.运行结果中把之前每行的均值累加结果都显示出来 我只想要输入的start,end之间所有行的和，也就是最后一行的结果才是我想要的，之前的代码显示的就只是一个结果，数据计算这部分的代码我并没有改动怎末显示的不一样呢
2.我用的是第一个方法，你提到的第二个方法在读取所有数据后，对数组进行计算，我想尝试一下，但脑子还是转不过来，数组也有行数列数，怎末对数组指定行之间进行计算呢，我看书上的例子都是对某一行进行操作


QQ截图20170131182927.png (22.54 KB, 下载次数: 372)

QQ截图20170131182927.png
作者: vvt    时间: 2017-1-31 19:31
1.我不清楚你现在的代码是怎样的？数据文件是怎样的？你做了什么改动。因此，你的第一个问题我没办法回答。
2.你可以从你需要的开始，循环到你结束的位置，进行求和（或其他操作）。并不一定必须从数组的第一个循环到最后一个。
作者: wxy    时间: 2017-1-31 21:20
这是现在的代码
[Fortran] 纯文本查看 复制代码
module typedef

    type tdd

        integer :: stat=0

        real ::  tmax , tmin ,  total , average

    end type

end module



PROGRAM MAIN

    USE typedef

    implicit none

    real::rain

    INTEGER,PARAMETER::temp=366

    type(tdd)::t(temp)! 保存气温值

    type(tdd):: total !日均温累加

    integer i , dy , start , end



    open(7,file="CN019067")

    read(7,*)

    write(*,*)"start:"

    read(*,*) start

    write(*,*)"end:"

    read(*,*) end



    do i=1,start-1

        read(7,*)

    enddo

    do i=start,end

        read(7,*)dy , t(i)%tmax , t(i)%tmin , rain



        t(i)%average=(t(i)%tmax+t(i)%tmin)/2    !日均温

        total%average=total%average+t(i)%average    !日均温累加=积温

    print *, dy,total%average    !

    enddo

    STOP

end PROGRAM

QQ截图20170131211931.png (10.44 KB, 下载次数: 424)

QQ截图20170131211931.png
作者: wxy    时间: 2017-1-31 21:21
本帖最后由 wxy 于 2017-1-31 21:33 编辑
vvt 发表于 2017-1-31 19:31
1.我不清楚你现在的代码是怎样的？数据文件是怎样的？你做了什么改动。因此，你的第一个问题我没办法回答。 ...

这是现在的代码和文件，一开始发的代码 比如1-50行求和，会只显示50行和的值，这个代码会显示每一行跟前几行的和，直到第50行，
[Fortran] 纯文本查看 复制代码
module typedef

    type tdd

        integer :: stat=0

        real ::  tmax , tmin ,  total , average

    end type

end module



PROGRAM MAIN

    USE typedef

    implicit none

    real::rain

    INTEGER,PARAMETER::temp=366

    type(tdd)::t(temp)! 保存气温值

    type(tdd):: total !日均温累加

    integer i , dy , start , end



    open(7,file="CN019067")

    read(7,*)

    write(*,*)"start:"

    read(*,*) start

    write(*,*)"end:"

    read(*,*) end



    do i=1,start-1

        read(7,*)

    enddo

    do i=start,end

        read(7,*)dy , t(i)%tmax , t(i)%tmin , rain



        t(i)%average=(t(i)%tmax+t(i)%tmin)/2    !日均温

        total%average=total%average+t(i)%average    !日均温累加=积温

    print *, dy,total%average    !

    enddo

    STOP

end PROGRAM

QQ截图20170131211931.png (10.44 KB, 下载次数: 436)

QQ截图20170131211931.png
作者: vvt    时间: 2017-1-31 21:52
本帖最后由 vvt 于 2017-1-31 21:54 编辑

把 print 放到循环外面不就可以了？
[Fortran] 纯文本查看 复制代码
  do i=start,end

    read(7,*)dy , t(i)%tmax , t(i)%tmin , rain

    t(i)%average=(t(i)%tmax+t(i)%tmin)/2    !日均温

    total%average=total%average+t(i)%average    !日均温累加=积温

  enddo

  print *, dy,total%average    !

如果你要全部365个都读完，然后再计算start到end的平均值，可以这样
[Fortran] 纯文本查看 复制代码
PROGRAM MAIN

  USE typedef

  implicit none

  real::rain

  INTEGER,PARAMETER::temp=366

  type(tdd)::t(temp)! 保存气温值

  type(tdd):: total !日均温累加

  integer i , dy , start , end



  open(7,file="CN019067")

  read(7,*)

  write(*,*)"start:"

  read(*,*) start

  write(*,*)"end:"

  read(*,*) end

  do i=1,temp

    read(7,*)dy , t(i)%tmax , t(i)%tmin , rain

    t(i)%average=(t(i)%tmax+t(i)%tmin)/2    !日均温

  end do

  total%average = 0

  do i=start,end

    total%average=total%average+t(i)%average    !日均温累加=积温

  enddo

  print *, dy,total%average    !

end PROGRAM

作者: wxy    时间: 2017-1-31 22:04
本帖最后由 wxy 于 2017-1-31 22:19 编辑
vvt 发表于 2017-1-31 21:52
把 print 放到循环外面不就可以了？

嗯嗯 太感谢了 然后就是请问一下 temp=366，366其实是行数 ，我想用获取行数的代码赋给temp，理论上是否是对的，temp=366是声明数组大小吧，获取行数代码放在后面赋值给temp 好像顺序反了，如果理论上可以 应该放在什么位置
[Fortran] 纯文本查看 复制代码
I=0

  do while (.true.)

    READ(8,*,end=100) line

    i=i+1

  end do

  100 CONTINUE

   m=i-4

  rewind(8)


作者: vvt    时间: 2017-1-31 23:46
本帖最后由 vvt 于 2017-1-31 23:48 编辑

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

  USE typedef

  implicit none

  real::rain

  type(tdd),allocatable::t(:)! 保存气温值

  type(tdd):: total !日均温累加

  integer i , m , dy , start , end

  character(len=3) :: line

  open(7,file="CN019067")

  I=0

  do

    READ(7,*,end=100) line

    i=i+1

  end do

  100 m = i - 1 !//你这个前面只多了一行

  rewind(7)

  allocate( t(m))  

  read(7,*)

  write(*,*)"start:"

  read(*,*) start

  write(*,*)"end:"

  read(*,*) end

  do i=1,m

    read(7,*)dy , t(i)%tmax , t(i)%tmin , rain

    t(i)%average=(t(i)%tmax+t(i)%tmin)/2    !日均温

  end do

  total%average = 0

  do i=start,end

    total%average=total%average+t(i)%average    !日均温累加=积温

  enddo

  print *, dy,total%average    !

end PROGRAM
作者: wxy    时间: 2017-2-1 08:55
vvt 发表于 2017-1-31 23:46
[mw_shl_code=fortran,true]PROGRAM MAIN
  USE typedef
  implicit none

谢谢，光照搬书上例子，把动态数组给忘了，还是要活学活用啊
作者: vvt    时间: 2017-2-1 10:02
你欠缺一点对程序的理解。
文件的行数，必须运行的时候才能判断，然后才能知道数组的大小。
运行以后才知道大小的数组，就得用动态数组。
有两个时间很重要：编译时，和运行时。

作者: wxy    时间: 2017-2-1 18:56
vvt 发表于 2017-2-1 10:02
你欠缺一点对程序的理解。
文件的行数，必须运行的时候才能判断，然后才能知道数组的大小。
运行以后才知道 ...

嗯嗯 明白了 谢谢
作者: wxy    时间: 2017-2-2 12:17
vvt 发表于 2017-2-1 10:02
你欠缺一点对程序的理解。
文件的行数，必须运行的时候才能判断，然后才能知道数组的大小。
运行以后才知道 ...

我有个关于批量处理的问题发了个新帖 麻烦帮忙看一下吧



欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2