Model of orocline dynamics