#set ye = 2010; set me = 06; set de = 15
#set ys = 2010; set ms = 06; set ds = 14
#set ye = 2010; set me = 07; set de = 05
#set ys = 2010; set ms = 07; set ds = 05
#set ye = 2010; set me = 11; set de = 13
#set ys = 2010; set ms = 11; set ds = 13
#set ye = 2011; set me = 05; set de = 01
#set ys = 2011; set ms = 05; set ds = 01
#set ye = 2011; set me = 07; set de = 16
#set ys = 2011; set ms = 07; set ds = 16
#set ye = 2012; set me = 08; set de = 26
#set ys = 2012; set ms = 08; set ds = 26
#set ye = 2012; set me = 08; set de = 29
#set ys = 2012; set ms = 08; set ds = 10
#set ys = 2012; set ms = 10; set ds = 06
#set ys = 2012; set ms = 08; set ds = 01
#set ys = 2013; set ms = 01; set ds = 01

#set ys = 1996; set ms = 01; set ds = 01
#set ye = 2013; set me = 02; set de = 09
# Tue Feb 19 12:28:37 EST 2013
#set ys = 2013; set ms = 02; set ds = 08
#set ye = 2013; set me = 02; set de = 17
# Mon Feb 25 10:48:48 EST 2013
#set ys = 2013; set ms = 02; set ds = 15
#set ye = 2013; set me = 02; set de = 23
# Thu Mar 28 11:17:12 EDT 2013
#set ys = 2013; set ms = 02; set ds = 23
#set ye = 2013; set me = 03; set de = 27

# CONVERTED to use sh_makeeqdef  TAH 130516
#set st = "2013/03/26" ; set en = "2013/05/14"
#set st = "2013/05/14" ; set en = "2013/06/04"
#set st = "2013/06/03" ; set en = "2013/09/13"  # Sun Sep 15 17:34:45 EDT 2013

#set st = "2010/06/03" ; set en = "2010/08/01"
#set st = "2013/09/13" ; set en = "2014/01/08"
#set st = "2014/01/01" ; set en = "2014/02/16"
#set st = "2014/02/15" ; set en = "2014/03/16"
#set st = "2014/04/13" ; set en = "2014/05/14"
#set st = "2014/05/13" ; set en = "2014/06/15"
#set st = "2014/06/13" ; set en = "2014/07/14"
set st = "2014/07/13" ; set en = "2014/08/12"

