Add pseudoeps and pseudodt to params to control pseudo time-stepping
when finding penalization field numerically.