SetMesh(128, 128, 1, 1e-9,1e-9,0.4e-9, 1, 1, 0)

Msat  =580e3
Aex   = 15e-12
enabledemag=false
alpha = 0.1
Ku1=0.59e6
anisU=vector(0,0,1)
Dind=0.0034089785

shiftregions=true

maxregion:=NREGION - 1
seed:=17
ext_makegrains(10e-9, maxregion, seed)

for i:=0; i<maxRegion; i++{
        for j:=i+1; j<maxRegion; j++{
                ext_ScaleExchange(i, j, 0.5)
        }
}


m =neelskyrmion(1, -1)
run(5e-11)
t=0.
minimize()

setsolver(5)
maxerr=0.001

ext_bubbleMz = -1.0
ext_centerbubble()

Xi=0.2
j=vector(5e12,2e12,0)

run(1e-10)
TOL:=1e-10
expectv("position", ext_bubblepos.average(), vector(-5.20721817016602e-08,-4.185146331787111e-08,0), TOL)

