Modify the initialization of H in the double arch geometry model

[DaehoHan]

The previous implementation gives weird PESs, since H at x=-Z is undefined.