27.10.2014 Views

title of the thesis - Department of Geology - Queen's University

title of the thesis - Department of Geology - Queen's University

title of the thesis - Department of Geology - Queen's University

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

yval= abs(z_syy(zi))<br />

xyval= abs(z_sxy(zi))<br />

subt=abs(xval-yval)<br />

s1 = 0.5*(xval+yval) + sqrt(subt*subt+4*xyval*xyval)<br />

s3 = 0.5*(xval+yval) - sqrt(subt*subt+4*xyval*xyval)<br />

ratio = s3/s1<br />

z_extra(zi)=ratio<br />

zi=z_next(zi)<br />

endloop<br />

bi=b_next(bi)<br />

endloop<br />

end<br />

ratios3s1<br />

plot hold z_extra fill block<br />

E.9.1 FISH Routine: ratio s3s1.fis<br />

def ratios3s1<br />

bi=block_head<br />

loop while bi#0<br />

zi=b_zone(bi)<br />

loop while zi#0<br />

xval= z_sxx(zi)<br />

yval= z_syy(zi)<br />

xyval= z_sxy(zi)<br />

s1 = 0.5*(xval+yval) + sqrt((xval-yval)*(xval-yval)+4*xyval*xyval)<br />

s3 = 0.5*(xval+yval) - sqrt((xval-yval)*(xval-yval)+4*xyval*xyval)<br />

ratio = s3/s1<br />

z_extra(zi)=ratio<br />

zi=z_next(zi)<br />

endloop<br />

bi=b_next(bi)<br />

endloop<br />

end<br />

ratios3s1<br />

plot fill z_extra block<br />

205

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!