fix for variable material properties in flux
Fixing initial condition.
boundary condition implementation
making output less fequent