#!/bin/csh

if ($#argv < 1) then
  echo ""
  echo "read_mle: neu_file mle_file"
  echo ""
  exit 1
endif


set neu_file = $1
set mle_file = $2
set noff = `grep offset $neu_file | wc -l`
set toff = `grep offset $neu_file | awk '{print $3}'`

echo -----------------------------------
echo $neu_file "               INPUT FILE"
echo $noff "                          NUMBER OF OFFSETS"
echo ----------------------------------------
echo "     Value     Sigma"
echo ------------ NORTH ---------------------

echo `grep "+NORT  INTER " $mle_file | awk '{print $4}'`\
     `grep "+NORT  INTER " $mle_file | awk '{print $6}'`\
     | awk '{printf"%10.4f%10.4f    INT_NORTH\n",$1,$2}'

echo `grep "+NORT  SLOPE " $mle_file | awk '{print $4}'`\
     `grep "+NORT  SLOPE " $mle_file | awk '{print $6}'`\
     | awk '{printf"%10.4f%10.4f    VEL_NORTH\n",$1,$2}'

echo `grep "+NORT  SIN" $mle_file | head -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    ANNUAL SIN NORT \n",$1,$2}'
echo `grep "+NORT  COS" $mle_file | head -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    ANNUAL COS NORT \n",$1,$2}'

echo `grep "+NORT  SIN" $mle_file | tail -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    SEMI-ANNUAL SIN NORT \n",$1,$2}'
echo `grep "+NORT  COS" $mle_file | tail -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    SEMI-ANNUAL COS NORT \n",$1,$2}'

set n = 1
while ($n <= $noff ) 
 echo `grep "+NORT  OFFSET " $mle_file | awk '{print $4,$6}' | head -n$n | tail -1` $toff[$n] $n \
 | awk '{printf"%10.4f%10.4f %10.4f    OFFSET %3s\n",$1,$2,$3,$4}'
@ n++
end

echo ------------ EAST ---------------------

echo `grep "+EAST  INTER " $mle_file | awk '{print $4}'`\
     `grep "+EAST  INTER " $mle_file | awk '{print $6}'`\
     | awk '{printf"%10.4f%10.4f    INT_EAST\n",$1,$2}'
                                                                                                                    
echo `grep "+EAST  SLOPE " $mle_file | awk '{print $4}'`\
     `grep "+EAST  SLOPE " $mle_file | awk '{print $6}'`\
     | awk '{printf"%10.4f%10.4f    VEL_EAST\n",$1,$2}'

echo `grep "+EAST  SIN" $mle_file | head -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    ANNUAL SIN EAST \n",$1,$2}'
echo `grep "+EAST  COS" $mle_file | head -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    ANNUAL COS EAST \n",$1,$2}'
                                                                                                                                           
echo `grep "+EAST  SIN" $mle_file | tail -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    SEMI-ANNUAL SIN EAST \n",$1,$2}'
echo `grep "+EAST  COS" $mle_file | tail -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    SEMI-ANNUAL COS EAST \n",$1,$2}'

set n = 1
while ($n <= $noff )
 echo `grep "+EAST  OFFSET " $mle_file | awk '{print $4,$6}' | head -n$n | tail -1` $toff[$n] $n \
 | awk '{printf"%10.4f%10.4f %10.4f    OFFSET %3s\n",$1,$2,$3,$4}'
@ n++
end

echo ------------ VERT ---------------------

echo `grep "+VERT  INTER " $mle_file | awk '{print $4}'`\
     `grep "+VERT  INTER " $mle_file | awk '{print $6}'`\
     | awk '{printf"%10.4f%10.4f    INT_VERT\n",$1,$2}'
                                                                                                                    
echo `grep "+VERT  SLOPE " $mle_file | awk '{print $4}'`\
     `grep "+VERT  SLOPE " $mle_file | awk '{print $6}'`\
     | awk '{printf"%10.4f%10.4f    VEL_VERT\n",$1,$2}'

echo `grep "+VERT  SIN" $mle_file | head -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    ANNUAL SIN VERT \n",$1,$2}'
echo `grep "+VERT  COS" $mle_file | head -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    ANNUAL COS VERT \n",$1,$2}'
                                                                                                                                           
echo `grep "+VERT  SIN" $mle_file | tail -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    SEMI-ANNUAL SIN VERT \n",$1,$2}'
echo `grep "+VERT  COS" $mle_file | tail -1` | awk '{print $4,$6}' | awk '{printf"%10.4f%10.4f    SEMI-ANNUAL COS VERT \n",$1,$2}'

set n = 1
while ($n <= $noff )
 echo `grep "+VERT  OFFSET " $mle_file | awk '{print $4,$6}' | head -n$n | tail -1` $toff[$n] $n \
 | awk '{printf"%10.4f%10.4f %10.4f    OFFSET %3s\n",$1,$2,$3,$4}'
@ n++
end

