______________________________________________________________________________ PR PROFILE is a module which operates on a map of an extended object such as a galaxy or an oversampled star, and models its surface brightness, central coordinates, inclination, and eccentricity as a function of distance from its center outward along the major axis. The output is written as a character file of type 1. Options exist for extending the profile to large distances away from the object by averaging some of the parameters as the solution procedes. One can extend the profile even further by fitting a polynomial to the run of parameters with radius. Interaction by the user with displays of the parameters as a function of radius enable one to take maximum advantage of this powerful option. ______________________________________________________________________________ Enter an option: 14-NOV-1991 14:44:52.75 OPTION=PR PROFILE, Version 2 Open the SAD file to be measured Enter the filename: (<CR> for default) FILE= GALAXY ....A default was taken here There is one map in the file Reading from the file 200 rows of 200 pixels in R4 format Enter a command: F = Coordinates, eccentricity, and inclination from file prepared by option CT M = Manual input of coordinates COMMAND=M ....Manual input was selected since no Grinnell was available Enter the column and row: (default = last position, fractional values allowed) POSITION=101 102 ....An educated guess Enter the estimated inclination: (default = 0.0 degrees) INCLINATION= 0.000000E+00 ....The default was taken Enter the estimated eccentricity: (default = 0.0) ECCENTRICITY= 0.000000E+00 ....The default was taken Enter the outer radius: RADIUS=20 ....Good signal-to-noise exists out to here Enter the sky level to be subtracted before integrating: (default = minimum) SKY= -999.999 ....The sky isn't known accurately at this time, so a default was taken. Enter the object's ID: (default = last entry) ID=30 Do you wish to add the results to an existing file? (default = N) OLD= F Enter the file name: (default = last entry) FILE_22=30 ....It is convenient to name the file the same as the ID. Later, all final results can be combined into one file. Enter the fractional increment: (default = 1.0) INCREMENT= 1.00000 ....The default is taken at the start Do you wish to solve for the center of each contour? (default = Y) SOLVE= T ....When the derivative of a contour is greater than or equal to zero, the solution has failed. This is because the signal-to-noise is to low or the fractional increment is too large. Derivative of contour 20 >= 0 Average of contours 5 through 10 Column Row Inclination Eccentricity 101.865 102.592 16.71 0.737 Derivative of contour 15 >= 0 Derivative of contour 19 >= 0 102.696 102.830 6.22 0.672 Derivative of contour 19 >= 0 102.736 102.971 4.42 0.656 Derivative of contour 19 >= 0 102.744 103.004 4.51 0.655 Derivative of contour 18 >= 0 102.743 103.011 4.59 0.654 102.743 103.012 4.64 0.654 102.743 103.011 4.66 0.654 Derivative of contour 20 >= 0 102.743 103.011 4.67 0.654 Derivative of contour 18 >= 0 102.743 103.010 4.67 0.654 Derivative of contour 19 >= 0 102.743 103.010 4.68 0.653 Derivative of contour 18 >= 0 Enter a choice: Y = Iterate 10 more times N = Do not iterate (default) CHOICE=Y ....We'll try again with a smaller fractional increment Enter the fractional increment: (default = 1.0) INCREMENT=.6 Do you wish to solve for the center of each contour? (default = Y) SOLVE= T Derivative of contour 18 >= 0 Average of contours 5 through 10 Column Row Inclination Eccentricity 102.743 103.010 4.68 0.653 102.743 103.010 4.68 0.653 Derivative of contour 19 >= 0 102.743 103.010 4.68 0.653 Derivative of contour 19 >= 0 Derivative of contour 20 >= 0 102.743 103.010 4.68 0.653 102.743 103.010 4.68 0.653 Derivative of contour 19 >= 0 102.743 103.010 4.68 0.653 Derivative of contour 19 >= 0 Derivative of contour 20 >= 0 102.743 103.010 4.68 0.653 Derivative of contour 20 >= 0 102.743 103.010 4.68 0.653 Derivative of contour 18 >= 0 102.743 103.010 4.68 0.653 102.743 103.010 4.68 0.653 Derivative of contour 19 >= 0 Enter a choice: Y = Iterate 10 more times N = Do not iterate (default) CHOICE= N ....It looks like the signal-to-noise is too high for contours 18 to 20, so we'll continue by averaging parameters in the solution form contours 17 on out. Do you wish to extend the profile? (default = N) EXTEND=Y Enter the inner and outer radii for the extension: RADII=17 30 ....We'll see if we can get out as far as contour 30 Enter a choice: Y = Iterate 10 more times N = Do not iterate (default) A = Iterate, but write single averages for Xc, Yc, inclination, and eccentricity CHOICE=Y Enter the fractional increment: (default = 1.0) INCREMENT= 1.00000 Do you wish to solve for the center of each contour? (default = Y) SOLVE= T Average of contours 17 through 22 Column Row Inclination Eccentricity 102.785 102.046 21.84 0.648 102.232 102.027 40.19 0.593 102.694 101.443 22.19 0.566 Derivative of contour 29 >= 0 Derivative of contour 30 >= 0 102.908 102.074 26.82 0.589 Derivative of contour 30 >= 0 102.543 102.182 27.31 0.485 102.393 102.082 15.36 0.513 Derivative of contour 28 >= 0 Derivative of contour 29 >= 0 Derivative of contour 30 >= 0 102.136 102.053 28.05 0.505 Derivative of contour 27 >= 0 Derivative of contour 28 >= 0 Derivative of contour 29 >= 0 Derivative of contour 30 >= 0 102.393 101.887 20.22 0.511 102.374 101.947 28.53 0.576 102.506 101.765 24.40 0.603 Enter a choice: Y = Iterate 10 more times N = Do not iterate (default) A = Iterate, but write single averages for Xc, Yc, inclination, and eccentricity CHOICE=Y Enter the fractional increment: (default = 1.0) INCREMENT=.5 ....We'll try decreasing this Do you wish to solve for the center of each contour? (default = Y) SOLVE= T Average of contours 17 through 22 Column Row Inclination Eccentricity 102.851 101.840 21.59 0.608 102.686 101.928 30.94 0.614 102.926 101.921 23.77 0.621 102.852 101.915 31.79 0.623 103.082 101.904 26.08 0.639 Derivative of contour 30 >= 0 102.980 101.918 27.85 0.646 Derivative of contour 28 >= 0 Derivative of contour 29 >= 0 Derivative of contour 30 >= 0 102.998 101.828 26.50 0.658 103.055 101.953 24.89 0.657 103.109 102.041 24.88 0.661 Derivative of contour 30 >= 0 103.065 102.097 24.94 0.670 ....The outer contours have too low a signal. Now we'll extend outward by fitting a polynomial to an input range of radii, and use those values to assign values to radii at large distances from the center. Enter a choice: Y = Iterate 10 more times N = Do not iterate (default) A = Iterate, but write single averages for Xc, Yc, inclination, and eccentricity CHOICE= N ....The program writes the fourier coefficients to the log file. These values do not appear on the screen. ID K Cos(4E) Sin(4E) Cos(2E) Sin(2E) Cos(E) Sin(E) Deriv 30 0 0.0000 0.0000 0.000 0.000 0.000 0.000 0.00 30 1 -2.1369 0.0000 0.000 0.000 0.000 0.000 -43.01 30 2 -0.9630 -4.5818 0.000 0.000 0.000 0.000 -107.08 30 3 -1.8790 1.2882 0.000 0.000 0.000 0.000 -107.36 30 4 -0.6090 2.2229 0.000 0.000 0.000 0.000 -72.68 30 5 2.7449 -0.4036 0.000 0.000 0.000 0.000 -56.20 30 6 0.3621 -0.8960 0.000 0.000 0.000 0.000 -39.37 30 7 0.5195 0.5935 0.000 0.000 0.000 0.000 -26.33 30 8 0.2962 1.2180 0.000 0.000 0.000 0.000 -17.33 30 9 -0.3984 2.1552 0.000 0.000 0.000 0.000 -12.24 30 10 1.2634 -0.1840 0.000 0.000 0.000 -0.001 -13.00 30 11 1.9597 0.5618 0.000 0.000 0.000 0.000 -6.98 30 12 -0.2242 0.1959 0.000 0.000 0.000 0.000 -4.29 30 13 -0.2622 -0.3169 0.000 0.000 0.000 0.000 -5.42 30 14 0.8146 -0.5501 0.000 0.000 0.000 0.000 -2.51 30 15 -0.2954 -0.9485 -0.001 0.000 0.000 -0.002 -3.00 30 16 -0.4682 -0.2498 -0.360 0.122 -0.182 0.495 -1.00 30 17 0.5661 -0.5060 -0.001 0.003 0.001 -0.001 -0.99 30 18 -1.0603 0.3686 0.656 0.035 -0.095 -0.797 -0.91 30 19 -0.1884 1.1182 -0.408 0.768 -0.722 0.791 -0.83 30 20 -1.5677 1.0231 -0.015 -0.033 -0.019 0.010 -0.76 30 21 -1.4268 1.4536 0.381 -0.224 -0.282 -0.228 -0.68 30 22 -0.0939 -0.1074 -0.311 -0.960 0.520 0.592 -0.61 30 23 -0.0932 -0.2462 0.122 -0.004 -0.040 -0.014 -0.54 30 24 0.0433 0.4375 0.352 0.099 -0.081 -0.174 -0.46 30 25 0.5516 0.3167 -0.249 0.783 -0.317 0.330 -0.39 30 26 -0.0505 0.2691 -0.140 0.086 0.359 -0.542 -0.31 30 27 -0.2233 -0.5427 -0.575 -1.617 0.239 -0.235 -0.24 30 28 0.3166 -0.0811 0.105 -1.202 -0.850 -0.673 -0.17 30 29 0.2280 0.7631 0.214 0.938 1.474 -0.396 -0.09 30 30 0.3042 -0.1018 0.738 -0.793 -0.550 0.387 -0.02 ....In order to select an input range, we must display the solution for X, Y, inclination, and eccentricity as a graph. Enter one or more commands: 1 = Displays average ADU vs. radius (default) 2 = Displays log(average ADU) vs. log(radius) 3 = Displays log(average ADU) vs. 1/4 power of radius 4 = Displays log(integrated average ADU) vs. log(radius) 5 = Displays X of center vs. radius 6 = Displays Y of center vs. radius 7 = Displays inclination vs. radius 8 = Displays eccentricity vs. radius 9 = Quits COMMANDS=5:8 Enter the ID of the object: (default = last entry) ID= 30 The minimum and maximum radii are 0.000 and 30.000. Enter the start radius and end radius for plotting: (default = entire range) RADII= 0.000000E+00 30.0000 Enter the plot devices: G = Graphics terminal (VT100/Retrographics and Visual 550, default) G1 = Graphics terminal (GraphOn) G2 = Graphics terminal (Codonics) G3 = Graphics terminal (Micro-term) F = File T = Tektronix 4662 flat bed plotter (1200 baud) DEVICES= G ....These same prompts are repeated for each of the four plots. . . . More objects to plot? (default = N) MORE= F Do you wish to smooth the profile parameters? (default = Y) SMOOTH= T Enter the degrees for smoothing the parameters Xc, Yc, inclination, and eccentricity: (default = 0 for all) Note: 0 = simple averaging, negative integer = no smoothing DEGREES= 0 0 0 0 ....We have selected simple averaging for each of the four parameters since no trend was visible in the plots. If any of the four outer output radii require an extrapolation, the largest of them will enforce that extrapolation on all four parameters. Enter the input and output pairs of radii for the smoothing of Xc: (default = previous entry range) RANGE=27 30 27 100 ....The input range is 27 to 30, and the output range is 27 to 100 The outer output radius will be assumed to have this value You can skip giving its value if you do not want a greater extension Enter the input and output pairs of radii for the smoothing of Yc: (default = previous entry range) RANGE= 27.0000 ....A default was taken 30.0000 27.0000 100.000 Enter the input and output pairs of radii for the smoothing of Inclination: (default = previous entry range) RANGE=26 30 26 100 Enter the input and output pairs of radii for the smoothing of Eccentricity: (default = previous entry range) RANGE=26 29 26 100 ....Contour number 30 was not as good for the eccentricity. Enter the sky level to be subtracted before integrating: (default = minimum) SKY= -999.999 ....We don't know the correct sky yet. Coefficients of smoothing functions Parameter Degree of coefficient Coefficient Xc 0 100.9878 Yc 0 102.9638 Inclination 0 32.7580 Eccentricity 0 0.6250 Enter one or more commands: 1 = Displays average ADU vs. radius (default) 2 = Displays log(average ADU) vs. log(radius) 3 = Displays log(average ADU) vs. 1/4 power of radius 4 = Displays log(integrated average ADU) vs. log(radius) 5 = Displays X of center vs. radius 6 = Displays Y of center vs. radius 7 = Displays inclination vs. radius 8 = Displays eccentricity vs. radius 9 = Quits COMMANDS= 1 ....Now we will look at the surface brightness distribution and gauge the sky level Enter the ID of the object: (default = last entry) ID= 30 The minimum and maximum radii are 0.000 and 100.000. Enter the start radius and end radius for plotting: (default = entire range) RADII= 0.000000E+00 100.000 Enter the plot devices: G = Graphics terminal (VT100/Retrographics and Visual 550, default) G1 = Graphics terminal (GraphOn) G2 = Graphics terminal (Codonics) G3 = Graphics terminal (Micro-term) F = File T = Tektronix 4662 flat bed plotter (1200 baud) DEVICES= G ....We use command R when the plot is displayed on a graphics terminal in order to magnify the plot so as to obtain an accurate sky level. This level is used later when writing the character file back to the map in order to reconstruct the galaxy. Here we find that the sky is 400.040 units. We also find that the data beyond a radius of 73 pixels along the major axis is not as good and so will be edited out of the character file. More objects to plot? (default = N) MORE= F Do you wish to smooth the profile parameters? (default = Y) SMOOTH=N ....We have already done this. ______________________________________________________________________________