用grads将nc格式数据转换为dat格式（z要改变）

yomi

本人的grads学的不好，处理数据遇到了瓶颈，希望能得到各位大神的指教。
我的数据是在noaa网站下载的，数据是nc格式的，在grads中打开后查看了一下ctl，ctl显示数据的z是1，t是1945（从1854年1月至2016年1月，时间间隔是一个月），现在想用fwrite命令将nc格式改为dat格式，并且想把月份放在z上（也就是希望转换之后的数据z层数是月份也就是12层，t是年份也就是原始数据对应的162年，请问在用fwrite命令的时候需要用循环还是怎样把数据转成我想要的dat格式？
在这里先谢过大神们了，求回复，求指教。

四叶草

楼主的想法理解起来还是清晰的！但是本人没有做过这类的提取，不过，回到grads数据排列方式，应该是可以实现的。两个循环嵌套，提取数据，之后自己写ctl，ctl描述中高度层描述为12个层次，时间描述为162个时次就行。楼主可以自己先尝试下。PS：楼主这样提取数据，是要固定年份，分析月份间的变化？还是？

yomi

四叶草 发表于 2016-3-19 08:54
楼主的想法理解起来还是清晰的！但是本人没有做过这类的提取，不过，回到grads数据排列方式，应该是可以实 ...

谢谢回复，我这样处理是因为要对月份的数据做滑动处理还要求每个月南海区域的区域平均，觉得这样把月份放在z上比较容易而且清楚些，我grads不是很好在网上查了nc转换dat但好像我这种的没看到，所以闷了几天也不知道怎么写，还想请问两个嵌套循环的话大概是怎样写呢？

四叶草

yomi 发表于 2016-3-19 09:43
谢谢回复，我这样处理是因为要对月份的数据做滑动处理还要求每个月南海区域的区域平均，觉得这样把月份放 ...

我的思路大致如下：
……
tt=1
while(tt<年数*12)
m=1
while(m<=12)
set t 'tt'
d var
tt=tt+1
m=m+1
endwhile
endwhile
……
以前没做过，不清楚对不对？楼主可以尝试下

臆小姐

楼主把nc文件处理成dat文件后面是要用别的软件计算的吧，那你在后面读文件的时候再把数据读成12*162的不是一样吗

yomi

四叶草 发表于 2016-3-19 10:28
我的思路大致如下：
……
tt=1

好的，太感谢你了，我试试

yomi

臆小姐 发表于 2016-3-19 10:36
楼主把nc文件处理成dat文件后面是要用别的软件计算的吧，那你在后面读文件的时候再把数据读成12*162的不是 ...

我现在没有用别的软件计算的打算，还是想用grads处理呢

yomi

四叶草 发表于 2016-3-19 10:28
我的思路大致如下：
……
tt=1

亲，还想请教你一下，你这个程序我不太明白，能不能解释一下呢？

四叶草

yomi 发表于 2016-3-20 13:58
亲，还想请教你一下，你这个程序我不太明白，能不能解释一下呢？

tt就是总的时次数，m是以12为周期进行循环控制。如果脚本没问题，内循环控制的是月份，外循环控制的是年份。