allow sim=NULL or obs=NULL; tested only with obs=NULL; implement color stratification