#!/bin/sh

cd $CCP4_SCR

cmtzsplit -stdin << eof
mtzin $CEXAM/rnase/rnase25F+F-.mtz
mtzout test1.mtz
 colin FNAT,SIGFNAT
 colout FP,SIGFP
mtzout test2.mtz
 colin FreeR_flag
 colout FreeR_flag
mtzout test3.mtz
 colin FPTNCD25(+),SIGFPTNCD25(+),FPTNCD25(-),SIGFPTNCD25(-)
 colout Fpl,SigFpl,Fmi,SigFmi
eof

[ $? != 0 ] && exit 1

cmtzjoin -stdin << eof
mtzout test4.mtz
mtzin test1.mtz
 colin FP,SIGFP
 colout FNAT,SIGFNAT
mtzin test2.mtz
 colin FreeR_flag
 colout FreeR_flag
mtzin test3.mtz
 colin Fpl,SigFpl,Fmi,SigFmi
 colout FPTNCD25(+),SIGFPTNCD25(+),FPTNCD25(-),SIGFPTNCD25(-)
eof