// geometry
	setgridsize(64, 128, 1)
	c := 3e-9
	setcellsize(c, c, c)

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

ext_makedefects(2, 3, 0.05)

anisU = vector(0, 0, 1)
Ku1.setRegion(2, 2e6)

//regions.save()
save(ku1)

// relax
	alpha = 3
	run(10e-9) 
	m_ := average(m)

save(m)
