lamp_mac/ATOMS/
gr_angle_r.pro
top gr_angle_r
result = gr_angle_r(w_in, angle_list, rcut_0, rcut_2)
******* ********** ** User callable ** The call is e.g. w1=gr_angle_r(w2,[angle_list],rcut_1,rcut_2) ** ** w2 is a workspace of fractional coords (i.e. between 0 and 1) read by xdatcar.pro ** side is side of box ** n_atomtypes is number of different atoms, ** atoms_list is a vector with n_atomtypes entries ** all variables read from CONTCAR ** angle_list is a vector length 3 of the 3 types of atoms defining the angle ** if the atoms are ordered O, Na, P, then the angle OPO would be [0,2,0] ** rcut_0 defines the radius between atoms 1 and 2 ** rcut_2 defines the radius between atoms 2 and 3 ** to limit the atoms for which angles are calculated ** rcut_0 and rcut_2 probably correspond to maxima in the respective g(r) ** ** method ********** ** loop thro frames and atoms ** choose a target atom and calculate histogram for each type of pair ** number of partials to calculate is n_atomtypes(n_atomtypes+1)/2 ** takes each frame and centre it ** then sort the atoms when ...
Parameters
- w_in
- angle_list
- rcut_0
- rcut_2
File attributes
Modification date: | Tue Oct 1 14:45:53 2013 |
Lines: | 147 |