// test relax with finite temperature


	setgridsize(128, 32, 1)
	setcellsize(500e-9/128, 125e-9/32, 3e-9)

	Msat  = 800e3
	Aex   = 13e-12
	alpha = 0.02
	m  = uniform(1, .1, 0)

	temp = 300

	relax()

	expect("T", temp.Average(), 300, 0)

