请教一个关于 fwrite 的问题

坏拖鞋_丑小鞋

初始数据是grib1格式，ctl描述文件：options yrev
ydef 91 linear 0.000000 1
xdef 181 linear 0.000000 1.000000
初始数据是y轴反向
现用fwrite语句将数据提取出来：
'set gxout fwrite'
   zz=1
while(zz<=9)
  'set x 1 181'
  'set y 1 91'
  'set z 'zz''
  'd hgtprs'
  zz=zz+1
endwhile

对于提取出来的数据而言，重新写的ctl文件里面不需要 options yrev 这一句
这是不是说明提取出来的数据，y维变成从南到北排列
如果用fortran处理的话，第一个格点表示最南，最西的那个点？

wangyanshuo

river 发表于 2016-10-9 18:07
无论哪个grads可用的资料，只要经过它自己的命令fwrite或者sdfwrite重新输出，那资料格式必定是按照grads ...

抱歉这么晚才顾得上回复您，我专门试验了下z层次排放顺序，如果x，y，z，t多维输出的话xyz其中必须至少固定其中1维才能一次输出，否则需要加循环，所以我固定了x，我的数据是从上到下即100hPa到1000hPa这样子排列，在ctl中也是如此说明的，输出前设z=1到17然后gs一次输出那么它的顺序也是按照我写的从100hPa到1000hPa，而不会翻转成1000hPa到100hPa。所以我想楼主的问题可能仅仅存在于x，y轴反向，z轴都会按照自己设定的排列。

river

原始ctl里面已经告诉grads这个资料是南北反向排列的，当你使用这个ctl画图或者提取资料的时候，grads就会自动把它反过来的，所以最终提取出来的资料都是按照grads的规则排列的，也就是由西向东，由南到北，由底层到高层的。所以你的理解是正确的

坏拖鞋_丑小鞋

river 发表于 2013-7-28 08:12

                               
登录/注册后可看大图

原始ctl里面已经告诉grads这个资料是南北反向排列的，当你使用这个ctl画图或者提取资料的时候，grads就会自 ...

好的，谢谢！

00吴广振

最近也在考虑相关问题

jeffggg

第一个格点就是最南，最西的那个点

wangyanshuo

river 发表于 2013-7-28 08:12
原始ctl里面已经告诉grads这个资料是南北反向排列的，当你使用这个ctl画图或者提取资料的时候，grads就会自 ...

您好，请问那如果ctl编写时候我直接ydef 73 linear 90 -2.5 的话，那么输出时候也会自动翻转成grads自己的排列方式吗？

river

wangyanshuo 发表于 2016-10-9 15:43
您好，请问那如果ctl编写时候我直接ydef 73 linear 90 -2.5 的话，那么输出时候也会自动翻转成grads自己 ...

你的想法是美好的，但是那么写可以吗。你自己是不是先要验证这样写对不对，能不能出图，才有后面的问题吧。如果这么些都是错的，那后面的问题没有任何意义

wangyanshuo

river 发表于 2016-10-9 15:59
你的想法是美好的，但是那么些可以吗。你自己是不是先要验证这样写对不对，能不能出图，才有后面的问题吧 ...

哦哦我还真没试过.......一般在matlab中就给翻转过来了，倒是等压面层次可以随意来....似乎等压面会按照自己写的顺序输出吗？...谢谢提醒！

river

wangyanshuo 发表于 2016-10-9 17:14
哦哦我还真没试过.......一般在matlab中就给翻转过来了，倒是等压面层次可以随意来....似乎等压面会按 ...

无论哪个grads可用的资料，只要经过它自己的命令fwrite或者sdfwrite重新输出，那资料格式必定是按照grads的规定来的