new
rest flac3D.sav
def showfirststrain
array ai(6)
zp=zone_head
loop while zp#null
dun=z_fsi(zp,ai)
e_th=ai(1)+ai(2)
e_tf=((ai(1)-ai(2))^2+(ai(4))^2)
z_extra(zp,1)=0.5*(e_th+squr(e_tf))
z_extra(zp,2)=0.5*(e_th-squr(e_tf))
zp=z_next(zp)
endloop
end
showfirststrain