% pxmas1 plots results of the interference fit for Dec. 25, 1996.
% RWB, June 20, 1998.
   nav=20;
   load ld98548a; % Dec. 25, with results of fit WITH the term in phi
		  % due to the velocity profile.
   lda=ld98548; % transfer structure of results.
   load ld98548b; % Dec. 25, WITHOUT effects of vel. profile.
   ldb=ld98548; % transfer second structure.
   %%%%%% Calculate normalized times.
   len1=length(lda.gmtsf);
   gmtsfn=(lda.gmtsf-lda.gmtsf(1))./(lda.gmtsf(len1)-lda.gmtsf(1));
   len2=length(lda.gmtus);
   gmtusn=(lda.gmtus-lda.gmtus(1))./(lda.gmtus(len2)-lda.gmtus(1));
   len3=length(lda.gmtfit);
   gmtfitn=(lda.gmtfit-lda.gmtfit(1))./(lda.gmtfit(len3)-lda.gmtfit(1));
   %%%%%% Average the rms array to cut down the noise.
   rmsaav=lda.rmsa;
   for i=1:len1-nav+1,rmsaav(i)=sum(rmsaav(i:i+nav-1))/nav;,end
   plot(gmtsfn,rmsaav,'black');
   hold on;
   plot(gmtusn,lda.hius,'-green');
   plot(gmtfitn,lda.rms+1,'blue');
   plot(gmtfitn,lda.vfit+1,'green');
   plot(gmtfitn,ldb.rms+10.5,'red');

%   plot([.25,.29],[.7,.7],'green')
   text(.3,.7,'current','color','green','fontsize',14);
   text(.3,3.5,'with velocity profile','color','blue','fontsize',14);
   text(.3,9.7,'DATA','color','black','fontsize',14);
   text(.3,13,'without velocity profile','color','red','fontsize',14);
   text(.3,14.5,'tide stage','color','green','fontsize',14);

   title('Interference Between Direct and Reflected Rays', ...
      'fontsize',20)
   xlabel('one lunar day, lunar noon at 0h Dec 25 1996 PST', ...
      'fontsize',20)


   hold off;

   print -djpeg -r85 ld98548.jpg