求助！关于写nc数据的ctll描述文件

笨北极熊笨

这是ctl描述文件，数据的说明在附件里：
DSET d:\data\u\1979u.nc     
TITLE  uwnd
undef -32787s
xdef 240 linear 0 1.5
ydef 121 linear -90 1.5
zdef 37 levels 1000 975 950 925 900 875 850 825 800  775 750 700 650 600 550 500 450
               400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1                                    
               
tdef 365 linear 12:00Z1Jan1979 24h
vars 1
u 37 99 uwnd
endvars
用xdfopen命令老是出错，希望哪个好心人帮我看看。

笨北极熊笨

Type: Analysis
Type of level: Pressure levels
Date: 1979-01-01 to 1979-12-31
Time: 12:00:00
Parameter: U velocity
Step: 0
Level: 1, 2, 3, 5, 7, 10, 20, 30, 50, 70, 100 to 250 by 25, 300 to 750 by 50, 775 to 1000 by 25
Grid: 1.5° x 1.5°

netcdf 1979u {
dimensions:
        longitude = 240 ;
        latitude = 121 ;
        levelist = 37 ;
        time = UNLIMITED ; // (365 currently)
variables:
        float longitude(longitude) ;
                longitude:units = "degrees_east" ;
                longitude:long_name = "longitude" ;
        float latitude(latitude) ;
                latitude:units = "degrees_north" ;
                latitude:long_name = "latitude" ;
        int levelist(levelist) ;
                levelist:units = "millibars" ;
                levelist:long_name = "pressure_level" ;
        int time(time) ;
                time:units = "hours since 1900-01-01 00:00:0.0" ;
                time:long_name = "time" ;
        short u(time, levelist, latitude, longitude) ;
                u:scale_factor = 0.00505391109093768 ;
                u:add_offset = 29.5350156689872 ;
                u:_FillValue = -32767s ;
                u:missing_value = -32767s ;
                u:units = "m s**-1" ;
                u:long_name = "U velocity" ;

// global attributes:
                :Conventions = "CF-1.0" ;
                :history = "2011-10-15 08:49:18 GMT by mars2netcdf-0.92" ;

shangbn305

这样写好像是不对，DEST是不是应该给定二进制原始数据的文件名呢，也就是.dat呢

xo2001

NC打开也需要写ctl么？直接用sdfopen 1979u.nc不行么？

笨北极熊笨

shangbn305 发表于 2013-3-8 10:02

                               
登录/注册后可看大图

这样写好像是不对，DEST是不是应该给定二进制原始数据的文件名呢，也就是.dat呢

我上动力论坛上看了，非标准的nc文件要写描述文件，现在在做毕业设计，但这个老是不对。俺们很忧愁。

笨北极熊笨

xo2001 发表于 2013-3-8 10:15

                               
登录/注册后可看大图

NC打开也需要写ctl么？直接用sdfopen 1979u.nc不行么？

我也这样认为的，但我们的研究生学长说要写呢，不会写。

xman

非标准nc写CTL确实麻烦。更离谱的是，我一同事用grads1.9能运行的gs文件，我的2.0版本居然运行不了~~~

笨北极熊笨

xman 发表于 2013-3-8 11:10

                               
登录/注册后可看大图

非标准nc写CTL确实麻烦。更离谱的是，我一同事用grads1.9能运行的gs文件，我的2.0版本居然运行不了~~~

那怎么写呢？有没有示范的例子？

river

我记得坛子里有人发过相关帖子的，你搜索一下吧

兰溪之水

http://bbs.06climate.com/forum.php?mod=viewthread&tid=6008