#!/bin/csh -f

set specName    = ft/test%03d.ft2
set simName     = sim/test%03d.ft2
set simDir      = sim
set difName     = dif/test%03d.ft2
set difDir      = dif
set inTabName   = peg_split28_pb
set auxTabName  = axt.tab
set outTabName  = nlin.tab
set errTabName  = err.nlin.tab
set noiseRMS    = 556814.187500
set specCount   = 2

set aRegSizeX  = 13
set dRegSizeX  = 13
set maxDX      = 0.0
set minXW      = 0.0
set maxXW      = 23.20
set simSizeX   = 1177

set aRegSizeY  = 18
set dRegSizeY  = 18
set maxDY      = 0.0
set minYW      = 0.0
set maxYW      = 34.67
set simSizeY   = 2048

set aRegSizeZ  = 2
set dRegSizeZ  = 2
set maxDZ      = 1.0
set minZW      = 0.0
set maxZW      = 4.00
set simSizeZ   = 2


seriesTab -in  $inTabName -list nlin.spec.list \
 -dx  $dRegSizeX -dy  $dRegSizeY \
 -adx $aRegSizeX -ady $aRegSizeY \
 -xzf 64 -yzf 72 -max \
 -out $auxTabName -adVar VOL -verb


cp $auxTabName $outTabName


if (!(-e $simDir)) then
   mkdir $simDir
endif

xyz2pipe -in $specName -verb \
| nmrPipe -fn SET -r 0.0 \
| pipe2xyz -to 0 -out $simName -ov

simSpecND -in $outTabName -list sim.spec.list \
          -mod   GAUSS1D  GAUSS1D  SCALE1D  \
          -w     $simSizeX  $simSizeY  $simSizeZ  \
          -apod None -verb


if (!(-e $difDir)) then
   mkdir $difDir
endif

addNMR -in1 $specName -in2 $simName -out $difName -sub