istwfk  = 1 in rf calculations