The first step is to check the navigation recorded in the *.gps file for each line of acquisition. I extracted the $GPRMC line from the CRP GPS files using an AWK script. (* refers to the line name, which for this day are L7F1, L8F1, F9L1, F10L1, L10F2, L11F1, L12F1, L13F1, L14F1, L15F1, L16F1, L17F1, L18F1, and L19F1).
AWK script awkit_tmp:
BEGIN {
FS=","
}
{
FS= ","
ARGC = 2
if ($1=="$GPRMC") {
gpstime = $2
hr = substr($2,1,2)
min = substr($2,3,2)
sec = substr($2,5,2)
latdeg = substr($4,1,2)
latmin = substr($4,3,6)
londeg = substr($6,2,2)
lonmin = substr($6,4,6)
declat = latdeg + (latmin/60)
declon = (londeg + (lonmin/60)) * -1
date = $10
day = substr($10,1,2)
month = substr($10,3,2)
year = substr($10,5,2)
printf("%8.6f, %8.6f, %s:%s:%s, %s, %s\n",declon, declat, hr, min, sec, date, ARGV[2])
}
}
This AWK script was initialized by "doawk" - shell script run under CYGWIN (UNIX like environment that runs under Windows):
files=`ls *.gps | cut -d. -f1`
for file in $files
do
awk -f awkhold $file.gps $file >> resnav_jd250.txt
done
