#!/bin/csh -f
unlimit
limit core 0

# from swir20:
# wormpep to end, using accnum as the ID
# this means trembl/inv.dat has to be part of the main data so
# wormpep can be compared to it

 ~pmr/nuswir/nuswir -b swir20 \
                -fSW ~pubseq/swiss38/swir20/new_seq.dat -K -P \
                -fSW ~pubseq/swiss38/swir20/upd_seq.dat -K -P \
                -fSW ~pubseq/swiss38/swir20/upd_ann.dat -K -P \
                -fSW ~pubseq/swiss38/seq.dat -K -P \
                -fTR ~pubseq/trembl12/arc.dat -S \
                -fTR ~pubseq/trembl12/fun.dat -S \
                -fTR ~pubseq/trembl12/hum.dat -S \
                -fTR ~pubseq/trembl12/mam.dat -S \
                -fTR ~pubseq/trembl12/mhc.dat -S \
                -fTR ~pubseq/trembl12/org.dat -S \
                -fTR ~pubseq/trembl12/phg.dat -S \
                -fTR ~pubseq/trembl12/pln.dat -S \
                -fTR ~pubseq/trembl12/pro.dat -S \
                -fTR ~pubseq/trembl12/rod.dat -S \
                -fTR ~pubseq/trembl12/unc.dat -S \
                -fTR ~pubseq/trembl12/vrl.dat -S \
                -fTR ~pubseq/trembl12/vrt.dat -S \
                -fTR ~pubseq/trembl12/inv.dat \
                -fWP ~pubseq/wormpep18/wormpep -w -I \
	-p12 -t1000 -N350 >&! swir20.out

# combine the main database with the stored TREMBL survivors

 cat swir20.fasta swir20.store >! swir20

# later, need to compare all these to tremblnew
# because the above only compares each trembl division to SW and
# we actually want to compare tremblnew to everything.
