Error reading in Map Length from res files

[PlantDr430]

Hello,

I know this project is old, but I was using it to follow the procedures in a recent paper that looked at recombination rates in fungi. So far this have been running smoothly, however, I noticed that in the log files for my LDhot runs my Map Length wasn't being reported correctly, and I don't know if this is affecting some of my results.

The head of my res_file from LDhat looks like this:

Loci	Mean_rho	Median	L95	U95
-1.000	191.55638	191.56499	175.19789	208.53172
0.01	0.00008	0.00006	0.00001	0.00024
0.02	0.00008	0.00006	0.00001	0.00024
0.021	0.00008	0.00006	0.00001	0.00024
0.022	0.00008	0.00006	0.00001	0.00024
0.029	0.00008	0.00006	0.00001	0.00024
0.092	0.00008	0.00006	0.00001	0.00024
0.106	0.00008	0.00006	0.00001	0.00024

By following your examples the Map Length should be 191.5 4Ner. However, when running LDhot the log file sets it as 0.192 4Ner:

Reading res: /data/wyka/recombination/Correct_run/LDhat_results/LDhat_results/alingment_103-1/sum_rates_res_kbp.txt
	Read 531 map positions
	Map Covers positions 0.01 - 29.87 kb
	Map Length = 0.192 4Ner

I found where this is occurring in your tools.cpp script, but I unfortunately don't know C++ and am unsure why it is dividing my Map Length by 1,000. I also know if this Map Length is being used elsewhere in the LDhot script and if it is using the correct Map Length for calculations or if it is using the incorrect version.

I have attached all files for reproduction:
locs_kbp.txt
sites.txt
sum_rates_res_kbp.txt
theta_0.01_lk_new_lk.txt

As well as the LDhot log file and results file:
theta_0.01_lk_ldhot.hotspots.txt
theta_0.01_lk_ldhot.log