#!/bin/csh -f
# divergenve of a wave field in 2d case


if ( $#argv != 3 ) then
echo "Usage:    Diverg ux uz div " 
exit(-1)
endif

echo " need the right order:  UX UZ"

set UX = $1  ; echo "esize=4">>$UX
set UZ = $2  ; echo "esize=4">>$UZ
set OUT = $3

echo "1.,-1" >jj
atoF <jj >2pt.h
echo "lag=0" >> 2pt.h
\rm jj

# first axis  is  z-axis
# duzdz
Conv filtin=2pt.h < $UZ > jdivz
# duxdx
Transp plane=12 < $UX | Conv filtin=2pt.h |Transp plane=12 >jdivx
# div = duzdz + duxdx
Add jdivx jdivz mode=sum > $OUT
# \rm jdivx jdivz 
