#!/bin/csh -f

#edit.bos.pl out1.bos/* > bos1.total
#edit edit.bos.pl

 foreach model ( bos1 bos2 )

    sort -k 3 -n $model.total > $model.total.sigma1
    echo 'sigma1 .15 .50 .85'
    head -n  683 $model.total.sigma1 | tail -1
    head -n 2276 $model.total.sigma1 | tail -1
    head -n 3868 $model.total.sigma1 | tail -1
    echo ' '

    sort -k 4 -n $model.total > $model.total.sigma2
    echo 'sigma2 .15 .50 .85'
    head -n  683 $model.total.sigma2 | tail -1
    head -n 2276 $model.total.sigma2 | tail -1
    head -n 3868 $model.total.sigma2 | tail -1
    echo ' '

    sort -k 5 -n $model.total > $model.total.sig
    echo 'kappa  .15 .50 .85'
    head -n  683 $model.total.sig   | tail -1
    head -n 2276 $model.total.sig   | tail -1
    head -n 3868 $model.total.sig   | tail -1
    echo ' '

    sort -k 6 -n $model.total > $model.total.kappa
    echo 'kappa  .15 .50 .85'
    head -n  683 $model.total.kappa | tail -1
    head -n 2276 $model.total.kappa | tail -1
    head -n 3868 $model.total.kappa | tail -1
    echo ' '

 end