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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 10129|回复: 5

[作图] lon和lat的单位问题!degrees_east等

[复制链接]

新浪微博达人勋

发表于 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


密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-3-17 12:53:28 | 显示全部楼层
;Tlevel&lat = lat
; Tlevel&lon = lon去掉注释
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-3-17 13:13:40 | 显示全部楼层
小鹿kong 发表于 2015-3-17 12:53
;Tlevel&lat = lat
; Tlevel&lon = lon去掉注释

去掉注释之后是另外一个错误 QQ图片20150317131041.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-3-18 09:43:04 | 显示全部楼层
copy_VarMeta(tem, geopClm0)错了
copy_VarMeta(VarFrom, VarTo),第一个有元数据(Meta)的,第二个是没有元数据的,你正好弄反了
Wrap结尾的函数结果会保留元数据,所以不用再copy一遍(当然原来的变量必须有元数据)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-3-18 10:41:15 | 显示全部楼层
小手一挥,分数到手。。。。

评分

参与人数 1金钱 -2 收起 理由
topmad -2 只保留一个给你...

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-3-18 11:59:48 | 显示全部楼层
井中月 发表于 2015-3-18 09:43
copy_VarMeta(tem, geopClm0)错了
copy_VarMeta(VarFrom, VarTo),第一个有元数据(Meta)的,第二个是没 ...

知道问题啦!好赞的
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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