NCL读取资料问题

fealbai

刚接触ncl
读取高度场资料：
Type: float
Total Size: 32376960 bytes
            8094240 values
Number of Dimensions: 3
Dimensions and sizes:   [time | 770] x [lat | 73] x [lon | 144]
Coordinates:
            time: [17067072..17628840]
            lat: [90..-90]
            lon: [ 0..357.5]
Number Of Attributes: 18
  level :       500
  long_name :   Monthly mean geopotential height
  actual_range :        ( -334.5161, 32321.1 )
  valid_range : ( -700, 35000 )
  units :       m
  add_offset :   0
  scale_factor :         1
  missing_value :       -9.96921e+36
  precision :   0
  least_significant_digit :     0
  GRIB_id :     7
  GRIB_name :   HGT
  var_desc :    Geopotential height
  dataset :     CDC Derived NCEP Reanalysis Products
  level_desc :  Multiple levels
  statistic :   Mean
  parent_stat : Other
  _FillValue :  -9.96921e+36
(0,0,0) 4982.193
(0,0,1) 4982.193
(0,0,2) 4982.193
(0,0,3) 4982.193
(0,0,4) 4982.193
(0,0,5) 4982.193
(0,0,6) 4982.193
(0,0,7) 4982.193
(0,0,8) 4982.193
(0,0,9) 4982.193
(0,0,10)        4982.193
(0,0,11)        4982.193
(0,0,12)        4982.193
(0,0,13)        4982.193
(0,0,14)        4982.193
(0,0,15)        4982.193
(0,0,16)        4982.193
(0,0,17)        4982.193
(0,0,18)        4982.193
(0,0,19)        4982.193
(0,0,20)        4982.193
(0,0,21)        4982.193
(0,0,22)        4982.193
(0,0,23)        4982.193
(0,0,24)        4982.193
值显示是正确的
但是读取风场 相对湿度等资料：
Type: short
Total Size: 16188480 bytes
            8094240 values
Number of Dimensions: 3
Dimensions and sizes:   [time | 770] x [lat | 73] x [lon | 144]
Coordinates:
            time: [17067072..17628840]
            lat: [90..-90]
            lon: [ 0..357.5]
Number Of Attributes: 17
  level :       500
  long_name :   Monthly mean u wind
  actual_range :        ( -66.37354, 124.3942 )
  unpacked_valid_range :        ( -125, 160 )
  units :       m/s
  add_offset :  202.65
  scale_factor :        0.01
  missing_value :       32766
  precision :   2
  least_significant_digit :     1
  var_desc :    u-wind
  dataset :     CDC Derived NCEP Reanalysis Products
  level_desc :  Multiple levels
  statistic :   Mean
  parent_stat : Other
  valid_range : ( -32765, -4265 )
  _FillValue :  32766
(0,0,0) -20393
(0,0,1) -20371
(0,0,2) -20349
(0,0,3) -20327
(0,0,4) -20305
(0,0,5) -20283
(0,0,6) -20260
(0,0,7) -20238
(0,0,8) -20216
(0,0,9) -20194
(0,0,10)        -20172
(0,0,11)        -20150
(0,0,12)        -20128
(0,0,13)        -20107
(0,0,14)        -20086
(0,0,15)        -20065
(0,0,16)        -20045
(0,0,17)        -20025
(0,0,18)        -20006
(0,0,19)        -19987
(0,0,20)        -19968
(0,0,21)        -19950
(0,0,22)        -19933
(0,0,23)        -19917
(0,0,24)        -19901
(0,0,25)        -19886
值都是这样的 问题出在哪？怎样解决
不胜感激

兰溪之水

u=f->u
u=u*scale_factor + add_offset

fealbai

兰溪之水 发表于 2012-4-3 19:33

                               
登录/注册后可看大图

u=f->u
u=u*scale_factor + add_offset

非常感谢！
我按照你写的编入脚本
load "/usr/local/ncl_ncarg-5.1.0/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "/usr/local/ncl_ncarg-5.1.0/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/local/ncl_ncarg-5.1.0/lib/ncarg/nclscripts/csm/contributed.ncl"
load "/usr/local/ncl_ncarg-5.1.0/lib/ncarg/nclscripts/csm/shea_util.ncl"

begin

        in1 =addfile("/scratch/student/bozk/uwnd.mon.mean.nc","r")
        in2 =addfile("/scratch/student/bozk/vwnd.mon.mean.nc","r")
       
        t  = in1->time
        lev= in1->level
        lat= in1->lat
        lon= in1->lon

        u  = in1->uwnd(:,5,:,:)
        u=u*scale_factor + add_offset

        v  = in2->vwnd(:,5,:,:)
        v=v*scale_factor + add_offset
end
这是报的错：
fatal:Variable (scale_factor) is undefined
fatal:Execute: Error occurred at or near line 16
这是哪的问题？
还有麻烦告诉我 加的命令是什么意思？

兰溪之水

u  = in1->uwnd(:,5,:,:)
U=u*u@scale_factor + u@add_offset

fealbai

兰溪之水 发表于 2012-4-3 20:30

                               
登录/注册后可看大图

u  = in1->uwnd(:,5,:,:)
U=u*u@scale_factor + u@add_offset

非常感谢！感激不尽，已经可以了，{:eb301:}

hzwjy

注意看变量类型
Type: float
Type: short

如果变量类型是非float的，可能就需要转换一下。
有相应的函数short2flt, byte2flt这些。

扬风三郎

正好出现问题，谢谢知道～～～

zhouderongneo

感谢楼主，感谢啦

hxyj

呵呵，好东西，刚刚被这给困惑了，谢谢。

空谷幽竹

兰溪之水 发表于 2012-4-3 20:30

                               
登录/注册后可看大图

u  = in1->uwnd(:,5,:,:)
U=u*u@scale_factor + u@add_offset

大神，为什么直接读出来不行，还要进行加减乘除哇