开启辅助访问 切换到窄版
搜索

Fortran Coder

 找回密码
 极速注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

[复制链接]

关于文件批量合并

[文件读写]
wxy 发表于 2017-4-22 14:23:42 查看161 回复7 阅读模式
本帖最后由 wxy 于 2017-4-22 14:25 编辑

我想将四个文件合并到一个新的文件中 ,出现图中的问题,通道号并没有冲突,不知道什麽原因,麻烦帮看一下
[Fortran] 纯文本查看 复制代码
program nc2dssat
    implicit none
    integer i ,m
    character,allocatable :: pre(:), rsd(:), tasmax(:), tasmin(:)
    character(len=512) :: filename_1, filename_2, filename_3, filename_4

    open(10,file="D:\miroc-esm-chem\historical\pr_v2\1951-1960\pre.dat") !降水
    open(11,file="D:\miroc-esm-chem\historical\rsds_v1\1951-1960\rsd.dat") !辐射
    open(12,file="D:\miroc-esm-chem\historical\tasmax_v1\1951-1960\tasmax.dat") !最高温
    open(13,file="D:\miroc-esm-chem\historical\tasmin_v1\1951-1960\tasmin.dat") !最低温

    do
        read(10,*)  filename_1
        read(11,*)  filename_2
        read(12,*)  filename_3
        read(13,*)  filename_4

        open(10,file=filename_1)
        open(11,file=filename_2)
        open(12,file=filename_3)
        open(13,file=filename_4)

        i=3652 !文件行数
        m=i
!===========================================================
!文件名配对,相同文件合并
!===========================================================
        if (filename_1==filename_2.and.filename_1==filename_3.and.filename_1==filename_4)then
            allocate(pre(m), rsd(m), tasmax(m), tasmin(m))

            open(14,file="D:\miroc-esm-chem\historical\test.txt")

            do i=1,m
                read(10,*)pre(i)
                read(11,*)rsd(i)
                read(12,*)tasmax(i)
                read(13,*)tasmin(i)
                write(14,"(5a,2x,6a,2x,6a,2x,5a)") rsd(i), tasmax(i), tasmin(i), pre(i)
            end do
            close(10)
            close(11)
            close(12)
            close(13)
            close(14)
            deallocate(rsd,tasmax,tasmin,pre)
        end if
    end do
end program

QQ截图20170422142058.png
回复

使用道具 举报

li913 发表于 2017-4-22 15:36:03
同一个文件被两次打开。
QQ截图20170422153520.png
fcode 发表于 2017-4-22 16:38:08
   open(10,file="D:\miroc-esm-chem\historical\pr_v2\1951-1960\pre.dat") !降水
   open(11,file="D:\miroc-esm-chem\historical\rsds_v1\1951-1960\rsd.dat") !辐射
   open(12,file="D:\miroc-esm-chem\historical\tasmax_v1\1951-1960\tasmax.dat") !最高温
   open(13,file="D:\miroc-esm-chem\historical\tasmin_v1\1951-1960\tasmin.dat") !最低温

这些,和

      open(10,file=filename_1)
      open(11,file=filename_2)
      open(12,file=filename_3)
      open(13,file=filename_4)

用了相同的通道号(10,11,12,13),并且同时打开了。

wxy  楼主 发表于 2017-4-22 18:38:50
本帖最后由 wxy 于 2017-4-22 19:21 编辑
fcode 发表于 2017-4-22 16:38
open(10,file="D:\miroc-esm-chem\historical\pr_v2\1951-1960\pre.dat") !降水
   open(11,file="D:\mi ...

改了通道号之后仍然是这个错误,现在还是同时打开吗
QQ截图20170422184012.png
fcode 发表于 2017-4-22 19:28:32
给出  pre.dat  , rsd.dat  , tasmax.dat  和 tasmin.dat 文件的内容
wxy  楼主 发表于 2017-4-22 19:36:35
fcode 发表于 2017-4-22 19:28
给出  pre.dat  , rsd.dat  , tasmax.dat  和 tasmin.dat 文件的内容

四个文件内容都是一样的
QQ截图20170422193747.png
wxy  楼主 发表于 2017-4-22 20:35:49
本帖最后由 wxy 于 2017-4-22 21:36 编辑
fcode 发表于 2017-4-22 19:28
给出  pre.dat  , rsd.dat  , tasmax.dat  和 tasmin.dat 文件的内容

刚才的问题解决了,麻烦了

wxy  楼主 发表于 2017-4-22 23:54:05
本帖最后由 wxy 于 2017-4-23 13:09 编辑
fcode 发表于 2017-4-22 19:28
给出  pre.dat  , rsd.dat  , tasmax.dat  和 tasmin.dat 文件的内容

还想请教一下,在批量读取五个dat文件列表时出现了问题,总是end of file,并新建一个filename_1,到filename_5的空白文件,不知道为什麽
[Fortran] 纯文本查看 复制代码
program nc2dssat
    implicit none
    integer i ,m
    character,allocatable :: pre(:), rsd(:), tasmax(:), tasmin(:)
    character(len=512) ::  filename_1, filename_2, filename_3, filename_4, filename_5, line
!==================================================================================
!文件名列表
!==================================================================================
    open(10,file="D:\miroc-esm-chem\historical\pr_v2\1951-1960\pre.dat") !降水
    open(11,file="D:\miroc-esm-chem\historical\rsds_v1\1951-1960\rsd.dat") !辐射
    open(12,file="D:\miroc-esm-chem\historical\tasmax_v1\1951-1960\tasmax.dat") !最高温
    open(13,file="D:\miroc-esm-chem\historical\tasmin_v1\1951-1960\tasmin.dat") !最低温
    open(14,file="D:\miroc-esm-chem\historical\test\historical.dat")!输入文件

    do
        read(10,*)  filename_1
        read(11,*)  filename_2
        read(12,*)  filename_3
        read(13,*)  filename_4
        read(14,*)  filename_5
!=====================================================================================
!批量打开文件
!=====================================================================================
        open(15,file="D:\miroc-esm-chem\historical\pr_v2\1951-1960\filename_1")
        open(16,file="D:\miroc-esm-chem\historical\rsds_v1\1951-1960\filename_2")
        open(17,file="D:\miroc-esm-chem\historical\tasmax_v1\1951-1960\filename_3")
        open(18,file="D:\miroc-esm-chem\historical\tasmin_v1\1951-1960\filename_4")
        open(19,file="D:\miroc-esm-chem\historical\test\filename_5")
!=====================================================================================
!文件行数
!=====================================================================================
        i = 0
        Do
          Read (9, *, End=100) line
          i = i + 1
        End Do
        100 m = i
        Rewind (9)

!===========================================================
!文件名配对,相同文件合并
!===========================================================
        if (filename_1==filename_2.and.filename_2==filename_3.and.filename_3==filename_4.and.filename_4==filename_5)then
            allocate(pre(m), rsd(m), tasmax(m), tasmin(m))

            do i=1,m
                read(15,*)pre(i)
                read(16,*)rsd(i)
                read(17,*)tasmax(i)
                read(18,*)tasmin(i)
                write(19,"(5a,2x,6a,2x,6a,2x,5a)") rsd(i), tasmax(i), tasmin(i), pre(i)
            end do
            deallocate(rsd,tasmax,tasmin,pre)
            close(15)
            close(16)
            close(17)
            close(18)
            close(19)
        end if
    end do
end program


QQ截图20170423003133.png
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则 长代码粘贴

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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