++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  +  SHELXC - Create input files for SHELXD and SHELXE - Version 2013/2  +
  +  Copyright (c) George M. Sheldrick 2003-13                           +
  +  I4_p9                           Started at 13:34:01 on 15 Dec 2020  +
  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

 SHELXC reads a filename stem (denoted here by 'xx') on the command line
 plus some instructions from 'standard input'. It writes some statistics to
 'standard output' and prepares the three files needed to run SHELXD and
 SHELXE. SHELXC can be called from a GUI by a command line such as:

 shelxc xx <t

 which would read the instructions from the file t, or (under most UNIX
 systems) by a simple shell script that includes the instructions, e.g.

 shelxc xx <<EOF
 CELL 49.70 57.90 74.17 90 90 90
 SPAG P212121
 SAD elastase.sca
 FIND 12
 <<EOF
 shelxd xx_fa
 shelxe xx xx_fa -s0.37 -m20 -h -b
 shelxe xx xx_fa -s0.37 -m20 -h -b -i

 which would also run shelxd to locate the sulfur atoms and shelxe (for
 both substructure enantiomers) to solve elastase by sulfur-SAD phasing.
 The reflection data may be in SHELX (.hkl), HKL2000 (.sca) or XDS
 XDS_ASCII.HKL format. Any names may be used for XDS reflection files,
 SHELXC recognises them by reading the first record.

 This script would read data from the .sca file and write the files
 xx.hkl (h,k,l,I,sig(I) in SHELX HKLF4 format for density modification
 by SHELXE or refinement with SHELXL) xx_fa.ins (cell, symmetry etc. for
 heavy atoms location  by SHELXD) and xx_fa.hkl (h,k,l,FA,sig(FA),alpha
 for both SHELXD and SHELXE). The starting phases for density
 modification are estimated as (alpha + heavy-atom phase in the simple
 approach used by SHELXE.

 For SIR or SIRAS, two input reflections files are specified by the
 keywords NAT and SIR or SIRA; for MAD at least two of the reflection files
 HREM, LREM, PEAK and INFL are required and NAT may also be given if higher
 resolution native data are available (e.g. SMet for SeMet MAD). Reflection
 data should be in SHELX .hkl or SCALEPACK .sca format; many other programs,
 including SCALA and XPREP, can output .sca format too. The keywords CELL,
 SPAG (space group) SPAG (space group) and FIND (number of heavy atoms) are
 always required, SFAC, MIND, NTRY, SHEL, ESEL and DSUL may be given and are
 written to the file xx_fa.ins for SHELXD. MAXM can be used to reserve
 memory in units of 1M reflections. For RIP phasing, NAT (or BEFORE) denotes
 the file before radiation damage and RIP (or AFTER) after radiation
 damage. For RIPAS the 'after' file must be called 'RIPA' and a keyword RIPW
 (default 0.6) gives the weight w to be assigned to the 'NAT' data in the
 estimation of the anomalous signal (a weight of 1-w is applied to the 'RIPA'
 data). DSCA (default 0.98) gives the factor to multiply the native data
 for SIR and SIRAS or the 'after' data for RIP after the data have been
 put on the same scale (this allows for the extra scattering power of the
 heavy atoms etc.); this can be critical for RIP phasing.

 ASCA (default 1.0) is a scale factor applied to the anomalous signal in a
 MAD experiment; to apply MAD to a small molecule, ASCA and DSCA should both
 be between 0 and 1, the best values have to be found by trial and error.
 Finally SMAD (without a number) sets the dispersive term to zero in a MAD
 experiment, equivalent to SAD using weighted mean anomalous differences
 from all the MAD datasets. This should always be tried whenever radiation
 damage is suspected.

 SHELXC also tests for and if necessary corrects the more common cases of
 inconsistent indexing when more than one dataset is involved. In addition,
 the mean value of |E^2-1| is calculated for each dataset to detect twinning.


 XC reads a filename stem (denoted here by 'xx') on the command line
 plus some instructions from 'standard input'. It writes some statistics to
 'standard output' and prepares the three files needed to run XM and
 XE. XC can be called from a GUI by a command line such as:

 xc xx <t

 which would read the instructions from the file t, or (under most UNIX
 systems) by a simple shell script that includes the instructions, e.g.

 xc xx <<EOF
 CELL 49.70 57.90 74.17 90 90 90
 SPAG P212121
 SAD elastase.sca
 FIND 12
 <<EOF
 xm xx_fa
 xe xx xx_fa -s0.37 -m20 -h -b
 xe xx xx_fa -s0.37 -m20 -h -b -i

 which would also run shelxd to locate the sulfur atoms and shelxe (for
 both substructure enantiomers) to solve elastase by sulfur-SAD phasing.
 The reflection data may be in SHELX (.hkl), HKL2000 (.sca) or XDS
 XDS_ASCII.HKL format. Any names may be used for XDS reflection files,
 XC recognises them by reading the first record.

 This script would read data from the .sca file and write the files
 xx.hkl (h,k,l,I,sig(I) in SHELX HKLF4 format for density modification
 by XE or refinement with XL) xx_fa.ins (cell, symmetry etc. for
 heavy atoms location  by XM) and xx_fa.hkl (h,k,l,FA,sig(FA),alpha
 for both XM and XE). The starting phases for density
 modification are estimated as (alpha + heavy-atom phase in the simple
 approach used by XE.

 For SIR or SIRAS, two input reflections files are specified by the
 keywords NAT and SIR or SIRA; for MAD at least two of the reflection files
 HREM, LREM, PEAK and INFL are required and NAT may also be given if higher
 resolution native data are available (e.g. SMet for SeMet MAD). Reflection
 data should be in SHELX .hkl or SCALEPACK .sca format; many other programs,
 including SCALA and XPREP, can output .sca format too. The keywords CELL,
 SPAG (space group) SPAG (space group) and FIND (number of heavy atoms) are
 always required, SFAC, MIND, NTRY, SHEL, ESEL and DSUL may be given and are
 written to the file xx_fa.ins for XM. MAXM can be used to reserve
 memory in units of 1M reflections. For RIP phasing, NAT (or BEFORE) denotes
 the file before radiation damage and RIP (or AFTER) after radiation
 damage. For RIPAS the 'after' file must be called 'RIPA' and a keyword RIPW
 (default 0.6) gives the weight w to be assigned to the 'NAT' data in the
 estimation of the anomalous signal (a weight of 1-w is applied to the 'RIPA'
 data). DSCA (default 0.98) gives the factor to multiply the native data
 for SIR and SIRAS or the 'after' data for RIP after the data have been
 put on the same scale (this allows for the extra scattering power of the
 heavy atoms etc.); this can be critical for RIP phasing.

 ASCA (default 1.0) is a scale factor applied to the anomalous signal in a
 MAD experiment; to apply MAD to a small molecule, ASCA and DSCA should both
 be between 0 and 1, the best values have to be found by trial and error.
 Finally SMAD (without a number) sets the dispersive term to zero in a MAD
 experiment, equivalent to SAD using weighted mean anomalous differences
 from all the MAD datasets. This should always be tried whenever radiation
 damage is suspected.

 XC also tests for and if necessary corrects the more common cases of
 inconsistent indexing when more than one dataset is involved. In addition,
 the mean value of |E^2-1| is calculated for each dataset to detect twinning.

 ==============================================================================

   145453 Reflections read from SAD file sad-p9_se_w2.sca

    21557 Unique reflections, highest resolution  1.744 Angstroms
    213.3 Friedel pairs used on average for local scaling

 Resl.   Inf.  3.87  3.07  2.68  2.44  2.26  2.13  2.02  1.93  1.86  1.80  1.74
 N(data)    2066  2024  1988  1948  2047  1911  2023  2030  1878  1805  1837
 Chi-sq     0.36  0.29  0.25  0.41  0.50  0.50  0.54  0.40  0.73  0.75  0.86
 <I/sig>    24.1  23.8  21.8  19.6  18.0  15.4  12.4  10.0   8.0   6.6   5.6
 %Complete  99.6 100.0 100.0 100.0 100.0 100.0 100.0 100.0  99.9  99.9  96.7
 <d"/sig>   2.00  1.75  1.90  1.81  1.69  1.42  1.25  1.09  1.05  1.03  0.95
 CC(1/2)    88.8  80.3  78.9  74.7  70.7  63.7  53.7  49.2  44.7  40.0  30.3

 For zero signal <d'/sig> and <d"/sig> should be about 0.80

    16206 Reflections written to file I4_p9_fa.hkl for input to SHELXD/E

    21557 Reflections written to file I4_p9.hkl for input to SHELXE

 Dataset  <|E^2-1|>  Max.resln.(A)
 SAD        0.767      1.744

 If <|E^2-1|> is less than about 0.68, it is likely that the structure
 is twinned or that F has been read in as intensity (specify -f!)

 File I4_p9_fa.ins written for input to SHELXD:

 TITL I4_p9_fa.ins SAD in I4
 CELL  0.98000  113.95  113.95   32.47   90.00   90.00   90.00
 LATT  -2
 SYMM -Y, X, Z
 SYMM -X, -Y, Z
 SYMM Y, -X, Z
 SFAC Se
 UNIT   160
 SHEL 999 2.2
 PATS
 FIND 5
 MIND -1.5 -0.1
 NTRY 1000
 SEED 1
 HKLF 3
 END 

  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  +  SHELXC for I4_p9 finished at 13:34:01 on 15 Dec 2020                +
  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++