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

LeoBao

'reinit'
'open d:\bylw\ncep_monmean\u.50.ctl'
'set gxout fwrite'
'set fwrite d:\bylw\ncep_monmean\uujp.grd'
k=1
while(k<=12)
'set x 1 144'
'set y 1 73'
'set z 'k''
'set t 1 '
'define b=ave(v,t+0,t=50)'
'set x 1 144'
'set y 1 73'
'set z ‘k’'
'set t 1 50'
'd v-b'
k=k+1
endwhile
'disable fwrite'
'reinit'
;
用上面的GrADS生成了一组数 x =144 ,y=73,z=12,t=50
在用fortran读数的时候，如果是这样读
read(1)(( a(i,j),i=1,144*73*12),t=50)
循环有没有错????

易小凯

你存取的顺序是x,y,t,z
最外层的循环是z

LeoBao

易小凯 发表于 2012-5-17 13:15

                               
登录/注册后可看大图

你存取的顺序是x,y,t,z
最外层的循环是z

想了很久，不知道怎么样用Grads把时间层调到最外围？能不能指导一下呢

做个霸气的木头

其实你这个数据生成的是3维：x=144、y=73、z=12，fortran读的时候设个三位数组就行了，外层循环是z，其次是y，最里层为x，ss（144,73,12）
要注意，你y读出来的比如ss（1,1,1）是第一层的纬度最低，经度最小的那个数据；

LeoBao

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

                               
登录/注册后可看大图

其实你这个数据生成的是3维：x=144、y=73、z=12，fortran读的时候设个三位数组就行了，外层循环是z，其次是 ...

我生成的数据是四维吧，fortran是老师给的多元回归，想把grads改一下，真心不知道怎么改

做个霸气的木头

LeoBao 发表于 2012-5-17 13:42

                               
登录/注册后可看大图

我生成的数据是四维吧，fortran是老师给的多元回归，想把grads改一下，真心不知道怎么改

你输出的是距平，是四维，我看错了，你想把时间输成外围是吗？

做个霸气的木头

可以把定义距平拉到外面像这样：
'reinit'
'open d:\bylw\ncep_monmean\u.50.ctl'
'set gxout fwrite'
'set fwrite d:\bylw\ncep_monmean\uujp.grd'
'set x 1 144'
'set y 1 73'
'set z 1 12'
'define b=ave(v,t+0,t=50)'
'set t 1 50'
'define c=v-b(t=1)'
'set t 1 50'
'set z 1 12'
'set x 1 144'
'set y 1 73'
'd c'
'disable fwrite'
'reinit'
;

易小凯

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

                               
登录/注册后可看大图

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

在输出的时候，还是一个场一个场的输出吧
我同学用grads连续输出的时候，发现有时会出错

易小凯

LeoBao 发表于 2012-5-17 13:17

                               
登录/注册后可看大图

想了很久，不知道怎么样用Grads把时间层调到最外围？能不能指导一下呢

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

做个霸气的木头

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

                               
登录/注册后可看大图

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

锟哥威武，今晚流行威武~