mirror of
https://github.com/thetechprepper/emcomm-tools-os-community.git
synced 2026-09-17 15:41:44 -04:00
thetechprepper
254 lines
7.7 KiB
Bash
Executable File
254 lines
7.7 KiB
Bash
Executable File
#!/bin/bash
|
|
#
|
|
# Author : Gaston Gonzalez
|
|
# Date : 23 May 2023
|
|
# Updated : 24 August 2025
|
|
# Purpose : Offline HF prediction using voacapl
|
|
#set -e
|
|
set -o pipefail
|
|
trap 'last_command=$current_command; current_command=$BASH_COMMAND' DEBUG
|
|
trap 'echo -e "${RED}\"${last_command}\" command failed with exit code $?.${NC}"' ERR
|
|
|
|
source /opt/emcomm-tools/bin/et-common
|
|
|
|
lookup_station() {
|
|
local type="$1" # call | grid | latlon
|
|
local value="$2" # callsign | grid | "lat,lon"
|
|
local json_file="$3"
|
|
|
|
case "$type" in
|
|
call)
|
|
curl -f -s "http://localhost:1981/api/license?callsign=${value}" > "$json_file"
|
|
if [[ $? -ne 0 ]]; then
|
|
echo -e "${RED}Location not found for callsign: ${value}.${NC}" >&2
|
|
exit 1
|
|
fi
|
|
;;
|
|
grid)
|
|
curl -f -s "http://localhost:1981/api/geo/grid?gridSquare=${value}" | jq '{lat: .position.lat, lon: .position.lon}' > "$json_file"
|
|
if [[ $? -ne 0 ]]; then
|
|
echo -e "${RED}Error converting grid to lat/lon using value: ${value}.${NC}" >&2
|
|
exit 1
|
|
fi
|
|
;;
|
|
latlon)
|
|
# Write directly to JSON for consistency
|
|
local lat lon
|
|
IFS=',' read -r lat lon <<< "$value"
|
|
if [[ -z $lat || -z $lon ]]; then
|
|
echo -e "${RED}Invalid lat,lon format: $value. Expected 'lat,lon'.${NC}" >&2
|
|
exit 1
|
|
fi
|
|
jq -n --arg lat "$lat" --arg lon "$lon" '{lat: ($lat|tonumber), lon: ($lon|tonumber)}' > "$json_file"
|
|
;;
|
|
*)
|
|
echo -e "${RED}Invalid station type: $type${NC}" >&2
|
|
exit 1
|
|
;;
|
|
esac
|
|
|
|
# Extract lat/lon regardless of type
|
|
local lat lon
|
|
lat=$(jq .lat "$json_file")
|
|
lon=$(jq .lon "$json_file")
|
|
|
|
echo "$lat $lon"
|
|
}
|
|
|
|
ET_SSN="${HOME}/.local/share/emcomm-tools/voacap/ssn.txt"
|
|
if [[ ! -e ${ET_SSN} ]]; then
|
|
echo -e "${RED}Sunspot number file not available:${WHITE}${ET_SSN}${NC}"
|
|
echo -e "${YELLOW}Try running the following to fetch the sunspot numbers (requires Internet):${NC}"
|
|
echo -e "${WHITE}cd ${HOME}/.local/share/emcomm-tools/voacap${NC}"
|
|
echo -e "${WHITE}./fetch-ssn.sh${NC}"
|
|
|
|
exit 1
|
|
fi
|
|
|
|
usage() {
|
|
echo "Usage: $(basename $0) [OPTIONS]"
|
|
echo
|
|
echo "Station options (choose one per TX and RX):"
|
|
echo " --tx-call CALLSIGN Transmitting station callsign"
|
|
echo " --tx-grid GRID Transmitting station Maidenhead grid"
|
|
echo " --tx-latlon LAT,LON Transmitting station coordinates in decimal degrees"
|
|
echo " --rx-call CALLSIGN Receiving station callsign"
|
|
echo " --rx-grid GRID Receiving station Maidenhead grid"
|
|
echo " --rx-latlon LAT,LON Receiving station coordinates in decimal degrees"
|
|
echo
|
|
echo "Other options:"
|
|
echo " -p POWER Output power [5|20|100|500|1500]"
|
|
echo " -m MODE Mode [AM|CW|JS8|SSB]"
|
|
echo
|
|
|
|
exit 1
|
|
}
|
|
|
|
tx_type=""
|
|
tx_value=""
|
|
rx_type=""
|
|
rx_value=""
|
|
power=""
|
|
mode=""
|
|
|
|
while [[ $# -gt 0 ]]; do
|
|
case "$1" in
|
|
--tx-call) tx_type="call"; tx_value="$2"; shift 2 ;;
|
|
--tx-grid) tx_type="grid"; tx_value="$2"; shift 2 ;;
|
|
--tx-latlon) tx_type="latlon"; tx_value="$2"; shift 2 ;;
|
|
--rx-call) rx_type="call"; rx_value="$2"; shift 2 ;;
|
|
--rx-grid) rx_type="grid"; rx_value="$2"; shift 2 ;;
|
|
--rx-latlon) rx_type="latlon"; rx_value="$2"; shift 2 ;;
|
|
-p) power="$2"; shift 2 ;;
|
|
-m) mode="$2"; shift 2 ;;
|
|
-*|--*) echo "Unknown option: $1"; usage ;;
|
|
*) echo "Unexpected argument: $1"; usage ;;
|
|
esac
|
|
done
|
|
|
|
# Validate required args
|
|
[[ -z $tx_type || -z $rx_type || -z $power || -z $mode ]] && usage
|
|
|
|
#echo "TX ($tx_type): $tx_value"
|
|
#echo "RX ($rx_type): $rx_value"
|
|
#echo "Power: $power"
|
|
#echo "Mode: $mode"
|
|
|
|
# Allow for case-insensitve match of the user-defined operating mode
|
|
MODE=$mode
|
|
case "${MODE,,}" in
|
|
js8)
|
|
MD="13.0"
|
|
;;
|
|
cw)
|
|
MD="24.0"
|
|
;;
|
|
ssb)
|
|
MD="38.0"
|
|
;;
|
|
am)
|
|
MD="49.0"
|
|
;;
|
|
*)
|
|
echo -e "${RED}Unknown mode: ${MODE}${NC}" >&2
|
|
usage
|
|
exit 1
|
|
;;
|
|
esac
|
|
|
|
ET_VOA_WORKING_DIR=$HOME/itshfbc/run
|
|
ET_VOA_REPORT=${ET_VOA_WORKING_DIR}/voacapl.txt
|
|
INP=${ET_VOA_WORKING_DIR}/voacapx.dat
|
|
OUT=${ET_VOA_WORKING_DIR}/voacapx.out
|
|
|
|
# Year = 2023
|
|
YEAR=$(echo $(date '+%Y'))
|
|
|
|
# Month (May) = 5.00
|
|
MONTH=$(date +'%m')
|
|
MONTH_FMT=$(date +'%-m.00')
|
|
|
|
#######################################################################
|
|
# TX Antenna
|
|
#######################################################################
|
|
|
|
if ! result=$(lookup_station "$tx_type" "$tx_value" tx-station.json); then
|
|
echo "Try specifying a grid or lat,lon for the transmitting station."
|
|
echo " --tx-grid GRID Transmitting station Maidenhead grid"
|
|
echo " --tx-latlon LAT,LON Transmitting station coordinates in decimal degrees"
|
|
exit 1
|
|
fi
|
|
read TL TK <<< "$result"
|
|
|
|
TL1=$( awk -v n1=$TL -v n2=90 -v n3=-90 'BEGIN {if (n1<n3 || n1>n2) printf ("%s", "a"); else printf ("%.2f", n1);}' )
|
|
|
|
# add North or South (N/S)
|
|
TLA=$( awk -v n1=$TL1 -v n2=0 'BEGIN {if (n1<n2) { n1=substr(n1,2); printf ("%6sS", n1); } else printf ("%6sN", n1);}' )
|
|
|
|
TK1=$( awk -v n1=$TK -v n2=180 -v n3=-180 'BEGIN {if (n1<n3 || n1>n2) printf ("%s", "a"); else printf ("%.2f", n1);}' )
|
|
|
|
# add East or West (E/W)
|
|
TLO=$( awk -v n1=$TK1 -v n2=0 'BEGIN {if (n1<n2) { n1=substr(n1,2); printf ("%7sW", n1); } else printf ("%7sE", n1);}' )
|
|
|
|
|
|
#######################################################################
|
|
# RX Antenna
|
|
#######################################################################
|
|
|
|
if ! result=$(lookup_station "$rx_type" "$rx_value" rx-station.json); then
|
|
echo "Try specifying a grid or lat,lon for the receiving station."
|
|
echo " --rx-grid GRID Receiving station Maidenhead grid"
|
|
echo " --rx-latlon LAT,LON Receiving station coordinates in decimal degrees"
|
|
exit 1
|
|
fi
|
|
read RL RK <<< "$result"
|
|
|
|
RL1=$( awk -v n1=$RL -v n2=90 -v n3=-90 'BEGIN {if (n1<n3 || n1>n2) printf ("%s", "a"); else printf ("%.2f", n1);}' )
|
|
|
|
# add North or South
|
|
RLA=$( awk -v n1=$RL1 -v n2=0 'BEGIN {if (n1<n2) { n1=substr(n1,2); printf ("%6sS", n1); } else printf ("%6sN", n1);}' )
|
|
|
|
RK1=$( awk -v n1=$RK -v n2=180 -v n3=-180 'BEGIN {if (n1<n3 || n1>n2) printf ("%s", "a"); else printf ("%.2f", n1);}' )
|
|
|
|
# add East or West
|
|
RLO=$( awk -v n1=$RK1 -v n2=0 'BEGIN {if (n1<n2) { n1=substr(n1,2); printf ("%7sW", n1); } else printf ("%7sE", n1);}' )
|
|
|
|
# Power settings
|
|
# Use 80% of user-defined power and express in killiwatts
|
|
# TODO: Make calculation dynamic.
|
|
PWR=$power
|
|
PW="0.0800"
|
|
echo "TX Power: ${PWR} watts"
|
|
if [ "$PWR" = "5" ]; then
|
|
PW="0.0040"
|
|
elif [ "$PWR" = "20" ]; then
|
|
PW="0.0160"
|
|
elif [ "$PWR" = "100" ]; then
|
|
PW="0.0800"
|
|
elif [ "$PWR" = "500" ]; then
|
|
PW="0.4000"
|
|
elif [ "$PWR" = "1500" ]; then
|
|
PW="1.2000"
|
|
else
|
|
echo "Unsupported power level."
|
|
usage
|
|
exit 1
|
|
fi
|
|
|
|
# Format for 117
|
|
# NOTE: The card requires <SSN>. (117.)
|
|
ssn=`grep "$YEAR $MONTH" ${ET_SSN} | awk '{print $5}' | cut -d"." -f1`
|
|
|
|
# read more about fine-tuning your input file:
|
|
# http://www.voacap.com/voacapw.html
|
|
# http://www.voacap.com/frequency.html
|
|
|
|
echo -e "\n"
|
|
|
|
cat << END | tee $INP
|
|
LINEMAX 999 number of lines-per-page
|
|
COEFFS CCIR
|
|
TIME 1 24 1 1
|
|
MONTH $YEAR $MONTH_FMT
|
|
SUNSPOT $ssn.
|
|
LABEL $TX$RX
|
|
CIRCUIT $TLA $TLO $RLA $RLO S 0
|
|
SYSTEM 1. 155. 3.00 90. $MD 3.00 0.10
|
|
FPROB 1.00 1.00 1.00 0.00
|
|
ANTENNA 1 1 2 30 0.000[samples/sample.00 ] 0.0 $PW
|
|
ANTENNA 2 2 2 30 0.000[samples/sample.00 ] 0.0 0.0000
|
|
FREQUENCY 3.60 5.30 7.1010.1014.1018.1021.1024.9028.20 0.00 0.00
|
|
METHOD 30 0
|
|
BOTLINES 8 12 21
|
|
TOPLINES 1 2 3 4 6
|
|
EXECUTE
|
|
QUIT
|
|
|
|
END
|
|
|
|
voacapl -s ~/itshfbc
|
|
echo "voacapl exited with $? using in=$INP and out=$OUT"
|
|
|
|
(tail -n +27 $OUT | head -n 5 && grep -e'- REL' $OUT && grep -e'- S DBW' $OUT) | ./rel.pl > ${ET_VOA_REPORT}
|
|
cat ${ET_VOA_REPORT}
|