dstemislimitswitch = .true.
dstemislimit = 1.D-4
