VOACAP Quick Guide Prelucrat
VOACAP Quick Guide Prelucrat
VOACAP Quick Guide Prelucrat
http://www.voacap.com/
What's New?
Current International Smoothed Sunspot Numbers, Ri, to be used in VOACAP, from NOAA/NGDC
(updated until 2018)
Current Space Weather: Solar Wind, Magnetogram, X-Rays, Protons
QTH Locator - find coordinates and Maidenhead locators on Google Maps
How to Run VOACAP under Apple Macintosh OS X? - a quick guide to making text-based
VOACAP predictions
Elevation-Angle Analyzer - generate statistical elevation-angle data for HFTA
o Statistical elevation-angle distribution - a case study
Graphical User Interface for VOACAP for Linux
Rockwell Collins offers VOACAP User's Guide on CD-ROM for free (shipping charges apply)
World Radio Broadcasting Stations - explore various BC transmitting sites on Google Maps
My blog in Finnish: VOACAP, ennusteet ja antennit
Sunrise and sunset times around the world
o Today's Sunrise and Sunset Times by DXCC country
o Today's Ionospheric Grayline by DXCC country
o Today's Sunrise and Sunset at Global HF BC Transmitting Sites
1. Downloading VOACAP
Download from Greg Hand's site, all new updates of the ITSHFBC package are available here
Alternative direct download, 8 Dec 2009 version (3.1 MB)
Download from the home of VOACAP, NTIA/ITS, this ITSHFBC page has not been updated for long
o The 8 Dec 2009 version
o The 26 Mar 2009 version
Download the Fortran source codes of the VOACAP Propagation Prediction Engine, in a special
directory as ZIP files.
VOACAP for Linux by J.A.Watson (M0DNS), an exciting project to make VOACAP available for
Linux users. Please contribute!
o Release 3 is out! Download the new VOACAP Linux GUI, Release 3!(1.5 MB gzipped file).
Now supports Area Coverage Maps with Jim Watson's VoaAreaPlot! This software is supplied
'as is' without any warranty, implied or otherwise. right-click and save the target. All users
of Release 1 are strongly encouraged to upgrade, due to the bugs in R1.
o VOACAP for Linux Software
The voacapl application: voacapl-0.5.7.tar.gz
The itshfbc directory: itshfbcFiles18Apr08.tar.gz
The python scripts to display the output
2. VOACAP - Backgrounds & Running the Prediction
General
VOACAP Overview
Screenshots of VOACAP
Modelling the real-life HF communication system into VOACAP
Setting VOACAP input parameters
Textual circuit prediction (Run > Circuit)
Graphical circuit prediction (Run > Graph)
Area coverage prediction (VOAAREA)
The calculation engine, VOACAPW
An introduction to VOACAPW.EXE
The FREQUENCY card in VOACAPX.DAT
1
Graphs & Colours
VOACAP & VOAAREA contours
VOACAP & VOAAREA plot colors
Y-axis scaling for VOACAP Time & Distance graphs
3. VOACAP - Understanding the Prediction
General
Ten Common Mistakes in Using VOACAP
Field Strength vs. Signal-to-Noise
Choosing the correct Sunspot Number
VOACAP Predictions at Low Frequencies
Comparing VOACAP, ICEPAC and REC533
Method 25: Ionospheric Parameters and Semi-thickness
MUF, SNR and REL
The MUF and SNR Distribution: Choosing the Best Frequency
Understanding Above-the-MUF Predictions
Calculating the MUFdays
Maintaining the Required Grade of Service: Calculating the Circuit Reliability for a Given Hour
Signal Power & Noise Power
Summing Signal Power from Multiple Modes
Calculating the Received RF Noise Power
Antennas & Radiation Angles
Method 15, or Transmitter & Receiver AntennaPatterns
VOACAP Angle Predictions
Minimum Angle in VOACAP
Predicting the Takeoff Angle Parameter TANGLE
Conversions
S DBW to S Meter Correlation Table
The Z Tables
4. VOACAP - Case Studies
VOACAP Case: Make your beacon predictions
VOAAREA Case: RUV Gufunes
5. VOACAP-related Papers
Using Propagation Predictions for HF DXing [2010] Visalia 2010 presentation by Dean Straw, N6BV.
Contest Antennas. DX or Domestic, What’s Your Pleasure? [2009] Visalia 2009 presentation by Dean
Straw, N6BV.
What’s the Deal About “NVIS”?, [2005] article courtesy by Dean Straw, N6BV
Another Way to View Propagation Predictions for DXing and Contesting [2006] Visalia 2006
presentation by Dean Straw, N6BV.
In Search of the Ultima Thule. Can VOACAP Help Find the Ultimate QTH? [2006] CCF/OHDXF
presentation by Jari Perkiömäki & George Lane.
Improved Guidelines for Automatic Link Establishment Operations [2005] by George Lane
o Graphs for the paper [2005]
Review of the High Frequency Ionospheric Communications Enhanced Profile Analysis & Circuit
(ICEPAC) Prediction Program [2005] by George Lane
Required Signal-to-Interference Ratios for Shortwave Broadcasting, [7.2 MB] by George Lane, 1997
Predicting Signal-to-Interference Probability in the High Frequency Band, [60 kB] by George Lane and
Greg Hand, 1997
Strategies Using Propagation Predictions for DXing and Contesting, by Dean Straw, N6BV
What is VOACAP Trying to Tell Me? by Dean Straw, N6BV
How Does VOACAP Compare With These Measurements, by Dean Straw, N6BV
6. NTIA/ITS HF software HELP Topics
General information on the NTIA/ITS HF Propagation Analysis Package
2
VOACAP. General information - Point-to-Point - Area Coverage - Signal-to-Interference.
Voice of America Coverage Analysis Program
o VOACAP Frequently Asked Questions (22-Feb-96 Version)
ICEPAC. General information - Point-to-Point - Area Coverage - Signal-to-Interference.
Ionospheric Communications Enhanced Profile Analysis and Circuit Prediction Program
REC533. General information - Point-to-Point - Area Coverage.
ITU-RS Recommendation 533. This has been updated to the ITU 1996 Version 3.
Data Files Frequently Asked Questions
Helpful documents on using the propagation models
o History (Part 1) of HF propagation models by Ray Rosich, 1978
o Introduction to the ICEPAC Technical Description
o Excerpts from Original ICEPAC User's Manual
o Corrections to Man-made Noise
o Note on MUFDAYs
o High Latitude Errors
o Corrected Signal-to-Noise Ratio Distribution
o A Comparison of Broadcast Quality to VOACAP using Methods 21, 22, and 30
o VOACAP Method 30. A Long Path/Short Path Smoothing Function
7. HamCAP User's Guide
HamCAP v1.5 User's Guide
HamCAP and Type 13 & 14 Antennas
8. Articles of interest
English
The Woody Woodpecker Story, Part 3 [PDF, 115 kB] by V. K. Lehtoranta, OH2LX
Jamming, or deliberate interference against radio broadcasting stations [PDF, 140 kB] by V. K.
Lehtoranta, OH2LX
www.radiojamming.info by Rimantas Pleikys
Finnish
Kokeiluja Faroksella - automaattinen kelimajakoiden kuuluvuudenmittaus vs. keliennusteet
Tee-se-itse täsmäennuste YLE Radio Finlandin lähetyksille
VOACAP ennustaa kuuluvuuden signaali-kohinasuhteen avulla
VOACAP taajuussuunnittelun työkaluna
Avaruussään vaikutuksesta HF-keliin
Avaruussään tarkkailua IonoProbella
Swedish
VOACAP förutspår hörbarheten med hjälp av signal-noiseförhållandet
Japanese
HamCAP 1.5 User's Guide [PDF, 49 kB]
9. Tools & downloadable files for use with VOACAP
Find coordinates and Maidenhead grid locators just by pointing and clicking
Browse the World with Google Maps, (start view: World)
Browse the World with Google Maps, (start view: Scandinavia)
World Radio Broadcasting Stations on Google Maps! new
QTH Locator (Google Maps), updated!
Downloadables
Software
MAKEVOA.EXE for creating VOACAP-compatible antenna models from HFTA, courtesy of R. Dean
Straw N6BV (zipped, 34 kB)
Elevation-Angle Analyzer for VOACAP, v1.0
Windows program (zipped, 830 kB) for conversion of on-line CTY.DAT to CTY.GEO [9 Dec 07]. Mac
and Linux versions available upon request!
Perl source code for conversion of WL_CTY.DAT to CTY.GEO
Coordinate files
3
Save the coordinate files (*.GEO) at C:\ITSHFBC\GEOCITY\
New country coordinate list from CTY.DAT by AD1C. [9 Dec 07]
GEO file for all Zones (WAZ). Courtesy of Peter Ackermann DL3NAA.
Latest Global HF Transmitter Site List. [5 Aug 06]
Latest country coordinate list by OH6BG. [30 Jul 06]
Latest NCDXF beacon coordinate list [30 Jul 06]
Finland city coordinate list
Antenna and configuration files
A collection of antennas for use with VOACAP is available. If you wish to contribute, please let me
know.
Fine-tuned VOACAP Point-to-Point parameter auto scaling settings, English. Save as
C:\ITSHFBC\DATABASE\Voacap.gph
Fine-tuned VOACAP Point-to-Point parameter auto scaling settings, Finnish. Tallenna nimellä
C:\ITSHFBC\DATABASE\Voacap.gph
Fine-tuned Y-axis scaling for VOACAP Time & Distance graphs. Save as
C:\ITSHFBC\DATABASE\voadist.gph
Fine-tuned colors for auto-contour shading for Point-to-Point plots. Save as
C:\ITSHFBC\DATABASE\Colors.gph
4
2 VOACAP - Backgrounds & Running the Prediction
2.1 General
VOACAP Overview
VOACAP is an improved and corrected version of IONCAP, retaining all of the theory as put forth by John
Lloyd, George Haydon, Donald Lucas and Larry Teters in the 1975-1985 time-frame with modifications
which were suggested/approved by George Lane, Donald Lucas,George Haydon and A. D. Spaulding (a
world authority on HF radio noise predictions).
Take a look at the VOACAP evolution chart (courtesy of George Lane, Lane Consultant).
Major improvements in efficiency, coding corrections and ease of understanding the IONCAP program were
made by Franklin Rhoads of the U.S. Navy Research Laboratory under the sponsorship of the Voice of
America (1985-1996). Many of the newer features in VOACAP and VOAAREA were designed and
implemented by Gregory Hand at the Institute for Telecommunication Sciences who created VOAAREA and
made many significant improvements to VOACAP.
Screenshots of VOACAP
VOACAP (for point-to-point circuit analysis) and VOAAREA (for area coverage analysis) are extremely
powerful tools for professional-grade HF circuit predictions. Contrary to the common belief, the VOACAP user
interface is often simpler than that of many other HF prediction software, and the result graphs and coverage
maps usually outperform those of the competing tools. There is no competition, however, as VOACAP is free!
Below you will catch a glimpse of the many exciting features VOACAP can offer.
2.1.2.1 The Input Screen
This is the input screen. Please note that there is only one screen, it cannot be much simpler than that.
This is the signal-to-noise (SNR) graph as a function of 2-30 MHz and time, for a circuit from the United
Nations in NYC to Cape Verde in Africa. Besides this graph, there are 19 additional graphs that describe
various properties of this particular circuit, such as signal power at receiver, and the optimum radiation angle of
the transmitter antenna.
6
2.1.2.3 User-defined Frequencies vs. Time Graphs
This is the 24-hour analysis graph of the given output parameter on user-defined frequencies from the UN to
Cape Verde, Africa. The user can choose from 23 output parameters; shown here is the parameter SDBW, i.e.
median signal power at receiver site. By looking at the graph, the user can easily understand the behaviour of
the lower frequencies and higher frequencies at the different times of the day on the given circuit.
This is the distance graph of the given output parameter on user-defined frequencies from the UN towards Cape
Verde, Africa at 01 UTC. The distance graph shows how the user-defined frequencies start to propagate away
from the transmitter towards the receiver site. The user can choose from 23 output parameters; shown here is
again the parameter SDBW, i.e. median signal power at different distances from the transmitter.
7
2.1.2.5 Coverage Maps
Area coverage maps are perhaps the most exciting part of VOACAP. The coverage maps are made with the
VOACAP companion program, VOAAREA. Coverage maps are produced by calculating a huge amount of
point-to-point circuits from the transmitter to the geographical area defined by the user. A 31x31 matrix, as
shown here, means that the coverage map was created from 961 point-to-point circuit calculations. When the
geographical area is large, a better resolution may be needed. VOAAREA allows a matrix of 361x361 (130,321
P-to-P circuits) to be calculated at maximum.
This particular coverage map shows the median Maximum Usable Frequency from Finland to the rest of the
world in November 2005 at 10 UTC. A total of 23 output parameters can be plotted on the map.
8
Modelling the real-life HF communication system into VOACAP
2.1.3.1 Overview of VOACAP input and output parameters
VOACAP operates in two visible modes, an input and output mode. The following two figures below give an
overview of parameters that need to be taken into account while constructing an HF communication system into
a working VOACAP model.
Below you will see the main screen with a set of reasonable input values for a high-power broadcasting station
with transmissions to Southern Africa:
As noted earlier, while using VOACAP or VOAAREA, be sure to select the CCIR coefficients. George Lane
explains,
"Do the URSI coefficients provide more accurate predictions? We do not know. What we do know is that
original IONCAP using the CCIR data base, corrects the predicted values by using the Transmission Loss
tables and determines the variability of the MUF values over the month using the F-tables. These are huge
corrections. They are based on taking the predicted ionospheric calculations using the CCIR data and
subtracting them from actual system performance data taken over the 1958 to 1963 time period. The differences
are collected into the two sets of tables I just mentioned. The corrections are applied to any path prediction you
make using IONCAP in order to predict signal power, SNR, reliability and required power gain. If you change
the ionospheric maps for only the F-region and do not recalibrate the predictions against measurement, you do
not know if you made an improvement or not. Chances are you have ruined the correlation."
These are the basic settings which you can initially set up as follows:
10
Parameter Recommendation
Method Select Method Number 20 (Complete System Performance). If you are
running circuits over 7000 km, use Method 30 (Short/Long Path Smoothing).
Year This is a user-specified value. It has no effect on calculations.
Coefficients Select CCIR (Oslo). This is the default for using VOACAP in long-term
predictions. Do not use URSI88 (Australian).
Time Select UT (Universal Coordinated Time). Local Mean Time does not mean
your local time but the local sun time at the receiver.
Transmitter Select the transmitter location from the pre-defined lists, sorted by city
(worldwide), by nation (worldwide) and by state (USA). You can also enter
the coordinates manually.
Receiver Select the receive location from the pre-defined lists, sorted by city
(worldwide), by nation (worldwide) and by state (USA). You can also enter
the coordinates manually.
Path Select Short or Long Path. The default is Short Path.
Freq(MHz) Enter the frequencies to be used for calculation. Use the zero (0) value to
remove a frequency from the set. A set of 11 user-specified frequencies can
be defined. Also three pre-defined frequency sets are available which can be
modified if needed.
Fprob Do not use all default values. Set foE to 1, foF1 to 1, foF2 to 1 and foEs to 0
(zero, instead of 0.7).
Groups
11
Parameter Recommendation
Months Enter the number of the month (e.g. December = 12.00). You can specify a total
of 10 months at a time. The format is MM.DD (month.day). VOACAP is not,
however, suitable for daily predictions. Therefore, do not specify any day since
if you do, the value of Coefficients will automatically be set to URSI88.
SSNs Enter the monthly smoothed sunspot number for the given month.
There are a number of sunspot resources available on the Net. However, be sure to use the values from the
following resource:
These are the current values for the year 2003 as of writing:
JAN FEB MAR APR MAY JUNE JULY AUG SEPT OCT NOV DEC
82.9 79.4 76.9 74.3 71.2 68.7 66.7 64.8 62.6 59.8 57.8 56.7
System
12
Parameter Recommendation
Man-made This is a critical setting. You can accept the default value of 145
noise (residential environment). If your receive location is in a quiet
environment, you can use 155 or higher. There are still other values
available, see the input help.
Min. take-off If you are a radio amateur or SWL, enter the take-off value of 3. For better
angle antenna installations, you can use the value of 1.
Req'd circuit Accept the default value of 90%.
reliability
Req'd SNR This is a critical setting. For a reasonable BC listening quality, use 67.
For a reasonable CW reception quality, use 24 (or 27) and for SSB, 45.
Multipath Accept the default value of 3 dB.
power
tolerance
Max. tolerable Accept the default value of 0.10 milliseconds.
time delay
Antennas
Parameter Recommendation
TxAnt= Select the transmitter antenna from the pre-defined lists. You must familiarise
yourself with the available antenna patterns using HFant program, included in
the VOACAP software package. This window is the place where you can
13
define different antennas for different frequency ranges. In the example above,
we have defined one antenna for the frequency range of 2 MHz (Min) to 30
MHz (Max).
MainBeam Enter the main beam of the antenna (in degrees).
at Rx Sets the main beam towards the receive location.
TxPower The transmitter power which is fed to the antenna. Usually, we need to assume
some power losses in transmission lines, etc. so you may set this value, for
instance, to 65% of the transmitter power.
Parameter Recommendation
Receive Select the receive antenna from the pre-defined lists. You must familiarise
Antenna= yourself with the available antenna patterns using HFant program, included in
the VOACAP software package. The default antenna SWWHIP.VOA
represents a typical shortwave whip antenna.
Receiver If it is a directional antenna, enter the main beam in degrees. For non-
bearing directional antennas, enter the value of 0.
at Tx Sets the main beam towards the transmitter location.
Gain This is valid only for isotropic antennas. The gain value is ignored if used
with other antenna types.
14
Textual circuit prediction (Run > Circuit)
2.1.5.1 VOACAP output parameters
Below you will find brief definitions for the 22 parameters in the final column:
Parameter Definition
MODE Ionospheric path for the most reliable mode (ie. the mode with the highest
reliability of meeting the REQ.SNR). For the Short Path Model, the number of
hops and mode type [E, F1, F2, Es layers] for the MRM are given, eg. 2F2. For
the Long Path Model (not to be confused with the long path circuit!), the mode
at the transmitter end and the mode at the receiver end are given, eg. F2F2 in
our example above.
1 2 3 4 5 6 7
123456789012345678901234567890123456789012345678901234567890123456789012345
COMMENT Any VOACAP default cards may be placed in the file: VOACAP.DEF
LINEMAX 55 number of lines-per-page
COEFFS CCIR
TIME 1 24 1 1
MONTH 200112.00
SUNSPOT 110.
LABEL YLE PORI JOHANNESBURG
CIRCUIT 61.47N 21.58E 26.25S 28.00E S 0
SYSTEM 1. 145. 3.00 90. 67.010.00 0.05
FPROB 1.00 1.00 1.00 0.00
ANTENNA 1 1 2 30 0.000[hfcc\HFBC_218.P15 ]160.0 350.0000
ANTENNA 2 2 2 30 0.000[default\SWWHIP.VOA ] 0.0 0.0000
FREQUENCY 6.07 7.20 9.7011.8513.7015.3517.7321.6525.89 0.00 0.00
METHOD 30 0
EXECUTE
QUIT
1.0 13.1 6.1 7.2 9.7 11.9 13.7 15.4 17.7 21.6 25.9 0.0 0.0 FREQ
F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 - - MODE
4.0 7.6 7.8 3.3 4.0 4.0 7.8 7.8 7.8 7.8 - - TANGLE
12.0 7.6 7.8 3.3 12.0 12.0 22.0 4.0 8.0 7.8 - - RANGLE
34.4 34.2 34.2 33.9 34.3 34.6 35.2 34.6 35.5 34.9 - - DELAY
360 290 297 314 332 387 466 402 574 458 - - V HITE
0.50 1.00 0.99 0.95 0.73 0.38 0.12 0.01 0.00 0.00 - - MUFday
142 146 143 142 138 145 157 198 268 289 - - LOSS
44 35 39 48 47 41 29 -5 -77 -96 - - DBU
-87 -90 -87 -84 -83 -90 -102 -143 -212 -233 - - S DBW
-167 -153 -156 -162 -165 -168 -169 -171 -173 -176 - - N DBW
80 63 69 78 83 78 68 28 -39 -58 - - SNR
13 16 12 10 11 16 26 66 133 138 - - RPWRG
0.74 0.24 0.57 0.74 0.78 0.70 0.51 0.03 0.00 0.00 - - REL
0.00 0.00 0.00 0.64 0.00 0.00 0.00 0.00 0.00 0.00 - - MPROB
0.35 0.19 0.28 0.34 0.38 0.33 0.23 0.03 0.00 0.00 - - S PRB
25.0 9.7 11.2 19.2 25.0 25.0 25.0 25.0 25.0 9.1 - - SIG LW
17.7 4.9 4.9 6.2 11.6 21.5 25.0 25.0 25.0 4.9 - - SIG UP
26.7 12.4 13.8 21.2 26.7 26.8 26.8 26.8 26.8 13.3 - - SNR LW
18.5 7.6 7.1 7.8 12.7 22.2 25.7 25.7 25.7 7.6 - - SNR UP
19.7 21.9 21.9 18.3 19.7 19.7 21.8 21.8 21.8 21.7 - - TGAIN
-0.9 -2.1 -1.9 -7.0 -0.9 -0.9 0.0 -5.4 -1.8 -1.9 - - RGAIN
54 51 55 57 56 51 41 1 -66 -71 - - SNRxx
There are 24 blocks by hour in this example prediction. One hour block, as seen above out of 24, contains 14
columns. The first column is the time ("1.0") in UT (Universal Coordinated Time). Hours are centered on the
hour. Thus, "1.0" represents the time from 0030 to 0130 UT centered at 0100 UT.
The second column (the data under "13.1", ie. 13.1 MHz) is the predicted median maximum usable frequency
(MUF) at that hour. Then follows the 11 user-specified frequency columns (the frequencies are in MHz). In our
example, only 9 frequencies have been specified, so two last columns out of 11 are empty. The final column
16
contains 22 parameters that have been calculated for these frequencies in the hour block. The first parameter
FREQ indicates the frequency (in MHz) used in calculation.
Our data example is again from the circuit of Pori (Finland) - Johannesburg (South Africa). Choose Graph
from the Run menu (abbreviated later in these web pages as Run > Graph). A new window will open,
containing all the available output parameters:
The graphical charts provide a nice overview of the parameters on the frequency/time scale. To view the exact
details of the values, click anywhere on the chart, the details appear below the menu bar:
While in a parameter output window, another parameter can be chosen by clicking the Parameters menu. The
frequency scale (on the left side) can be changed from the Scale menu. If you wish the window to be pasted
into another program, click the to Clipboard menu. To print the window, choose from the Print menu.
Some of the most prominent parameters to look at are SNR, REL andSNRxx. Below are those parameters from
our circuit. The thick black line that runs through the chart is the median MUF.
17
2.1.6.1 SNR (aka SNR50)
The SNR values (in dB-Hz) that can be maintained 50% of days (15 days) in the month.
18
2.1.6.3 SNRxx (aka SNR90)
The SNR values (in dB-Hz) that can be maintained 90% of days (27 days) in the month.
From a point-to-point circuit analysis we can proceed to an area coverage analysis. It sometimes happens that a
frequency that is most suitable for a certain point-to-point circuit, is not the best possible for providing good
service to a specific area.
19
Let us take our example circuit of Pori (Finland) to Johannesburg (South Africa) at 1 UTC in December 2001.
We have had an analysis of this circuit in the article Understanding Above-the-MUF Predictions. It shows that
the best frequency would be 11.9 MHz. Now let us see how that frequency would work for the African
continent.
The following screen reflects the basic settings of VOAAREA - a companion program to VOACAP - needed
for running the coverage analysis:
We are interested in the circuit reliability in Africa, therefore we choose REL from Parameters. You may
choose several parameters at a time. As some locations in our target area are 10,000 km away, we choose
theMethod 30, Short/Long Path Model Smoothing. Be sure to select CCIR (Oslo) coefficients (not URSI88).
Also, keep the Fprob values as shown above. The following screen shots explain how to set some other input
values shown above:
2.1.7.1 Plot Center
2.1.7.2 Groups
20
2.1.7.3 System
When the calculations are done, choose Run > Plot results. A window will appear, choose File > Open
TEMP.VG1:
21
To plot the map, just choose Plot to > Window. You can also plot or print the map to the printer. To view the
remaining maps, choose File > Open. The file names are of format TEMP.VG? where the question mark is a
number between 1 to 7.
Furthermore, in this window you can choose several output parameters for viewing if you did not do so in the
main screen. Just select them by clicking the Parameters button.
If you wish to change the colours or values to be plotted, click theContours button. The following window will
appear:
Then click the REL button. The following window will appear:
22
You can change the values and colours here. Note that a value of -999 is used to fill unused lines. Be sure to
select a proper colour for that value too! If you change any of the default values, click the Values assigned
are:button. The text AUTO scaled will change to USER defined. Otherwise, your settings will be ignored.
23
2.2 The calculation engine, VOACAPW
An introduction to VOACAPW.EXE
The FREQUENCY card in VOACAPX.DAT
An introduction to VOACAPW.EXE
VOACAPW.EXE is the program module in the VOACAP program package that is responsible for VOACAP
propagation calculations. The graphical VOACAP user interface program is mainly responsible for collecting
the user input and putting it into a format that can be fed into the VOACAPW.EXE program. Further, the
graphical VOACAP suite of programs take care of presenting the output that comes as the result of
VOACAPW calculations.
VOACAPW.EXE can be found in the C:\itshfbc\bin_win\ directory (provided the user has accepted the default
VOACAP installation options). This program can be invoked in a variety of ways, depending on the type of
calculations we wish to run.
Below is the summary of the syntax that comes with the program:
|------------------------------------------------------------------------------|
| The propagation modules VOACAPW.EXE, ICEPACW.EXE, or REC533W.EXE cannot be |
| executed individually unless you have created the proper formated data input |
| file. You can execute: |
| POINTWIN.EXE (point-to-point interface) |
| AREAWIN.EXE (area coverage interface) |
| SI_WIN.EXE (signal-to-interference interface) |
| AREAINV.EXE (inverse area coverage interface) |
| HFANTWIN.EXE (antenna interface) |
| individually from this directory. |
|------------------------------------------------------------------------------|
| If you create your own input files in the proper formats, you can use the |
| examples below as to how to execute the various propagation modules. |
| The examples below are given for VOACAP. To make it run for ICEPAC or |
| REC533, just substitute the desired model for VOACAP |
|------------------------------------------------------------------------------|
| Execute with: |
| VOACAPW.EXE [S] directory inputfile outputfile (P-P circuit) |
| VOACAPW.EXE [S] directory inputfile outputfile (P-P graph) |
| VOACAPW.EXE [S] directory VOACAPD.DAT VOACAPD.OUT (P-P distance) |
| VOACAPW.EXE [S] directory VOACAPT.DAT VOACAPT.OUT (P-P time) |
| VOACAPW.EXE [S] directory BATCH (P-P batch) |
| VOACAPW.EXE [S] directory BATCH deckname.dat (P-P NEW batch) |
| VOACAPW.EXE [S] directory BATCH deckname.dat nam.out (P-P NEW batch) |
| VOACAPW.EXE [S] directory AREA CALC VOAAREAW.CIR (AREA batch) |
| VOACAPW.EXE [S] directory AREA method pathname (AREA single) |
| VOACAPW.EXE [S] directory INV CALC VOAAREAW.CIR (AREA INVERSE batch) |
| VOACAPW.EXE [S] directory INV method pathname (AREA INVERSE single) |
| VOACAPW.EXE [S] directory CIRAF pathname (CIRAF single) |
| where: |
| [S] = SILENT, then no messages are written to output window |
| [S] means either put in SILENT or nothing. |
| directory = full pathname to the install directory (e.g. c:\ITSHFBC) |
| method = (CALC/SCREEN/PRINT) |
| pathname = pathname below directory\AREADATA\ of input data file |
| (e.g. DEFAULT\DEFAULT.VOA) |
| deckname.dat = a proper formated input "card" deck. |
| |
| The input and output files ifor point-to-point runs will be in the |
| ..\RUN directory. |
| |
| You will probably need to also include the full pathname to the EXE file. |
| Thus, using the default installation directory, the P-P circuit would be: |
| |
24
| c:\itshfbc\bin_win\voacapw.exe c:\itshfbc voacapx.dat voacapx.out |
| |
| Note: Not all options will work for REC533 (e.g. distance, time, & INVERSE) |
|------------------------------------------------------------------------------|
| The only options you probably want to use are: |
| VOACAPW.EXE [S] directory inputfile outputfile (P-P circuit) |
| VOACAPW.EXE [S] directory AREA method pathname (AREA single) |
| VOACAPW.EXE [S] directory BATCH deckname.dat nam.out (P-P NEW batch) |
|------------------------------------------------------------------------------|
| |
| |
|------------------------------------------------------------------------------|
| To execute the world map plotting program from a command prompt: |
| |
| worldwin.exe run_directory data_directory filename noplot |
| where: |
| run_directory = RUN directory (e.g. C:\ITSHFBC\RUN) |
| data_directory= data sub-directory (e.g. AREADATA) |
| filename - file in ..\AREADATA\ directory of the |
| form subdir\????????.mGn (Grid file) |
| where:subdir=sub-directory under ..\AREADATA |
| m = V = VOACAP |
| m = I = ICEPAC |
| m = R = REC533 |
| m = C = COVMAP |
| n = group number of plot [1-9] |
| = 'voaareax.da?' means read for filename |
| = 'areapltx.dat' means read for filename |
| noplot = ' ' = do worldwin plot with wait |
| = 'n' = do worldwin plot with no wait |
| = 'b' = do worldwin plot with no wait,force B/W |
| = 'x' = do not do worldwin |
| = 'p' = do worldwin to Printer only |
| = 's' = do worldwin to Screen only |
| = 'm' = map only, no contours |
| Will create ????????.mxn (contour file) |
| where: x = index [1-9] of parameters contoured |
| |
| Thus, it can be exeuted with: |
| c:\itshfbc\bin_win\worldwin c:\itshfbc\run areadata default\default.ig1 |
| |
|------------------------------------------------------------------------------|
|------------------------------------------------------------------------------|
To run the calculations in silent mode, enter the following command in the C:\itshfbc\bin_win directory:
To run the calculations in normal mode (when the calculations appear in the output window), use the following
command:
As we are running VOACAPW.EXE from the command prompt, we will need to make sure that the input file
is in the format accepted by the program. Take a look at the following the example VOACAPX.DAT file
25
generated by the graphical VOACAP program after you choose the Circuit command from the Run menu for
the circuit of Pori (Finland) to Johannesburg (South Africa):
COMMENT Any VOACAP default cards may be placed in the file: VOACAP.DEF
LINEMAX 55 number of lines-per-page
COEFFS CCIR
TIME 1 24 1 1
MONTH 2002 7.00
SUNSPOT 96.
LABEL PORI JOHANNESBURG
CIRCUIT 61.48N 21.78E 26.25S 28.00E S 0
SYSTEM 1. 155. 3.00 90. 67.0 3.00 0.10
FPROB 1.00 1.00 1.00 0.00
ANTENNA 1 1 2 30 0.000[hfcc\HFBC_218.P15 ]175.0 350.0000
ANTENNA 2 2 2 30 0.000[default\SWWHIP.VOA ] 0.0 0.0000
FREQUENCY 6.07 7.20 9.7011.8513.7015.3517.7321.6525.89 0.00 0.00
METHOD 30 0
EXECUTE
QUIT
The definition lines for the entire year could be as follows (the year 2002 with the predicted smoothed sunspot
numbers):
MONTH 2002 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.0010.0011.0012.00
SUNSPOT 110. 109. 106. 103. 101. 98. 96. 93. 90. 87. 85. 81.
2. You can define a total of up to 20 antennas: ten transmit and ten receive antennas. If you were using the
graphical VOACAP program, you could define a total of five antennas, four transmit antennas and one receive
antenna. In our example above, we had one transmit antenna and one receive antenna:
Let us assume we have an amateur radio station with a good selection of antennas at our disposal. Our
definition of antennas for VOACAPW.EXE could be as follows:
26
(3) The minimum frequency for the antenna (in MHz).
(4) The maximum frequency for the antenna (in MHz).
(5) The design frequency. However, for isotropic antennas, it is the gain
(as in the case of the CCIR.000 antenna). The design frequency of zero
in our example means that the design frequency = operating frequency.
(6) The path to the antenna definition file under the C:\itshfbc\antennas
directory.
(7) The main beam of the antenna.
(8) For transmit antennas: the power (in kW). NOTE! You can define a
different power level for each of the transmit antennas!
For receive antennas: the gain (in dBi, applicable only to
isotropic antennas).
3. You can run several Methods in one go. Otherwise, you would run only one Method per one file. In our
example circuit, we defined Method 30 (Short/Long Path Smoothing):
METHOD 30 0
(1) (2)
Explanation:
If we wanted to run both Method 30 (Short/Long Path Smoothing) and Method 9 (HPF-MUF-FOT Graph) for
our circuit, we would write the following lines:
METHOD 30 0
EXECUTE
METHOD 9 0
EXECUTE
Now VOACAPW.EXE first executes Method 30 for the input data entered above, and then Method 9 for the
same input data. If you take a look at the result file (VOACAPX.OUT), you will see that the graphs for our
newly-specified 12 months follows the Method 30 output data.
If we implement all our modifications so far, the input file looks like this:
COMMENT Any VOACAP default cards may be placed in the file: VOACAP.DEF
LINEMAX 55 number of lines-per-page
COEFFS CCIR
TIME 1 24 1 1
MONTH 2002 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.0010.0011.0012.00
SUNSPOT 110. 109. 106. 103. 101. 98. 96. 93. 90. 87. 85. 81.
LABEL PORI JOHANNESBURG
CIRCUIT 61.48N 21.78E 26.25S 28.00E S 0
SYSTEM 1. 155. 3.00 90. 67.0 3.00 0.10
FPROB 1.00 1.00 1.00 0.00
ANTENNA 1 1 2 11 12.000[default\CCIR.000 ] 0.0 0.1000
ANTENNA 1 2 11 15 0.000[omat\3elhw.ANT ]175.0 0.1000
ANTENNA 1 3 15 22 0.000[omat\3elone.ANT ]175.0 0.1000
ANTENNA 1 4 22 24 0.000[omat\3eloq.ANT ]175.0 0.1000
ANTENNA 1 5 24 27 0.000[omat\3elo3q.ANT ]175.0 0.1000
ANTENNA 1 6 27 30 0.000[omat\3eltwo.ANT ]175.0 0.1000
ANTENNA 2 7 2 11 6.000[default\CCIR.000 ] 0.0 6.0000
ANTENNA 2 8 11 15 0.000[omat\3elhw.ANT ]357.0
ANTENNA 2 9 15 22 0.000[omat\3elone.ANT ]357.0
ANTENNA 2 10 22 30 0.000[omat\3eloq.ANT ]357.0
FREQUENCY 6.07 7.20 9.7011.8513.7015.3517.7321.6525.89 0.00 0.00
METHOD 30 0
EXECUTE
METHOD 9 0
EXECUTE
27
QUIT
Please note that running this file produces a 300-page result file!
COMMENT Any VOACAP default cards may be placed in the file: VOACAP.DEF
LINEMAX 55 number of lines-per-page
COEFFS CCIR
TIME 1 24 1 1
MONTH 2002 1.00 2.00 3.00 4.00 5.00 6.00 7.00 8.00 9.0010.0011.0012.00
SUNSPOT 110. 109. 106. 103. 101. 98. 96. 93. 90. 87. 85. 81.
LABEL PORI JOHANNESBURG
CIRCUIT 61.48N 21.78E 26.25S 28.00E S 0
SYSTEM 1. 155. 3.00 90. 67.0 3.00 0.10
FPROB 1.00 1.00 1.00 0.00
ANTENNA 1 1 2 11 12.000[default\CCIR.000 ] 0.0 0.1000
ANTENNA 1 2 11 15 0.000[omat\3elhw.ANT ]175.0 0.1000
ANTENNA 1 3 15 22 0.000[omat\3elone.ANT ]175.0 0.1000
ANTENNA 1 4 22 24 0.000[omat\3eloq.ANT ]175.0 0.1000
ANTENNA 1 5 24 27 0.000[omat\3elo3q.ANT ]175.0 0.1000
ANTENNA 1 6 27 30 0.000[omat\3eltwo.ANT ]175.0 0.1000
ANTENNA 2 7 2 11 6.000[default\CCIR.000 ] 0.0 6.0000
ANTENNA 2 8 11 15 0.000[omat\3elhw.ANT ]357.0
ANTENNA 2 9 15 22 0.000[omat\3elone.ANT ]357.0
ANTENNA 2 10 22 30 0.000[omat\3eloq.ANT ]357.0
FREQUENCY 6.07 7.20 9.7011.8513.7015.3517.7321.6525.89 0.00 0.00
METHOD 30 0
EXECUTE
FREQUENCY 2.00 3.70 7.0510.0714.1818.1021.2324.9528.30 0.00 0.00
EXECUTE
QUIT
Please note that while running different circuits in one go, you must use a fixed antenna configuration for all
the circuits! In our example we have chosen an omnidirectional 17-dBi antenna for testing purposes.
There are a few things you can do with the FREQUENCY card (line) that are not available when using the
VOACAP GUI for creating the data input file.
Below is the normal structure of VOACAPX.DAT (as created by the VOACAP GUI) when I want to know
about the path from Vaasa (Finland) to New York (USA) on 7.05, 10.12 and 14.05 MHz:
COMMENT Any VOACAP default cards may be placed in the file: VOACAP.DEF
LINEMAX 55 number of lines-per-page
COEFFS CCIR
TIME 1 24 1 1
MONTH 2003 7.00
SUNSPOT 63.
LABEL VAASA NEW YORK
CIRCUIT 63.10N 21.60E 40.72N 74.00W S 0
SYSTEM 1. 155. 3.00 90. 24.0 3.00 0.10
FPROB 1.00 1.00 1.00 0.00
ANTENNA 1 1 2 30 0.000[default\CCIR.000 ] 0.0 0.1000
ANTENNA 2 2 2 30 0.000[default\SWWHIP.VOA ] 0.0 0.0000
FREQUENCY 7.0510.1214.05 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00
METHOD 30 0
EXECUTE
QUIT
First of all, it is important to note that the placement of decimal points for frequencies on the FREQUENCY
line is crucial since they are in fixed positions. You can mess the frequencies up big time if you are not careful.
This is probably one of the reasons that there is an interface program that takes care of things like this.
Typically, the FREQUENCY line accepts up to 11 user defined frequencies. However, there are other ways of
using the FREQUENCY line.
29
2.2.2.1 Calculating from 2 to 30 MHz
If you want to calculate the prediction for the whole HF range on the given path from 2.0 MHz to 30.0 MHz in
1.0 MHz steps, you write as follows:
FREQUENCY -1
1.0 14.3 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 10.0 0.0 0.0 FREQ
2F2 4 E 4F2 3F2 3F2 3F2 3F2 3F2 3F2 2F2 - - MODE
7.7 3.0 21.5 13.8 11.0 10.3 10.0 10.1 10.4 3.1 - - TANGLE
22.5 21.3 23.9 22.8 22.4 22.3 22.2 22.2 22.3 21.8 - - DELAY
434 92 380 362 304 289 284 286 292 291 - - V HITE
0.50 1.00 1.00 1.00 1.00 1.00 1.00 0.99 0.94 0.99 - - MUFday
171 203 183 168 158 156 155 154 154 157 - - LOSS
-18 -61 -36 -22 -14 -7 -4 -6 -5 -1 - - DBU
-151 -181 -154 -141 -135 -133 -132 -132 -133 -132 - - S DBW
-171 -144 -149 -152 -154 -156 -158 -159 -161 -163 - - N DBW
20 -37 -5 11 19 23 26 27 28 30 - - SNR
30 81 48 32 24 19 17 17 17 16 - - RPWRG
0.42 0.00 0.00 0.08 0.27 0.46 0.55 0.58 0.60 0.64 - - REL
0.00 0.00 0.00 0.02 0.12 0.24 0.32 0.42 0.00 0.59 - - MPROB
0.20 0.00 0.04 0.11 0.17 0.21 0.24 0.27 0.28 0.27 - - S PRB
25.0 16.4 16.4 16.4 16.4 16.4 16.8 18.5 20.1 21.0 - - SIG LW
25.0 8.3 8.3 8.3 8.3 8.3 8.3 8.8 9.1 8.8 - - SIG UP
25.9 19.8 19.2 18.8 18.4 18.2 18.3 19.7 21.2 22.0 - - SNR LW
25.4 13.3 12.6 12.0 11.6 11.1 10.8 10.9 11.0 10.5 - - SNR UP
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 - - TGAIN
-2.0 -7.5 0.0 -0.5 -1.1 -1.2 -1.3 -1.3 -1.2 -7.4 - - RGAIN
-6 -57 -24 -8 0 5 7 7 7 8 - - SNRxx
The MUF is 14.3 MHz at 01 UTC, followed by frequencies starting from 2.0 MHz. The output file consists of
many pages...
As you can see, the first frequency on the line is zero. What that will do is the output for frequencies:
30
XMTR 2-30 + 0.0 dBi[default\CCIR.000 ] Az= 0.0 OFFaz=295.6 0.100kW
RCVR 2-30 2-D Table [default\SWWHIP.VOA ] Az= 0.0 OFFaz= 32.6
3 MHz NOISE = -155.0 dBW REQ. REL = 90% REQ. SNR = 24.0 dB
MULTIPATH POWER TOLERANCE = 3.0 dB MULTIPATH DELAY TOLERANCE = 0.100 ms
1.0 14.3 11.9 7.1 10.1 14.1 0.0 0.0 0.0 0.0 0.0 0.0 0.0 FREQ
2F2 2F2 3F2 2F2 2F2 - - - - - - - MODE
7.7 3.7 10.0 3.1 6.0 - - - - - - - TANGLE
22.5 21.9 22.2 21.8 22.2 - - - - - - - DELAY
434 310 284 291 379 - - - - - - - V HITE
0.50 0.90 1.00 0.99 0.56 - - - - - - - MUFday
171 156 155 157 164 - - - - - - - LOSS
-18 -1 -4 -1 -10 - - - - - - - DBU
-151 -136 -132 -132 -143 - - - - - - - S DBW
-171 -166 -158 -163 -170 - - - - - - - N DBW
20 30 26 30 27 - - - - - - - SNR
30 19 16 16 23 - - - - - - - RPWRG
0.42 0.62 0.55 0.65 0.56 - - - - - - - REL
0.00 0.00 0.33 0.59 0.00 - - - - - - - MPROB
0.20 0.29 0.24 0.28 0.26 - - - - - - - S PRB
25.0 24.1 16.8 21.1 25.0 - - - - - - - SIG LW
25.0 9.0 8.3 8.9 20.0 - - - - - - - SIG UP
25.9 24.9 18.3 22.0 25.9 - - - - - - - SNR LW
25.4 10.3 10.8 10.6 20.4 - - - - - - - SNR UP
0.0 0.0 0.0 0.0 0.0 - - - - - - - TGAIN
-2.0 -6.0 -1.3 -7.4 -3.2 - - - - - - - RGAIN
-6 5 8 8 1 - - - - - - - SNRxx
The MUF is 14.3 MHz, the FOT is 11.9 MHz and the rest of the frequencies are my 7.05, 10.12 and 14.05
MHz.
FREQUENCY 0
or
FREQUENCY 0 -1
1.0 14.3 2.0 3.7 5.5 5.5 7.2 8.9 10.7 12.4 14.1 15.9 17.6 FREQ
2F2 4 E 3F2 3F2 3F2 3F2 3F2 2F2 2F2 2F2 2F2 2F2 MODE
7.7 3.0 15.2 10.6 10.5 10.0 10.4 3.2 4.0 6.5 7.7 7.7 TANGLE
22.5 21.3 23.1 22.3 22.3 22.2 22.3 21.8 21.9 22.3 22.5 22.5 DELAY
434 92 394 295 294 284 292 296 319 396 434 434 V HITE
0.50 1.00 1.00 1.00 1.00 1.00 0.95 0.97 0.84 0.54 0.27 0.10 MUFday
171 203 173 157 157 155 154 157 156 165 187 215 LOSS
-18 -61 -25 -11 -10 -7 -5 -1 -2 -12 -34 -61 DBU
31
-151 -181 -144 -134 -134 -132 -133 -133 -136 -144 -167 -195 S DBW
-171 -144 -151 -155 -155 -158 -161 -164 -167 -170 -174 -177 N DBW
20 -37 7 21 21 26 28 31 31 26 7 -18 SNR
30 81 36 22 21 17 17 16 19 24 44 68 RPWRG
0.42 0.00 0.04 0.35 0.37 0.54 0.60 0.65 0.63 0.54 0.19 0.02 REL
0.00 0.00 0.00 0.16 0.17 0.34 0.00 0.58 0.00 0.00 0.00 0.00 MPROB
0.20 0.00 0.09 0.19 0.20 0.25 0.28 0.28 0.29 0.24 0.11 0.03 S PRB
25.0 16.4 16.4 16.4 16.4 17.6 20.0 21.5 25.0 25.0 25.0 25.0 SIG LW
25.0 8.3 8.3 8.3 8.3 8.4 9.1 9.7 9.7 21.8 25.0 25.0 SIG UP
25.9 19.8 18.9 18.3 18.3 19.0 21.0 22.4 25.8 25.9 26.1 26.3 SNR LW
25.4 13.3 12.2 11.4 11.3 10.8 10.9 11.2 10.9 22.2 25.3 25.3 SNR UP
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 TGAIN
-2.0 -7.5 -0.4 -1.2 -1.2 -1.3 -1.2 -7.1 -5.4 -2.8 -2.0 -2.0 RGAIN
-6 -57 -12 2 3 7 7 8 5 0 -20 -44 SNRxx
Here, the MUF is again 14.3 MHz, the first frequency to the right of the MUF is 2.0 MHz and the last is 17.6
MHz which is the HPF (MUFday = 0.10) at 01 UTC. The HPF is obviously different from hour to hour.
The remaining nine frequencies are "computed on a non-linear scale based on the E-MUF and F-MUF values".
This comes from the original IONCAP and can be elaborated by those more knowledgeable in the IONCAP
theory. :-)
Thanks to Greg Hand for explaining the FREQUENCY parameters.
2.3.1.1 Syntax
xfact, icode, num, [ con(num) ]
Where:
xfact = multiplying factor
The max value of num = 6 can only be used in the case of icode = 0.
32
** in reality, VOACAPW.EXE accepts 10 user-defined values. These need to be defined in the data input card.
Explanations:
Explanations:
Explanations:
The steps are based on a complicated algorithm, but it essentially finds 6 contours within the min and max. In
this case it would find <6,12,18,24,30,>30.
Explanations:
33
Contours start at this value and go up in steps of 10 [= CON(1)] for 6 contours. The upper limit is 200 [=
CON(2)].
icode=3 is opposite of icode=1, and is used for values that are really negative, but are calculated as positive
values (such as LOSS & RPWRG).
COLORS GPH Colors for auto contour shading for Point-to-Point plots
COLORS CON Colors for auto contour shading for Area Coverage plots
Colors.gph - original
8 7 1 2 3 15 Color indexes
1 4 6 56 58 59 Black/White shade codes
dark grey red gren yell whit
grey
For point-to-point graphs we could use the following color scale (originally suggested by Dan AC6LA
at http://www.qsl.net/ac6la/mapcontrol.html) [from highest to lowest value]:
For the color abbreviations, see file COLORS.WIN in the same DATABASE directory. Subsequently, you
would change the first line of COLORS.GPH to:
7 14 10 3 9 1 Color indexes
34
We will now change the Area Coverage colors accordingly:
Colors.con - original
15 7 6 14 2 3 9
whit grey cyan pale grn yell orng
blue
Let us use the color scale suggested by Dan AC6LA (from highest to lowest value):
Subsequently, you would change the two first lines of COLORS.CON to:
7 5 14 10 3 9 1
grey magn llbl lgrn yell orng red
Now these 7 colors are always automatically selected for Area Coverage plots, so no laborous manual color
definitions while plotting are needed from our part.
And while using the VOAAREA program, my settings under the LAYERS button are:
Grid = Grey
Countries = Blue (or Black)
Zones = Ignore
Cities = Ignore
MainBeam = Red
Contours = Shading, no contours
35
The Distance Plot calculates the variation of output parameter values with regards to distance at the user-
specified hour.
By default, most of the output parameter values are automatically scaled. There are, however, some parameters
which could benefit from a fixed Y-scale. The Y-axis scaling info is found inC:\itshfbc\database\voadist.gph.
2.3.3.1 Syntax
The variable names are:
xfact,icode,zmin,zmax,zlab,zinc
Example 1
SDBW = Median signal power at receiver (dBW)
Signal Power at Receiver (dBW)
1. 0 -160. -50. 10. 2.
Explanation:
Example 2
SNR = Median signal-to-noise ratio (dB)
Signal-to-Noise (Median) (dB)
1. 0 0. 100. 10. 2.
Explanation:
Set the minimum value to 0. This is, in fact, too low. A minimum value of 20 would do.
Set the maximum value to 100. For big BC stations, this can be too low a max value...
Set the number label at every 10th increment
Set the Y-axis tick mark at every 2nd increment
Example 3
SNRxx = Signal-to-Noise ratio (dB) at Req. Rel.
Signal-to-Noise at Req. Rel. (dB)
1. 0 0. 100. 10. 2
Explanation:
Set the minimum value to 0. This is, in fact, too low. A minimum value of 20 would do.
Set the maximum value to 100. Now, this can be too high a value. However, now at least, the SNR and SNRxx
scales are the same.
Set the number label at every 10th increment
Set the Y-axis tick mark at every 2nd increment
36
3 VOACAP - Understanding the Prediction
3.1 General
Ten Common Mistakes in Using VOACAP
This story is supposed to be taken in a light note although these are the potential risk areas that can turn
your prediction upside down :-)Some of the issues covered here seem to be rather commonplace according to
my searches using Google... however, it's another story how serious these all really are: some are very serious
and some are not. You may have compiled altogether different selections, and therefore, it would be nice to
hear about them :-) Enjoy!
1. The use of the URSI coefficients. - This seems to be a more common mistake than you might expect, and it
may even happen without your knowledge! Use the CCIR Coefficients. See also Mistake 5.
2. The choice of a wrong Smoothed Sunspot Number for the month. - The resource that maintains the
(predicted) SSN figures based on the Lincoln-McNish smoothing function to be used with VOACAP is located
presently at the address below.
These are the sunspot numbers used in the database reduction for the maps used in IONCAP and now
VOACAP.
ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUNSPOT_NUMBERS/sunspot.predict
I don't know where the hfradio.org site, for example, takes its numbers from but they are not from this resource.
3. The use of an unsuitable Required SNR (REQ.SNR) value for the transmission mode. - This, too, is
difficult for a number of people - and of course it varies case by case, depending on the grade of service you
wish to achieve. However, you can immediately see if there is something wrong with these values, say 73 for
CW ;-) Perhaps a common start value for CW is 24/27, 40 for SSB and 67 for BC. And, yes, the values are in
dB*Hz.
4. The use of a wrong transmitting power. - The transmitting power to use is the power at the antenna
feedpoint, not at the back of your transmitter :-) Therefore, a rule of thumb could be that you use 70% of the
transmitter output power as the power at the feedpoint. I gladly admit I break this rule on a regular basis ;-) And
you really have to pay attention and enter all values in kW, so 1 watt is 0.001 kW.
5. Entering a Day in the Month/SSN Groups window. - This is where you can make an unconscious mistake
(see Mistake 1). In Groups, there should only be the number of the month (e.g. "7.00" or "8.00" and so on) and
its corresponding SSN. If you enter a day within the month (say "7.29" where 7 is July and 29 is the day), you
will ruin your prediction :-) Using the day automatically forces VOACAP to use the URSI Coefficients for
calculations. And it does so without any warnings!
6. Setting a wrong Required Reliability (REQ.REL.) value. - This can be somewhat a matter of taste but I
would suggest always to use the value of 90%. Of course, if you know what you are doing, then it's OK :) I can
figure out at least one advantage with 90%: you can determine the best operating frequency just by looking for
the highest SNRxx value, right?
7. Setting too much/little man-made noise. - This is again difficult. The default value of 145 is a good starting
point. Here in Finland, if you go to the countryside, you can regularly enter 164 (remote) since the only noise
37
comes from within your receiver - or from your head ;-) In the cities, it is a different story. In the suburbs, I
guess 155 is, at least in this country, suitable.
8. Setting far too optimistic a Minimum Angle. - As discussed inanother article, the default value of 0.1
degrees is usually too good to be true in the context of amateur radio. There are rare exceptions, though. Use
3.0 degrees and it will be good. This is especially true if you use isotropic antennas; 0.1 degrees is absolutely
out of question!
9. Selecting a wrong antenna. - You have to investigate your antenna and its characteristics with the HFANT
software, no exceptions :-) You cannot be always certain that the azimuth and elevation patterns are as you
expect. Otherwise, you may choose an antenna which has nothing to do with the case in hand - so the
prediction will be in danger.
10. The use of a wrong Method. - Most of us needs only a few Methods. I use regularly Method 30,
sometimes Method 25 and Method 21. If I calculate FOT-MUF-HPF curves, then I use Method 9. A good
starting point is Method 30, really. For the other Methods, you will need to know what you are doing since only
part of the available Methods were tested and verified during the development of VOACAP. For instance,
Method 9 was not among them but Methods 30, 25 and 21 were.
George Lane: Think of an HF circuit as a telephone link. On one end a modulated signal is introduced which
contains information (voice, data, etc). The modulated signal is amplified by the transmitter and raised to RF.
The RF is then introduced to the transmit antenna which must efficiently accept this power and launch it into a
window in space.
This window is defined by the azimuthal spread of the great circle route (usually plus or minus about 2 degrees
for point-to-point circuits or by the width of the broadcast target area). The vertical height of the window is
determined by the ionospheric modes which will support that frequency to the intended target area. The trip
through the ionosphere attenuates the signal power until it is delivered to the receive location as a power flux
density (or as some people say, a field strength).
At this point we must consider the receive antenna pattern. This pattern should be such that it only accepts
waves arriving within the ionospheric window defined by azimuth angle and arrival angles of the ionospheric
modes (hops). However, the receive antenna accepts some energy at all angles. Any receive antenna over
perfect earth will have an average or integrated gain of 3 dB (it receives both the direct wave and also the
ground reflected wave - hence twice the power or 3 dB). RF noise power is arriving at the receive site from all
directions (atmospheric, man-made and galactic). This noise power is picked up by the receive antenna with 3
dB of gain in the lossless case. The receiver accepts both the signal and all of the noise power within the
bandpass filters of the receiver. The ability of the person or device to extract information from the receive
signal is dependent on the signal power being sufficiently above the noise floor such that intelligence is
detectable.
Field strength or signal power is only half of the story. A poorly designed receive antenna will deny reception
of a signal just as though the transmitter had been shut off. A receive antenna in a high noise environment will
mask the signal such that no intelligence can be detected.
Early work prior to World War II in broadcasting only used the field strength as a measure of communications
success. However, since the 1940's radio engineers have used Signal-to-Noise ratio as the preferred measure of
communications success. SNR is directly related to intelligence transfer. Also at great expense very accurate
RF noise prediction models have been developed for use in computing the SNR.
38
Choosing the correct Sunspot Number
I am a user of VOACAP or VOACAP-derived software (for instance HamCAP, MultiProp, WinCAP Wizard 3 or
ACE-HF). I do not know which of the available sunspot numbers (SSNs) I should use in any of these
programs mentioned. Many people and major HF propagation sites suggest using daily or monthly sunspot
numbers from the SEC (Space Environment Center,http://sec.noaa.gov/ftpdir/weekly/Predict.txt) or from
SIDC, Belgium (http://sidc.oma.be).
George Lane: In short, use the smoothed International Sunspot Number (SSN) values in VOACAP -- or in
any other interface program that uses VOACAP as their calculation engine. Smoothed SSNs can be found
from the National Geophysical Data Center
at:ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUNSPOT_NUMBERS/sunspot.predict. They maintain the
(predicted) SSN figures based on the Lincoln-McNish smoothing function. These are the sunspot numbers used
in the database reduction for the worldwide ionospheric maps used in IONCAP and now VOACAP.
The variation of "daily solar flux" over the days of the month is taken into account by VOACAP
through the use of the MUF-DAY tables. The trick is to get the circuit MUF predicted correctly using the
smoothed International Sunspot Number from the National Geophysical Data Center. This also requires using
the CCIR coefficients (same epoch as the F-days tables) and NOT the URSI-88 coefficients (a totally different
solar cycle for the foF2 mixed with the old foF1 and foE). Also, the day of the month should not be
specified, as this will force the use of the URSI-88 coefficients.
Method 26 will show the FOT-MUF-HPF values by hour for the month. That is where you can see the expected
variation of the MOFs due to daily variations in the solar flux. I wish we had changed the designation of FOT
to 90% of the days and the HPF to 10% of the days, as they represent the lower and upper decile range of the
expect MOF values over the days of the month at that hour.
To be more elaborate, the solar energy reaching the ionosphere is dependent on the number and size of the
sunspots. The effects take two to three days to substantially change the ionization in the ionosphere. The net
effect is that on a given HF link the Maximum Observed Frequency (MOF) for the path at a given hour will
change over the days of the month. An old rule of thumb was that the MOFs would vary around the
monthly median value (MUF) for that hour by ± 15%. In the early 1960's the scientists at ITS in Boulder
CO collected a fairly substantial data base over a full solar cycle during the IGY (1958-1963) period.
An investigation of the distribution of daily values of Maximum Observed Frequency (the MUF at a given hour
and day in the month) about their monthly median was carried out. Three points in the distribution over the
days of the month were considered; values of the daily MUF exceeded 0.90, 0.50 and 0.10. These points are
now referred to as the optimum working frequency (Frequence Optimum de Travail or FOT), the Maximum
Usable Frequency (MUF) and the Highest Probable Frequency (HPF). [Davis, R. M. and N. Groome (1964),
Variations of the 3000-km MUF in time and space (private communication)]
Data used in this study was derived from measurements at 13 stations representing a range of geomagnetic
latitudes from 71 degrees S to 88 degrees N. The variation in frequency about the MUF was represented in
tables of upper and lower decile values for low, medium and high Sunspot Number. Each table showed values
for a given season, local time in hours 00, 04, 08, etc., and each 10 degree of geographic latitude from 10
degrees to 80 degrees, north or south.
The study indicated that the distribution of MUFs is wider at night than in the daytime and wider at low
latitudes than high latitudes in the daytime. Again in daytime the distribution is wider in summer than winter,
except at high latitudes where the reverse is true.
The sunspot number dependence is weaker, but in daytime the difference between two ratios seems to increase
39
with sunspot number at latitudes higher than 40 degrees and to decrease with increasing sunspot number at
latitudes below 40 degrees. The distribution of frequency variation was mostly a function of the foF2 and not
the M(3000)F2 (the factor used to convert the vertical incidence critical frequency to the oblique path MUF);
therefore, the distributions are assumed valid for any oblique path. [Davis, R. and N. Groome (1965), The
Effect of Auroral Zone Absorption on High Frequency System Loss, (private communication)] [Lucas, D.L.
and G. W. Haydon, [1966] Predicting Statistical Performance Indexes for High Frequency Ionospheric
Telecommunications Systems, ESSA Technical Report IER 1-ITSA 1.]
In 1966, tables of the FOT-MUF-HPF distribution as a function of local time at the transmitter, season,
smoothed sunspot number (SSN) and geomagnetic latitude (of the transmit site) were incorporated in the first
widely used HF ionospheric radio performance prediction program, ITSA-1 [Lucas and Haydon, 1966]. These
are often referred to as the F-Days tables. The value of F-DAY was loosely defined as the probability of
ionospheric support (i.e. the fraction of the days that the operating frequency is below the MUF).
The development of IONCAP changed these definitions, as the F-DAYS values were only used for the lowest
order mode for a circuit-hour. This is also called the MUF mode and is very carefully computed using
convergence of up to 5 iterations. Modified distributions are assumed for the MUFs about each of the higher
order modes (up to 20) depending on the layer; Es, E, F1 and F2 [Lloyd, et al., 1978] [Lloyd, J. L., G. W.
Haydon, D. L. Lucas and L. R. Teters, [1978] Estimating the Performance of Telecommunication Systems
Using the Ionospheric Transmission Channel; Volume I: Techniques for Analyzing Ionospheric Effects Upon
HF Systems {DRAFT}, US Army CEEIA Technical Report EMEO-PED-79-7].
The authors of IONCAP often stated that F-DAYS are not used in IONCAP. However, in VOACAP the "F-
Day" factor given for each of the user defined operating frequencies was changed to show the value actually
being used for the most reliable mode (mrm) for that frequency. This new factor is named MUF-DAY in
VOACAP [Lane, 2001].
My paper at the Ionospheric Effects Symposium (Alexandria VA, May 3-5, 2005) happens to be about
using VOACAP for frequency assignment for ALE links. Rockwell Collins provided me with nearly 22 days
of data on a long path from France to Cedar Rapids (over 7,000 km). I had a real headache as the data was
taken in July 2004 and the International Smoothed Sunspot number kept changing each month until Feb 2005.
It takes 6 months for the smoothing function to settle in on a final value!
Collins had kept track of the frequency having the highest SNR at each hour for a scan of 20 frequencies. My
contention was that ALE systems should follow the circuit MUF as that is where the signal is the
strongest. A linear correlation analysis of the measured median best SNR frequency versus VOACAP
predicted MUF provided a correlation coefficient (r2) of 86%. That is amazingly high considering the fact that
the MUF mode changed from 2F2 during the day to 3F2 at night on that path.
I also looked at the correlation between the measured upper decile and lower decile of the frequencies being
used per hour with the predicted HPF and the FOT. There was good agreement but not as good as at the MUF.
Hope this is helpful in what seems like a never ending discussion of sunspot numbers. By the way, ICEPAC
was supposed to be a daily prediction program but it seems that work on the program ended long before it was
ready to be released. But that is the topic of my second paper.
George Lane: Poor George Haydon and John Lloyd (the creators of IONCAP) pulled their hair out when I said
IONCAP had to go down to 2 MHz. I knew there are real life situations where 2 MHz is a needed frequency
such as in northern Germany, Norway, Alaska, etc. in the winter time. Army radio links were always failing
and it wasn't until we got chirp sounders on those links that we found the MUF was going below 3 MHz at
40
night on short paths. George Haydon said there was very little data below 4 MHz but there was some for short
paths that did go down to 2 MHz. So they modeled a fit to those cases. Risky, but it has proven to give good
results for Near Vertical Incidence Skywave (NVIS) situations.
Radio amateurs and SWLs (who like to DX on 1.8 MHz band or pick up signals from very distant MW stations
at night) are finding VOACAP does not predict such contacts at 2 MHz whereas they know that they are
making such contacts. What is happening is that VOACAP is a ray hop model for normal HF skywave
propagation and it does not consider the 'wave guide' sort of propagation that occurs with the residual E layer
on all dark paths for the medium wave band.
John Wang at the FCC has just had his worldwide medium wave prediction model adopted by the ITU-R. I
think John Wang's model is excellent and it should be used for medium wave coverage predictions. I would not
try mixing it with VOACAP unless you run both models at 1.8 MHz and select the higher signal power.
VOACAP is quite good for NVIS and Wang's model is good for long distance MW paths in full darkness.
Since ITS delivered IONCAP to the Army in 1978, I have fought to keep IONCAP theory unmodified.
VOACAP is a cleaned up and corrected version of IONCAP and it retains all of the theory as put forth by
Lloyd, Haydon and Lucas. I don't recommend taking portions of Wang's model to be put into VOACAP. The
reason is that both programs are semi-empirical in that the theory has been forced to fit certain measured data
but not the same set. From a purely statistical point of view, one should not mix two dissimilar models without
redoing the correlation analysis.
George Lane: Let me first address the time discontinuity in the ITS programs VOACAP, ICEPAC and
REC533. There are several large global maps of ionospheric data used in these programs. That data was
collected during the IGY (the International Geophysical Year, late 50's to early 60's). The data was reduced to
tables which were valid for 2 or 3 month periods at the even hour. These maps are the critical frequency for the
E, F1, F2, M-3000, F-days distribution and the EXCESS SYSTEM LOSS TABLES (more about them later).
Of equal importance to the signal power calculation is the Noise Power. Atmospheric radio noise was mapped
in 3 month groups and 4 hour time blocks. The original data on an hourly basis, as collected, is lost.
I remember sitting in with the old Signal Corps Radio Propagation Agency engineers as they discussed whether
it is realistic to make computer predictions at each month from tables where the monthly dependence was lost
by averaging. I bet their eyes would pop out of their heads if they knew Johnny-come-latelies are now making
daily predictions with these averaged models! Certainly, you can force the programs to give a smooth transition
across the days of the month but it is just a trick which has little bearing on reality. If you have made HF
measurements for as many years as I have, you know that the day to day variations are not smooth across the
days of the month at a given hour. However, we can quite accurately compute the distribution of the variation
such that we can say only 3 days at this hour during the month should be this bad. But we have no idea when
they will occur during the month.
41
The fix was to create a table of differences at the median, upper and lower decile for predicted versus measured
signal power. This was done for one epoch of ionospheric data so that high and low sunspot data was included.
This difference table became the EXCESS SYSTEM LOSS TABLE (transmission loss table in IONCAP). It is
the fudge factor which brings the prediction signal power down to the level of the measured signal power as a
function of path length, geomagnetic latitude, hour, month and sunspot level.
John Lloyd, creator of IONCAP, points out very vividly that if you change any of the data base tables to a
different epoch, you must recompute the Excess System Loss correction maps. Since the actual data no longer
exists, you have no way of doing so. This is why you may get very poor predictions running VOACAP with the
URSI coefficients. It could happen that they are really more accurate but it would be a fluke. From a
statistician's point of view, you can have no confidence in the performance predictions when using the URSI
coefficients which are from a different epoch than the rest of the maps.
I can hear someone saying let's use CCIR Data Base D or D-1 or D-2. I attended meetings for the development
of Rec.533. I became intrigued when I discovered that they had subtracted 8 dB from the median values in
Excess System Loss tables. They only used the residuals of the table less the 8 dB. CCIR never did use any of
the decile tables, just the median values. Then they ran REC-533 and subtracted the signal powers from those
contained in Data Base D-1. They found the median difference was 9.2 dB. So the net difference was 9.2 - 8.0
or 1.2 dB! Rather amazing since they only use the foF2 maps and compute F-1 and E layer from them.
Someone noticed that some of the circuits in DataBase D-1 were greatly different than predictions. So they
started removing these offending circuits from the data base. I became very curious because they always
seemed to fall in the region of the 1F2 to 2F2 transition. In talking with the German who had processed most of
this data, I learned that he had removed the actual antenna pattern gain from the measured signal power by
assuming a fix layer height for the F2 layer. Therefore, the normalized data in D-1 was extremely dependent on
the nomogram used to assess the takeoff and arrival angle of the signal. This meant that CCIR was forcing
Rec.533 to fit the nomogram with a fixed layer height!
About the time I retired, CCIR was attempting to 'de-normalize' D-1. Anyhow I have never been too keen on
Rec. 533. I especially dislike the use of median values and then applying some variation tables at the end of the
prediction. PC's have enough memory and are fast enough today that you don't need all of the simplifications
which were needed in the 1980's when most of the Rec.533 development work was done. The intent for 'little
252', HFBC-84 and then Recommendation 533 was to develop a very fast running computer model which could
analyze the full International Broadcast frequency schedule on a seasonal basis using inexpensive desk top
computers. As such, the development was very forward thinking and successful. There was no way of knowing
at that time how fast computer technology would overtake events.
The reason VOACAP was more 'correct' than ICEPAC is that Method 30 in VOACAP uses a smoothing
function between the short path (ray hop) and long (forward scatter) models as a function of distance between
7000 and 10,000 km. ICEPAC uses the abrupt transition between short path and long path at 10,000 km. It
turns out that in the ray hop model on that path between Germany and the USA finds an ionospheric control
point with a tremendous signal loss whereas the long path model computes a weak but detectable signal
42
assuming a scatter propagation mechanism. Therefore, ICEPAC found no signal and VOACAP found a weak
scatter signal which was actually there.
The following refer to each of the control points (in this case 5 of them). The number of control points varies
with the distance between TX & RX. "k" will be the control point.
The Sporadic E values in this example are 0.0 because the FPROB variable ("0.0" for Es) removes the sporadic
E layer.
What is the term "semi-thickness"? I see it being used in the explanations above.
43
George Lane: It is an old term dating back to the 1940's. Back then they used a set of templates which they
would slide over the ionogram based on an assumed parabolic shape to the ionospheric layer of interest. They
would find the curve which had the best fit when the template was aligned with the height of the maximum
electron density for the layer. The parabolic curve with the best fit was numbered with the semi-thickness (half
thickness) of the parabola. It was a term needed to determine the actual ray path as it passed through the various
layers.
The ionogram plots show electron density as a function of height. It was recognized early on that these plots
tended to show the layers as having a parabolic shape. If you consult a geometry text, you will find that there is
a locus point for a parabola which is centered on the line where the tangent to the line is vertical. This is the
height of the maximum electron density. Each parabola has a width factor associated with the locus point. This
is where the line tangent becomes horizontal. By curve-fitting to the ionograms, the scaler was able to find the
semi-thickness of the layer based on simple geometry. These parameters are needed in order to determine the
amount of bending which occurs as the ray at a particular frequency penetrates the 3 layers. Newer prediction
models such as PropLab actually compute the ray path for a particular electron density profile. IONCAP was
the first program to attempt this in a quasi-ray trace approach.
In other words, the MUF is the frequency for which ionospheric support is predicted on 50% of the days of the
month, ie. 15 days out of 30 days. So on a given day communications may or may not succeed on the frequency
marked as the MUF.
To ensure a good communication link between two locations, the operating frequency is typically chosen below
the predicted MUF. It is often claimed that the optimal operating frequency lies somewhere between 80-90% of
the MUF (e.g. if the MUF is 10 MHz, the optimal frequency would be around 8-9 MHz). However, in
VOACAP it is the predicted SNR distribution using Complete System Performance methods (e.g. Methods 20,
21, 22 or 30) that determines which frequencies provide an acceptable grade of service.
44
the user-specified frequencies are always the most reliable mode. For a more detailed discussion,
seeCalculating MUFdays.
Below are the four SNR output parameters needed for analysis:
1.0 13.1 6.1 7.2 9.7 11.9 13.7 15.4 17.7 21.6 25.9 0.0 0.0 FREQ
F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 - - MODE
80 63 69 78 83 78 68 28 -39 -58 - - SNR
26.7 12.4 13.8 21.2 26.7 26.8 26.8 26.8 26.8 13.3 - - SNR LW
18.5 7.6 7.1 7.8 12.7 22.2 25.7 25.7 25.7 7.6 - - SNR UP
54 51 55 57 56 51 41 1 -66 -71 - - SNRxx
The SNR indicates the dB-Hz value that can be maintained on 50% of the days (ie. on 15 days) in the month.
In our example above on 11.9 MHz, the SNR value is 83 (dB-Hz).
The SNRxx (ie. SNR90, provided the REQ.REL. is 90%) indicates the dB-Hz value that can be maintained on
90% of the days (ie. on 27 days) in the month. In our example above on 11.9 MHz, the SNRxx value is 56 (dB-
Hz). This can be calculated as SNR - SNR LW (or 83 - 27 = 56 in our example).
And finally, the SNR10 (calculated as SNR + SNR UP) is the dB-Hz value that can be maintained on 10% of
the days (ie. on 3 days) in the month. In our example above on 11.9 MHz, the SNR10 value is appr. 96 (dB-
Hz).
The two most prominent parameters to consider in search of the best frequency are the SNR and SNR90
values. As a rule of thumb, look for the highest SNR value and the highest SNR90 value. Let us assume
that the required SNR we wish to maintain in our circuit is 67 (not a good but still a reasonable listening quality
in international broadcasting). We will see that the SNRxx is below 67 at all our frequencies which means none
of them cannot maintain that grade of service on 27 days out of 30 days. Then we will need to look for the
highest SNR. Of our frequencies, the best would be 11.9 MHz with the SNR value of 83.
Conclusion
In conclusion, 11.9 MHz is the best candidate for the operating frequency at 01 UTC
during that month. 11.9 MHz is also below the predicted MUF of 13.1 MHz for that mode.
1.0 13.1 6.1 7.2 9.7 11.9 13.7 15.4 17.7 21.6 25.9 0.0 0.0 FREQ
F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 - - MODE
80 63 69 78 83 78 68 28 -39 -58 - - SNR
13 16 12 10 11 16 26 66 133 138 - - RPWRG
0.74 0.24 0.57 0.74 0.78 0.70 0.51 0.03 0.00 0.00 - - REL
26.7 12.4 13.8 21.2 26.7 26.8 26.8 26.8 26.8 13.3 - - SNR LW
18.5 7.6 7.1 7.8 12.7 22.2 25.7 25.7 25.7 7.6 - - SNR UP
54 51 55 57 56 51 41 1 -66 -71 - - SNRxx
45
The RPWRG is related to the SNR90 and REQ.SNR. In our example above, the REQ.SNR was set to 67.
The RPWRG (the required power gain) parameter tells us how many desibels are needed in the
communication system to achieve the SNR90 value of 67. It is calculated as REQ.SNR - SNRxx (or 67 - 56 =
11 on 11.9 MHz). As the value of the RPWRG is positive in our example, it means that many desibels are
needed for our system; if the value had been negative, that many desibels would have been in excess (ie.
unnecessary) to achieve the required SNR for 27 days out of 30 days.
This parameter relates to the (communication) system design. In our example on 11.9 MHz, we should consider
what measures we could take to add the necessary 11 desibels to the system: doubling the transmitting power
would give us 3 desibels, using a more powerful transmitter antenna could give us a few desibels more, and at
the receiving end we could choose, say, a 3-element Yagi instead of the whip antenna which would still
contribute some more desibels.
The REL is related to the SNR and REQ.SNR, and is defined as a circuit reliability factor. It tells us the
percentage of days in the month when the SNR value will equal to or exceed the REQ.SNR. The SNRxx tells
us which SNR value can be achieved on 90% of the days (27 days) in the month. If the SNRxx would have
been 67, then the value of REL had been 0.90 (or 90%, which is the REQ.REL. we have specified) and the
RPWRG would have been zero (0).
Conclusion
The REL value of 0.78 on 11.9 MHz suggests that the required SNR of 67 can be achieved
on 78% of days in the month. To translate the percentage value to the number of days, take
a look at the Z Tables. We will see that 78% equals to 23 days.
For a more detailed discussion on how the REL value is calculated, see the article Maintaining the Required
Grade of Service.
The MUF is the middle value of the daily variations of MOF at that hour so it might only occur once in the
month! We really need to think about the MOF (hourly median maximum observed frequency on a single day).
For frequencies below the MOF we have spectral reflection of the wave from the ionospheric layer. At the
MOF itself, the low ray and the high ray add for 3 dB enhancement (this is called the junction frequency). As
we go above the MOF, the ray penetrates deeper into the ionosphere. Here in the upper regions of the
ionosphere there exists blobs of ionized gases. As the ray passes through a blob, it is slightly bent. After
passing through enough blobs it is reflected back toward the earth. This is called the Above-the-MUF mode.
Just above the MOF the loss is not too much more than spectral reflection. But as the operating frequency is
higher than the MOF, the losses go up exponentially as there are fewer and fewer blobs with sufficient electron
density to bend the ray.
Eventually as the frequency goes higher and higher all of the energy penetrates the ionosphere and is lost. In
VOACAP the Above-the-MUF loss is limited to 25 dB. Personally, I think it is too low and probably should be
allowed to go to 40 to 50 dB.
46
So what is VOACAP doing? It is computing the distribution of the MOFs over the days of the month.
MUFDAY tells you the fraction of the days for which the operating frequency is below the MUF. For those
days which are Above-the-MUF, Above-the-MUF losses are included. You will see the SNR LW getting larger
and larger as the operating frequency approaches the MUF. As you go above the MUF, SNR LW approaches
27 dB which is a function of that 25 dB limit, I mentioned before.
If you only need a few days per month to meet your REQ SNR, then you may well be able to work frequencies
at or above the MUF. But if you require REL = 90%, then SNR LW comes into play and the SNRxx is much
lower.
1.0 13.1 6.1 7.2 9.7 11.9 13.7 15.4 17.7 21.6 25.9 0.0 0.0 FREQ
F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 - - MODE
4.0 7.6 7.8 3.3 4.0 4.0 7.8 7.8 7.8 7.8 - - TANGLE
12.0 7.6 7.8 3.3 12.0 12.0 22.0 4.0 8.0 7.8 - - RANGLE
34.4 34.2 34.2 33.9 34.3 34.6 35.2 34.6 35.5 34.9 - - DELAY
360 290 297 314 332 387 466 402 574 458 - - V HITE
0.50 1.00 0.99 0.95 0.73 0.38 0.12 0.01 0.00 0.00 - - MUFday
142 146 143 142 138 145 157 198 268 289 - - LOSS
44 35 39 48 47 41 29 -5 -77 -96 - - DBU
-87 -90 -87 -84 -83 -90 -102 -143 -212 -233 - - S DBW
-167 -153 -156 -162 -165 -168 -169 -171 -173 -176 - - N DBW
80 63 69 78 83 78 68 28 -39 -58 - - SNR
13 16 12 10 11 16 26 66 133 138 - - RPWRG
0.74 0.24 0.57 0.74 0.78 0.70 0.51 0.03 0.00 0.00 - - REL
0.00 0.00 0.00 0.64 0.00 0.00 0.00 0.00 0.00 0.00 - - MPROB
0.35 0.19 0.28 0.34 0.38 0.33 0.23 0.03 0.00 0.00 - - S PRB
25.0 9.7 11.2 19.2 25.0 25.0 25.0 25.0 25.0 9.1 - - SIG LW
17.7 4.9 4.9 6.2 11.6 21.5 25.0 25.0 25.0 4.9 - - SIG UP
26.7 12.4 13.8 21.2 26.7 26.8 26.8 26.8 26.8 13.3 - - SNR LW
18.5 7.6 7.1 7.8 12.7 22.2 25.7 25.7 25.7 7.6 - - SNR UP
19.7 21.9 21.9 18.3 19.7 19.7 21.8 21.8 21.8 21.7 - - TGAIN
-0.9 -2.1 -1.9 -7.0 -0.9 -0.9 0.0 -5.4 -1.8 -1.9 - - RGAIN
54 51 55 57 56 51 41 1 -66 -71 - - SNRxx
George Lane: The MUF = 13.1 MHz at that hour. So a broadcast frequency of 13.7 is only a little bit above the
MUF. At least 38% of the days will be below the MUF (MUFday = 0.38) and reception will be great. On 62%
of the days reception will be not-so-good to terrible... So on almost half of the days of the month propagation
will be just great. That is why the median SNR should be high.
So then, how about 11.9 MHz? 13.7 MHz is only 70% reliable (REL = 0.70) whereas 11.9 is 78% reliable
(REL = 0.78). At the median, the Signal-to-noise ratio (SNR50) at 11.9 MHz is 5 dB higher than at 13.7 MHz.
Clearly there is about 1 S-unit advantage to using 11.9 MHz band over that of the 13.7 MHz band on at least
half of the days of the month. Also on the bad days of the month, 11.9 will work better than 13.7 MHz so that is
another advantage.
I fear at 15.4 MHz, you are getting into the noise of the program. By that I mean, VOACAP will give
predictions even when it has no idea what is going to happen. One would normally stay away from using this
frequency for this reason:
47
SNR10 = 94 dB-Hz
SNR50 = 68 dB-Hz
SNR90 = 41 dB-Hz
If we believe the prediction, then VOACAP is saying that 80% of the days of the month (bounded by SNR90
and SNR10) will have a SNR somewhere between 41 and 94 dB-Hz. That is a spread of 53 dB or a factor of
200,000! If the program could talk, it would tell you that it doesn't have any idea what is going to happen on
15.4 MHz. At best you can say there is a 12% (MUFday = 0.12) chance of good service and 88% chance it will
be less than good.
There is one more variable that one should look at when dealing with frequencies way above the MUF. That is
the Signal Power at the receiver, S PWR (aka S DBW). Often you will find that the power is actually at or
below the receiver threshold. In that case the receiver thermal noise becomes controlling and the grade of
service is much less than the SNR would indicate. In our example, 15.4 MHz has a predicted S PWR of -102
dBW. If we assume the receiver input impedance is 50 Ohms and that the transmission lineloss is minimal, then
the median voltage at the input of the receiver is -102 + 137 or 35 dB relative to 1 microvolt, or 56 uV (see [1]
below; see also: S DBW to S Meter Correlation Table). For a good receiver this is an adequate input voltage.
But if we impose a desired reliability of 90%, the signal power reaching the receiver is -102 dBW - 25 dB (SIG
LW) = -127 dBW. Using the previous logic, then the voltage at the receiver under ideal conditions is only 3 uV.
For most receive systems which have antenna efficiency, transmission line and impedance mismatch losses,
this signal voltage does not overcome thermal noise even in a good commercial receiver.
If you look at S PWR at 17.7 MHz, you will find that the voltage at the input of the receiver has dropped to 0.5
uV which would require a very expensive receiver with good sensitivity and little impedance mismatch loss
between the antenna and the receiver.
Conclusions
VOACAP makes the assumption that the external RF noise power is controlling. At high
frequencies where the signal power and the noise power become very low, this assumption may
fail. If we only look at the SNR, we may be fooled into thinking the grade of service will be good
when in fact the receiver noise is greater than the signal plus external noise power.
The authors of IONCAP (the parent program for VOACAP) always warned us that we
should only believe the predictions when the reliability was approximately 90% or higher.
In your example, if I were the frequency planner for that broadcast, I would select 11.9
MHz as the best frequency. But I still would expect that on some days transmissions on
15.4 MHz will come booming in. The chances that transmissions on 17.7 MHz will be
heard on any day of the month depends greatly on the receiver and its antenna.
So the conclusion is that reception on frequencies which are above the predicted MUF are
certainly possible and in some cases reception will be possible on more than 15 days of the
month at that frequency.
3.2.3.1 [1] Equations to Calculate The Median Voltage at the Input of the
Receiver
Ohm's Law:
E = I*R or I = E/R
Power formula:
P = EI
48
P is in Watts, E is in Volts and I is in Amps and R is in Ohms.
Then:
P = E * [E/R] = (E^2)/R
Convert to dB:
10 Log E^2 = 20 Log E = 10 Log [P*R]
Assume R = 50 Ohms:
E(dBuV) - 120 = 10 Log P + 10 Log (50 Ohms) = P(dBW) + 17
E(dBuV) = P(dBW) + 17 + 120 = P(dBW) + 137
Note
VOACAP computes the signal power [S DBW] which is equivalent to the P(dBW) in the above equations. DO
NOT confuse the predicted field strength in the VOACAP output, DBU (dB relative to 1 microVolt per meter), with
the voltage at the receiver as given by dBuV (dB relative to 1 microVolt) in the above equations.
First we need to run Method 9, HPF-MUF-FOT graph, for our circuit. We are interested in the frequencies at
01 UT (= GMT). Below are the MUF, FOT and HPF values from Method 9:
The HPF, FOT and MUF represent the Highest Possible Frequency, Optimum Working Frequency (or
Frequency of Optimum Traffic) and Median Maximum Usable Frequency, respectively, for the hourly MOF
(median maximum observed frequency) distribution. At the HPF, communication can be supported on 10% of
the days (3 days) in the month; at the FOT, 90% of the days (27 days) in the month and, as discussed above, at
the MUF, 50% of the days (15 days) in the month.
A Rule of Thumb
The practical meaning of FOT and HPF [As obtained from Method 9 in VOACAP] is that they determine the frequency
range within which we should find80% of the measured MOFs.
Let us choose two frequencies 9.7 and 13.7 MHz for our calculations. 9.7 MHz is below the predicted MUF so
we choose the formula 1:
Then for 13.7 MHz, above the predicted MUF, we choose the formula 2:
To confirm these values, below you see that VOACAP has calculated the same values for those two frequencies
(0.95 = 95% and 0.38 = 38%):
1.0 13.1 6.1 7.2 9.7 11.9 13.7 15.4 17.7 21.6 25.9 0.0 0.0 FREQ
F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 F2F2 - - MODE
0.50 1.00 0.99 0.95 0.73 0.38 0.12 0.01 0.00 0.00 - - MUFday
50
Maintaining the Required Grade of Service: Calculating the
Circuit Reliability for a Given Hour
Let us assume you wish to know on how many days in a month a certain grade of service can be
maintained on a given hour between two locations. The grade of service is expressed in VOACAP by the
parameter REQ.SNR (the required signal-to-noise ratio in dB-Hz). This is a critical setting that must be
considered carefully. The correct value to select is dependent on the transmission mode, and the reception
quality we like to maintain. As a rule of thumb, a reasonable listening quality in the AM mode can be
maintained by using the REQ.SNR value of 67, for CW the REQ.SNR can be, say, 27 and for SSB, 45.
Remember always to set REQ.REL. to 90%.
When your required SNR (REQ.SNR) is equal to or less than the predicted SNR, use the following:
When your required SNR (REQ.SNR) is greater than the predicted SNR, use the following:
SWR uses a 100-watt transmitter with a horizontal half-wave dipole (beamed at 150 degrees) at 8 meters above
the ground. The receive antenna in Copenhagen is the default shortwave whip antenna, SWWHIP.
This is what we get as a result from VOACAP by running the circuit (Method 20):
18.0 6.8 6.2 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 FREQ
1F2 1F2 - - - - - - - - - - MODE
35.9 29.6 - - - - - - - - - - TANGLE
4.2 3.9 - - - - - - - - - - DELAY
392 308 - - - - - - - - - - V HITE
0.50 0.76 - - - - - - - - - - MUFday
123 119 - - - - - - - - - - LOSS
20 23 - - - - - - - - - - DBU
-104 -101 - - - - - - - - - - S DBW
-158 -157 - - - - - - - - - - N DBW
54 56 - - - - - - - - - - SNR
36 30 - - - - - - - - - - RPWRG
51
0.06 0.04 - - - - - - - - - - REL
0.00 0.00 - - - - - - - - - - MPROB
0.11 0.12 - - - - - - - - - - S PRB
20.3 17.1 - - - - - - - - - - SIG LW
9.4 5.4 - - - - - - - - - - SIG UP
22.2 19.4 - - - - - - - - - - SNR LW
10.9 7.8 - - - - - - - - - - SNR UP
-2.5 -4.6 - - - - - - - - - - TGAIN
-0.7 -0.3 - - - - - - - - - - RGAIN
31 37 - - - - - - - - - - SNRxx
Please note that for the circuit calculation we have assumed some power losses so that only 70% of the
transmitting power goes to the antenna. That explains the power of 0.070kW in the above example. The same
"70% rule" applies to our second example below.
As we can see, the required SNR (67) is greater than the predicted SNR (56) for that hour. The poor
performance of this circuit is mainly due to the low transmitting power. Furthermore, the broadcast frequency
(6.17 MHz, or 6.2 as given above) is too close to the MUF (6.8 MHz). The MUF is defined as the median
maximum usable frequency for a given path, month, SSN and hour. On each day of the month at this hour,
there is a maximum observed frequency (MOF) for this mode (1F2). The median of this distribution is called
the MUF for that mode. In other words, the MUF is the frequency for which ionospheric support is predicted on
50% of the days of the month.
The required SNR is greater than the predicted SNR, so we need to use the formula 2:
Now, look up the z value of 1.805 in Table 2. The percentage value will be around 3.6% which corresponds to
just about 1 day. By the way, if the power to the antenna would have been 100 watts instead of 70 watts, we
could have got 7% (2 days) as the result. Looking back now to the circuit results, you will see that the value
of REL, the Circuit Reliability Factor, is 0.04, ie. 4%. This is the same answer VOACAP gives us without
any laborous calculations from our part!
The Answer
If the required SNR is 67, then the circuit reliability is predicted to be about 3.6% (REL 0.04 ie. 4%). That means
you can maintain a reasonable grade of broadcasting service between Virrat and Copenhagen only on one day out
of 30 days in December 2001 on the frequency of 6170 kHz at 1800 UTC.
52
2. Case YLE Radio Finland
In December 2001, YLE Radio Finland transmits from Pori (Finland) to Africa on 15.520 MHz with a
transmitter power of 500 kW between 17-18 UTC. Calculate the number of days in a month to achieve a 67-dB
circuit reliability at these hours in Johannesburg (South Africa).
This circuit is almost 10,000 kms so we are using Method 30, Short/Long Path Smoothing:
17.0 20.4 15.5 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 FREQ
F2F2 F2F2 - - - - - - - - - - MODE
7.0 7.2 - - - - - - - - - - TANGLE
14.0 7.2 - - - - - - - - - - RANGLE
35.0 34.1 - - - - - - - - - - DELAY
462 281 - - - - - - - - - - V HITE
0.50 0.87 - - - - - - - - - - MUFday
143 140 - - - - - - - - - - LOSS
47 49 - - - - - - - - - - DBU
-87 -85 - - - - - - - - - - S DBW
-172 -167 - - - - - - - - - - N DBW
85 83 - - - - - - - - - - SNR
9 2 - - - - - - - - - - RPWRG
0.81 0.88 - - - - - - - - - - REL
0.00 0.00 - - - - - - - - - - MPROB
0.40 0.47 - - - - - - - - - - S PRB
25.0 15.1 - - - - - - - - - - SIG LW
20.4 5.5 - - - - - - - - - - SIG UP
26.7 17.5 - - - - - - - - - - SNR LW
21.2 7.2 - - - - - - - - - - SNR UP
21.8 21.9 - - - - - - - - - - TGAIN
-0.5 -2.4 - - - - - - - - - - RGAIN
58 65 - - - - - - - - - - SNRxx
Now the REQ.SNR is lower than the predicted SNR (83) so we'll apply
The z value of 1.17 (now from Table 1) translates approximately to 88%, which corresponds to 26 days. Again,
we will note that VOACAP has already calculated the same value for us in REL (0.88 = 88%).
18.0 17.4 15.5 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 FREQ
F2F2 F2F2 - - - - - - - - - - MODE
4.0 4.0 - - - - - - - - - - TANGLE
12.0 12.0 - - - - - - - - - - RANGLE
34.4 34.2 - - - - - - - - - - DELAY
350 317 - - - - - - - - - - V HITE
0.50 0.84 - - - - - - - - - - MUFday
143 137 - - - - - - - - - - LOSS
45 51 - - - - - - - - - - DBU
-87 -81 - - - - - - - - - - S DBW
-170 -168 - - - - - - - - - - N DBW
83 87 - - - - - - - - - - SNR
11 7 - - - - - - - - - - RPWRG
0.77 0.83 - - - - - - - - - - REL
53
0.00 0.00 - - - - - - - - - - MPROB
0.37 0.42 - - - - - - - - - - S PRB
25.0 25.0 - - - - - - - - - - SIG LW
20.9 9.5 - - - - - - - - - - SIG UP
26.7 26.6 - - - - - - - - - - SNR LW
21.6 10.8 - - - - - - - - - - SNR UP
19.7 19.7 - - - - - - - - - - TGAIN
-0.9 -0.9 - - - - - - - - - - RGAIN
56 60 - - - - - - - - - - SNRxx
Here, z = (87 - 67) / (ABS(26.6) / 1.28) = 0.962, which is approximately 83% or 25 days (see Table 1). Note
again the REL of 0.83 (83%) given by VOACAP.
The Answer
YLE Radio Finland is predicted to maintain a reasonable grade of broadcasting service to Johannesburg (South
Africa) on 26 days out of 30 days in December 2001 on the frequency of 15.520 MHz at 1700 UTC, and on 25 days
out of 30 days at 1800 UTC.
George Lane: It helps to see what is happening when looking at all modes output (Method 25) of VOACAP or
ICEPAC. VOACAP is attempting to compute the hourly median signal power values for the 30 days in the
month at that hour and SSN. It does this for up to 21 different possible modes [ e.g. 2F2, 3F2, 4F2, 2F1, 3F1,
4F1, 4E, 5E and 6E for the high ray and the low ray modes plus the 4Es, 5Es and 6Es for a total of 21 ]. For
each of the possible modes, the program computes the monthly median and the power levels which should be
exceeded on 27 days of the month and on 3 days per month (lower and upper decile). Again remember, these
are predictions of the hourly medians at a given hour over the days of the month and not what is happening
during a single hour of one day.
It is assumed that the presence of multiple modes will increase the available power over the days of the month
compared to what any individual mode might contribute by itself. The calculation process is very simple. The
signal power in dBW is converted to Watts and the values are added to obtain the combined signal power. The
sum is then computed as dBW. This process is done for the median signal power and the upper and lower
decile signal powers to give the monthly distribution. Usually, the contribution of secondary modes to the
monthly distribution of hourly median signal powers is not much but may be 4 to 6 dB. This happens usually
when the antenna patterns support a secondary mode better than the primary mode (a bad design situation but
very common).
In practical terms, the program is attempting to cover what really happens. I have had the opportunity to be
making signal power measurements on paths where there was an oblique ionospheric sounder. Some days the
1F2 mode would be stronger but on other days when the 2F2 MUF would rise sufficiently, that mode would be
stronger. When we look at the distribution of the signal powers at that hour over the 30 days of the month, the
2F2 mode occasionally provided a stronger signal than the 1F2. This means that 2F2 mode improved the
performance of the circuit than if we only had the 1F2 mode. So the presence of other modes which become
stronger than the most reliable mode on some days, will enhance the signal power available than if we had only
the single mode by itself.
54
I think people get confused about the prediction process. They think that the mode contributions are being
added for a single day in the month which doesn't make much sense. We all know that multiple modes at a
given instant can cause harmful multipath which does NOT improve the circuit performance. What the program
is doing is saying that on the days when the primary mode is weak, is there another mode which may come into
play? If so, then the monthly distribution of signal powers is higher than if no other mode is available.
Another problem is that many people think in terms of field strength measurements where the antenna for field
strength meter combines all of the incoming signals into one field strength reading. VOACAP is based on
signal power measurements at the receiver and not field strength measurement. VOACAP computes the power
of each arriving ray path using the appropriate antenna patterns for that takeoff or arrival angle. The summing
process of the individual mode contributions to obtain a combined power distribution for an hour over the days
of the month is forced to fit actual measurements by applying the Transmission Loss Table (Table 7) in the
IONCAP theory manual.
There is no theoretical justification for adding the signal powers of the various modes to obtain the combined
signal power distribution. But there is good agreement with measurement which is what this program is all
about. If you look at Method 25 output, you will see that the contribution of the other modes is not very much
over that of the most reliable mode by itself. But these are the small adjustments in the prediction process that
has made the IONCAP family of programs better than the older or simpler models.
By the way, the IONCAP family of programs does have a multipath probability calculation which is supposed
to give an estimate of the probability that the presence of other modes will cause serious multipath conditions.
Sadly, this calculation is in error. However, funding ran out before I could get the corrections into VOACAP.
George Lane: Yes, the gain of your receive antenna is indeed ignored in the calculation of the noise power.
The noise power is computed from the CCIR radio noise model where the noise power is assumed to arrive
from all directions equally. The use of the theoretical short lossless whip (1/32nd of a wavelength) was to give
slightly higher noise for the lower arrival angles. The main assumption here is that any antenna over ground has
an average gain factor of 2 since all power directed into the ground plane is reflected back into the upper half
space. Thus if the noise arrives from all angles equally, any antenna will produce approximately as much noise
power at the receiver as the short lossless whip.
However, if the antenna is loaded and has a resistive loss, then that loss must be removed from the noise power
as well as the received signal power. Rhombics and terminated vee antennas are examples where the efficiency
factor of the antenna must be subtracted from the CCIR noise power. [Note: the power gain values in the
Method 13 and 14 tables include the efficiency loss] Short vertical antennas such as whips and inverted-L
antennas also have an efficiency loss which is attributed to resistive losses in the ground plane. You can see
these values by running a short fixed length monopole using Method 14. The efficiency values are printed at
the bottom of the gain table.
Antennas such as horizontal dipoles, yagi, log periodic antennas are assumed to have no efficiency losses. By
the way, this is a bad assumption when the dipole element is closer than 1/16th wavelength over the ground. At
that point, ground losses begin to show up in the input impedance of the antenna and are a real resistive loss!
What about the SWWHIP receive antenna that produces less noise power than an isotrope? In the SWWHIP
file we can see that an efficiency factor of -4.8 dB has been applied.
55
George Lane: Basically, you need to integrate the power gain pattern of the antenna in the half space over the
ground. This is often called the 'average' gain of the antenna. If the gain factor is 2, then the efficiency is 0 dB.
In the case of my shortwave whip model (based on actual measurements of a 15 ft. whip antenna over a ground
plane with 15 radials of 15 ft. each at 2 MHz in 1943 at Ohio State University [classified SECRET at the
time]), the average gain is 0.331 rather than 2 (10Log(0.331) = -4.8 dB). This value is subtracted from the
received noise power, only. The gain pattern of the antenna is applied to the signal power at the appropriate
arrival angles.
George Lane: It is not necessary to run Method 15 in order to obtain valid circuit predictions. However, it is
necessary to input the antenna data correctly in order to obtain valid circuit predictions. It is very easy to
generate a wrong antenna pattern and it is very difficult to determine that the antenna pattern is wrong when
you are just looking at the circuit predictions. Therefore, I recommend that users always run Method 15 when
beginning an analysis to make sure the pattern for the transmit and receive antennas look reasonable for the
antenna you are trying to model.
The pattern is shown as a table of gain values in dB relative to an isotrope in free space for the great circle route
azimuth defined by the geographic coordinates you specified for the transmit and receive locations. The pattern
is shown by elevation angles with 0 degrees being along the horizon and 90 degrees being straight up (y-axis)
and by frequencies from 2 to 30 MHz (x-axis). For example, if you are modeling a horizontal dipole antenna at
one quarter wave length above ground, you would expect to see the maximum gain of 5 to 8 dBi at the high
elevation angles and negative gain values at the low elevation angles. Also the pattern should be relatively
constant from frequency to frequency across the table.
George Lane: When you input the antenna data, you must specify an antenna azimuth or select "at" the
transmit or receive site. If you specify an azimuth, then the program will determine the difference between the
great circle route azimuth and the one you specified. Then the pattern will be computed for that many degrees
from the mainbeam azimuth for that antenna. For some antennas which are omnidirectional, one can ignore the
off-azimuth. The isotrope, whip and discone are some examples of omni-directional antennas.
The antennas do have a real azimuth and that should be entered into VOACAP. Remember that VOACAP uses
azimuths relative to True North. Often the azimuth of the antenna as built is specified in terms of Magnetic
North as obtained from a compass reading. You will need to correct magnetic azimuths to true north azimuths
using the correct declination angle for that geographic location.
It is difficult to visualize the antenna patterns, as Method 15 gives only textual, not graphical, results.
George Lane: When you want to generate an antenna pattern for the antenna of interest, you should use
HFANT. With HFANT you can see the pattern plotted out on a graph. When you have obtained patterns that
you feel represent the antenna in question, be sure to name the antenna and save it in your antenna folder. Then
56
when you run VOACAP, you can pull that antenna design from your antenna file. I still recommend that on the
first run for that circuit, you run Method 15 just to make sure the right patterns are being called.
I remember being called by an individual who was getting very poor predictions on a path that should have
worked very well with a horizontal dipole antenna. After much work, I noticed the dipole was nearly 90
degrees off-azimuth! It seems he thought the mainbeam azimuth of the dipole was the direction of the wire and
not broadside to it. If he had run Method 15, he would have soon discovered that the dipole was not working
the way it should.
I have calculated MUF, SNR, etc. without calculating the antenna pattern with Method 15. Is this fatal to
the results?
George Lane: NO. The MUF is dependent on the circuit, time of day, month and sunspot number and is
independent of the antenna, transmitter power, etc. But, SNR, REL, RPWRG, Sig PWR, and DBU are
dependent on the system parameters which include the antenna patterns.
George Lane: You are absolutely correct about the takeoff and arrival angle calculations in VOACAP. They
are fixed values for simple ray hops. Exotic modes such as N and M are ignored, too. The reason is that there
was no measured data which would provide a statistical representation of the angle distribution.
Around 1990, Dr. Leo McNamara invited me to see a demonstration of the SkyLoc system which had been set
up at the FCC HFDF site here in Columbia, MD. This system provided full angle of arrival data which was
sampled at a very rapid rate. What amazed me was that the data collected over 5 minutes or so was showing
strongest signals at arrival angles spanning 10 to 30 degrees and with a azimuthal spread of plus or minus about
2 to 3 degrees. We collected quite a few samples for different stations and then ran VOACAP. It was
interesting that the model predicted the arrival angle within a few degrees for the actual angle observed with the
highest number of occurrences. But it was very obvious that many angles were being used over a short period
of time.
I was invited back to the FCC HFDF site in Columbia after I retired from VOA. It seems that they built a
home-brew design for a phased array which obtains full angle of arrival data in real time for any signal you
want to look at... not only from Columbia MD but all of the 13 other remote FCC DF sites.
Yes, takeoff and arrival angle predictions are WOEfully inadequate in VOACAP and all other models that I
know of. I suspect that the minute-to-minute fading we observe in signal power is a direct result of ever
changing angles. I think the effect is very large too... probably 8 to 16 dB. It was my suspicion that the great
range in takeoff and arrival angles is why theory and measurement never agreed.
What I saw with the real-time SkyLoc HFDF system is that in a few minutes, an entire envelope of arrival
angles were measured. These arrival angles may be clustered around a single median or about several medians
for multiple ray hop modes. What would be nice is to have enough worldwide data to develop the statistical
window for the takeoff/arrival angles. The current FCC HFDF network might provide that data for a person
interested in exploring this issue.
When I was a Lt. in the Army, I was tutored in HF system design by Tim Shaw of the old Radio Propagation
Agency. I was responsible for the design of the antennas being installed in SE Asia and Tim cautioned me not
to use antennas whose beam widths were less than 10 degrees in either the vertical or horizontal plane. He cited
cases where they had used antenna with beam widths as narrow as 5 degrees only to find that they experienced
deep fading which he attributed to instantaneous variations in the arrival angle. Nearly 30 years later, I watched
57
in real time as the SkyLoc system recorded the distribution of arrival angles. How I wished Tim was still alive
to see that! Time after time we would see plus or minus 4 to 8 degrees around the central angle per mode.
George Lane: According to Donald Lucas, using 0.1 degrees gives you a better calculation of the frequency
dependence of the ionosphere (i.e. the MUF for the lowest order mode for the circuit). John Lloyd told me the
value of 3 degrees was better from a practical standpoint. It depends on the actual horizon clearance at both the
transmit and receive site. If you are using a realistic radiation pattern taking into account any horizon
obstructions, then the low angles below 3 degrees will disappear anyhow. The reason for this is that gain at
angles below 3 degrees for antennas over real ground is very low and also because John Lloyd programmed
IONCAP (now VOACAP) to favor angles of 3 degrees or higher.
So Lucas is right (use 0.1 degrees for minimum angle) when the radiation patterns for the antenna have very
low gains at near-horizon angles. But if you don't know the low angle gain and/or are using isotropes, then
Lloyd is right (use 3.0 degrees for minimum angle). If you see many cases where TANGLE is below 3 degrees,
you should suspect you have used erroneous input somewhere.
Is there something inherent to VOACAP to make the model always pick the lower takeoff and arrival
angles?
George Lane: The simple answer is: Yes, the model likes the lower angles because the path length is less and
the transmission losses are less.
The more complicated answer: Yes, VOACAP can miss the most reliable mode which is at a higher angle.
How can this be? We have to remember that VOACAP only considers 3 modes for each ionospheric layer. So
if VOACAP starts with 2F2 mode, then only the 3F2 mode and the 4F2 mode will be considered for the most
reliable mode. But let's say that the antenna patterns have their maximum gain at higher angles, will the 5F2
mode be considered. No, it will not.
What can we do to prevent VOACAP from finding the 'best' mode for the antennas we are using? The answer is
fairly simple. If both the transmit and receive antenna have positive gain in dBi at angles below 3 degrees and
there are no horizon obstructions, then one sets the minimum angle to 0.1 degrees in the input to VOACAP. If
the transmit and receive antennas have patterns that combine to give maximum gain at angles above 3 degrees,
then set the minimum angle to 3.0 degrees.
Let's say I have a site which has an obstruction at the horizon of 9 degrees, should I use a minimum angle of 9
degrees? No! When you specify minimum angles above 3 degrees, you can make the program reduce the
MUF. This causes the program to make all kinds of stupid predictions.Many people think that the minimum
angle is used to account for horizon obstructions and that is not true. The program lets you decide whether
you can use angles of less than 3 degrees but should not be used to block higher angles above 3 degrees as that
modifies the ionosphere. The best solution is to modify the antenna pattern for horizon obstructions using a
program such as NEC.
58
So I set the minimum angle to 3.0 degrees in VOACAP. Now for these antennas with higher takeoff angles, I
find the program selects the 5F2 mode at 21 degrees as the most reliable mode. Why did this happen? Before
the program had found the 2F2 mode at 2.4 degrees as the lowest order mode. Then it applied the 3 hop rule
and looked at the 3F2 and the 4F2 modes and stopped. When I used the 3 degree minimum angle, the 2F2 mode
was blocked, so VOACAP started with the 3F2 mode at 9.4 degrees and then found the 4F2 mode at 15.2
degrees and the 5F2 mode at 20.7 degrees. The 5F2 mode has a median SNR of 23 dB*Hz. This is getting close
to a SNR that can be used for reception of manual CW between trained operators. Also the MUF remained at
16.4 MHz for the 2F2 mode.
I then set the minimum angle to 9.4 degrees to notch out the 3F2 mode. This caused VOACAP to reduce the
MUF from 16 MHz to 13 MHz. That is not good as we have now modified the ionosphere with an obstruction
angle. THAT IS A NO-NO!
Conclusion
For most amateur applications, one should use a minimum angle of 3.0 degrees. The default value in
VOACAP is 0.1 degrees. The reason for this is that the Voice of America spends a lot of money to build
antennas that do produce considerable gain at angles below 3 degrees so we changed the IONCAP default
from 3 degrees to 0.1 degrees for that reason.
A note for VOA: If you use one of your variable beam arrays in the high angle mode, you should change the
minimum angle from 0.1 to 3 degrees!
A word of caution about isotropic patterns: Never use an isotrope with the minimum angle set to 0.1
degrees. The reason is that a gain of 0 dBi at 0.1 degrees is almost infinite gain from an antenna over a ground
plane. VOACAP will lock on to that impossible gain and make predictions for modes that are really below the
horizon. The program will ignore higher angle modes which are above the three hops which the program will
consider.
George Lane: For path lengths less than 7000 km, the angle prediction takes into account the gain of both the
TX and RX antennas. Remember that VOACAP develops a reflectrix table of possible takeoff and arrival
angles which will be supported on that path hour and operating frequency. The program finds the most reliable
mode and its associated angle based on combining the path loss for the possible angles and the combined gain
of the TX and RX antenna at each possible angle. For longer paths (over 7000 km) where the long-path
calculation model (Method 21) can become the dominant propagation mechanism, the best takeoff angle and
best arrival angle are computed and printed out based on the ionospheric conditions between 1000 and 2000 km
of each end of the circuit.
So, from a modelling point of view, is it better to use an isotropic RX antenna with +6 dBi gain, for instance,
when we actually do not know all RX antennas, or should we just choose an antenna that could be perhaps
most likely?
George Lane: Good question. There are several factors. First we don't know what type of antenna is being
used at the distant station. The directivity pattern is probably similar to a dipole but we do not know how well
the antenna performs at the low angles nor do we know if the antenna is pointed at our location. It is the similar
problem when you are dealing with ground-to-air HF communications. I have found engineers spend great
sums of money to obtain accurate models of the aircraft antenna radiation pattern. But in practice, you don't
know the orientation of the aircraft with respect to the ground station. So when in doubt about the distant
station radiation pattern, the isotropic pattern is a reasonable first assumption. Any antenna placed over earth
produces an average gain of 3 dBi (half of the power is radiated into the hemisphere and the other half is
59
reflected from the ground for a net gain of 3 dB over that of an isotrope in free space). Therefore, I would
recommend using up to a 3 dBi omni-directional pattern (isotrope +3 dBi).
One must be careful with isotropes, however. Zero dBi or +3 dBi at 0.1 degrees elevation angle is a tremendous
amount of gain compared to what is possible over real earth (-20 to -30 dBi is more realistic). With such gain at
low angles, one will force VOACAP into erroneously predicting low angle modes with a very high signal
power. To minimize the effect of this high gain at low angles, one should set the Minimum Angle at 3.0
degrees. Even a gain of 3 dBi at 3 degrees above the horizon requires a very good antenna.
The second part of your question is whether or not there is a better antenna assumption for the distant station. If
you are only dealing with a point-to-point circuit, one could assume, for example, a half wave dipole at 30 feet
in height as a typical antenna. However, if you are using VOAAREA, one must specify a bearing angle for the
receive antenna. With the assumed dipole, this bearing will be used at each grid point for the receive dipole.
That means the radiation pattern of the assumed receive antenna will vary from one grid point to the next. One
could use Inverse VOAAREA which would allow one to point the dipole at the grid points toward your receive
location. But then again, one is faced with the problem of what bearing you should use for your antenna. I still
prefer VOAAREA using an isotrope with up to 3 dB additional gain and a Minimum Angle set at 3.0 degrees
for the general case when you don't know the details of the distant stations.
3.5 Conversions
S DBW to S Meter Correlation Table
The Z Tables
There are a few assumptions: the receiver input impedance is 50 Ohms and the transmission lineloss is
minimal. Also we assume that the S meter is correctly calibrated. In practice, S meters are mostly imprecise and
can only give relative readings.
60
The Z Tables
3.5.2.1 Converting Z values to Percentage Values and Corresponding
Number of Days
Use the two tables below to convert the value of z to the probability (%) value and to the corresponding number
of days.
for converting a percentage value to days, just find the value. The corresponding number of days is in
the right column.
for calculating the REL, use Table 1 if the REQ.SNR is lower or equals to the predicted SNR. Use
Table 2, if the REQ.SNR is greater than the predicted SNR.
for calculating the MUFday, use Table 1 if the FREQ is equal to or below the predicted MUF. Use
table 2 if the FREQ is above the predicted MUF.
Table 1. Table 2.
z % Days/month z % Days/month
61
2.10 98.21 29.46 2.10 1.79 0.54
2.15 98.42 29.53 2.15 1.58 0.47
2.20 98.61 29.58 2.20 1.39 0.42
2.25 98.78 29.63 2.25 1.22 0.37
2.30 98.93 29.68 2.30 1.07 0.32
2.35 99.06 29.72 2.35 0.94 0.28
2.40 99.18 29.75 2.40 0.82 0.25
2.45 99.29 29.79 2.45 0.71 0.21
2.50 99.38 29.81 2.50 0.62 0.19
2.55 99.46 29.84 2.55 0.54 0.16
2.60 99.53 29.86 2.60 0.47 0.14
2.65 99.60 29.88 2.65 0.40 0.12
2.70 99.65 29.90 2.70 0.35 0.10
2.75 99.70 29.91 2.75 0.30 0.09
2.80 99.74 29.92 2.80 0.26 0.08
2.85 99.78 29.93 2.85 0.22 0.07
2.90 99.81 29.94 2.90 0.19 0.06
2.95 99.84 29.95 2.95 0.16 0.05
3.00 99.87 29.96 3.00 0.13 0.04
3.05 99.89 29.97 3.05 0.11 0.03
3.10 99.90 29.97 3.10 0.10 0.03
3.15 99.92 29.98 3.15 0.08 0.02
3.20 99.93 29.98 3.20 0.07 0.02
3.25 99.94 29.98 3.25 0.06 0.02
3.30 99.95 29.99 3.30 0.05 0.01
3.35 99.96 29.99 3.35 0.04 0.01
3.40 99.97 29.99 3.40 0.03 0.01
These automatic beacon stations transmit every three minutes on the given frequency, 24 hours a day. A
transmission consists of the callsign sent in Morse code at 22 words per minute followed by four one-second
dashes. The callsign and the first dash are sent at 100 watts. The remaining dashes are sent at 10 watts, 1 watt
and 0.1 watts.
The different colours will indicate the reliability (%) to achieve the required SNR for communication. The best
reliability is shown in white and the worst in dark gray. Look at the colour legend for values.
Now when you are all set, try running the circuit analyses from various NCDXF beacon locations to your
receive location, and you will get a glimpse of best times for openings to all corners of the world.
63
Start by reading the RNW (Radio Nederland) article on
RUV,http://www.rnw.nl/realradio/features/html/iceland020111.html. Look for the following data:
transmitter coordinates
transmitter power
times & frequencies of transmission
transmitter antenna
antenna beams
As a result, the VOAAREA main input screen should look like this:
64
3.5.4.3 Run the prediction
Choose Run > Calculate > Save{temp.VOA}/Calculate.
To plot the map after calculations, choose Run > Plot results.
Open the first file to be plotted: File > open TEMP.VG1.
If you do not wish to change the default colours or values, just choose Plot to > Window.
To view the 18 UT map, first Exit the current map, and choose File > Open > temp.vg2.
Voila!
This is approximately what you should see (I have changed some of the default colours for my maps):
65