Three steps left

1. main for initialize
2. write viscosity stencil
3. change in FGH