ncl计算假等值线相当位温总是显示没有图形

没有，如果

f=addfile("d:\zh.nc","r")

tm=short2flt(f->t(92,{850},:,:))
rh=short2flt(f->r(92,{850},:,:))
p=850.

wks=gsn_open_wks("png","theta-850-24-08")

es=6.112*exp(17.67*(tm-273.15)/(tm-29.65))
q=rh*(0.62197*es/(p-es))/100.
e=p*q/(0.62197+q)+10^(-10)
tlcl=55.0+2840.0/(3.5*log(tm)-log(e)-4.805)
theta=tm*(1000./p)^(0.2854*(1.0-0.28*q))      
eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))

copy_VarCoords(tm,eqt) ; assign coordinates
;eqt!0        ="time"   
;eqt&time =tm&time   
;eqt!1        ="level"  
;eqt&level=tm&level
;eqt!0       ="latitude"   
;eqt&latitude =tm&latitude   
;eqt!1       ="longitude"  
;eqt&longitude=tm&longitude


res=True
res@mpMinLatF             =  17          ; Asia limits
res@mpMaxLatF             =  55
res@mpMinLonF             =  72
res@mpMaxLonF             = 136


res@mpLimitMode       = "Corners"
res@mpLeftCornerLatF  = 15
res@mpRightCornerLatF = 50
res@mpLeftCornerLonF  = 90
res@mpRightCornerLonF = 130

;res@mpLambertParallel1F = 10   
;res@mpLambertParallel2F = 6   
;res@mpLambertMeridianF  =104.06
res@mpDataBaseVersion   ="MediumRes"
res@mpDataSetName       ="Earth..4"
res@mpOutlineOn         =True
res@mpFillOn            =False
res@mpOutlineSpecifiers =(/"China:states","Taiwan"/)

res@pmTickMarkDisplayMode="Always"
res@tmXTOn               =False
res@tmYROn               =False
res@gsnLeftString        = ""
res@gsnRightString       = ""

res@cnLineLabelsOn       =True

plot=gsn_csm_contour_map(wks,eqt,res)



运行结果：
warning:ContourPlotInitialize: scalar field is constant; no contour lines will appear; use cnConstFEnableFill to enable fill

pandasp1213

一步步检查计算的变量，看看哪一步出问题了

没有，如果

pandasp1213 发表于 2018-8-25 09:37
一步步检查计算的变量，看看哪一步出问题了

因为是初学又是自己自学，关于ncl变量的使用还不是很清楚；冒昧请教是不是需要对变量的维数进行定义，是否是计算过程中变量的维数不匹配。

pandasp1213

没有，如果 发表于 2018-8-29 15:04
因为是初学又是自己自学，关于ncl变量的使用还不是很清楚；冒昧请教是不是需要对变量的维数进行定义，是 ...

如果维数不匹配，运行会报错的。输出看一下每一步计算的变量数值有没有问题

没有，如果

pandasp1213 发表于 2018-8-30 07:08
如果维数不匹配，运行会报错的。输出看一下每一步计算的变量数值有没有问题

Variable: es
Type: float
Total Size: 260640 bytes
            65160 values
Number of Dimensions: 2
Dimensions and sizes:   [lat | 181] x [lon | 360]
Coordinates:
            lat: [90..-90]
            lon: [ 0..359]
Number Of Attributes: 1
  _FillValue :  -9.99e+33
查看计算变量显示相关信息只有缺省值，是否是代表变量没有正常进行计算得到相关数值？

pandasp1213

可以用printMinMax函数简单看一下，如果是缺省值，那就有问题了

willshordan

真的很有用的帖子