提问：局地温度诊断方程为什么第二项计算出来符号不对，想请教各位前辈

Jyy1108

我推导的p坐标系下局地温度诊断方程：
其中我想计算的第二项垂直非绝热加热：(Γd-Γ)ω，画出的图结果感觉是差了一个负号，想请教一下各位前辈

hhhel

没有代码没有图，不知道你的问题到底是啥。计算的时候垂直速度是omega吗，垂直速度和垂直方向单位是Pa吗

Jyy1108

hhhel 发表于 2024-3-30 10:13
没有代码没有图，不知道你的问题到底是啥。计算的时候垂直速度是omega吗，垂直速度和垂直方向单位是Pa吗

不好意思不好意思，就是感觉结果是相反的；然后垂直速度确实也是Omega；这是我的代码
f1 = addfile("E:\data\NCEP1\air.mon.mean.nc","r")
;print(f1)

f2 = addfile("E:\data\NCEP1\omega.mon.mean.nc","r")
;print(f2)

time = f1->time
it_s = 196101
it_e = 202012
YYYY = cd_calendar(time, -1)
rec_s = ind(it_s.eq.YYYY)
rec_e = ind(it_e.eq.YYYY)

air = f1->air(rec_s:rec_e,{1000:100},:,:)
;printVarSummary(air)
omega = f2->omega(rec_s:rec_e,:,:,:)
;printVarSummary(omega)

;;;;JJA;;;;;
air_JJA = month_to_season(air, "JJA")
;printVarSummary(air_JJA)
;air_JJA = air_JJA+273.15
omega_JJA = month_to_season(omega,"JJA")
;printVarSummary(omega_JJA)

p = air_JJA&level
;print(p)
;printVarSummary(p)

dtdp = center_finite_diff_n(air_JJA, p, False, 0, 1)
;printVarSummary(dtdp)

R = 287
cp = 1004
RT = R*air_JJA
;printVarSummary(RT)

P_conform = conform_dims(dimsizes(RT), p, 1)

rd = RT/(P_conform*cp)
;printVarSummary(rd)

adia = (rd-dtdp)*omega_JJA
copy_VarMeta(air_JJA, adia)
;printVarSummary(adia)

adv_850 = adia(:,{850},:,:)
printVarSummary(adv_850)

hhhel

垂直气压p应该转化为Pa？

Jyy1108

hhhel 发表于 2024-3-31 10:17
垂直气压p应该转化为Pa？

好滴，谢谢前辈，我试试