#!/bin/sh  
# \
exec oagtclsh "$0" "$@"

set auto_path [linsert $auto_path 0  /usr/local/oag/apps/lib/$env(HOST_ARCH)]
set auto_path [linsert $auto_path 0 /usr/local/oag/lib_patch/$env(HOST_ARCH)]
APSStandardSetup

# bucnh length
set lengthList "35 40 50 60 70 80 90 100"
set fileList ""
foreach length $lengthList {
    set lengthf [format %03.0lf $length]
    exec resistiveWallPower Al.power.${lengthf}ps -fullHeight=6.3e-3 \
        -explicitsum \
        -resistivityFile=Al.resistivity.Powell_1959 \
        -current=200 -bunches=48 -sigmat=$length \
        -length=4.45 -temperatureRange=beg=20.0,end=80.0,number=13
    lappend fileList Al.power.${lengthf}ps
}

eval exec sddscombine $fileList Al.power.variousSigmaT -overWrite

exec sddssort Al.power.variousSigmaT -pipe=out \
    -col=Temperature,decr \
    | sddsregroup -pipe=in Al.power.variousTemp \
    -newpara=Temperature -newcol=SigmaT

catch {eval exec rm [glob Al.power.variousSigmaT?* Al.power.variousTemp?*]}

exec sddsplot3 Al.power.variousSigmaT -col=Temperature,* Al.power.variousSigmaT \
    -sep=nameindex -split=pages -group=nameindex -grap=line,vary -order=spectral \
    "-leg=para=SigmaT,format=%.0f ps" &

exec sddsplot3 Al.power.variousTemp -col=SigmaT,* Al.power.variousTemp \
    -sep=nameindex -split=pages -group=nameindex -grap=line,vary  -order=spectral \
    "-leg=para=Temperature,format=%.1f K" &

exit
