Fortran Coder

标题: 求30年的5-10月的逐日降水平均值 [打印本页]

作者: 陌曦 时间: 2020-12-3 15:56
标题: 求30年的5-10月的逐日降水平均值
有1981-2010年共30年的5-10月（184天）的逐日降水的30个二进制（bin）文件，经度为0-360，共240个格点，纬度为-90-90，共121个格点。想输出5-10月份（184天）逐日降水平均值的一个二进制文件，不知道该怎么求，求指教！

作者: 风平老涡 时间: 2020-12-3 23:51
需要知道原始数据格式，二进制只是存储格式。

作者: 陌曦 时间: 2020-12-4 19:19

风平老涡发表于 2020-12-3 23:51
需要知道原始数据格式，二进制只是存储格式。

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

Program main

Integer:: nx=240, ny=121,nt=31,i,j,m,k

real surp(240,121,31,30),dp(240,121,31)

character::  filename, year



do k=1981,2010  

   write(year,'(i4)') k

   filename=year//'.bin'

   open(11,file='F:\cut\surp-\'//filename,form='binary')

   do m=1,nt

     do j=1,ny

           do i= 1,nx

       read(12)surp(i,j,m,k-1980)

          enddo                                                        

          enddo

   enddo

   Close(11)

enddo





do i= 1,nx

do j= 1,ny

do m= 1 ,nt

   ss=0.0

   do k = 1,30

      ss= ss+surp(i,j,m,k)

   enddo

   dp(i,j,m)=ss/30.0

enddo

enddo

enddo

open(12,file='F:\surp.interim.daily.climato.bin',form='binary')

write(12)(((dp(i,j,m),i=1,240),j=1,121),m=1,31)        

End Program        main

我试图一下输出5-10月的平均，显示写出溢出，我把之前的二进制数据按月拆开后，想改成先只输出5月份的数据还是报错C:\Users\cjq\Desktop

作者: 陌曦 时间: 2020-12-4 19:20

风平老涡发表于 2020-12-3 23:51
需要知道原始数据格式，二进制只是存储格式。

运行结果显示run-time error F6600: WRITE(internal)
- internal file overflow
有什么办法，不知道该怎么搞

作者: Transpose 时间: 2020-12-4 19:57

陌曦发表于 2020-12-4 19:20
运行结果显示run-time error F6600: WRITE(internal)
- internal file overflow
有什么办法，不知道该怎 ...

character(len=20)

没有给定的话，字符串长度默认为1，不够拿来存储4位字符

作者: 陌曦 时间: 2020-12-4 20:31

Transpose 发表于 2020-12-4 19:57
character(len=20)

没有给定的话，字符串长度默认为1，不够拿来存储4位字符

天啊谢谢！！

作者: 陌曦 时间: 2020-12-5 15:48

陌曦发表于 2020-12-4 20:31
天啊谢谢！！

请问我在读取上面出来的数据时？

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

open(11,file='F:\ERA_Interim_1.5\cut5-10\surp.interim.daily.climato.bin',FORM='UNFORMATTED',ACCESS='DIRECT',action='read',RECL=NX*121)

LOOPDAY: DO K=1,nt

  print *,'Begin time = ',k

        HT  =-9.99E33        

        ffG =-9.99e33       

        tg  =-9.99e33        

        ug  =-9.99e33        

        wg  =-9.99e33        

        ft  =-9.99e33        

        fv  =-9.99e33        

        fw  =-9.99e33       

  read (11,rec=k)((orog1(i,j),i=1,nx),j=1,121)

  pres(:,1:119)=orog1(:,2:120)

  if (yrev) then

    pres(:,1:119)=pres(:,119:1:-1)

  endif

出现这种报错是为何，该怎么办呢？
Begin time = 1
run-time error F6508: READ(F:\ERA_Interim_1.5\cut5-10\surp.interim.daily.climato.bin)
- too many bytes read from unformatted record

作者: necrohan 时间: 2020-12-8 08:22
要找数据格式的说明才能读入，或者上传部分数据，或者有大神可以帮你分析

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)