#!/bin/csh -f

set SNR = 10.0
set nx = 512
set ny = 3072
set master = `ls ../../*PRM | head -1`

set PRF = `grep PRF $master |awk -F"=" '{print $2}'`
set SC_vel = `grep SC_vel $master|awk -F"=" '{print $2}'`
set earth_radius = `grep earth_radius $master|awk -F"=" '{print $2}'`
set SC_height = `grep SC_height $master|awk -F"=" '{print $2}'`
set rng_samp_rate = `grep rng_samp_rate $master | head -1 | awk '{print $3}'`

set ground_vel = `echo $SC_vel $earth_radius $SC_height | awk '{print $1/sqrt(1+$3/$2)} '`
echo "ground velocity: " $ground_vel
set azi_size = `echo $ground_vel $PRF | awk '{printf "%.6f",$1/$2}' `
set rng_size = `echo $rng_samp_rate | awk '{printf "%.6f", 299792458.0/$1/2}' `
echo "azi pixel size", $azi_size
echo "rng pixel size", $rng_size

awk '{if ($5 > '$SNR' ) print $1,$3,$4,$5}' freq_xcorr.dat > azi.dat
awk '{if ($5 > '$SNR' ) print $1,$3,$2,$5}' freq_xcorr.dat > rng.dat

set xmin = `gmt gmtinfo azi.dat -C |awk '{print $1}'`
set xmax = `gmt gmtinfo azi.dat -C |awk '{print $2}'`
set ymin = `gmt gmtinfo azi.dat -C |awk '{print $3}'`
set ymax = `gmt gmtinfo azi.dat -C |awk '{print $4}'`

set xinc = `echo $xmax $xmin $nx |awk '{printf "%d", ($1-$2)/($3-1)}'`
set yinc = `echo $ymax $ymin $ny |awk '{printf "%d", ($1-$2)/($3-1)}'`

echo "xinc:"$xinc"  yinc:"$yinc

gmt blockmedian azi.dat -R$xmin/$xmax/$ymin/$ymax -I$xinc/$yinc -Wi | awk '{print $1, $2, $3}' > azi_b.dat
gmt blockmedian rng.dat -R$xmin/$xmax/$ymin/$ymax -I$xinc/$yinc -Wi | awk '{print $1, $2, $3}' > rng_b.dat

gmt xyz2grd azi_b.dat -R$xmin/$xmax/$ymin/$ymax  -I$xinc/$yinc -Gaoff.grd
gmt grdmath aoff.grd $azi_size MUL = azi_offset.grd
gmt xyz2grd rng_b.dat -R$xmin/$xmax/$ymin/$ymax  -I$xinc/$yinc -Groff.grd
gmt grdmath roff.grd $rng_size MUL = rng_offset.grd