set posf = 'Week_All.Pos'
if( ! -e $posf ) then 
     grep 'NEU Reference position' ../COM/ts_com/*pbo*final_nam08.pos | awk '{printf(" %9.6f %9.6f %s_GPS \n",$6,$5,substr($0,15,4) ) }' >! $posf
     awk '{print substr($3,1,4)}' $posf >! $posf:r.sites
     echo Created  $posf and  $posf:r.sites
endif
set MM = 3.5
sh_makeeqdef -start $st -end $en -f $posf -m $MM


exit


echo Range $ys $ms $ds to $ye $ms $ds
set MM = 3.5  # Set minimum magnitude to return from catalog (See grep below also)
#wget -O eqf.out "http://eqint.cr.usgs.gov/neic/cgi-bin/epic/epic.cgi?SEARCHMETHOD=2&CLAT=0.0&CLON=0.0&CRAD=0.0&FILEFORMAT=4&SEARCHRANGE=HH&SLAT2=85&SLAT1=19&SLON2=-30&SLON1=-180&SUBMIT=Submit+Search&SYEAR=${ys}&SMONTH=${ms}&SDAY=${ds}&EYEAR=${ye}&EMONTH=${me}+&EDAY=${de}&LMAG=${MM}&UMAG=9&NDEP1=&NDEP2=&IO1=&IO2="
 wget -O eqf.out "http://neic.usgs.gov/cgi-bin/epic/epic.cgi?SEARCHMETHOD=2&CLAT=0.0&CLON=0.0&CRAD=0.0&FILEFORMAT=4&SEARCHRANGE=HH&SLAT2=85&SLAT1=06&SLON2=-30&SLON1=-180&SUBMIT=Submit+Search&SYEAR=${ys}&SMONTH=${ms}&SDAY=${ds}&EYEAR=${ye}&EMONTH=${me}+&EDAY=${de}&LMAG=${MM}&UMAG=9&NDEP1=&NDEP2=&IO1=&IO2="
#
# Now see how many large EQ's we have
grep ' PDE' eqf.out | awk -v MM=${MM} '{if ( $9 > MM ) {{s=s+1} {cs = 3.e-8*13^($9)} {cd = 2.5e-3*5^($9)}{printf("* EQ_DEF M %3.1f\n eq_def %2.2d  %8.3f %8.3f %8.1f 8 %4d %2d %2d %2d %2d %9.4f\n eq_rename %2.2d\n eq_coseis %2.2d  0.001 0.001 0.001 %10.3f %10.3f %10.3f\n",$9,s,$6,$7,cd+8,$2,$3,$4,substr($5,1,2),substr($5,3,2),cs,s,s,cs,cs,cs)}}}' >! eq.def
#
\mv eqf.out eq_${ye}${me}${de}.out
\cp eq.def  eq_${ye}${me}${de}.def
echo "* EQDEFS for $ys $ms $ds to $ye $me $de Generated `date`" >! eq_${ye}${me}${de}.def
echo "Created eq_${ye}${me}${de}.out and eq_${ye}${me}${de}.def"
#
# Now see how many of these earthquakes may have displaced sites
#
grep 'eq_def' eq.def | awk '{print $2, $3,$4,$5,$7,$8,$9,$10,$11,$12}' >&! eqo.dat
#
# Now loop over earthqakes
set posf = 'Week_All.Pos'
#grep '^ ' ../../control/PBO_all.pos | awk '{print " ",$1,$2,$13}' >! $posf
# MOD TAH: 121006 Generate site list from POS files
grep 'NEU Reference position' ../COM/ts_com/*pbo*final_nam08.pos | awk '{printf(" %9.6f %9.6f %s_GPS \n",$6,$5,substr($0,15,4) ) }' >! $posf
echo "Testing `wc -l $posf` site locations"
set num = `wc -l eqo.dat | awk '{print $1}'`
set line = 0
set snum = 0

while ( $line < $num )
   @ line = $line + 1
   set eqa = `tail -n +${line} eqo.dat | head  -1`
#   echo "Testing EQ "$eqa[1]" Dist "$eqa[4] " POS "$eqa[2] $eqa[3] " Date "$eqa[5] $eqa[6] $eqa[7]" Time " $eqa[8] $eqa[9] " CoS " $eqa[10]
   set nn = `grep '^ ' $posf | awk -v lat=$eqa[2] -v lng=$eqa[3] -v d=$eqa[4] -v eq=$eqa[1] '{cd=sin($2/57.296)*sin(lat/57.296)+cos($2/57.296)*cos(lat/57.296)*cos((lng-$1)/57.296)} {dist = atan2(sqrt(1-cd^2),cd)*6378.15} {if ( (sin($2/57.296)*sin(lat/57.296)+cos($2/57.296)*cos(lat/57.296)*cos((lng-$1)/57.296)) > cos(d/6378.153)) {print " EQ ",eq,$3,dist,d}}' | wc -l`
   if( $nn > 0 ) then
      @ snum ++
#      echo "SEQ Earthquake # $snum" 
      echo "*---------------------------------------------------------------------" >> eq_${ye}${me}${de}.def
      echo "SEQ Earthquake # $snum" >> eq_${ye}${me}${de}.def
      grep '^ ' $posf | awk -v lat=$eqa[2] -v lng=$eqa[3] -v d=$eqa[4] -v eq=$eqa[1] -v cs=$eqa[10] '{cd=sin($2/57.296)*sin(lat/57.296)+cos($2/57.296)*cos(lat/57.296)*cos((lng-$1)/57.296)} {dist = atan2(sqrt(1-cd^2),cd)*6378.15} {cm=cs*(8.0/dist)^2} {if ( (sin($2/57.296)*sin(lat/57.296)+cos($2/57.296)*cos(lat/57.296)*cos((lng-$1)/57.296)) > cos(d/6378.153)) {printf("* EQ %2.2d %s %10.2f %10.2f CoS %10.1f mm\n",eq,$3,dist,d,cm*1000.)}}' >> eq_${ye}${me}${de}.def
      set lns = `echo $line | awk '{print ($1-1)*4+1}'`
      tail -n +$lns eq.def | head -4>> eq_${ye}${me}${de}.def
   endif
end

