#!/bin/bash
data=${1-ex1.dat}		# data is assumed to have x-y-err

# fit model and print model
#
# polynomial
trend1d -W $data -Fxmr -N2r > $data.fit # will produce x-model-reduced

# fourier
#trend1d -W $data -Fxmr -Nf6r > $data.fit # will produce x-model-reduced

ofile=$data.ps

reg=`minmax $data -I1/0.1`
#reg=-R-1/13/-2/2

echo $0: working on data $data
echo $0: determined region $reg
proj=-JX7/3
# plot data and errors
psxy -Y4 $reg $proj $data -Sc0.02 -Ey -Gred -Ba5f.5:"x":/a1f.1:"y"::."data and model fit":Wesn -K  > $ofile
# and model
psxy $data.fit $reg $proj -O -W10,darkorange -K >> $ofile
#
# and below detrended model
gawk '{print($1,$3)}' $data.fit | \
    psxy -Y-4.5 $reg $proj  -Sc0.02  -Gred -Ba5f.5:"x":/a1f.1:"y"::."detrended data":WeSn -O >> $ofile





echo $0: output in $ofile
