求助！定义变量时候总是出错

zhangdee

下面这段是定义的ncl代码：

p_plane       = wrf_user_intrp3d(pf, pf,"v",plane, 90.,False)

lon_plane     = wrf_user_intrp2d(lon,plane,90.,False)

v_plane!0     ="lev"

v_plane&lev   =p_plane(:,0)

lev@long_name ="pressure"

lev@units     ="hPa"

lev@positive  ="down"

v_plane!1     ="lon"

v_plane&lon   = lon_plane

这是报错信息fatal:Variable (lev) is undefined, can not assign attribute (long_name)

这是变量输出的信息：

Variable: p_plane

Type: float

Total Size: 19200 bytes

            4800 values

Number of Dimensions: 2

Dimensions and sizes:        [Vertical | 96] x [Horizontal | 50]

Coordinates:

Number Of Attributes: 3

  _FillValue :        9.96921e+36

  description :         

  Orientation :        Cross-Section: (0,25) to (49,25) ; center=(25,25) ; angle=90

Variable: lon_plane

Type: float

Total Size: 200 bytes

            50 values

Number of Dimensions: 1

Dimensions and sizes:        [Horizontal | 50]

Coordinates:

Number Of Attributes: 3

  _FillValue :        9.96921e+36

  units :        degree_east

  description :        LONGITUDE, WEST IS NEGATIVE

Variable: v_plane

Type: float

Total Size: 19200 bytes

            4800 values

Number of Dimensions: 2

Dimensions and sizes:        [Vertical | 96] x [Horizontal | 50]

Coordinates:

Number Of Attributes: 4

  _FillValue :        9.96921e+36

  units :        m s-1

  description :        y-wind component

  Orientation :        Cross-Section: (0,25) to (49,25) ; center=(25,25) ; angle=90

想问有没有人知道错误在哪里

guoguohh

按照书上给的步骤是：
先命名v_plane!0   ="lev"
再设定lev的数值lev_pts  =p_plane(:,0)
给数值赋单位lev_pts@units    ="hPa"
最后把数值和单位赋值给v_plane&lev  =lev_pts

tSSkill

可能是lev数组没定义，建议这样试试
lev=p_plane(:,0)
v_plane&lev   =lev