sensible heat and latent flux output from satmedmfvdifq.F are out of range after first time step.

[edwardstrobach]

Good evening,

I wrote out a diagnostic for the gabls case to examine the heat(i) and evap(i) fields in satmedmfvdifq.F. I wanted to examine the time evolution of these fields for different experiments since this field is used to calculate the convective velocity scale and buoyancy production. I'm not sure what's going on with these fields but I'm seeing the following in the output:

Dataset('output.nc').variables['heat_2d'][:]
array([[ 1.20925240e-01],
[ 9.96920997e+36],
[ 9.96920997e+36],
[ 9.96920997e+36],
[ 9.96920997e+36],
[ 9.96920997e+36],
[ 9.96920997e+36],

Note, the buoyancy production and convective velocity fields look okay. This gives me the impression that the heat(i) and evap(i) fields are not truly out of range.