Currently the model engine assumes main output variables are represented as m/timestep.
Each timestep, this response variable is converted to m^3/s based on the feature's area, and then given to the nexus (here).
A serious unit problem exists if the main output variable in the BMI model is NOT provided in a units of m. The model engine should probably attempt to convert the response variable it gets here to m for internal use at this point, and if the unit is unconvertible in this case, it should likely terminate at this point with a prominent error.