HKL2000で処理できなかったフレーム数、反射数の多いデータをSCALAで処理するために
エラー
HKL2000のscalepack16m/scalepackvirus/scalepackribo/scalepackmf6500などをもってしても10000マイのデータは処理できません。 かといって、MOSFLM->SCALAやってみると
Scala: *** Too many batches, increase MAXBAT & recompile ***
というエラーが出ます・・・。
あぁ、こんなにデータとるんじゃなかった、などと言う前に・・・。
注意：以下の話は、CCP4のソースファイルをダウンロードしてビルドしている場合にのみ有効です。
scala_/parameter.fh
% cd $CCP4/src/ % vi scala_/parameter.fh
以下のパラメータを変更。
MAXBAT
バッチ番号（フレーム数の最大）
maxbat=5000 -> maxbat=10000
MAXRFL
最大反射数というのは分かるが、何の最大反射数かは不明。多分、等価反射やと思う
maxrfl = 400 -> maxrfl=4000
lib/src/ccplib.f
% cd $CCP4/lib/src/ % vi ccplib.f
として以下のパラメータ。
KPRI
テーブルサイズ・・・、って何？まぁ、直しとこう。 みんなのひろば、また大御所先生からのご指摘でこのKPRIは素数であるべき、とのこと。好みの数字に近い素数を選択して入れてください
KPRI=5003 -> KPRI=10007
3箇所あるよ。
再コンパイル
ccplibを使っているプログラムも再構築したいので（特にMTZファイル操作関連）
% cd $CCP4/ % make clean % make
とした。
SCALA 単体を再コンパイル
上記作業でいけると思いますが、SCALAのコンパイルのみやり直す時は
% cd $CCP4/src/ % make scala
でいけます。
