Electric field
disp v
coor x,y
shap %1 %2
gaus %3
load 0.0
coef tn
mate t0 alfa 300;0.6e8;
stif
$c6 sigma=t0*alfa/tn
dist=+[v/x;v/x]*sigma+[v/y;v/y]*sigma
end
Temperature field
disp u
coor x,y
shap %1 %2
gaus %3
coef un vn
mate almd t0 alfa 400;300;0.6e8;
stif
$c6 sigma=t0*alfa/vn
$cv vx={vn/x}
$cv vy={vn/y}
$c6 v2=vx**2+vy**2
dist=+[u/x;u/x]*almd+[u/y;u/y]*almd
load=+[u]*v2*sigma
end