*heading
Testing basic EM capability - Expected results: A_r + i A_i
Where A_r = A_i = -y e_1 + x e_2
*preprint,model=yes,history=yes
***restart,write,freq=100
*parameter
cond_fct= 1.e-8
cond_cu = 58.0e6
conduct = cond_fct*cond_cu
** mesh related
nelsx = 10
nbl   = 1
nbr   = nbl+nelsx
ntl   = nbl+nelsx*nbr
ntr   = ntl+nbr-1
**
nblt  = ntr+1
nbrt  = nblt+nelsx
ntlt  = nblt+nelsx*nbr
ntrt  = ntlt+nbr-1
**
nblpx = nbl+1
nblpy = nbl+nbr
nbrmx = nbr-1
nbrpy = nbr+nbr
ntrmx = ntr-1
ntrmy = ntr-nbr
ntlpx = ntl+1
ntlmy = ntl-nbr
**
n11 = nbl
n12 = n11+1
n13 = n12+nbr
n14 = n13-1
n15 = nblt
n16 = n15+1
n17 = n16+nbr
n18 = n17-1
**
x1 = 0.0
y1 = 0.0
x2 = 10.0
y2 = 0.0
x3 = 10.0
y3 = 10.0
x4 = 0.0
y4 = 10.0
**
*node,nset=nbl
<nbl>, <x1>,<y1>
*node,nset=nbr
<nbr>, <x2>,<y2>
*nfill,nset=nbot
nbl,nbr,<nelsx>,1
****
*node,nset=nblt
<nblt>, <x1>,<y1>, 2.d0
*node,nset=nbrt
<nbrt>, <x2>,<y2>, 2.d0
*nfill,nset=nbott
nblt,nbrt,<nelsx>,1
****
*node,nset=ntl
<ntl>, <x4>,<y4>
*node,nset=ntr
<ntr>, <x3>,<y3>
*nfill,nset=ntop
ntl,ntr,<nelsx>,1
****
*node,nset=ntlt
<ntlt>, <x4>,<y4>, 2.d0
*node,nset=ntrt
<ntrt>, <x3>,<y3>, 2.d0
*nfill,nset=ntopt
ntlt,ntrt,<nelsx>,1
*nfill,nset=nall
nbot, ntop, <nelsx>,<nbr>
nbott,ntopt,<nelsx>,<nbr>
*element,type=emc3d8
1, <n11>,<n12>,<n13>,<n14>, <n15>, <n16>, <n17>, <n18>
*elgen,elset=emc3d8
1, <nelsx>,1,1, <nelsx>,<nbr>,<nelsx>
*elset,elset=one
emc3d8
*elset,elset=bcx0,generate
1,91,10
*elset,elset=bcy0,generate
1,10,1
*elset,elset=bcx1,generate
10,100,10
*elset,elset=bcy1,generate
91,100,1
*surface,name=bcxy0
bcx0,s6
bcy0,s3
*surface,name=bcx1
bcx1,s4
*surface,name=bcy1
bcy1,s5
**
*orientation,name=oricyl,system=cylindrical
0.,0.,0., 0.,0.,1.
3,
*solid section,elset=one,material=mone
*material,name=mone
*density
1.
*dielectric
8.854187817e-12
*electrical conductivity
<conduct>
*magnetic permeability
1.2566370614e-6
**
*step
*electromagnetic,low frequency, time harmonic
60.d0
*d em potential
bcxy0,s,mvp,,  0.d0, 1.,1.,1.,0.d0, 1.,1.,1.
bcx1, s,mvp,, 10.d0, 0.,1.,0.,10.d0, 0.,10.,0.
bcy1, s,mvp,, 10.d0,-1.,0.,0.,10.d0,-1.,0.,0.
*decurrent
one,cjnu, 1.0, 1.,1.,1., oricyl
*decurrent,imaginary
one,cjnu, 1.0, 1.,1.,1., oricyl
*output,field
*element output
EME,EMB
*end step