PROGRAM ARRSNP

To group together equivalent measurements and give some statistics
on their degree of equivalence, indices may be submultiples of integers.

Description:

Arranges observations of elements of polarisation matrices according to the matrix indices, the reflectionindices, and the symmetry given by the crystal data. Averages of repeated measurements of the same item, and of groups of equivalent items are made.

Input:

The crystal data file must contain:

Symmetry information on S cards
I cards which select the different options. The words on I cards which are recognised are: DTYP,MULH,MULK,MULL,PRIN,REJ,RNUM,WGHT
DTYP is mandatory and the number which follows indicates the format of the data only 11 and 12 for SNP data

DTYP = 11 IGOR type data (rhistory crpd outout) ext .pol
DTYP = 12 MAD data from logfiles (rlogfile crpd) ext .pol

All the other words are optional:
Any or all of MULH, MULK and MULL each followed by the number by which

the corresponding index h,k,or \$\ell\$ should be multiplied to obtain an integer.

REJ if given and followed by the value 1 a list of reflections to

be rejected is expected.

RNUM if given and followed by the value 1 the record numbers will

be stored with the reflections in the output file as well as in the listing file.

WGHT The number which follows selects the weighting scheme

= 0 (default) Unit weights for individual measurements, and

the mean of equivalents weighted according to their esd's.

= 1 Statisitical weights used for individual measurements, and

esd's used for weighting for mean of equivalents.

= 2 Unit weights used throughout.

A data file of the type indicated by DTYP above. Its default extension is .pol
If REJ=1 a file containing the record numbers of reflections to be rejected. The numbers must be in increasing order except that a negative number means reject all reflections with record numbers between the previous positive number given and that whose record number is the absolute value just read (inclusive). The default extension for this file-name is .rej

Output:

A listing file: ARRSNP.lis which reports how the crystal data file was interpreted and lists the data which have been arranged. In this output all reflections equivalent to one another are grouped together and averaged.
Within the groups multiple measurements of the same reflection are placed together. These are arranged with h changing most rapidly, k next and l most slowly.
The esd's of the means are calculated both from the esd's of the values being averaged, and from the weighted deviations of these values from the mean. The larger of the two values is printed and used for further calculation.
Some statistics on the degree of equivalence between supposedly equivalent measurements are also given in the listing file and on the terminal.
A file with extension .arr destined to be read by further processing programs contains all the data arranged as above but with no averaging. These data can subsequently be read by function INPUAR or equivalent code. One call to INPUAR fills in the COMMON /REFS/ with all the data for one group of equivalent reflections. (COMMON /REFS/K(3,2),LL(48,2),R(500,2),SCALE(2),INC,II,FF(3,2) )
The final index IS in REFS is to enable data to be read from more than one source, for merging etc. LL(I,IS) contains the number of entries in R for the Ith reflection in the group. The entries in R for each reflection are:

II+1 to II+3 : its indices hkl
II+4 to II+LL: groups of "value" "esd" ("other data")

or groups of "record number" "value" "esd" ("other data") if RNUM 1 was given.

II starts at 1 for the first reflection in the group and should be incremented by LL(I,IS)+2 after obtaining each subgroup I of data (the extra two slots provide space in which to store the means of repeated measurements). If LL(J+1,IS)=3 the Jth reflection is the last in the group and the next three entries in R are the indices hkl with which to label the group of equivalents. See subroutine STATIS for an example of using the data stored in REFS.

Notes:

AVERPA reads and processes .arr files containing SNP data

Running the program:

On running the program the user is asked for the name of the crystal data file, then if REJ=1, for the file containing the rejection list.
The names of one or more data files of the type specified by DTYP are then given interactively. After the last file has been read the program reports the number of records sorted, and asks for a name for the output file of sorted data.

Calls:

AINOUT ARPRIN ASUNIT CENTRE ERRMES GMREV GMZER IICD3 INDFIX INPUTI JFIX MESS NEWLIN NOPFIL OPSYM PLREAD PREFIN REJECT ROTSYM SYMOP SYMUNI WTMEAN

Common blocks used:

/ARRDAT/ to use INC NREF IBUF NUMVAL ICHNGE

/IOUNIT/ to use LPT ITO LUNI

/NSYM/ to use NOPC

/SCRACH/ to use all members

*** ARRSNP new by PJB C4.2 October 2007 ***

Data Collection and Reduction . . . . . . . Main Program

---

---

---

P. Jane Brown

Institut Laue Langevin,
Grenoble, FRANCE