## GARANT MACROS: appPeaks - appends expected peaks of desired spectrum
##
## USAGE: appPeaks <spectrum> <read.seq> <read.prot> <append.peaks>
##
##  Peaks expected are generated for the spectrum of type <spectrum> and the
## sequence <read.seq>. The frequencies are taken from the atom list <read.prot>
## and from the assigned peaks in <append.peaks>. All expected peaks that are
## not present in <append.peaks> are appended to this peak list. The color of
## the appended peaks is set to 5. Inconsistent assignments are reported on
## standard output.
##
## EXAMPLE: appPeaks HCCH24 tendmutDIANA.seq tendmut.prot hcch.peaks
##
## SEE ALSO: spectra types

parameter spec seq prot peaks

# if defined call macro $p1 else issue command spectrum $p1
if (macro('$spec')) then
   $spec			# call macro with default setup
else
   spectrum $spec
end if

load sequence $seq 		# sequence

make coherences			# generate all possible coherences
make couplings			# generate couplings
make noe			# generate noe's
make peaks			# generate expected peaks

peaklist $spec $peaks $prot 0.0
check inconsistencies

load frequencies $prot
append expected peaks $spec 0 5
write peaks $spec $peaks $prot

quit