runbulgerot [parameter=value]
This program calculates the rotation curve for an oblate spheroidal bulge. It assumes a projected surface density that follows a Sersic profile. The central surface magnitude M_0, characteristic radius R_0 and Sersic index n_sers must be given. Also, the inclination and intrinsic axis ratio of the bulge must be given. Finally, the distance and the Galactic foreground extinction of the galaxy is needed.
% runbulgerot out=sersic1.tab % rotcurves rotcur 0,1,1,5,1 sersec1.tab radii=0:2:0.01 mode=omega yrange=0,400
The original version of bulgerot had a hardcoded set of input parameters
(see BENCHMARK below), but the modified version for NEMO (and runbulgerot)
reads them from stdin in a certain order, viz. 
% bulgerot << EOF_BULGEROT bulgerot.dat # log sersic:n=1 # galaxy 15.0d0 1.0d0 1.0d0 # mu0, r0, n 0.0d0 1.0d0 10.0d6 0.0d0 # inc, q, dist, ar 0.0d0 0.01d0 250 # rstart, rstep, nradii EOF_BULGEROT
#    radius     intensity      density         mass        velocity   
#    (kpc)      (kg/m^2)       (kg/m^3)       (Msun)        (km/s)     
    2.490000  0.224088D-20  0.827590D-40  0.323728D+09  0.235621D+02
    2.500000  0.182313D-20  0.671960D-40  0.323728D+09  0.235148D+02
@inferno:   219.102u 0.000s 3:39.28 99.9%     gcc      3.4.6  i7-3820 @ 3.60GHz
@inferno:   167.281u 0.000s 2:47.42 99.9%     gfortran 4.4.7  i7-3820  @ 3.60GHz
@nemo2:      77.328u 0.004s 1:17.50 99.7%     gfortran 4.7.3  i7-3630QM @ 2.40GHz
@xps13:      43.85u  0.00s  0:43.85 100%      gfortran 9.3.0  i5-1135G7 @ 2.40GHz
https://ui.adsabs.harvard.edu/abs/2008MNRAS.385.1359N (Noordermeer, E. 2008, MNRAS, 385, 1359)
NEMO/usr/noordermeer/bulgerot/ bulgerot.f runbulgerot.c
2007? V1 original fortran release as bulgerot EN 11-sep-13 V1.1 NEMO interface PJT