（已解决）用GrADS输出的数据fortran读取的问题

LeoBao

易小凯 发表于 2012-5-17 20:05

                               
登录/注册后可看大图

有人已经给你回复了
你试试。要是有问题的话，用循环把高度和时间进行固定，一个场一个场地输出。

问题解决了，谢谢
'reinit'
'open d:\bylw\ncep_monmean\u.50.ctl'
'set gxout fwrite'
'set fwrite d:\bylw\ncep_monmean\uujp.grd'
it=1
while(it<=50)
k=1
while(k<=17)
'set x 1 144'
'set y 1 73'
'set z 'k''
'set t 'it''
'define b=ave(v,t=1,t=50)'
'set x 1 144'
'set y 1 73'
'set z 'k''
'set t 'it''
'd v-b'
k=k+1
endwhile
it=it+1
endwhile
'disable fwrite'
'reinit'
;

LeoBao

做个霸气的木头 发表于 2012-5-17 13:56

                               
登录/注册后可看大图

可以把定义距平拉到外面像这样：
'reinit'
'open d:\bylw\ncep_monmean\u.50.ctl'

问题解决了，谢谢
'reinit'
'open d:\bylw\ncep_monmean\u.50.ctl'
'set gxout fwrite'
'set fwrite d:\bylw\ncep_monmean\uujp.grd'
it=1
while(it<=50)
k=1
while(k<=17)
'set x 1 144'
'set y 1 73'
'set z 'k''
'set t 'it''
'define b=ave(v,t=1,t=50)'
'set x 1 144'
'set y 1 73'
'set z 'k''
'set t 'it''
'd v-b'
k=k+1
endwhile
it=it+1
endwhile
'disable fwrite'
'reinit'
;

做个霸气的木头

LeoBao 发表于 2012-5-17 22:55  问题解决了，谢谢  'reinit'  'open d:\bylw\ncep_monmean\u.50.ctl'

恩恩，解决就好！

饼酱

正好需要！感谢！

张掖

这样写出的数据是二进制格式的，怎么写成nc格式的，可以sdfopen

张掖

做个霸气的木头 发表于 2012-5-17 13:56
可以把定义距平拉到外面像这样：
'reinit'
'open d:\bylw\ncep_monmean\u.50.ctl'

这样报错。too many varying dimensions

张掖

LeoBao 发表于 2012-5-17 22:50
问题解决了，谢谢
'reinit'
'open d:\bylw\ncep_monmean\u.50.ctl'

你这样行吗？我试验的，结果都是第一层的数据，lev不管怎么变化，绘制的图都是一样的

dafengche

顶楼主，虽然我也不知道答案，但在学习中……

Lufeng22

如果把这封装一下搞个函数就更好了

no_passion

感谢楼主分享，很有用