Python3 RF Exposure Calculator Code

Select Operating Mode, Antenna Type, and length of Time transmitting.
Tables of Calculations will be based on 5, 10, 25, 50, 75, 100, 250, 500, 750 1000, and 1500 Watts, for the Bands 160m through 23cm, will be generated.

ARRL_Based_RF-Exposure-Calculator_All-Bands_5-10-25-50-75-100-250-500-750-1000-1200-1500W.03.py

'''
 Created: 24.06.2022
Modified: 03.12.2023

Code written with python 3.12.0
Code tested with python 3.9.13, 3.10.11, 3.11.6, 3.12.0, and 3.13.0a2

This Code accepts Type of Mode, Type of Antenna, and Number of Minutes Transmitting within a 30 Minute Contact period. It then calculates the max Density and Safe distance for Controlled and Uncontrolled RF Exposure; at the Powers of 5, 10, 25, 50, 75, 100, 250, 500, 750, 1000, 1200, and 1500 Watts; for the Bands 160m through 23cm.
'''

# ------------------------------------------------------------------------------
import os                           # Built in Python Module / Library.   
import math as m                    # Built in Python Module / Library.
import sys                          # Built in Python Module / Library.

# ------------------------------------------------------------------------------
def Clear_Screen():
 try:
  print('\033[2J\033[3J\033[1;1H')
 except:
  os.system('cls') # Clear a 'Windows' Terminal Window.
# ------------------------------------------------------------------------------
def ResizeWindow(lRows, lRowWidth):
 sys.stdout.write("\x1b[8;{rows};{cols}t".format(rows=lRows, cols=lRowWidth))
# ------------------------------------------------------------------------------
def TimeAvgPercent(tx, rx, interval):
 # interval is either 6min or 30min per spec in a controlled vs. uncontrolled environment respectively

cycle = tx + rx
 remainder = interval % cycle
 txcompletecycles = m.floor( interval / cycle) * tx #// tx minutes for complete cycles
 
 if(tx>=interval):
  return 1.0
 elif(cycle>=interval):
  return (tx / interval)	 #// edge case where 1 tx+rx cycle is larger than interval
 elif(tx>remainder):
  return ((txcompletecycles + remainder) / interval);
 else:
  return ((txcompletecycles+tx) / interval);
# ------------------------------------------------------------------------------

entryErrorMessage='Try again, with a valid Entry.'
groundRef=0.64 # If Ground Adjustment was not considered groundRef=25

tWidth=132

dashLine=('-'*tWidth)
ResizeWindow(45, tWidth)
Clear_Screen()

modeList='A: Conversational CW (duty cycle=40%)'+'\n'+'B: Conversational SSB, heavy speech processing (duty cycle=50%)'+'\n'+'C: Conversational SSB, no speech processing (duty cycle=20%)'+'\n'+'D: FM                (duty cycle=100%)'+'\n'+'E: AM                (duty cycle=100%)'+'\n'+'F: FSK/RTTY          (duty cycle=100%)'+'\n'+'G: AFSK SSB          (duty cycle=100%)'+'\n'+'H: Carrier Always On, i.e., for Tuning Up      (duty cycle=100%)'+'\n'+'I: For all others, or if unknown, assume worst case scenario (duty cycle=100%)'

antennaList='A: 1/2 wave Dipole    2.15 dBi'+'\n'+'B: 10 Element Yagi   15.1 dBi'+'\n'+'C:  2 Element Yagi    5.9 dBi'+'\n'+'D:  3 Element Yagi    8.1 dBi'+'\n'+'E:  4 Element Yagi    9.1 dBi'+'\n'+'F:  5 element Yagi   10.1 dBi'+'\n'+'G:  6 Element Yagi   11.1 dBi'+'\n'+'H:  8 Element Yagi   13.1 dBi'+'\n'+'I: Delta Loop         5.2 dBi'+'\n'+'J: Four Square        5.2 dBi'+'\n'+'K: G5RV               1.0 dBi'+'\n'+'L: Hex Beam           5.0 dBi'+'\n'+'M: Moxon              6.0 dBi'+'\n'+'N: 1/4 Wave Vertical  1.5 dBi'+'\n'+'O: Windom (OCD)       2.0 dBi'

