[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[IDV #NVU-263991]: Error when calculating temperature tendency
- Subject: [IDV #NVU-263991]: Error when calculating temperature tendency
- Date: Tue, 31 May 2022 14:27:44 -0600
> Hi,
>
> I'm trying to calculate this:
>
> def TEST_LHomega(t, V, omega, theta, theta1, theta2):
> T_over_theta = quo(t, theta)
> dtheta = sub(theta1, theta2)
> p1 = extractPressureFromNWPGrid(theta1)
> p2 = extractPressureFromNWPGrid(theta2)
> dp = sub(p1, p2)
> dthetadp = quo(dtheta,dp)
>
> return mul(mul(T_over_theta, dthetadp),omega)
>
> where: t=temperature (K), theta = potential temperature (K), V=true wind
> vector (m/s), omega=vertical pressure velocity (Pa/s).
>
> The final units should be K/s but I'm getting K*m/(kg*s^2), it seems that
> the units from the omega field are not reported correctly. How can I fix it?
>
> Thank you,
>
> Paula
>
If only the unit of omega field incorrected, you can use the newUnit(..) to set
the
correct unit. Example: newUnit(grid, "topo", "m")
Yuan
Ticket Details
===================
Ticket ID: NVU-263991
Department: Support IDV
Priority: High
Status: Open
===================
NOTE: All email exchanges with Unidata User Support are recorded in the Unidata
inquiry tracking system and then made publicly available through the web. If
you do not want to have your interactions made available in this way, you must
let us know in each email you send to us.