爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8615|回复: 5

[作图] 麻烦熟练工帮忙看下NCL文件问题出在哪里

[复制链接]

新浪微博达人勋

发表于 2012-7-26 10:26:16 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
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
data = addfile("air.mon.mean.nc","r")
T = data->air(336:696,:,:,:)
Tavgtime = dim_avg_Wrap(T(level|:,lat|:,lon|:,time|:))
Tavgtimelon = dim_avg_Wrap(Tavgtime(level|:,lat|:,lon|:))
Tavgtimelon!0 = "lev"
Tavgtimelon!1 = "lat"
result = new(73,integer)
lat=fspan(-90,90,73)
lat@units="degree_north"
result!0 = "lat"
result&lat = lat
do ilat = 0,72
  do ilev = 1,15
    if ((Tavgtimelon(ilev,ilat) < Tavgtimelon(ilev-1,ilat)) .and.(Tavgtimelon(ilev,ilat) < Tavgtimelon(ilev+1,ilat))) then
          result(ilat) = ilev
        end if
  end do
end do
wks = gsn_open_wks("pdf","tropopouse")
res = True
plot = gsn_csm_xy(wks,result&lat,result,res)
end
出来的图是一条水平直线,麻烦熟练工们帮忙看看问题所在
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-7-26 17:56:54 | 显示全部楼层
不懂想做啥,帮顶等楼下。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-7-26 22:34:47 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-7-27 08:49:08 | 显示全部楼层
realmadridhz 发表于 2012-7-26 22:34
想通过IF语句那边,通过温度递减率的变化得到对流层顶所在层次

Tavgtimelon(ilev,ilat) < Tavgtimelon(ilev-1,ilat)
错在这里,NCL语法里的“小于”应该是 “ .LT. ” 和FORTRAN类似。
“<” 这个符号是另外的用法,可以理解为“取小值”
剩下的就不用说了吧。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-7-27 09:42:12 | 显示全部楼层
hzwjy 发表于 2012-7-27 08:49
Tavgtimelon(ilev,ilat) < Tavgtimelon(ilev-1,ilat)
错在这里,NCL语法里的“小于”应该是 “ .LT. ”  ...

恩,谢谢了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-7-27 09:49:30 | 显示全部楼层
真棒啊{:soso_e163:}{:soso_e163:}
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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