Units of output of budget routines (related to multiplication with density)

[afmoene]

Dear all,

Not a major issue, but rather a cause of confusion. In the computation of budget terms at each timestep each term (as well as the total TKE and SFS TKE) is multiplied with the appropriate density. This multiplication is, as far as I can see, related to the introduction of the anelastic approximation.
However, after the computation of the mean, at the end of the averaging period, the density is not divided out. This implies that the units of the terms stored in the output files are not m^2/s^3 but kg/m^-1/s^3 (NB: the units given now in the NetCDF are wrong anyway: m/s^2).

The question is: is the solution simply a matter of dividing the current output by a the temporal and spatial mean rhobh(k)? Or is there more to it?

Best, Arnold

dales/src/modbudget.f90

Line 158 in ff38f8d

call ncinfo(ncname( 1,:),'tker','Resolved TKE','m/s^2','tt')