++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  +  SHELXC - Create input files for SHELXD and SHELXE - Version 2013/2  +
  +  Copyright (c) George M. Sheldrick 2003-13                           +
  +  P41212_fibro                    Started at 18:30:57 on 19 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.

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

    10868 Reflections read from NAT file nat-fibro-nat.sca

    10868 Unique reflections, highest resolution  1.605 Angstroms

 Resl.   Inf.  3.58  2.84  2.48  2.26  2.09  1.97  1.87  1.79  1.72  1.66  1.61
 N(data)    1074  1026  1004   960  1029   942   999   935   968   981   950
 <I/sig>    53.3  53.3  48.3  41.7  37.2  28.4  21.5  16.5  12.7   6.5   4.6
 %Complete  94.3  99.1  99.5  99.7  99.6  99.5  99.8  99.6  99.5  99.3  92.1

    10868 Reflections written to file P41212_fibro.hkl for input to SHELXE

   221819 Reflections read from SAD file sad-fibro-sad.sca

     5667 Unique reflections, highest resolution  1.997 Angstroms
    194.2 Friedel pairs used on average for local scaling

 Resl.   Inf.  4.45  3.53  3.08  2.80  2.60  2.45  2.33  2.22  2.14  2.06  2.00
 N(data)     597   542   532   505   512   492   487   540   466   529   465
 Chi-sq     0.89  1.60  3.75  2.33  2.27  1.93  2.21  1.95  1.69  6.60  1.27
 <I/sig>   155.2 145.9 118.1  98.0  79.5  63.7  49.5  36.7  28.8  19.1  11.8
 %Complete  95.7  96.6  96.2  96.7  97.5  97.0  97.0  96.6  96.1  95.7  92.8
 <d"/sig>   5.03  2.80  1.95  2.10  1.80  1.63  1.43  1.21  1.08  1.03  1.06
 CC(1/2)    96.2  81.8  63.8  52.7  53.9  38.6  36.3  17.8  10.8  11.5  11.7

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

     4044 Reflections written to file P41212_fibro_fa.hkl for input to SHELXD/E

 Dataset  <|E^2-1|>  Max.resln.(A)
 NAT        0.717      1.605
 SAD        0.735      1.997

 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 P41212_fibro_fa.ins written for input to SHELXD:

 TITL P41212_fibro_fa.ins SAD in P41212
 CELL  1.54000   37.86   37.86  108.42   90.00   90.00   90.00
 LATT  -1
 SYMM 1/2-Y, 1/2+X, 1/4+Z
 SYMM -X, -Y, 1/2+Z
 SYMM 1/2+Y, 1/2-X, 3/4+Z
 SYMM 1/2-X, 1/2+Y, 1/4-Z
 SYMM Y, X, -Z
 SYMM 1/2+X, 1/2-Y, 3/4-Z
 SYMM -Y, -X, 1/2-Z
 SFAC S
 UNIT   320
 SHEL 999 2.5
 PATS
 FIND 10
 MIND -3.5
 DSUL 4
 NTRY 10000
 SEED 1
 HKLF 3
 END 

  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  +  SHELXC for P41212_fibro finished at 18:30:58 on 19 Dec 2020         +
  ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++