# MIF 2.1

set PI [expr {4*atan(1.)}]
set MU0 [expr {4*$PI*1e-7}]

set Hx -24.6
set Hy   4.3
set Hz   0.0

Specify Oxs_BoxAtlas:atlas {
  xrange {0 500E-9}
  yrange {0 125E-9}
  zrange {0 3e-9}
}

Specify Oxs_RectangularMesh:mesh {
  cellsize {1e-9 1e-9 1e-9}
  atlas :atlas
}

Specify Oxs_UniformExchange {
  A  13E-12
}

Specify Oxs_Demag {}

Specify Oxs_FixedZeeman [subst {
   multiplier [expr 0.001/$MU0]
   comment {Convert field values from mT to A/m}
   field {$Hx $Hy $Hz}
}]

Specify Oxs_RungeKuttaEvolve:evolver {
  do_precess 1
  gamma_LL 2.21e5
  alpha 0.02
  max_timestep 0.2e-12
}

Specify Oxs_TimeDriver {
  evolver :evolver
  mesh :mesh
  stopping_time 1e-12
  stage_iteration_limit {1 {0 5000} :expand:}
  stage_count 5001
  Ms 8e5
  m0 { Oxs_FileVectorField {
    atlas :atlas
    norm  1
    file stdprob4-start.omf
  } }
}