lon和lat的单位问题！degrees_east等

lavenderSC361 · 发表于 2015-3-17 12:27:44

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

用NCL编程计算两层之间的平均温度，程序如下：
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;************************************************
begin
;************************************************
  g=9.80665
  R=287
  ln=log(500./200.)
  item=g/(R*ln)
; print(ln)
;************************************************
; open netCDF file
  diri="/cygdrive/J/NCEP/height/"
  fili="hgt.mon.mean.nc"
f= addfile(diri+fili, "r")
geop  =f->hgt(84:695,5:9,{10:45},{80:150})    ;195501-200512 500hPa-200hPa
; printVarSummary(geop)
time=f->time
lev=f->level
lon=f->lon
lat=f->lat
  ;lat@units = "degrees_north"
  ;lon@units = "degrees_east"
geopClm0 = clmMonTLLL( geop )  ;气候平均
;printVarSummary(geopClm0)

;*******计算两层之间平均温度*************

  tem=(geopClm0(:,4,:,:)-geopClm0(:,0,:,:))*item-273.15
copy_VarMeta(tem, geopClm0)
missing=geop@missing_value

printVarSummary(tem)
;*********计算区域距平*************
Tlevel=dim_rmvmean_n_Wrap(tem(4,:,:),0)
copy_VarMeta(Tlevel,geop)
; print(Tlevel)
  Tlevel@missing_value=  missing
  Tlevel@_FillValue=  missing
  Tlevel!0 ="lat"
  Tlevel!1 ="lon"
  ；Tlevel&lat = lat
； Tlevel&lon = lon
  lat@units = "degrees_north"
  lon@units = "degrees_east"
; print(tem@missing_value)

printVarSummary(Tlevel)

计算可以得出结果，但是总是提示没有加上单位，错误如下 QQ图片20150317122356.png