#!/usr/bin/env tclsh

# Copyright (C) Navraj S. Pannu 2007, Leiden University
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
#
#

set anomalous       [lindex $argv 0]
set mtz2scacommand  [lindex $argv 1]
set inputmtz        [lindex $argv 2]
set output          [lindex $argv 3]
set ip              [lindex $argv 4]
set sigip           [lindex $argv 5]

if { $anomalous == "0" } { 

    set mtz_script "XTAL xtal1\n DNAME data1 \n COLUmn I=$ip SI=$sigip\nOUTS\nOUTPut $output\nEND\n"

    set command "$mtz2scacommand HKLIN $inputmtz << \"$mtz_script\""
    catch {eval exec $command } output 
    puts $output

} else {
    set im        [lindex $argv 6]
    set sigim     [lindex $argv 7]

    set mtz_script "XTAL xtal1\n DNAME data1 \n COLUmn I+=$ip SI+=$sigip I-=$im SI-=$sigim\nOUTS\nOUTPut $output\nEND\n"

    set command "$mtz2scacommand HKLIN $inputmtz << \"$mtz_script\""
    catch {eval exec $command } output 
    puts $output
}