Debugged model code works?, save multipanel plots

Debug intensity profile computation in Sean's model code.
Remove unused code, move code for better flow, add comments.
Save circulated version of multipanel plots (with last tweaks to
various plots).  First draft of guide to model physics and code
implementation.