用一个ctl批量处理多个nc文件

旺旺小小酥 · 发表于 2017-4-16 17:20:45

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

有20年数据的nc文件，现在也写好了可以批量处理的ctl，但是不知道为什么在出图的时候一直显示Notice: Automatic Grid Interpolation Taking Place。Cannot contour grid - all undefined values，我一开始觉得应该是经纬度的问题，但是设置了好几遍还是出不来图，之前试都可以出来等值线图的。我贴出来我的ctl文件和我在linux系统中运行grads的情况，希望各位大神提点一下，感激不尽！

yqy1204 · 发表于 2017-4-23 10:05:32

我想提取1985-2015年的土壤湿度，数据是逐月的，但gs运行总是有问题，楼主求教啊

'reinit'
y=1985

while(y<=2015)

'set gxout fwrite'
'set fwrite E:\shuju\content\'%y%m'.dat'
'sdfopen E:\shuju\GLDAS\GLDAS_NOAH10_M.A%y%m.020.nc'

   while(m<=12)
      'set z 't1
      'set lon 70 140'
      'set lat 15 55'
      'd Soil moisture content'
      m=m+1
   endwhile

'reinit'
y=y+1
endwhile
'reinit'

旺旺小小酥 · 发表于 2017-4-19 11:25:15

问题已解决，是在linux系统中连接数据文件时不小心把软连接指向的文件改了。

yqy1204 · 发表于 2017-4-19 11:38:49

旺旺小小酥发表于 2017-4-19 11:25
问题已解决，是在linux系统中连接数据文件时不小心把软连接指向的文件改了。

您好，我想问一下20年数据的nc文件，写批量处理的ctl就可以画图了吗，可否发一份ctl文件给我，求大神

旺旺小小酥 · 发表于 2017-4-19 11:46:07

论坛里有专门的帖子说ctl的http://bbs.06climate.com/forum.php?mod=viewthread&tid=6737
另外我改好的ctl，已经在帖子里啦。

旺旺小小酥 · 发表于 2017-4-23 20:11:20

原谅我也不是很精通，光看gs文件，也看不太出来，你可以把你报错的信息也贴出来吗？
另外，这种情况我是用一个ctl处理多个nc文件，在通过define variation=ave(variation,t+初始时间，t+结束时间，时间间隔）这个命令来实现的，关于处理多个nc文件的ctl我已经贴出来了

yqy1204 · 发表于 2017-4-24 23:40:05

旺旺小小酥发表于 2017-4-23 20:11
原谅我也不是很精通，光看gs文件，也看不太出来，你可以把你报错的信息也贴出来吗？
另外，这种情况我是用 ...

可否留个联系方式

旺旺小小酥 · 发表于 2017-4-25 10:30:42

可以，加我qq吧，472578644

xiaojiu · 发表于 2017-7-4 22:47:51

感谢分享，我也遇到这个问题了

adaajiu · 发表于 2017-10-29 16:37:19

学习了，谢谢

		自动登录	找回密码
密码			立即注册

[图形美化] 用一个ctl批量处理多个nc文件

登录后查看更多精彩内容~

相关帖子