请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3825|回复: 2

[脚本编辑] ERA Interim等熵面数据怎么写ctl?

[复制链接]

新浪微博达人勋

发表于 2016-4-13 22:41:22 | 显示全部楼层 |阅读模式

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

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

x
下载了ERA Interim等熵面数据,是不规则的nc资料,用ncdump后写了ctl,但打开后提示出错,请各位指点!
ctl如下:
dset D:\netcdf\is201305.nc
title  is
undef 1.e30
dtype netcdf
xdef  321 linear  60.0  0.125
ydef  241 linear  30.0  0.125
zdef  15 levels  850 700 600 530 475 430 395 370 350 330 315 300 285 275 265
tdef  124 linear 00Z01may2013  6hr      
VARS   7
pres    15 99 Pressure(Pa)
pv      15 99 Potential vorticity(K m**2 kg**-1 s**-1)
q       15 99 Specific humidity(kg kg**-1)
vo      15 99 Vorticity (relative)(s**-1)
d       15 99 Divergence(s**-1)
u       15 99 U component of wind(m s**-1)
v       15 99 V component of wind(m s**-1)
ENDVARS

文件头如下:
netcdf is201305 {
dimensions:
        longitude = 321 ;
        latitude = 241 ;
        level = 15 ;
        time = UNLIMITED ; // (124 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 level(level) ;
                level:long_name = "levelist" ;
        int time(time) ;
                time:units = "hours since 1900-01-01 00:00:0.0" ;
                time:long_name = "time" ;
                time:calendar = "gregorian" ;
        short pres(time, level, latitude, longitude) ;
                pres:scale_factor = 1.63915574157891 ;
                pres:add_offset = 54680.2877219339 ;
                pres:_FillValue = -32767s ;
                pres:missing_value = -32767s ;
                pres:units = "Pa" ;
                pres:long_name = "Pressure" ;
        short pv(time, level, latitude, longitude) ;
                pv:scale_factor = 7.90401599271671e-009 ;
                pv:add_offset = 0.000207577008696949 ;
                pv:_FillValue = -32767s ;
                pv:missing_value = -32767s ;
                pv:units = "K m**2 kg**-1 s**-1" ;
                pv:long_name = "Potential vorticity" ;
        short q(time, level, latitude, longitude) ;
                q:scale_factor = 3.42865958892848e-007 ;
                q:add_offset = 0.0112343460090831 ;
                q:_FillValue = -32767s ;
                q:missing_value = -32767s ;
                q:units = "kg kg**-1" ;
                q:long_name = "Specific humidity" ;
                q:standard_name = "specific_humidity" ;
        short vo(time, level, latitude, longitude) ;
                vo:scale_factor = 1.10947607275713e-008 ;
                vo:add_offset = 7.38325657864937e-005 ;
                vo:_FillValue = -32767s ;
                vo:missing_value = -32767s ;
                vo:units = "s**-1" ;
                vo:long_name = "Vorticity (relative)" ;
                vo:standard_name = "atmosphere_relative_vorticity" ;
        short d(time, level, latitude, longitude) ;
                d:scale_factor = 6.73893680605735e-009 ;
                d:add_offset = -1.43547884087534e-005 ;
                d:_FillValue = -327press67s ;
                d:missing_value = -32767s ;
                d:units = "s**-1" ;
                d:long_name = "Divergence" ;
                d:standard_name = "divergence_of_wind" ;
        short u(time, level, latitude, longitude) ;
                u:scale_factor = 0.00148065684688775 ;
                u:add_offset = 23.9375044912177 ;
                u:_FillValue = -32767s ;
                u:missing_value = -32767s ;
                u:units = "m s**-1" ;
                u:long_name = "U component of wind" ;
                u:standard_name = "eastward_wind" ;
        short v(time, level, latitude, longitude) ;
                v:scale_factor = 0.00152109774294892 ;
                v:add_offset = -0.233609669965224 ;
                v:_FillValue = -32767s ;
                v:missing_value = -32767s ;
                v:units = "m s**-1" ;
                v:long_name = "V component of wind" ;
                v:standard_name = "northward_wind" ;
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-3-22 10:46:43 | 显示全部楼层
请问楼主是怎么下载这个等熵面资料的呀,没有找到。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-12-15 21:24:22 | 显示全部楼层
本帖最后由 爱喝柠檬汁 于 2018-12-15 21:25 编辑

我也遇到同样的问题,请问楼主解决了吗o(╥﹏╥)o
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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