ERA Interim等熵面数据怎么写ctl？

xry_2005

下载了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" ;

大南瓜

请问楼主是怎么下载这个等熵面资料的呀，没有找到。。。

爱喝柠檬汁

我也遇到同样的问题，请问楼主解决了吗o(╥﹏╥)o