Fortran读取grd文件

勇敢的悲伤

绯宜物语 发表于 2014-8-15 09:48
这样是吧？

'sdfopen f:\nc\1.nc'

大致差不多吧，你i可以试一试。。。

lysx

定义完变量var后，建议加上set x 1 721; set y 1 361。然后再d该变量，不然用fortran读完后，你会发现数据错位了。

cenhui99

这个问题我正好也要呢

绯宜物语

lysx 发表于 2014-8-15 10:57
定义完变量var后，建议加上set x 1 721; set y 1 361。然后再d该变量，不然用fortran读完后，你会发现数据 ...

嗯嗯，好的。谢谢，，

绯宜物语

cenhui99 发表于 2014-8-15 10:58
这个问题我正好也要呢

嗯嗯，一起努力学吧……

绯宜物语

勇敢的悲伤 发表于 2014-8-15 10:27
大致差不多吧，你i可以试一试。。。

好的，谢谢……

勇敢的悲伤

绯宜物语 发表于 2014-8-15 14:31
好的，谢谢……

不客气，大家一起学习嘛~

绯宜物语

勇敢的悲伤 发表于 2014-8-15 15:52
不客气，大家一起学习嘛~

写了这个程序，不知道怎么改，可以帮我看看吗？
program main
integer i,t
real var(i),year(t)
open(1,file='f:\1gc.grd',status='old',form='binary')
read(1) (var(i),i=1,1320)
do t=1,110
year(t)=(var(t*12-11)+var(t*12-10)+var(t*12-9)+var(t*12-8)+var(t*12-7)+var(t*12-6)+var(t*12-5)+var(t*12-4)+var(t*12-3)+var(t*12-2)+var(t*12-1)+var(t*12))/12
enddo
open(2,file='f:\1gc.txt',status='new')
write(*,*) 't2m=', year
close(1)
close(2)
end

紫菱

这个问题不算复杂，为什么要用fortran呢，可以直接在grads解决所有问题的。

绯宜物语

紫菱 发表于 2014-8-20 15:33
这个问题不算复杂，为什么要用fortran呢，可以直接在grads解决所有问题的。

grads和Fortran都是刚开始学，不太会……我是要画逐年平均和10年平均的图，可是用grads不出图……'sdfopen f:\nc\1.nc'
'enable print f:\nian2.gmf'
'print'
'set x 1 720'
'set y 1 361'
'set z 1 1'
'set t 1 1320'
'define var=tloop(aave(t2m,x=1,x=720,y=1,y=361))'
i = 1
a=1
b=12
while(i <= 110)
a=i*12-11
b=i*12
'set t 'a' 'b''
'define var1=ave(var,t='a',t='b')'
endwhilei=i+1
'set x 1'
'set y 1'
'set z 1 1'
'set cmark 0'
'set ccolor 2'
'd var1'
'disable print';
可是运行好久都不出来。