The LAMP Application

Generated by IDLdoc

single page | use frames     summary     class     fields     routine details     file attributes

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