nclh画的era5假相当位温剖面图为什么是这个样子啊啊啊啊

shanglina

t=f[:]->T_GDS0_ISBL(:,9:31,{25},{70:100})
        rh=f[:]->R_GDS0_ISBL(:,9:31,{25},{70:100})
        t_mean=dim_avg_n_Wrap(t,0)
       rh_mean=dim_avg_n_Wrap(rh,0)
f1=addfile("pl7.1.00.grib","r")
p=f1->lv_ISBL0(9:31)
printVarSummary(t_mean)      
tlcl = tlcl_rh_bolton( t_mean,rh_mean,(/1,1/))
printMinMax(tlcl,0)
printMinMax(p,0)
p0:=conform(t_mean,p,0)
printVarSummary(p0)
w=mixhum_ptrh(p0,t_mean,rh_mean,1)
copy_VarCoords(rh_mean,tlcl)
copy_VarCoords(rh_mean,w)
printMinMax(w,0)
tlcl1=pot_temp_equiv_tlcl(p0,t_mean,tlcl,w,(/0,1,0,1/))
tlcl1=smth9_Wrap(tlcl1,0.50,0.25,True)
错误在哪里啊

houtian

你认为应该是什么样子？

眼里有星辰

这怎么了？

shanglina

houtian 发表于 2023-6-22 05:01
你认为应该是什么样子？

感觉这个等值线不该这么平直

shanglina

眼里有星辰 发表于 2023-6-24 15:03
这怎么了？

感觉应该有不稳定层结，不该如此稳定

眼里有星辰

shanglina 发表于 2023-6-25 20:29
感觉应该有不稳定层结，不该如此稳定

找个文献复现一下，可以复现就没问题，不能的话找找数据脚本的问题

shanglina

眼里有星辰 发表于 2023-6-27 20:36
找个文献复现一下，可以复现就没问题，不能的话找找数据脚本的问题

好的！！谢谢！

kongdechao

请问问题解决了吗

shanglina

kongdechao 发表于 2024-1-4 14:40
请问问题解决了吗

没有画这个剖面图了