Scaling with SCALA (or better: aimless): Difference between revisions

m
no edit summary
mNo edit summary
Line 2: Line 2:




An example script to convert XDS_ASCII.HKL into a [http://www.ccp4.ac.uk CCP4] multi-record MTZ file is:
An example script (contributed by Clemens Vonrhein) to convert XDS_ASCII.HKL into a [http://www.ccp4.ac.uk CCP4] multi-record MTZ file is:


  #!/bin/csh -f
  #!/bin/csh -f
Line 41: Line 41:
== Using pointless ==
== Using pointless ==


Another way of going from a XDS formatted reflection file to a multi-record MTZ file (suitable for SCALA and TRUNCATE, but not directly for structure solution) is via the [[pointless]] program
A better way of going from a XDS formatted reflection file to a multi-record MTZ file (suitable for SCALA and TRUNCATE, but not directly for structure solution which requires single-record MTZ files) is via the [[pointless]] program
(see also "[ftp://ftp.ccp4.ac.uk/ccp4/6.0.2/prerelease/runningpointlessandscala.pdf Running Pointless and Scala]"):
(see also "[ftp://ftp.ccp4.ac.uk/ccp4/6.0.2/prerelease/runningpointlessandscala.pdf Running Pointless and Scala]"). This way is better (according to PRE), but pointless has only been made available recently (and one needs the latest version of the program). It's very easy to do:


  pointless -copy xdsin XDS_ASCII.HKL hklout XDS_ASCII.mtz
  pointless -copy xdsin XDS_ASCII.HKL hklout XDS_ASCII.mtz
Line 60: Line 60:
Note: another program to work with those multi-record MTZ files is [http://www.ccp4.ac.uk/dist/html/reindex.html REINDEX].
Note: another program to work with those multi-record MTZ files is [http://www.ccp4.ac.uk/dist/html/reindex.html REINDEX].


== Using xdsconv to get the MTZ file directly (i.e. circumventing SCALA) ==


The XDSCONV way of obtaining a MTZ file would be just:
== See also ==
INPUT_FILE=XDS_ASCII.HKL
INCLUDE_RESOLUTION_RANGE=50 1  ! optional
OUTPUT_FILE=temp.hkl  CCP4    FRIEDEL'S_LAW=FALSE
(XDSCONV picks up space group and cell parameters from XDS_ASCII.HKL.) After running [[XDSCONV]] one would then run
f2mtz HKLOUT temp.mtz<F2MTZ.INP
cad HKLIN1 temp.mtz HKLOUT junk_xdsconv.mtz<<EOF
LABIN FILE 1 E1=FP E2=SIGFP E3=DANO E4=SIGDANO E5=ISYM
LABOUT FILE 1 E1=FP E2=SIGFP E3=DANO_sulf E4=SIGDANO_sulf E5=ISYM_sulf
END
EOF


ISYM column is important if you want to run SHARP afterwards.
[[XDSCONV]] to get the MTZ file directly (i.e. circumventing SCALA)
2,652

edits