Pentru Vlad

Adauga la comentarii

Convert tone-GPa

#!/usr/bin/python
#transform psi to tones and MPa
from math import *
r = 5.0
tones = input('Enter tones value: ')
MPa = (tones * 9.8 * 10**3) / (pi * r**2)
print 'MPa = ', MPa

Reciprocal space for TEM

#!/usr/bin/python
#Angle between two planes in reciprocal space
from math import *
from MatPy.interact import *
deg=pi/180
h1 = input('h1 = ')
k1 = input('k1 = ')
l1 = input('l1 = ')
h2 = input('h2 = ')
k2 = input('k2 = ')
l2 = input('l2 = ')
print('Crystallograpgic system:\n')
print('1 - Cubic\n')
print('2 - Tetragonal\n')
print('3 - Ortorombic\n')
print('4 - Hexagonal\n')
sys = input('1, 2, 3 or 4 => ')
if sys == 1:
a = input('a = ')
g_rec=[1/a**2 0 0 0 1/a**2 0 0 0 1/a**2]
elseif sys == 2 
a = input('a = ')
c = input('c = ')
g_rec=[1/a**2 0 0 0 1/a**2 0 0 0 1/c**2]
elseif sys == 3
a = input('a = ')
b = input('b = ')
c = input('c = ')
g_rec=[1/a**2 0 0 0 1/b**2 0 0 0 1/c**2]
elseif sys == 4
a = input('a = ')
c = input('c = ')
g_rec=[4/(3*a**2) 2/(3*a**2) 0 2/(3*a**2) 4/(3*a**2) 0 0 0 1/c**2]
end

p1=[h1 k1 l1]
p2=[h2 k2 l2]
cos_teta=(p1*g_rec*p2'/(sqrt(p1*g_rec*p1')*sqrt(p2*g_rec*p2')))
teta=rad2deg(acos(cos_teta))
zone = cross(p1,p2)

Lumea blogului

Oare cum si-ar inchipui un pescarus lumea lui pierduta in labirinturile blogului. O fi ca intr-o cusca, o fi ca un zbor sub cerul nemarginit?