I realized that I hadn't added code for writing/reading SCALAR variables
in system_io.C, so I just checked in some changes that fix that. In the
process, I noticed that it looks like the if(read_additional_data) block
in System::read_parallel_data overwrites the solution vector when it
should actually be writing into the extra vectors in the system, right?
The fix is just to change solution to pos->second...