# MAS CSA spectrum of a single nuclei
# Uses the csa specific method

# 1.8 sec on a 166 MHz Pentium II w/Linux

proc main {} {
    global par

    set sfrq 400
    set iso 5
    set aniso 13

    csainit 128 128 500 rep168
    csaspec [expr $aniso*$sfrq] 0.25

    set lst {}
    foreach i [csareturnint [expr $iso*$sfrq]  31.88] {
        lappend lst [concat [lrange $i 1 2] 100 0]
    }
    set f [fcreate -np 8192 -sw 20000 -type spe]
    faddpeaks $f 0 $lst
    fsave $f $par(name).spe
    funload $f
}