爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4401|回复: 3

[求助] fortran批处理读取文件出现invalid real问题

[复制链接]

新浪微博达人勋

发表于 2018-4-19 21:02:01 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
program main
implicit none
integer,parameter::dy=365,nfile=16
real,parameter::a=28.5,b=23.5,c=33.5
real,dimension(dy,nfile)::var
real,dimension(dy)::u,v
integer x,y,l
character*200 filename
open(11,file='D:\grads\Year\1987\input.txt')
do x=1,nfile,1
read(11,*) filename
open(12,file='D:\grads\Year\1987\1987P01.dat',form='binary')
read(12,*) ((var(y,x),y = 1,dy),x=1,nfile)
! print *,var
enddo
close(11)
close(12)

open(3,file='D:\grads\TXT\1989\u.txt',action='write',form='formatted',status='new')
open(4,file='D:\grads\TXT\1989\v.txt',action='write',form='formatted',status='new')  
do y=1,nfile
    do x=1,dy
       if (var(x,y) == -9.96921e+36)then
          l=x
         var(l,y)=0
     endif
    enddo      
enddo

do x=1,dy
       u(x) = 0.5*(var(x,12) + var(x,13) - var(x,4) - var(x,5))
       v(x) = 0.25/cos(a) * (var(x,5) + 2*var(x,9) + var(x,13) - var(x,4) - 2*var(x,8) - var(x,12))
  write(3,*) u(x)
  write(4,*) v(x)
enddo
close(3)
close(4)
end program main

各位叔叔阿姨哥哥姐姐和小伙伴们,上述是本人在批量处理文件时的程序,input.txt内为所有data数据的路径+文件名,一共16行,即每一年有16个data数据。读取过程中出现了一个问题,如下图

注:data文件转化成txt文件后为一列real型数据,无特殊字符和描述文件说明出现,写了ctl配合用nc也可以直接画出图。
想了好几天,问了同学也没想明白问题出在哪里?拜托各位大神帮帮忙,小生感激不尽。

提示出现的问题

提示出现的问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-5-3 16:45:38 | 显示全部楼层
没有人帮帮我麽
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-5-3 16:45:44 | 显示全部楼层
{:5_275:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2018-5-3 19:37:56 | 显示全部楼层

我也想学习一下,发一份给我吧,谢啦!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

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