quantheory / mg2-python-driver Goto Github PK

Should the flux at the ground be zero? If so, the numerical flux should be set to match this instead of being based on the k=nlev cell quantity.

The use of "dum1" for cloud ice and "dum" for cloud liquid are meant to incorporate evaporation of these elements outside of a cloud (and should enter back into the equations elsewhere), but cause a discontinuity in the numerical flux across cell interfaces.

A better approach would be to treat this dynamic as a sink and incorporate it accordingly.