¥Write a separate dynamical code for determining boundary evolution (about 4
zones deep)
¥Solve magnetic induction, and
continuity equations
¥Substitute ILCT-derived flow field for
the momentum equation
(this provides the data driving)
¥Fully couple the coronal MHD code to
the boundary code