allBands=[1.8, 2.0, 3.5, 4.0, 5.332, 5.405, 7.0, 7.3, 10.1, 10.15, 14.0, 14.35, 18.068, 18.168, 21.0, 21.45, 24.890, 24.990, 28.0, 29.7, 50.0, 54.0, 144.0, 148.0, 222.0, 225.0, 420.0, 450.0, 902.0, 928.0, 1240.0, 1300.0]

tWattsList=[5, 10, 25, 50, 75, 100, 250, 500, 750, 1000, 1200, 1500]

print()
print(modeList)
mode=(input('Respective Letter: ')).upper()
try:
 if(mode.isnumeric()):raise
except:
 print(entryErrorMessage)
 exit()
 
textSpacer=''
if(mode=='A'):
 modeDC=0.4
 tMode='CW'
elif(mode=='B'):
 modeDC=0.5
 tMode='SSB'
elif(mode=='C'):
 modeDC=0.2
 tMode='SSB'
elif(mode=='D'):
 modeDC=0.2
 tMode='FM'
elif(mode=='E'):
 modeDC=1.0
 tMode='AM'
elif(mode=='F'):
 modeDC=1.0
 tMode='FSK/RTTY'
elif(mode=='G'):
 modeDC=1.0
 tMode='AFSK SSB'
elif(mode=='H'):
 modeDC=1.0
 tMode='Carrier only'
 textSpacer='   '
elif(mode=='I'):
 modeDC=1.0
 tMode='Worst Case'
 textSpacer='  '
else:
 print(entryErrorMessage)
 exit()

print()
print(antennaList)
antenna=(input('Respective Letter: ')).upper()
try:
 if(antenna.isnumeric()):raise
except:
 print(entryErrorMessage)
 exit()

if(antenna=='A'):
 antennaGain=2.15
 tAntenna='1/2 wave Dipole'
elif(antenna=='B'):
 antennaGain=15.1
 tAntenna='10 Element Yagi'
elif(antenna=='C'):
 antennaGain=5.9
 tAntenna='2 Element Yagi'
elif(antenna=='D'):
 antennaGain=8.1
 tAntenna='3 Element Yagi'
elif(antenna=='E'):
 antennaGain=9.1
 tAntenna='4 Element Yagi'
elif(antenna=='F'):
 antennaGain=10.1
 tAntenna='5 Element Yagi'
elif(antenna=='G'):
 antennaGain=11.1
 tAntenna='6 Element Yagi'
elif(antenna=='H'):
 antennaGain=13.1
 tAntenna='8 Element Yagi'
elif(antenna=='I'):
 antennaGain=5.2
 tAntenna='Delta Loop'
elif(antenna=='J'):
 antennaGain=5.2
 tAntenna='Four Square'
elif(antenna=='K'):
 antennaGain=1.0
 tAntenna='G5RV'
elif(antenna=='L'):
 antennaGain=5.0
 tAntenna='Hex Beam'
elif(antenna=='M'):
 antennaGain=6.0
 tAntenna='Moxon'
elif(antenna=='N'):
 antennaGain=1.5
 tAntenna='1/4 Wave Vertical'
elif(antenna=='O'):
 antennaGain=2.0
 tAntenna='Windom'
else:
 print(entryErrorMessage)
 exit()

print()
try:
 txTime=input('Transmit duty cycle: (time Tx-ing, in Minutes [30 max.]): ')
 txTime=float(txTime)
 #rxTime=input('Transmit duty cycle: (time Rx-ing, in Minutes [30 max.]): ')

rxTime=30-txTime
 print('Transmit duty cycle: (time Rx-ing, in Minutes [30 max.]): '+str(int(rxTime)))
 tResult=input("Press '<return>' to continue.")
 rxTime=float(rxTime)

