Fix flux units of output data

cbegeman commented 5 years ago

Contains two changes to data output method:

we modify units of resolved fluxes so they are consistent with user-defined input and output units. Both resolved and SGS fluxes are computed as rho_0 * w'X' * output_conversion_factor.
scalarflux_output_conversion also needs to take units 1/rho_0 as any flux input to diffusion_s is assumed to take the form rho_0 * w'X
For MOST method 'mcphee', we output T,S at the top boundary in profiles. This requires excluding the nzt point from the post-facto calculation of resolved w'T' and w'S'

xylar commented 5 years ago

I don't really feel qualified to review this one. Feel free to suggest someone else if you want a review.

cbegeman commented 5 years ago

@qingli411 @lconlon Would either of you like to review this PR?

lanl / palm_lanl