每十年平均的纬向平均气温曲线

绯宜物语

lqouc 发表于 2014-8-5 15:38
理论上是互相不影响的，你先输出数据，然后disablefwrite之后再display一遍数据，然后画图就应该可以了。

我也不知道为什么这样，还是得分开.'sdfopen f:\nc\1.nc'
'set gxout fwrite'
'set fwrite f:\1new.grb'
'set x 1 720'
'set y 1 360'
'set z 1 1'
'set t 1 1320'
'define var1=ave(t2m,x=1,x=720)'
'define var2=ave(var1,t=1,t=1320)'
'set ccolor 2'
'set cmark 0'
'd var2'
'disable fwrite'
'enable print f:\t2m0.gmf'
'print'
'disable print'

绯宜物语

lqouc 发表于 2014-8-5 15:38
理论上是互相不影响的，你先输出数据，然后disablefwrite之后再display一遍数据，然后画图就应该可以了。

我想画每十年平均的纬向平均气温曲线，写了下面的gs文件：
'sdfopen f:\nc\1.nc'
'set x 1 721'
'set y 1 361'
'set z 1 1'
'set t 1 1320'
'define t2mave=ave(t2m,x=1,x=721)'
'define var1=ave(t2mave,t=1,t=120)'
'define var2=ave(t2mave,t=121,t=240)'
'define var3=ave(t2mave,t=241,t=360)'
'define var4=ave(t2mave,t=361,t=480)'
'define var5=ave(t2mave,t=481,t=600)'
'define var6=ave(t2mave,t=601,t=720)'
'define var7=ave(t2mave,t=721,t=840)'
'define var8=ave(t2mave,t=841,t=960)'
'define var9=ave(t2mave,t=961,t=1080)'
'define var10=ave(t2mave,t=1081,t=1200)'
'define var11=ave(t2mave,t=1201,t=1320)'
'set cmark 0'
'set ccolor 1'
'd var1'
'set cmark 0'
'set ccolor 2'
'd var2'
'set cmark 0'
'set ccolor 3'
'd var3'
'set cmark 0'
'set ccolor 4'
'd var4'
'set cmark 0'
'set ccolor 5'
'd var5'
'set cmark 0'
'set ccolor 6'
'd var6'
'set cmark 0'
'set ccolor 7'
'd var7'
'set cmark 0'
'set ccolor 8'
'd var8'
'set cmark 0'
'set ccolor 9'
'd var9'
'set cmark 0'
'set ccolor 10'
'd var10'
'set cmark 0'
'set ccolor 11'
'd var11'
;
可是画出来的图， 和之前画的110年平均的一样……这是怎么回事呀？可以帮我看看吗？谢谢……

绯宜物语

同时想用这个数据计算全球平均温度随年的变化
'sdfopen f:\nc\1.nc'
'set t 1 1320'
'define var=aave(t2m,x=1,x=721,y=1,y=361)'
t = 1
while(t < 1320)
'set t 't' 't+11''
'd ave(var,t='t',t='t+11')'
t=t+11  
endwhile
;
可是画不出来图……

绯宜物语

绯宜物语 发表于 2014-8-7 15:45
同时想用这个数据计算全球平均温度随年的变化
'sdfopen f:\nc\1.nc'
'set t 1 1320'

'sdfopen f:\nc\1.nc'
'set x 1'
'set y 1'
'set z 1 1'
'set t 1 1320'
'define var=aave(t2m,x=1,x=721,y=1,y=361)'
i = 1
while(i <= 110)
'set t 'i'*12-11,'i'*12'
'define var1=ave(var,t='i'*12-11,t='i'*12)'
'd var1'
endwhile
;
又写了这样的……可是运行了好久没有反应……

lqouc

嗯，才看到。
第一个程序在计算完纬向平均后加一个set t1
第二个程序要用tloop，用法自己查一下

绯宜物语

lqouc 发表于 2014-8-9 13:41
嗯，才看到。
第一个程序在计算完纬向平均后加一个set t1
第二个程序要用tloop，用法自己查一下

嗯呢，谢谢……第一个做出来了。第二个的意思是这个吗：
'sdfopen f:\nc\1.nc'
'set t 1 1320'
'define var=aave(t2m,x=1,x=721,y=1,y=361)'
i = 1
while(i <= 110)
'set t 'i'*12-11,'i'*12'
'define var1=tloop（ave(var,t='i'*12-11,t='i'*12)）'
'd var1'
endwhile
；
其实不太理解tloop到底是什么意思……

绯宜物语

lqouc 发表于 2014-8-9 13:41
嗯，才看到。
第一个程序在计算完纬向平均后加一个set t1
第二个程序要用tloop，用法自己查一下

我打算把数据输出后再做，写的Fortran程序不对，可以帮我看看嘛？program main
implicit none
integer t
real var
open(1,file='f:\1gc.grd',status='old',form='binary')
do t=1,1320
read(1) var
enddo
open(2,file='f:\1gc.txt',status='new')
write(*,*) 't2m=', var
close(1)
close(2)
end

输出的就只有一个数。

这个提示有错误，但我不太会改。
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
可以帮我看看嘛？实在弄不出来了。

绯宜物语

lqouc 发表于 2014-8-9 13:41
嗯，才看到。
第一个程序在计算完纬向平均后加一个set t1
第二个程序要用tloop，用法自己查一下

版主……我要输出1900-2009年逐年的冬季平均，写了下面的gs文件，可是第一年的数据没有输出，这个要怎么设置呀？帮我看看吧……非常感谢……

'sdfopen f:\nc\1.nc'
'set gxout fwrite'
'set fwrite f:\1wi.grd'
'set x 1 720'
'set y 1 361'
'set z 1 1'
while(n = 1)
    'define var4=ave(t2m,t='n',t='n+1')'
    'd var4'
endwhile
n=12
while(n<=1308)
'define var4=ave(t2m,t='n',t='n+2')'
n=n+12
'd var4'
endwhile
'disable fwrite'
;

绯宜物语

lqouc 发表于 2014-8-9 13:41
嗯，才看到。
第一个程序在计算完纬向平均后加一个set t1
第二个程序要用tloop，用法自己查一下

版主……我要输出1900-2009年逐年的冬季平均，写了下面的gs文件，可是第一年的数据没有输出，这个要怎么设置呀？帮我看看吧……非常感谢……

'sdfopen f:\nc\1.nc'
'set gxout fwrite'
'set fwrite f:\1wi.grd'
'set x 1 720'
'set y 1 361'
'set z 1 1'
while(n = 1)
    'define var4=ave(t2m,t='n',t='n+1')'
    'd var4'
endwhile
n=12
while(n<=1308)
'define var4=ave(t2m,t='n',t='n+2')'
n=n+12
'd var4'
endwhile
'disable fwrite'
;

lqouc

绯宜物语 发表于 2014-8-25 20:15
版主……我要输出1900-2009年逐年的冬季平均，写了下面的gs文件，可是第一年的数据没有输出，这个要怎么 ...

我看没啥问题，你可以第一年不用循环啊