#!/bin/bash

# Extend a file having (x, By) to have uniform field as a function of y
# This doesn't satisfy Maxwell's equations but gives the idea

# We'll just use 5 replications of the data 
# We'll also assume 3 GeV momentum and no vertical deflection
sddscombine BvsX.sdds BvsX.sdds BvsX.sdds BvsX.sdds BvsX.sdds  -pipe=out \
    | sddsprocess -pipe "-define=col,y,i_page 2 - 0.01 *,units=m" \
    "-define=parameter,H,3e9 c_mks /,units=T*m" \
    "-define=col,xp,By H /" \
    "-process=xp,max,%sMax" \
    "-define=col,xpFactor,xp xpMax /" \
    "-define=col,ypFactor,0" \
    | sddscombine -pipe -merge \
    | sddssort -pipe=in -column=y,incr -column=x,incr deflectionMap.sdds

sddscontour -shade -xyz=x,y,xpFactor deflectionMap.sdds