if((txTime>30) or (rxTime>30)): raise
except:
 print(entryErrorMessage)
 exit()
print()

topStringRow01='   Freq.   Power   '+textSpacer+'Mode        Antenna      Antenna   TX  RX        CMD           CMSD              UMD            UMSD'
topStringRow02='    MHz      W                              '+textSpacer+'Gain dBi  Min Min    (mw/cm2)      mtrs | ft         (mw/cm2)       mtrs | ft'

bottomStringRow01=' CMD: Controlled Maximum Density. The Amateur Radio Operator is the Controlled.'
bottomStringRow02='CMSD: Controlled Minimum Safe Distance.'
bottomStringRow03=' UMD: UnControlled Maximum Density. Anyone other than the Amateur Radio Operator is the UnControlled.'
bottomStringRow04='UMSD: UnControlled Minimum Safe Distance.'

Clear_Screen()

for watts in tWattsList:

print(topStringRow01)
 print(topStringRow02)
 print(dashLine)

watts=float(watts)

pwr1=1000.0*watts*modeDC*TimeAvgPercent(txTime, rxTime, 6)
 eirp1=pwr1*m.pow(10,antennaGain/10.0)
 numerator1=groundRef*eirp1

pwr2=1000.0*watts*modeDC*TimeAvgPercent(txTime, rxTime, 30)
 eirp2=pwr2*m.pow(10,antennaGain/10.0)
 numerator2=groundRef*eirp2
 
 for freq in allBands:
  if (freq < 1.34):
   CMaxDensity = 100.0
   UMaxDensity = 100.0
  elif (freq < 3):
   CMaxDensity = 100.0
   UMaxDensity = 180.0 / (m.pow(freq, 2.0))
  elif (freq < 30):
   CMaxDensity = 900.0 / (m.pow(freq, 2.0))
   UMaxDensity = 180.0 / (m.pow(freq, 2.0))
  elif(freq < 300.0):
   CMaxDensity = 1.0
   UMaxDensity = 0.2
  elif (freq < 1500.0):
   CMaxDensity = freq / 300.0
   UMaxDensity = freq / 1500.0
  elif(freq < 100000.0):
   CMaxDensity = 5.0
   UMaxDensity = 1.0
  else:
   input("The FCC does not have exposure limits above 100 GHz")
   exit()
  
  denominator=CMaxDensity*m.pi
  CMinDistance=m.sqrt(numerator1 / denominator)
  CMinDistance=CMinDistance/30.48
  CMinSafeDistance=CMinDistance*0.3048
  
  denominator=UMaxDensity*m.pi
  UMinDistance=m.sqrt(numerator2 / denominator)
  UMinDistance=UMinDistance/30.48
  UMinSafeDistance=UMinDistance*0.3048

print(('{:5.3f}'.format(freq)).rjust(8, ' ')+'   '+(str(int(watts))).rjust(4, ' ')+'  '+'{:^8s}'.format(tMode)+'  '+'{:^17s}'.format(tAntenna)+'  '+'{:^4s}'.format(str(antennaGain)).rjust(5, ' ')+'   '+(str(int(txTime))).rjust(2, ' ')+'  '+(str(int(rxTime))).rjust(2, ' ')+' '+'   '+('{:5.3f}'.format(CMaxDensity)).rjust(8, ' ')+'    '+('{:.3f}'.format(CMinSafeDistance)).rjust(7, ' ')+' | '+('{:5.3f}'.format(CMinDistance)).rjust(7, ' ')+'    '+('{:5.3f}'.format(UMaxDensity)).rjust(8, ' ')+'    '+('{:5.3f}'.format(UMinSafeDistance)).rjust(7, ' ')+' | '+('{:5.3f}'.format(UMinDistance)).rjust(7, ' '))
 
 print()
 print(dashLine)

print(bottomStringRow01)
print(bottomStringRow02)
print(bottomStringRow03)
print(bottomStringRow04)
print(dashLine)

