modularized some things

subfunctions return fields with results
demand shock added by shifting simulation distribution around