Altogether there are 4 psfiles: -rw-r--r-- 1 hsalo tati 664779 Sep 19 17:43 spiral_md500_2.ps -rw-r--r-- 1 hsalo tati 664779 Sep 19 17:43 spiral_md500_1.ps -rw-r--r-- 1 hsalo tati 664779 Sep 19 17:43 spiral_md333_2.ps -rw-r--r-- 1 hsalo tati 664779 Sep 19 17:43 spiral_md333_1.ps -- 2 different time intervals for Mdisk=0.333 and 0.5 isolated runs. The decomposions are for m=2 logarithmic spirals, using particles within 50 arcsecs from the center. In each plot the MIDDLE frame shows the amplitude vs. time and p. The UPPERMOST frame is the same using the inclination of the spiral with respect to a trailing tangential direction, (tan(inclination)=m/p). So p=0 corresponds to inclination = 90 degrees p=30 corresponds to inclination of 3.8 degrees (trailing) p=-30 to -3.8 degrees = 176.2 degrees (leading) etc. The LOWERMOST panel is the mean amplitude (over the time interval) vs. p. I have used equally spaced p values (dp=0.5) in this plot, so the near-radial spirals are rather poorly sampled (p=0.5 corresponds to 75 degrees, p=1 to 63, p=1.5 to 53 etc). So, this is just to get an idea. From the stored data I can produce similar plots for any m, time-interval, radius-interval, with desired values of p. So we can have a look at that during your visit. I hope I have not forgot any essential info.