The Code above is pure and simple Python.
For colored Text, by installing the 'termcolor' Module, use the Code below.

ARRL_Based_RF-Exposure-Calculator_w|Termcolor_All-Bands_5-10-25-50-75-100-250-500-750-1000-1200-1500W.03.py

''' ''' Created: 24.06.2022 Modified: 03.12.2023 Code written with python 3.12.0 Code tested with python 3.9.13, 3.10.11, 3.11.6, 3.12.0, and 3.13.0a2 This Code accepts Type of Mode, Type of Antenna, and Number of Minutes Transmitting within a 30 Minute Contact period. It then calculates the max Density and Safe distance for Controlled and Uncontrolled RF Exposure; at the Powers of 5, 10, 25, 50, 75, 100, 250, 500, 750, 1000, 1200, and 1500 Watts; for the Bands 160m through 23cm. ''' # ------------------------------------------------------------------------------ import os # Built in Python Module / Library. import math as m # Built in Python Module / Library. import sys # Built in Python Module / Library. from termcolor import colored as cd # Must be downloaded and installed via 'pip3 install termcolor'. # ------------------------------------------------------------------------------ def Clear_Screen(): try: print('\033[2J\033[3J\033[1;1H') except: os.system('cls') # Clear a 'Windows' Terminal Window. # ------------------------------------------------------------------------------ def ResizeWindow(lRows, lRowWidth): sys.stdout.write("\x1b[8;{rows};{cols}t".format(rows=lRows, cols=lRowWidth)) # ------------------------------------------------------------------------------ def TimeAvgPercent(tx, rx, interval): # interval is either 6min or 30min per spec in a controlled vs. uncontrolled environment respectively cycle = tx + rx remainder = interval % cycle txcompletecycles = m.floor( interval / cycle) * tx #// tx minutes for complete cycles if(tx>=interval): return 1.0 elif(cycle>=interval): return (tx / interval) #// edge case where 1 tx+rx cycle is larger than interval elif(tx>remainder): return ((txcompletecycles + remainder) / interval); else: return ((txcompletecycles+tx) / interval); # ------------------------------------------------------------------------------ entryErrorMessage='Try again, with a valid Entry.' groundRef=0.64 # If Ground Adjustment was not considered groundRef=25 tWidth=132 magentaDashLine=cd('-'*tWidth, 'magenta') ResizeWindow(45, tWidth) Clear_Screen() modeList='A: Conversational CW (duty cycle=40%)'+'\n'+'B: Conversational SSB, heavy speech processing (duty cycle=50%)'+'\n'+'C: Conversational SSB, no speech processing (duty cycle=20%)'+'\n'+'D: FM (duty cycle=100%)'+'\n'+'E: AM (duty cycle=100%)'+'\n'+'F: FSK/RTTY (duty cycle=100%)'+'\n'+'G: AFSK SSB (duty cycle=100%)'+'\n'+'H: Carrier Always On, i.e., for Tuning Up (duty cycle=100%)'+'\n'+'I: For all others, or if unknown, assume worst case scenario (duty cycle=100%)' antennaList='A: 1/2 wave Dipole 2.15 dBi'+'\n'+'B: 10 Element Yagi 15.1 dBi'+'\n'+'C: 2 Element Yagi 5.9 dBi'+'\n'+'D: 3 Element Yagi 8.1 dBi'+'\n'+'E: 4 Element Yagi 9.1 dBi'+'\n'+'F: 5 element Yagi 10.1 dBi'+'\n'+'G: 6 Element Yagi 11.1 dBi'+'\n'+'H: 8 Element Yagi 13.1 dBi'+'\n'+'I: Delta Loop 5.2 dBi'+'\n'+'J: Four Square 5.2 dBi'+'\n'+'K: G5RV 1.0 dBi'+'\n'+'L: Hex Beam 5.0 dBi'+'\n'+'M: Moxon 6.0 dBi'+'\n'+'N: 1/4 Wave Vertical 1.5 dBi'+'\n'+'O: Windom (OCD) 2.0 dBi' allBands=[1.8, 2.0, 3.5, 4.0, 5.332, 5.405, 7.0, 7.3, 10.1, 10.15, 14.0, 14.35, 18.068, 18.168, 21.0, 21.45, 24.890, 24.990, 28.0, 29.7, 50.0, 54.0, 144.0, 148.0, 222.0, 225.0, 420.0, 450.0, 902.0, 928.0, 1240.0, 1300.0] tWattsList=[5, 10, 25, 50, 75, 100, 250, 500, 750, 1000, 1200, 1500] print() print(modeList) mode=(input(cd('Respective Letter: ', 'cyan'))).upper() try: if(mode.isnumeric()):raise except: print(entryErrorMessage) exit() textSpacer='' if(mode=='A'): modeDC=0.4 tMode='CW' elif(mode=='B'): modeDC=0.5 tMode='SSB' elif(mode=='C'): modeDC=0.2 tMode='SSB' elif(mode=='D'): modeDC=0.2 tMode='FM' elif(mode=='E'): modeDC=1.0 tMode='AM' elif(mode=='F'): modeDC=1.0 tMode='FSK/RTTY' elif(mode=='G'): modeDC=1.0 tMode='AFSK SSB' elif(mode=='H'): modeDC=1.0 tMode='Carrier only' textSpacer=' ' elif(mode=='I'): modeDC=1.0 tMode='Worst Case' textSpacer=' ' else: print(entryErrorMessage) exit() print() print(antennaList) antenna=(input(cd('Respective Letter: ', 'cyan'))).upper() try: if(antenna.isnumeric()):raise except: print(entryErrorMessage) exit() if(antenna=='A'): antennaGain=2.15 tAntenna='1/2 wave Dipole' elif(antenna=='B'): antennaGain=15.1 tAntenna='10 Element Yagi' elif(antenna=='C'): antennaGain=5.9 tAntenna='2 Element Yagi' elif(antenna=='D'): antennaGain=8.1 tAntenna='3 Element Yagi' elif(antenna=='E'): antennaGain=9.1 tAntenna='4 Element Yagi' elif(antenna=='F'): antennaGain=10.1 tAntenna='5 Element Yagi' elif(antenna=='G'): antennaGain=11.1 tAntenna='6 Element Yagi' elif(antenna=='H'): antennaGain=13.1 tAntenna='8 Element Yagi' elif(antenna=='I'): antennaGain=5.2 tAntenna='Delta Loop' elif(antenna=='J'): antennaGain=5.2 tAntenna='Four Square' elif(antenna=='K'): antennaGain=1.0 tAntenna='G5RV' elif(antenna=='L'): antennaGain=5.0 tAntenna='Hex Beam' elif(antenna=='M'): antennaGain=6.0 tAntenna='Moxon' elif(antenna=='N'): antennaGain=1.5 tAntenna='1/4 Wave Vertical' elif(antenna=='O'): antennaGain=2.0 tAntenna='Windom' else: print(entryErrorMessage) exit() print() try: txTime=input('Transmit duty cycle: (time Tx-ing, in Minutes [30 max.]): ') txTime=float(txTime) #rxTime=input('Transmit duty cycle: (time Rx-ing, in Minutes [30 max.]): ') rxTime=30-txTime print('Transmit duty cycle: (time Rx-ing, in Minutes [30 max.]): '+str(int(rxTime))) tResult=input(cd('<return>', 'cyan', attrs=['blink'])) rxTime=float(rxTime) if((txTime>30) or (rxTime>30)): raise except: print(entryErrorMessage) exit() print() topStringRow01=cd(' Freq. Power '+textSpacer+'Mode Antenna Antenna TX RX CMD CMSD UMD UMSD', 'green') topStringRow02=cd(' MHz W '+textSpacer+'Gain dBi Min Min (mw/cm2) mtrs | ft (mw/cm2) mtrs | ft', 'green') bottomStringRow01=cd(' CMD: Controlled Maximum Density. The Amateur Radio Operator is the Controlled.', 'cyan') bottomStringRow02=cd('CMSD: Controlled Minimum Safe Distance.', 'cyan') bottomStringRow03=cd(' UMD: UnControlled Maximum Density. Anyone other than the Amateur Radio Operator is the UnControlled.', 'cyan') bottomStringRow04=cd('UMSD: UnControlled Minimum Safe Distance.', 'cyan') Clear_Screen() for watts in tWattsList: print(topStringRow01) print(topStringRow02) print(magentaDashLine) watts=float(watts) pwr1=1000.0*watts*modeDC*TimeAvgPercent(txTime, rxTime, 6) eirp1=pwr1*m.pow(10,antennaGain/10.0) numerator1=groundRef*eirp1 pwr2=1000.0*watts*modeDC*TimeAvgPercent(txTime, rxTime, 30) eirp2=pwr2*m.pow(10,antennaGain/10.0) numerator2=groundRef*eirp2 for freq in allBands: if (freq < 1.34): CMaxDensity = 100.0 UMaxDensity = 100.0 elif (freq < 3): CMaxDensity = 100.0 UMaxDensity = 180.0 / (m.pow(freq, 2.0)) elif (freq < 30): CMaxDensity = 900.0 / (m.pow(freq, 2.0)) UMaxDensity = 180.0 / (m.pow(freq, 2.0)) elif(freq < 300.0): CMaxDensity = 1.0 UMaxDensity = 0.2 elif (freq < 1500.0): CMaxDensity = freq / 300.0 UMaxDensity = freq / 1500.0 elif(freq < 100000.0): CMaxDensity = 5.0 UMaxDensity = 1.0 else: input("The FCC does not have exposure limits above 100 GHz") exit() denominator=CMaxDensity*m.pi CMinDistance=m.sqrt(numerator1 / denominator) CMinDistance=CMinDistance/30.48 CMinSafeDistance=CMinDistance*0.3048 denominator=UMaxDensity*m.pi UMinDistance=m.sqrt(numerator2 / denominator) UMinDistance=UMinDistance/30.48 UMinSafeDistance=UMinDistance*0.3048 print(('{:5.3f}'.format(freq)).rjust(8, ' ')+' '+(str(int(watts))).rjust(4, ' ')+' '+'{:^8s}'.format(tMode)+' '+'{:^17s}'.format(tAntenna)+' '+'{:^4s}'.format(str(antennaGain)).rjust(5, ' ')+' '+(str(int(txTime))).rjust(2, ' ')+' '+(str(int(rxTime))).rjust(2, ' ')+' '+' '+('{:5.3f}'.format(CMaxDensity)).rjust(8, ' ')+' '+('{:.3f}'.format(CMinSafeDistance)).rjust(7, ' ')+' | '+('{:5.3f}'.format(CMinDistance)).rjust(7, ' ')+' '+('{:5.3f}'.format(UMaxDensity)).rjust(8, ' ')+' '+('{:5.3f}'.format(UMinSafeDistance)).rjust(7, ' ')+' | '+('{:5.3f}'.format(UMinDistance)).rjust(7, ' ')) print() print(magentaDashLine) print(bottomStringRow01) print(bottomStringRow02) print(bottomStringRow03) print(bottomStringRow04) print(magentaDashLine)

Pure and simple Python Code:
Download - Contextual Click and save, 'ARRL_Based_RF-Exposure-Calculator_All-Bands_5-10-25-50-75-100-250-500-750-1000-1200-1500W.03.py' Source code here.

Python Code using 'termcolor' Module:
Download - Contextual Click and save, 'ARRL_Based_RF-Exposure-Calculator_w|Termcolor_All-Bands_5-10-25-50-75-100-250-500-750-1000-1200-1500W.03.py' Source code here.

Return to KF3G