1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33
|
#include <clipper/clipper.h>
#include <clipper/clipper-ccp4.h>
//#include <clipper/clipper-ccp4.h>
#include <iostream>
using namespace clipper;
using namespace clipper::data32;
int main(int argc, char** argv)
{
CCP4MTZfile file;
HKL_info mydata;
HKL_data<F_sigF> myfsig(mydata);
file.open_read(argv[1]);
file.import_hkl_info( mydata, false );
file.import_hkl_data( myfsig, argv[2] );
Spacegroup sp = file.spacegroup();
file.close_read();
file.open_append(argv[1],"junk1.mtz");
file.export_hkl_data( myfsig, "*/*/anom");
file.close_append();
file.open_write("junk2.mtz");
file.export_hkl_info( mydata );
file.export_hkl_data( myfsig, "*/*/anom");
file.close_write();
sp.debug();
}
|