Bugfix to upgraded boundary conditions: apparently the solver does not like to have NAN values even if it does not have to evolve them; changed the '#ifdef CHECK' values in the 2nd order routines to be 1.e60 instead of NAN.