Refactored, and finished debugging of passing of ds/dx, ds/dy from CISM -> Albany in CISM-Albany
run.  With this capability, CISM-Albany results match well with those calculated in the CISM
glissade dycore.  This resolves the question of why there were discrepancies b/w
CISM-Albany and Glissade :).

Also modified felix_driver.cpp such that "Body Force" type is set internally in CISM-Albany
run.  If pointers to ds/dx, ds/dy are passed from CISM->Albany that are non-null,
"Body Force" is set to "FO Surface Grad Provided"; otherwise, it is set to "FO INTERP SURF GRAD"
(meaning ds/dx and ds/dy are calculated in Albany from the surface height field, as before).
Merge branch 'master' of github.com:gahansen/Albany