ncl用era5数据画假相当位温出现问题，请各位大神帮忙看下

zaozedi

各位大神，我用erad5数据画的假相当位温，出现好多小圈，设置的等值线338以上的也没画出来，换了时次画338以上也没画出来。但是用grads可以画出来呢，也没有这么多小圈，不知道问题出现在哪了
f     = addfile("/cygdrive/E/30.nc","r")
        lon = f->longitude
  lat = f->latitude
  tk=short2flt(f->t(10,{850},:,:))
        rh=short2flt(f->r(10,{850},:,:))
  printVarSummary(rh)

        prs = 850
  tc = tk - 273.16
  copy_VarMeta(tk, tc)
  tc@units = "degC"
  es =6.112*exp(17.67*tc/(tc+243.5))
  qs=0.622*es/(prs-0.378*es)
  q=rh*qs/100
  e=prs*q/(0.622+q)
  tlcl=55.0+2840.0/(3.5*log(tk)-log(e)-4.805)
  theta=tk*((1000/prs)^(0.2854*(1.0-0.28*q)))
  se=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))
  se=smth9_Wrap(se,0.50,0.25,True)
  copy_VarMeta(rh,se)
        wks = gsn_open_wks("png","/cygdrive/E/se")

  res = True
  res@cnLinesOn = True
  res@gsnAddCyclic = False
  ;res@cnFillOn = True
  ;res@cnLevelSelectionMode = "ManualLevels"
  ;res@cnMinLevelValF = 320
  ;res@cnMaxLevelValF = 420
  ;res@cnLevelSpacingF = 4

  res@cnLevelSelectionMode  = "ExplicitLevels"
  res@cnLevels              =(/320,324,328,330,334,338,342,346,350/)
  res@cnLevelSpacingF=2  
  res@mpMinLonF = 100
  res@mpMaxLonF = 120
  res@mpMinLatF = 20
  res@mpMaxLatF = 40
  res@mpFillOn = False
  res@gsnLeftString = ""
        plot = gsn_csm_contour_map(wks,se,res);水平分布图
end

cmh201128

你好，请问你解决了吗？是什么原因造成的呢？

zaozedi

cmh201128 发表于 2023-11-8 20:26
你好，请问你解决了吗？是什么原因造成的呢？

没有，我间隔调大调小都是这样，不知道啥原因

zaozedi

zaozedi 发表于 2023-11-14 11:11
没有，我间隔调大调小都是这样，不知道啥原因

你也遇到这样的问题吗

Eudora99

LZ，你好，请问为什么不用算子@scale_factor 和@add_offset对原始数据进行换算了呢？

lc2016

Eudora99 发表于 2023-12-15 21:47
LZ，你好，请问为什么不用算子@scale_factor 和@add_offset对原始数据进行换算了呢？

short2flt函数

lc2016

你这俩图的范围不一样，同一个范围应该是一样的