Hi Morris, Here is the monotonic spline version of INIT1.f. All my changes to INIT1.f are commented with notes that begin with "c rayr:". The changes boil down to this: if we have at least two obs for a quadrant, don't subtract the Holland profile from the observations and just fit the observations with a monotonic spline, starting with the r=rmax,v=vmax point (not the r=0, v=0 central point). It may be that this will work fine with only 1 observation- I didn't investigate this case. The header routine for the monotonic spline is called akim2.f (a misnomer, I'm afraid) and it calls pchim and pchfe, that are both contained in the code monospline.f. The routine xerror.f is included because it is summoned from the monospline routines when something goes wrong. I got the monotonic spline stuff from Sergei Frolov, whom I believe scavenged it off of the web somewhere. You will need to add akim2.f, monospline.f and xerror.f to the axisymet makefile to make this work. Good luck and feel free to call me (401-788-9031) if it doesn't work! -Ray