#!/bin/csh

set mosflm_output_file="../Mtz/fps_2_1_143.mtz"
set resolution="20.0 2.0"

##################################################
#SORTMTZ
sortmtz hklout ../Mtz/L1_sort.mtz<< EOF_sortmtz  > L1_sortmtz.log
H K L M/ISYM BATCH I SIGI FRACTIONCALC RATDELSD IPR SIGIPR
${mosflm_output_file}
EOF_sortmtz


#SCALA
scala \
HKLIN ../Mtz/L1_sort.mtz\
HKLOUT ../Mtz/L1_scala.mtz\
<< eof-scala >& L1_scala.log
!RUN 1 BATCH 1 TO 143
!RESOLUTION 20.00 2.00
RESOLUTION ${resolution}
CYCLES 10
SDCORRECTION  1.25  0.03
EXCLUDE SDMIN 3.00  
REJECT SCALE 6.00 6.00 REJECT
REJECT MERGE 6.00 6.00 KEEP
SCALES BATCH BROTATION SPACING 5.00 
!RESTORE
!INITIAL
INTENSITIES PROFILE PARTIALS 
ANOMALOUS ON
FINAL PARTIALS
OUTPUT AVERAGE
eof-scala

truncate:
truncate hklin ../Mtz/L1_scala\
         hklout ../Mtz/L1_trunc << eof_t >>! L1_truncate.log
truncate yes
ranges 10
labout F=FP SIGF=SIGFP
eof_t