Problem: there doesn't seem to be an output variable to display the vertical temperature profile of the water column, which is a crucial diagnostic variable for assessing the lake's energy balance (and would also be crucial for biogeochemistry if someone wanted to add that to the model).
Proposed solution: create an output variable OUT_LAKE_TEMP that is an array of lake layer temperatures. Set it equal to VIC's internal lake_var->temp array. This would need to be patterned after OUT_SOIL_TEMP or some other similar variable, in terms of which files to change, how to handle the multiple dimensions, etc.
LAKES = TRUE
OUT_LAKE_TEMP
that is an array of lake layer temperatures. Set it equal to VIC's internallake_var->temp
array. This would need to be patterned afterOUT_SOIL_TEMP
or some other similar variable, in terms of which files to change, how to handle the multiple dimensions, etc.