jsv5D
;... Maxshr.xfm   
; Calculates lines for shear strain plot
; ...Define hyperbolic functions. The ones
; supplied with sigma plot do not seem to
; working correctly
ssinh(x)=(exp(x)-exp(-x))/2
ccosh(x)=(exp(x)+exp(-x))/2
;... Set range for Rxz
col(1) = data(1,8,.1)
Rxz = col(1)
Erxz=ln(col(1))
;... Calculate Maximum shear strain
col(2) = ssinh(Erxz)
;... Calculate minimum octahedral strain
col(3)=0.6666*sqrt(2)*ssinh(Erxz)
;... Calculate maximum octahedral strain
col(4)=0.6666*sqrt(ccosh(Erxz)^2+ccosh(Erxz)-2)
;... Write shear strain tics used on right axis
col(5)= data(0,4,.5)
;... Angular shear strain for left axis tics
col(6)= arctan(col(5))