#!/bin/csh -f

scalepack << eof-scale

resolution 200.0 2.50            [Resolution limits for all data to be scaled]
                                 [optional]

number of zones 10               [Number of resolution shells                ]
estimated error 0.02 0.03 0.04   [Estimated error for each resolution shell  ]
 0.05 0.06 0.08 0.11 0.14 0.17 
 0.22 

error scale factor 1.3           [Multiplicative factor applied to input     ]
                                 [sigmas                                     ]
[scale restrain 0.02]            [this command restrains scale factors       ]
                                 [between consecutive batches by appling     ]
                                 [expectation that the scale factors deviate ]
                                 [by 2% or so between consecutive batches.   ]

[B restrain 0.5]                 [this command restrains B factors           ]
                                 [between consecutive batches by appling     ]
                                 [expectation that the B factors deviate     ]
                                 [by 0.5 A**2 (angstrom square) or so between]
                                 [consecutive batches.                       ]
space group p4
output file gf5.sca              [output file containing h,k,l,I,sigmaI  in  ]
                                 [the case when anomalous flag is absent, or ]
                                 [h,k,l,I+,sigmaI+,I-,sigmaI-, with anomalous]
                                 [flag on                                    ]
[scale anomalous]
anomalous                        [anomalous flag on                          ]

reference batch 1                [reference batch in scaling procedure. In   ]
                                 [older versions use 'reference film' command]
                                 [instead                                    ]

add partials  1 to 60            [only partials between 1-60 will be used ]

ignore overloads                 [all reflection with saturated, (overloaded)]
                                 [ pixels, will be ignored                   ]
postrefine 10                    [Number of cycles for postrefinement, if    ]
                                 [you set 'postrefine 0', program will skip  ]
                                 [postrefinement procedure                   ]
fit crystal a* 1 to 60           [cell will be refined on the basis of the   ]
                                 [data present in the batches 1-60 only.     ]
                                 [The common value of a* will be fitted for  ]
                                 [all batches in the range 1-60              ]
[fit crystal b* 1 to 60]
fit crystal c* 1 to 60
[fit crystal alpha 1 to 60]
[fit crystal beta  1 to 60]
[fit crystal gamma 1 to 60]

fit batch rotx 1 to 60           [crystal orientation of a spindle axis angle]
                                 [will be refined for each batch separately. ]
                                 [rotz can not be fitted in postrefinement   ]
                                 [for each batch separately                  ]
fit batch roty 1 to 60
fit crystal mosaicity 1 to 60    [one mosaicity value will be fitted for all ]
                                 [batches                                    ]

rejection probability 1.e-4      [expected fraction of outliers in the data  ]


write rejection file             [outliers will be written into 'reject' file] 
@reject
format denzo_ip                  [format of the input intensity data, default]
                                 [output format from denzo. In general, you  ]
                                 [can have many different formats in one     ]
                                 [scale job, useful when data were collected ]
                                 [on different detectors                     ]
sector 1 to 60
FILE 1 './gf5###.x'              
                                 [batch numbers will start at 1 for file     ]
                                 [lyso001.x and will increase to 60 for file ]
                                 [lyso060.x. All *.x files are in current    ]
                                 [directory                                  ]

eof-scale