求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
需要知道原始数据格式，二进制只是存储格式。

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

查看代码

打印?

01	Program main

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

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

04	character::  filename, year

05

06	do k=1981,2010

07	write(year,'(i4)') k

08	filename=year//'.bin'

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

10

do m=1,nt

11

do j=1,ny

12	do i= 1,nx

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

14

enddo

15

enddo

16

enddo

17

Close(11)

18

enddo

19

20

21	do i= 1,nx

22	do j= 1,ny

23	do m= 1 ,nt

24

ss=0.0

25	do k = 1,30

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

27

enddo

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

29

enddo

30

enddo

31

enddo

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

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

34	End Program main

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

陌曦

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

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

Transpose

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

character(len=20)

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

陌曦

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

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

天啊  谢谢！！

陌曦

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

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

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

查看代码

打印?

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

02	LOOPDAY: DO K=1,nt

03	print *,'Begin time = ',k

04	HT  =-9.99E33

05	ffG =-9.99e33

06	tg  =-9.99e33

07	ug  =-9.99e33

08	wg  =-9.99e33

09	ft  =-9.99e33

10	fv  =-9.99e33

11	fw  =-9.99e33

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

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

14	if (yrev) then

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

16

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

要找数据格式的说明才能读入，或者上传部分数据，或者有大神可以帮你分析