Set ctlim = 0.96 in Nygaard_2022. which corresponds to limiting the term in...

Set ctlim = 0.96 in Nygaard_2022. which corresponds to limiting the term in the sqrt in epsilon to 3 as done in the Ørsted matlab code