(important) correct TOURSMNAR for single beta1, sigma1, beta2sp = 0, etc
dim(param) = 2q + 5 isntead of 3q + 5