- neumann boundary for the implicit methods ...
- output stuff ... SEC, HR ...
- check ALPHA with Jean ... why is it constant this makes no sense somehow 
- fix bug when no boundary nodes are specified but boundary flags are true ...
