I am trying to figure out how to get the correct variables to calculate energy and moisture budget in WRF in a radiative-convective equilibrium (RCE) setting. For moisture budget at RCE, I used the equation:

LH = P_e

where LH = Surface latent heat flux (W m-2), so basically evaporation; and Precipitation is the sum of cumulus and grid-scale (Large-scale) precipitation (P, in mm/day), converted to W m-2 using the formula P_e = Lv * P / 86400, where P is in mm/day, and Lv is the latent heat of evaporation.

This seems to work and for my RCE setting, it largely holds that LH = P_e

However, I am having problem getting the energy budget right. There seems to be many variables in WRF and I am not sure which ones I should use.

The equation I am currently using is:

Radiative cooling, Qc = SW + LW

Heating, Qh = LH + HFX

where

SW = SWDNT + SWUPB - SWUPT - SWDNB

LW = LWDNT + LWUPB - LWUPT - LWDNB

(SWDNT is instantaneous downwelling shortwave flux at top, SWUPB is instantaneous upwelling shortwave flux at bottom, SWUPT is upwelling SW flux at top, SWDNB is downwelling SW flux at bottom, and so on)

LH = Surface latent heat flux

HFX = Sensible heat flux ("Upward heat flux at the surface")

In principle, at RCE we should get Qc = Qh, but I am getting a discrepancy of around 10 W m-2 between Qc and Qh, with more heating than cooling. This makes me wonder if I am using the wrong variables to calculate my energy budget, or that my model has not reached RCE (although I doubt this is the case as the moisture budget seems to hold?).

Would anyone be able to advice?

Thank you