处理ncep日资料内存不够

chan_hsu

chuliziliao-ncep - 副本.f90 (1.24 KB, 下载次数: 0)

2015-4-8 15:26 上传

点击文件名下载附件
请问大家 fortran处理ncep日资料内存不够，警告为“total image size 461840384 exceeds max (268435456); image may not run” 该怎么处理？

1. program hgt500
   
2. implicit none
   
3. integer,Parameter::x=144,y=73,M=53 !1961-2013年
   
4. integer no_leap(12),leap(12)
   
5. integer i,j,dy,mo,yr,nyr
   
6. real h(x,y,31,12,M)
   
7. data leap/31,29,31,30,31,30,31,31,30,31,30,31/
   
8. data no_leap/31,28,31,30,31,30,31,31,30,31,30,31/
   
9. 
   
10. open(10,file='hgt.grd',form='binary')
    
11. do yr=1,M   
    
12.         nyr=yr+1960
    
13.                 do mo=1,12
    
14.                         !闰年
    
15.                         if(mod(nyr,4)==0) then
    
16.                                 do dy=1,leap(mo)
    
17.                                         do j=1,y   
    
18.                                                 do i=1,x
    
19.                                                         read(10) h(i,j,dy,mo,yr)
    
20.                                 enddo;enddo;enddo
    
21. !非闰年
    
22.                         else  
    
23.                                 do dy=1,no_leap(mo)
    
24.                                         do j=1,y   
    
25.                                                 do i=1,x
    
26.                                                         read(10) h(i,j,dy,mo,yr)
    
27.                                 enddo;enddo;enddo
    
28.                         endif
    
29.                 enddo;enddo
    
30. close(10)
    
31. print*,'read ok'
    
32. 
    
33. open(30,file='hgt_may-aug.grd',form='binary')
    
34.         do yr=1,M   
    
35.                 nyr=yr+1960
    
36.                         do mo=5,8
    
37.                                 !闰年
    
38.                                 if(mod(nyr,4)==0) then
    
39.                                         do dy=1,leap(mo)
    
40.                                                 do j=1,y   
    
41.                                                         do i=1,x
    
42.                                                                 write(30) h(i,j,dy,mo,yr)
    
43.                                         enddo;enddo;enddo
    
44.                                 !非闰年
    
45.                                 else  
    
46.                                         do dy=1,no_leap(mo)
    
47.                                                 do j=1,y   
    
48.                                                         do i=1,x
    
49.                                                                 write(30) h(i,j,dy,mo,yr)
    
50.                                         enddo;enddo;enddo
    
51.                                 endif
    
52.            enddo;enddo
    
53. 
    
54. end program

复制代码

lqouc

warning又不是报错，无视掉直接运行，如果不行就调大堆栈（根据你的编译器，具体方法论坛搜下），再不行就吧数据分开算喽。

言深深

你这种可以：
1，忽视warning看看；
2，如果无效，降低纬度，你这里甚至都可以不用数组都行。一边read，一边write

baiduwenku

在win7 64位下编译运行，64位往往支持的内存更多一些。
IVF2011+VS2010+WIN7 64的环境。