/*
 * Simple calculations of half wavelengths of ham bands.
 *
 * Original soure code by Mike Markowski AB3AP
 * Thu Jun 28 07:01:26 EDT 2012

 * Source code modified by IU7RAL lenght in meters (instead of feet) - BandPlan IARU R1 adoption
 * 04/05/2023
*/

#include <stdio.h>

main() {
	printHalfwaves();
}

/*
 * Print ranges of half wavelengths for ecah ham band.
 */
printHalfwaves() {

	/* rw(1830., 1850.);       /* 160m */
	/* rw(3500., 3800.);       /* 80m */
	/* rw(5351.5, 5366.5);      /* 60m */
	rw(7000., 7200.);       /* 40m */
	rw(10100., 10150.);     /* 30m */
	rw(14000., 14350.);     /* 20m */
	rw(18069., 18168.);     /* 17m */
	rw(21000., 21450.);     /* 15m */
	rw(24890., 24990.);     /* 12m */
	rw(28000., 29700.);     /* 10m */
	
}

/*
 * For a given frequency range, calculate the half wavelength range and print
 * it.  In addition, print up to 4th multiples of each range up to the length
 * of 160m half wavelength.
 *
 * Comments are also printed out, assuming that the output will saved to a file,
 * and that file used by gnuplot for plotting.
 */
rw(double min_kHz, double max_kHz) {
	
	double lambda0_ft, lambda1_ft;
	double loFreq_MHz = 7; /* qui mettere la frequenza pił bassa in MHz della banda pił lunga che si desidera ascoltare */
	double lambdaMax_ft = 2 * 142.6464 / loFreq_MHz; /* Max wavelength in band. */
	int n;

	double qtr_ft = 142.6464 / loFreq_MHz / 2;
	printf("# %.3f to %.3f kHz, too short for %f MHz\n",
		min_kHz, max_kHz, loFreq_MHz);
	printf("%.3f 0\n%.3f 1\n%.3f 1\n%.3f 0\n\n",
		0., 0+(1e-3), qtr_ft, qtr_ft+(1e-3));
	
	n = 1; /* First multiple. */
	do {
		lambda0_ft = n * 142.6464 / (max_kHz * 1e-3);
		lambda1_ft = n * 142.6464 / (min_kHz * 1e-3);
		/* Print in format gnuplot expects. */
		printf("# %.3f to %.3f kHz, multiple %d\n",
			min_kHz, max_kHz, n);
		printf("%.3f 0\n%.3f 1\n%.3f 1\n%.3f 0\n\n",
			lambda0_ft-(1e-3), lambda0_ft,
			lambda1_ft, lambda1_ft+(1e-3));
		/* Prepare for next multiple. */
		n++;
	/* Change '5' in next line to max number of multiples to calculate. */
	} while (lambda1_ft < lambdaMax_ft && n < 5);
}