#!/bin/csh
set nonomatch

elegant chrom1.ele 

sddsxref chrom1.twi chrom1.fin -pipe=out -leave=* -transfer=parameter,MAL.DP -rename=parameter,MAL.DP=delta \
    | sddsprocess -print=column,Rootname,%s\#%ld,ElementName,ElementOccurence -pipe=in chrom1.tmp

sddsvslopes chrom1.tmp -pipe=out -columns=betax,betay,alphax,alphay,etax,etaxp -independentVariable=delta \
| sddsxref -pipe chrom1.tmp -take=s \
| sddsprocess -pipe=in chrom1.fit \
    "-define=col,bx,betaxSlope betaxIntercept /" \
    "-define=col,by,betaySlope betayIntercept /" \
    "-define=col,ax,alphaxSlope alphaxIntercept bx * -" \
    "-define=col,ay,alphaySlope alphayIntercept by * -" \
    "-define=col,Wx,ax sqr bx sqr + sqrt" \
    "-define=col,Wy,ay sqr by sqr + sqrt" 

\rm chrom1.tmp

sddsplot -column=s,W? -sep chrom1.fit -title=





  
     
