#!/bin/bash

#
# sample the LOS tiles
#

rm los e n u std
echo "" > los

grdmath sample.grd 1 EQ 0 NAN = mask.1.grd
grdmath sample.grd 2 EQ 0 NAN = mask.2.grd
grdmath sample.grd 3 EQ 0 NAN = mask.3.grd
grdmath sample.grd 4 EQ 0 NAN = mask.4.grd
grdmath sample.grd 5 GE 0 NAN = mask.5.grd
grdmath sample.grd 5 GE 0 NAN = mask.6.grd
grdmath sample.grd 5 GE 0 NAN = mask.7.grd

cat numbers | while read lon lat num
do
cat list | while read a 
do
losgrd="maskmodel_"$a".ll.grd"
egrd="lle_"$a".grd"
ngrd="lln_"$a".grd"
ugrd="llu_"$a".grd"
stdgrd="std_"$a".ll.grd"

echo $a
echo $num

# mask
grdsample manmask.grd `grdinfo -I- $losgrd` `grdinfo -I $losgrd` -Gmask1.grd -F
grdsample mask.$num.grd `grdinfo -I- $losgrd` `grdinfo -I $losgrd` -Gmask2.grd

# blockmedian
# los
echo los
grdmath mask1.grd $losgrd MUL mask2.grd MUL = tmp.grd 
grd2xyz tmp.grd -bo | blockmedian -I$lon/$lat `grdinfo -I- tmp.grd` -bi > lostmp 
cat los lostmp >> tmp 
mv tmp los
wc -l los

# std
echo std
grdmath 1 $losgrd ISNAN SUB 0 NAN = mask3.grd
grdcut $stdgrd `grdinfo -I- $losgrd` -Gstdtmp.grd
grdmath mask1.grd stdtmp.grd MUL mask2.grd MUL mask3.grd MUL = tmp.grd
grd2xyz tmp.grd -bo | blockmedian -I$lon/$lat `grdinfo -I- tmp.grd` -bi >> std
wc -l std

awk '{print $1,$2}' lostmp > ll
grdtrack ll -G$egrd `grdinfo -I- $egrd` >> e
grdtrack ll -G$ngrd `grdinfo -I- $ngrd` >> n
grdtrack ll -G$ugrd `grdinfo -I- $ugrd` >> u

done
done

rm tmp.grd mask1.grd mask2.grd etmp.grd ntmp.grd utmp.grd stdtmp.grd mask3.grd
rm ll