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

realmadridhz

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
出来的图是一条水平直线，麻烦熟练工们帮忙看看问题所在

hzwjy

不懂想做啥，帮顶等楼下。

realmadridhz

hzwjy 发表于 2012-7-26 17:56

                               
登录/注册后可看大图

不懂想做啥，帮顶等楼下。

想通过IF语句那边，通过温度递减率的变化得到对流层顶所在层次

hzwjy

realmadridhz 发表于 2012-7-26 22:34

                               
登录/注册后可看大图

想通过IF语句那边，通过温度递减率的变化得到对流层顶所在层次

Tavgtimelon(ilev,ilat) < Tavgtimelon(ilev-1,ilat)
错在这里，NCL语法里的“小于”应该是 “ .LT. ” 和FORTRAN类似。
“<” 这个符号是另外的用法，可以理解为“取小值”
剩下的就不用说了吧。

realmadridhz

hzwjy 发表于 2012-7-27 08:49

                               
登录/注册后可看大图

Tavgtimelon(ilev,ilat) < Tavgtimelon(ilev-1,ilat)
错在这里，NCL语法里的“小于”应该是 “ .LT. ”  ...

恩，谢谢了

houtianyu1989

真棒啊{:soso_e163:}{:soso_e163:}