Use the multi_usno_runs script to calculate star/gal source densities at Ra,Dec sky positions. The local file SGV.values can be used to adjust the stellarity index and V magnitude cuts.

 

% star_dens  
Usage: star_dens filename file.out
arg1 = filename for list of [RAsex,DECsex,Name]
arg2 = Name of output file 

% cat list.1 
19:40:55.5817 +26:50:32.590 NGC6815
19:01:22.2533 +11:36:11.352 NGC6738
12:30:22.2533 +31:36:11.352 NGalPole1
12:30:22.2533 +21:36:11.352 NGalPole2

% cat SGV.values 
0.0 11.5 19.5 

% star_dens list.1 Results  

% cat Results  
19:40:55.5817   +26:50:32.590   1.7104149       NGC6815         21.74375
19:01:22.2533   +11:36:11.352   2.6158173       NGC6738         14.69318
12:30:22.2533   +31:36:11.352   84.1621399      NGalPole1       0.66765
12:30:22.2533   +21:36:11.352   82.8211365      NGalPole2       0.97960