# The four stage scripts of the line 31-81 flow, concatenated (Claude Opus 5.5, 2026-10-17).

######## flow1.sh
#!/bin/ksh
# Stage 1: preprocessing and CDP sort for NPRA 31-81 (Claude Opus 5.5)
export PATH=$HOME/seismic/su/bin:/usr/local/bin:$PATH CWPROOT=$HOME/seismic/su
cd $HOME/seismic/work
t() { echo "--- $1"; shift; /usr/bin/time ksh -c "$*" 2>&1 | egrep "Real|User|System" | tr "\n" " "; echo; }
t kill_dead  "sukill key=trid a=2 min=2 max=2 < geom.su > k.su 2>/dev/null || cp geom.su k.su"
t gain       "sugain tpow=2 < k.su > g.su"
t mute       "sumute key=offset xmute=-5225,-1500,0,1500,5225 tmute=0.62,0.20,0.02,0.20,0.62 < g.su > m.su"
t decon      "supef minlag=0.004 maxlag=0.16 pnoise=0.001 mincorr=0.2 maxcorr=2.5 < m.su > d.su"
t filter     "sufilter f=8,12,60,80 < d.su > f.su"
t sort       "susort cdp offset < f.su > cdp.su"
ls -l *.su

######## flow2.sh
#!/bin/ksh
# Stage 2: NMO + stack, normal and reversed spread (Claude Opus 5.5)
export PATH=$HOME/seismic/su/bin:/usr/local/bin:$PATH CWPROOT=$HOME/seismic/su
cd $HOME/seismic/work
t() { echo "--- $1"; shift; /usr/bin/time ksh -c "$*" 2>&1 | egrep "Real|User|System" | tr "\n" " "; echo; }
t nmo        "sunmo par=nmo.par smute=1.6 < cdp.su > nmo.su"
t stack      "sustack key=cdp normpow=1 < nmo.su > stack.su"
t rev_geom   "geom3181 reverse=1 < f.su 2>/dev/null | susort cdp offset > cdp_rev.su"
t rev_stack  "sunmo par=nmo.par smute=1.6 < cdp_rev.su | sustack key=cdp normpow=1 > stack_rev.su"
ls -l stack*.su

######## flow3.sh
#!/bin/ksh
# Stage 3: trace editing + balance, restack (Claude Opus 5.5)
export PATH=$HOME/seismic/su/bin:/usr/local/bin:$PATH CWPROOT=$HOME/seismic/su
cd $HOME/seismic/work
t() { echo "--- $1"; shift; /usr/bin/time ksh -c "$*" 2>&1 | egrep "Real|User|System" | tr "\n" " "; echo; }
t edit       "suwind key=ep reject=110 < f.su | suwind key=trid min=1 max=1 > e.su"
t balance    "sugain pbal=1 < e.su > b.su"
t sort       "susort cdp offset < b.su > cdpb.su"
t nmo_stack  "sunmo par=nmo.par smute=1.6 < cdpb.su | sustack key=cdp normpow=1 > stack2.su"
ls -l stack2.su

######## flow4.sh
#!/bin/ksh
# Stage 4: post-stack migration, three methods (Claude Opus 5.5)
export PATH=$HOME/seismic/su/bin:/usr/local/bin:$PATH CWPROOT=$HOME/seismic/su
cd $HOME/seismic/work
t() { echo "--- $1"; shift; /usr/bin/time ksh -c "$*" 2>&1 | egrep "Real|User|System" | tr "\n" " "; echo; }
TM=0.0,0.45,0.75,1.05,1.35,1.65,1.95,2.25,2.55,3.2,4.0,6.0
VM=7750,7750,8050,8550,9200,9950,10350,10750,12150,13750,15000,17000
t stolt   "sustolt cdpmin=1 cdpmax=536 dxcdp=55 tmig=$TM vmig=$VM smig=1 < stack2.su > mig_stolt.su"
t phase   "sumigps dx=55 tmig=$TM vmig=$VM < stack2.su > mig_ps.su"
t kirch   "sumigtk dxcdp=55 tmig=$TM vmig=$VM < stack2.su > mig_tk.su"
ls -l mig_*.su