jsv5D
;..OctNadai.xfm   4/25/94
;... function calculates contours of conventional
; octahedral shear strain for Nadai strain plot
;... contour values should be placed in col(25)
n=size(col(25))  ;number of contours
;... calculate Ed for each contour value
eta=2.25*col(25)^2 +1
col(26)=0.40825* ln(eta+sqrt(eta^2-1))
theta1=data(-30,30,2)
theta2=data(-90,90,2)
;... calculate Syprime, Rxy pairs for contours
for i=1 to n do
  ;... location of next pair of columns
  cell(27,1)=27+2*i-1
  ;... calculate Sy prime
  col(cell(27,1))=exp(0.8165*cell(26,i)*sin(theta1))
  ;... calculate Rxz
  col(cell(27,1)+1)=exp(1.4142*cell(26,i)*cos(theta1))
  ;... location of next pair of columns
  cell(27,1)=27+2*n+(2*i-1)
  ;... calculate Sy prime
  col(cell(27,1))=exp(0.8165*cell(26,i)*sin(theta2))
  ;... calculate Rxz
  col(cell(27,1)+1)=exp(1.4142*cell(26,i)*cos(theta2))
end for