爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 23584|回复: 11

[分享资料] 用Grads读netcdf的时候出现了个问题 请教ctl该怎么写

[复制链接]

新浪微博达人勋

发表于 2012-7-11 01:54:05 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
由于刚刚接触grads 很多东西不太熟悉 还请各位帮忙
研究室别人用GMT做的海表温度(SST)的nc文件
我只是想用grads将这个netcdf转换成fortran能读的binary
刚一尝试使用就遇到了问题 搜了相关内容也还没太搞清
貌似用xdfopen那个办法对我而言比较可行
但是具体的ctl文件该如何写 还请留下宝贵意见

这是问题部分
ga-> sdfopen SST_20120709_funkabay.nc
Scanning self-describing file:  SST_20120709_funkabay.nc
gadsdf: SDF file has no discernable X coordinate.
  To open this file with GrADS, use a descriptor file with an XDEF entry.
  Documentation is at http://iges.org/grads/gadoc/SDFdescriptorfile.html

用ncdump看了这个文件的信息
[puma@puma ~/share]$ ncdump -h SST_20120709_funkabay.nc
netcdf SST_20120709_funkabay {
dimensions:
     x = 555 ;
     y = 277 ;
variables:
     double x(x) ;
          x:long_name = "x" ;
          x:actual_range = 139., 144. ;
     double y(y) ;
          y:long_name = "y" ;
          y:actual_range = 40.5, 43. ;
     float z(y, x) ;
          z:long_name = "z" ;
          z:_FillValue = NaNf ;
          z:actual_range = -999., 29.2049999237061 ;

// global attributes:
          :Conventions = "COARDS/CF-1.0" ;
          :title = "/home/navi/netcdf/funkabay/sst/SST_20120709_funkabay.nc" ;
          :history = "xyz2grd -R139/144/40.5/43 -I0.00902/0.00905 /home/navi/modis_cod/funkabay/5days_binned/sst_20120709_funkabay.flat -G/home/navi/netcdf/funkabay/sst/SST_20120709_funkabay.nc -ZTLf" ;
          :GMT_version = "4.5.7 [64-bit]" ;
}
如果哪位可以帮忙 我实在感之不尽

另外后面的由netcdf转成binary那一步骤
请问是不是用lats4d 来完成 那-format 是不是选择  grads_grib 就可以

第一次发帖 请各位多多包涵


密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2012-7-11 08:55:28 | 显示全部楼层
搜一下兰溪的那个帖子吧  http://bbs.06climate.com/forum.php?mod=viewthread&tid=6008
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-7-11 10:52:51 | 显示全部楼层
很简单 我给你一个类似的例子 你研究一下
DSET F:\ncData\test\121700200.cty
title laps_cloud_type
UNDEF 1.0E37
XDEF x 251 LINEAR 51.0578613   0.2000000
YDEF y 211 LINEAR 7.3617783      0.2000000
ZDEF z 41         LINEAR 1 1
TDEF record 1 LINEAR 18jun2012 1dy
VARS 1
cty 41 99 cloud type
ENDVARS

我解释一下
以XDEF为例
XDEF x 251 LINEAR 51.0578613   0.2000000
x表示nc文件的x维。 此行是将nc中的x维定义成XDEF,后面几个你应该懂吧
如果没有Z和Time维,那就删掉这两个的定义试试。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-7-11 10:53:12 | 显示全部楼层
打开的时候
使用  xdfopen *.ctl
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-7-11 20:28:23 | 显示全部楼层

非常感谢你的帮助 不过我写了下面这样的一个ctl依然报错
z应该就是我要的 海表温度了 不过does not have any non-coordinate variables.
这是哪里出问题了呢 希望能得到您的帮忙

DSET /home/puma/share/SST_20120709_funkabay.nc
TITLE modis_daily_sst
UNDEF -999
XDEF x 555 LINEAR 139 0.0108
YDEF y 277 LINEAR 40.5 0.012635
VARS 1
z 1 0 SST
ENDVARS


ga-> xdfopen SST.ctl
Scanning Descriptor File:  SST.ctl
SDF file has no discernable time coordinate -- using default values.
gadsdf: SDF file does not have any non-coordinate variables.
SDF Descriptor file SST.ctl was not successfully opened & parsed.
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-7-11 20:45:06 | 显示全部楼层
mofangbao 发表于 2012-7-11 08:55
搜一下兰溪的那个帖子吧  http://bbs.06climate.com/forum.php?mod=viewthread&tid=6008

谢谢 麻烦能帮我看下接下来发生的情况是为什么吗?

Scanning Descriptor File:  SST.ctl
SDF file has no discernable time coordinate -- using default values.
gadsdf: SDF file does not have any non-coordinate variables.
SDF Descriptor file SST.ctl was not successfully opened & parsed.

怎么读不了我的variable呢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-7-11 22:02:50 | 显示全部楼层
仔细看了下这个
http://www.iges.org/grads/gadoc/descriptorfile.html
终于搞明白了哈 非常感谢以上两位的帮助
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-8-18 15:14:07 | 显示全部楼层
Pumaxun 发表于 2012-7-11 22:02
仔细看了下这个
http://www.iges.org/grads/gadoc/descriptorfile.html
终于搞明白了哈 非常感谢以上两位 ...

我也出现和你一样的问题,请问怎样解决的?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-8-23 10:26:19 | 显示全部楼层
我也想弄明白
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-10-31 18:22:17 | 显示全部楼层
马克一记啊啊啊啊啊啊啊啊啊啊啊啊
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表