NCL画假相当位温图

大仓鼠

天玖Y 发表于 2019-4-22 20:03
在最后画图的plot里面设置经度就行

您好,谢谢您的指导,还有个问题,请问第三行time  = fi->initial_time0_hours表示什么,初始时刻吗,我下载的grib2数据里没有这个变量

憧睿诠煦

发现ncl有计算相当位温的函数，可以看这个链接
https://www.ncl.ucar.edu/Document/Functions/Contributed/pot_temp_equiv_tlcl.shtml

syt

许市 发表于 2019-3-16 10:09
哪个函数?

我用的这个函数pot_temp_equiv_tlcl，如下：

  tlcl    = tlcl_rh_bolton (t, rh, (/1,1/))           ; degK, degK
  p       = lev;*0.01
  printMinMax(tlcl,0)
  printMinMax(p,0)
  w    = mixhum_ptrh(conform(t,p,1), t ,rh,-1 ) ;  mixing ratio (g/kg)
  copy_VarCoords(rh,tlcl)
  copy_VarCoords(rh,w)
  printVarSummary(w)
  printMinMax(w,0)
  thetae_tlcl = pot_temp_equiv_tlcl(conform(t,lev,1), t, tlcl, w, (/0,1,0,1/)) ; iounits ==>  hPa, degK, g/kg, degK
  printVarSummary(thetae_tlcl)  

Bert

syt 发表于 2020-7-19 12:08
我用的这个函数pot_temp_equiv_tlcl，如下：

  tlcl    = tlcl_rh_bolton (t, rh, (/1,1/))          ...

你好，你能提供一下用ncl画假相当位温图完整脚本吗？

许市

syt 发表于 2020-7-19 12:08
我用的这个函数pot_temp_equiv_tlcl，如下：

  tlcl    = tlcl_rh_bolton (t, rh, (/1,1/))          ...

哇，谢谢啦

syt

Bert 发表于 2020-8-25 20:28
你好，你能提供一下用ncl画假相当位温图完整脚本吗？

画的话，就和其他剖面图一样了。脚本比较乱，不好找……

Bert

syt 发表于 2020-9-6 11:56
画的话，就和其他剖面图一样了。脚本比较乱，不好找……

你好，能否提供一个其他的剖面图脚本，我想在就是ncl画剖面图，画不对

woshiwzm

syt 发表于 2020-7-19 12:08
我用的这个函数pot_temp_equiv_tlcl，如下：

  tlcl    = tlcl_rh_bolton (t, rh, (/1,1/))          ...

你好，我也想请问能提供全部脚本吗？

versatile

syt 发表于 2020-7-19 12:08
我用的这个函数pot_temp_equiv_tlcl，如下：

  tlcl    = tlcl_rh_bolton (t, rh, (/1,1/))          ...

你好，我想问一下你这里面的t和p是不是分别为tk 和pressure呀？然后lev也是pressure？

versatile

syt 发表于 2020-7-19 12:08
我用的这个函数pot_temp_equiv_tlcl，如下：

  tlcl    = tlcl_rh_bolton (t, rh, (/1,1/))          ...

还有就是我想问问这里为啥要用conform呀，tk和p的维数本来不就是相同的吗