I would suggest using our project SIMuRG ( https://simurg.iszf.irk.ru/; Yasyukevich et al., 2020)
Yasyukevich et al. SIMuRG: System for Ionosphere Monitoring and Research from GNSS // GPS solutions V. 24, 69, 2020. https://link.springer.com/article/10.1007%2Fs10291-020-00983-2.
You can use GPS-TEC analysis application Developed by Dr. Gopi Krishna Seemala.
Brief features of this application - GPS data analysis
Ability to batch process the input files (RINEX format etc) for example: all files of the month, year, all stations and all files in directory.
Gets ephemeris from IGS navigation file, has the ability to download the navigation file automatically if connected to internet; unless it finds the file in the same directory as data.
Calculates TEC from the observation data of GPS Rinex, Novatel (reads only ID43 records), SCINDA (.scn files) and Leica (for leica formats, file name convention is yet to fix, works with our convention names now).
Process cycle slips in phase data
Read satellite biases from DCB IGS code files, if not available calculates them.
Calculates the receiver bias
Calculates the inter-channel biases for different satellites in the receiver.
Plots the vertical TEC values on screen and writes ascii output files (*.CMN & *.STD) in the same directory of data file.
The quite reliable source could be the RINEX GPS-TEC software as others already suggested to you. Just check the version you will obtain regarding automatic download of DCB files.
You may try GPS-TEC (Ver. 2.9.5, https://seemala.blogspot.com/). The software can automatically download DCBs to your local directory as long as you have an internet connection. However, you may not be able to process observation data in Rinex 3.03 currently.
You may also try GAMIT/GLOBK( Ver. 10.71, http://geoweb.mit.edu/gg/).