diff tools/rdock/data/pmf/smoothed/smoother.m @ 0:bc03dbb6eb37 draft

planemo upload commit 781926e52355f7805db8d9a4ccafeff397b19aa4-dirty

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tools/rdock/data/pmf/smoothed/smoother.m	Mon Aug 29 03:38:13 2016 -0400
@@ -0,0 +1,2310 @@
+## natural-cubic-spline interpolation
+## usage: yspline = spline(x,y,xspline)
+## example: 
+##	x = 0:10; y = sin(x);
+##	xspline = 0:0.1:10; yspline = spline(x,y,xspline);
+##	##plot(x,y,"+",xspline,yspline);
+## Given the vectors x and y, which tabulate a function, with 
+## x(1) < x(2) < x(3) <... or x(1) > x(2) > x(3) >..., and given
+## the vector xspline, this function returns a natural-cubic-spline 
+## interpolated vector yspline. 
+## author: Zdenek Remes, May 22, 1999
+
+
+function ynew = spline(x,y,xnew)
+[x,index]=sort(x);
+y=y(index);
+n=length(y);
+y2(1)=0.0;
+y2(n)=0.0;
+u(1)=0.0;
+for i=2:n-1
+  sig=(x(i)-x(i-1))/(x(i+1)-x(i-1));
+  p=sig*y2(i-1)+2.0;
+  y2(i)=(sig-1.0)/p;
+  u(i)=(y(i+1)-y(i))/(x(i+1)-x(i))-(y(i)-y(i-1))/(x(i)-x(i-1));
+  u(i)=(6.0*u(i)/(x(i+1)-x(i-1))-sig*u(i-1))/p;
+endfor;
+k=n-1;
+while (k >= 1)
+  y2(k)=y2(k)*y2(k+1)+u(k);
+  k--;
+endwhile;
+
+i1=1; in=length(xnew);
+
+#if (xnew(1) < x(1)) 
+#  error("spline: bad xspline");
+#endif;
+#if (xnew(in) > x(n)) 
+#  error("spline: bad xspline");
+#endif;
+  
+if (xnew(1) == x(1)) 
+  ynew(1)=y(1);
+  i1=2;
+endif;
+if (xnew(in) == x(n)) 
+  ynew(in)=y(n);
+  in=in-1; 
+endif;
+  
+
+for i=i1:in  
+  khi=n;
+  klo=1;
+  while ((khi-klo) > 1)
+    k=floor((khi+klo)/2);
+    if (x(k) > xnew(i)) 
+      khi=k;
+    else
+      klo=k;  
+    endif;
+  endwhile;
+  h=x(khi)-x(klo);
+  a=(x(khi)-xnew(i))/h;
+  b=(xnew(i)-x(klo))/h;
+  ynew(i)=a*y(klo)+b*y(khi)+((a^3-a)*y2(klo)+(b^3-b)*y2(khi))*(h*h)/6.0;
+endfor;
+endfunction;
+
+## Savitzky-Golay smoothing filter
+## usage: [xsavgol,ysavgol]=savgol(x,y,nl,nr,m)
+## example: x=0:0.01:3;y1=sin(x.^3);y=y1+(rand(1,301)-0.5)/3;
+##	    [xsavgol,ysavgol]=savgol(x,y,10,10,2);
+##	    ##plot(x,y,"+",xsavgol,ysavgol,x,y1)
+## Given vectors x, y containing a tabulated data y=f(x) with 
+## equally spaced x's this function calculates smoothed data
+## ysavgol=g(xsavgol) by Savitzky-Golay smoothing filter.
+## nl is the number of leftward (past) data points used, while
+## nr is the number of rightward (future) data points, making
+## the total number of data points used nl+nr+1. m is the order
+## of the smoothing polynomial, also equal to the highest
+## conserved moment; usual values are m=2 or m=4.
+## The idea of Savitzky-Golay filtering is to smooth the
+## underlying data y=f(x) within the moving window not by a
+## constant (whose estimate is the average), but by a poly-
+## nomial of higher order. Thus for a point y(i) the function
+## savgol fits by a least-squares method a polynomial to
+## points y(i-nl), ..., y(i+nr) in the moving window, and
+## then set g(i-nl+1) to the value of that polynomial at
+## position x(i). 
+## Zdenek Remes, Mai 22, 1999
+	
+function [xnew,ynew]=savgol(x,y,nl, nr, M)
+
+    if max(diff(x,2))>100*eps 
+	error("The x's must be equally spaced.")
+    endif
+    for i=-nl:nr
+	for j=0:M
+	    A(i+nl+1,j+1)=i^j;
+	endfor
+    endfor
+    AA=inv(A'*A);
+    for i=-nl:nr
+	cc=0;
+	for m=0:M
+	    cc=cc+AA(1,m+1)*i^m;
+	endfor
+    c(i+nl+1)=cc;
+    endfor
+    
+    nx=length(x);
+    for i=nl:nx-nr-1
+	yy=0;
+	for j=-nl:nr
+	    yy=yy+c(j+nl+1)*y(i+j+1);
+	endfor
+        xnew(i-nl+1)=x(i+1);
+	ynew(i-nl+1)=yy;
+    endfor    
+endfunction
+
+## x=0:0.01:3;
+## y1=sin(x.^3);
+## y=y1+(rand(1,301)-0.5)/3;
+## [xsavgol,ysavgol]=savgol(x,y,10,10,2);
+## plot(x,y,"+",xsavgol,ysavgol,x,y1);
+## pause;
+zero=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ];
+##load -force CPNA.dat;
+##[xs,ys]=savgol(x,y,5,5,4);
+##save -ascii CPNA.sm xs ys;
+
+load -force CFBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFBr.sm xs ys;
+printf("CFBr.sm\n");
+load -force CFC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFC0.sm xs ys;
+printf("CFC0.sm\n");
+load -force CFC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFC3.sm xs ys;
+printf("CFC3.sm\n");
+load -force CFCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFCF.sm xs ys;
+printf("CFCF.sm\n");
+load -force CFCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFCL.sm xs ys;
+printf("CFCL.sm\n");
+load -force CFCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFCN.sm xs ys;
+printf("CFCN.sm\n");
+load -force CFCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFCO.sm xs ys;
+printf("CFCO.sm\n");
+load -force CFCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFCP.sm xs ys;
+printf("CFCP.sm\n");
+load -force CFCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFCW.sm xs ys;
+printf("CFCW.sm\n");
+load -force CFF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFF_.sm xs ys;
+printf("CFF_.sm\n");
+load -force CFFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFFe.sm xs ys;
+printf("CFFe.sm\n");
+load -force CFHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFHL.sm xs ys;
+printf("CFHL.sm\n");
+load -force CFMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFMg.sm xs ys;
+printf("CFMg.sm\n");
+load -force CFMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFMn.sm xs ys;
+printf("CFMn.sm\n");
+load -force CFN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFN0.sm xs ys;
+printf("CFN0.sm\n");
+load -force CFNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFNA.sm xs ys;
+printf("CFNA.sm\n");
+load -force CFNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFNC.sm xs ys;
+printf("CFNC.sm\n");
+load -force CFND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFND.sm xs ys;
+printf("CFND.sm\n");
+load -force CFNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFNP.sm xs ys;
+printf("CFNP.sm\n");
+load -force CFNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFNR.sm xs ys;
+printf("CFNR.sm\n");
+load -force CFNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFNS.sm xs ys;
+printf("CFNS.sm\n");
+load -force CFOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFOA.sm xs ys;
+printf("CFOA.sm\n");
+load -force CFOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFOC.sm xs ys;
+printf("CFOC.sm\n");
+load -force CFOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFOD.sm xs ys;
+printf("CFOD.sm\n");
+load -force CFOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFOE.sm xs ys;
+printf("CFOE.sm\n");
+load -force CFOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFOR.sm xs ys;
+printf("CFOR.sm\n");
+load -force CFOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFOS.sm xs ys;
+printf("CFOS.sm\n");
+load -force CFP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFP_.sm xs ys;
+printf("CFP_.sm\n");
+load -force CFSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFSA.sm xs ys;
+printf("CFSA.sm\n");
+load -force CFSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFSD.sm xs ys;
+printf("CFSD.sm\n");
+load -force CFV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFV_.sm xs ys;
+printf("CFV_.sm\n");
+load -force CFZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFZn.sm xs ys;
+printf("CFZn.sm\n");
+load -force CFcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFcF.sm xs ys;
+printf("CFcF.sm\n");
+load -force CFcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CFcP.sm xs ys;
+printf("CFcP.sm\n");
+load -force CNBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNBr.sm xs ys;
+printf("CNBr.sm\n");
+load -force CNC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNC0.sm xs ys;
+printf("CNC0.sm\n");
+load -force CNC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNC3.sm xs ys;
+printf("CNC3.sm\n");
+load -force CNCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNCF.sm xs ys;
+printf("CNCF.sm\n");
+load -force CNCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNCL.sm xs ys;
+printf("CNCL.sm\n");
+load -force CNCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNCN.sm xs ys;
+printf("CNCN.sm\n");
+load -force CNCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNCO.sm xs ys;
+printf("CNCO.sm\n");
+load -force CNCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNCP.sm xs ys;
+printf("CNCP.sm\n");
+load -force CNCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNCW.sm xs ys;
+printf("CNCW.sm\n");
+load -force CNF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNF_.sm xs ys;
+printf("CNF_.sm\n");
+load -force CNFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNFe.sm xs ys;
+printf("CNFe.sm\n");
+load -force CNHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNHL.sm xs ys;
+printf("CNHL.sm\n");
+load -force CNMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNMg.sm xs ys;
+printf("CNMg.sm\n");
+load -force CNMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNMn.sm xs ys;
+printf("CNMn.sm\n");
+load -force CNN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNN0.sm xs ys;
+printf("CNN0.sm\n");
+load -force CNNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNNA.sm xs ys;
+printf("CNNA.sm\n");
+load -force CNNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNNC.sm xs ys;
+printf("CNNC.sm\n");
+load -force CNND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNND.sm xs ys;
+printf("CNND.sm\n");
+load -force CNNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNNP.sm xs ys;
+printf("CNNP.sm\n");
+load -force CNNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNNR.sm xs ys;
+printf("CNNR.sm\n");
+load -force CNNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNNS.sm xs ys;
+printf("CNNS.sm\n");
+load -force CNOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNOA.sm xs ys;
+printf("CNOA.sm\n");
+load -force CNOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNOC.sm xs ys;
+printf("CNOC.sm\n");
+load -force CNOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNOD.sm xs ys;
+printf("CNOD.sm\n");
+load -force CNOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNOE.sm xs ys;
+printf("CNOE.sm\n");
+load -force CNOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNOR.sm xs ys;
+printf("CNOR.sm\n");
+load -force CNOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNOS.sm xs ys;
+printf("CNOS.sm\n");
+load -force CNP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNP_.sm xs ys;
+printf("CNP_.sm\n");
+load -force CNSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNSA.sm xs ys;
+printf("CNSA.sm\n");
+load -force CNSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNSD.sm xs ys;
+printf("CNSD.sm\n");
+load -force CNV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNV_.sm xs ys;
+printf("CNV_.sm\n");
+load -force CNZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNZn.sm xs ys;
+printf("CNZn.sm\n");
+load -force CNcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNcF.sm xs ys;
+printf("CNcF.sm\n");
+load -force CNcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CNcP.sm xs ys;
+printf("CNcP.sm\n");
+load -force COBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COBr.sm xs ys;
+printf("COBr.sm\n");
+load -force COC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COC0.sm xs ys;
+printf("COC0.sm\n");
+load -force COC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COC3.sm xs ys;
+printf("COC3.sm\n");
+load -force COCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COCF.sm xs ys;
+printf("COCF.sm\n");
+load -force COCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COCL.sm xs ys;
+printf("COCL.sm\n");
+load -force COCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COCN.sm xs ys;
+printf("COCN.sm\n");
+load -force COCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COCO.sm xs ys;
+printf("COCO.sm\n");
+load -force COCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COCP.sm xs ys;
+printf("COCP.sm\n");
+load -force COCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COCW.sm xs ys;
+printf("COCW.sm\n");
+load -force COF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COF_.sm xs ys;
+printf("COF_.sm\n");
+load -force COFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COFe.sm xs ys;
+printf("COFe.sm\n");
+load -force COHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COHL.sm xs ys;
+printf("COHL.sm\n");
+load -force COMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COMg.sm xs ys;
+printf("COMg.sm\n");
+load -force COMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COMn.sm xs ys;
+printf("COMn.sm\n");
+load -force CON0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CON0.sm xs ys;
+printf("CON0.sm\n");
+load -force CONA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CONA.sm xs ys;
+printf("CONA.sm\n");
+load -force CONC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CONC.sm xs ys;
+printf("CONC.sm\n");
+load -force COND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COND.sm xs ys;
+printf("COND.sm\n");
+load -force CONP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CONP.sm xs ys;
+printf("CONP.sm\n");
+load -force CONR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CONR.sm xs ys;
+printf("CONR.sm\n");
+load -force CONS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CONS.sm xs ys;
+printf("CONS.sm\n");
+load -force COOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COOA.sm xs ys;
+printf("COOA.sm\n");
+load -force COOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COOC.sm xs ys;
+printf("COOC.sm\n");
+load -force COOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COOD.sm xs ys;
+printf("COOD.sm\n");
+load -force COOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COOE.sm xs ys;
+printf("COOE.sm\n");
+load -force COOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COOR.sm xs ys;
+printf("COOR.sm\n");
+load -force COOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COOS.sm xs ys;
+printf("COOS.sm\n");
+load -force COP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COP_.sm xs ys;
+printf("COP_.sm\n");
+load -force COSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COSA.sm xs ys;
+printf("COSA.sm\n");
+load -force COSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COSD.sm xs ys;
+printf("COSD.sm\n");
+load -force COV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COV_.sm xs ys;
+printf("COV_.sm\n");
+load -force COZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COZn.sm xs ys;
+printf("COZn.sm\n");
+load -force COcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COcF.sm xs ys;
+printf("COcF.sm\n");
+load -force COcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii COcP.sm xs ys;
+printf("COcP.sm\n");
+load -force CPBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPBr.sm xs ys;
+printf("CPBr.sm\n");
+load -force CPC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPC0.sm xs ys;
+printf("CPC0.sm\n");
+load -force CPC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPC3.sm xs ys;
+printf("CPC3.sm\n");
+load -force CPCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPCF.sm xs ys;
+printf("CPCF.sm\n");
+load -force CPCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPCL.sm xs ys;
+printf("CPCL.sm\n");
+load -force CPCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPCN.sm xs ys;
+printf("CPCN.sm\n");
+load -force CPCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPCO.sm xs ys;
+printf("CPCO.sm\n");
+load -force CPCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPCP.sm xs ys;
+printf("CPCP.sm\n");
+load -force CPCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPCW.sm xs ys;
+printf("CPCW.sm\n");
+load -force CPF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPF_.sm xs ys;
+printf("CPF_.sm\n");
+load -force CPFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPFe.sm xs ys;
+printf("CPFe.sm\n");
+load -force CPHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPHL.sm xs ys;
+printf("CPHL.sm\n");
+load -force CPMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPMg.sm xs ys;
+printf("CPMg.sm\n");
+load -force CPMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPMn.sm xs ys;
+printf("CPMn.sm\n");
+load -force CPN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPN0.sm xs ys;
+printf("CPN0.sm\n");
+load -force CPNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPNA.sm xs ys;
+printf("CPNA.sm\n");
+load -force CPNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPNC.sm xs ys;
+printf("CPNC.sm\n");
+load -force CPND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPND.sm xs ys;
+printf("CPND.sm\n");
+load -force CPNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPNP.sm xs ys;
+printf("CPNP.sm\n");
+load -force CPNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPNR.sm xs ys;
+printf("CPNR.sm\n");
+load -force CPNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPNS.sm xs ys;
+printf("CPNS.sm\n");
+load -force CPOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPOA.sm xs ys;
+printf("CPOA.sm\n");
+load -force CPOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPOC.sm xs ys;
+printf("CPOC.sm\n");
+load -force CPOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPOD.sm xs ys;
+printf("CPOD.sm\n");
+load -force CPOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPOE.sm xs ys;
+printf("CPOE.sm\n");
+load -force CPOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPOR.sm xs ys;
+printf("CPOR.sm\n");
+load -force CPOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPOS.sm xs ys;
+printf("CPOS.sm\n");
+load -force CPP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPP_.sm xs ys;
+printf("CPP_.sm\n");
+load -force CPSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPSA.sm xs ys;
+printf("CPSA.sm\n");
+load -force CPSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPSD.sm xs ys;
+printf("CPSD.sm\n");
+load -force CPV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPV_.sm xs ys;
+printf("CPV_.sm\n");
+load -force CPZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPZn.sm xs ys;
+printf("CPZn.sm\n");
+load -force CPcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPcF.sm xs ys;
+printf("CPcF.sm\n");
+load -force CPcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii CPcP.sm xs ys;
+printf("CPcP.sm\n");
+load -force HHBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHBr.sm xs ys;
+printf("HHBr.sm\n");
+load -force HHC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHC0.sm xs ys;
+printf("HHC0.sm\n");
+load -force HHC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHC3.sm xs ys;
+printf("HHC3.sm\n");
+load -force HHCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHCF.sm xs ys;
+printf("HHCF.sm\n");
+load -force HHCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHCL.sm xs ys;
+printf("HHCL.sm\n");
+load -force HHCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHCN.sm xs ys;
+printf("HHCN.sm\n");
+load -force HHCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHCO.sm xs ys;
+printf("HHCO.sm\n");
+load -force HHCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHCP.sm xs ys;
+printf("HHCP.sm\n");
+load -force HHCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHCW.sm xs ys;
+printf("HHCW.sm\n");
+load -force HHF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHF_.sm xs ys;
+printf("HHF_.sm\n");
+load -force HHFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHFe.sm xs ys;
+printf("HHFe.sm\n");
+load -force HHHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHHL.sm xs ys;
+printf("HHHL.sm\n");
+load -force HHMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHMg.sm xs ys;
+printf("HHMg.sm\n");
+load -force HHMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHMn.sm xs ys;
+printf("HHMn.sm\n");
+load -force HHN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHN0.sm xs ys;
+printf("HHN0.sm\n");
+load -force HHNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHNA.sm xs ys;
+printf("HHNA.sm\n");
+load -force HHNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHNC.sm xs ys;
+printf("HHNC.sm\n");
+load -force HHND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHND.sm xs ys;
+printf("HHND.sm\n");
+load -force HHNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHNP.sm xs ys;
+printf("HHNP.sm\n");
+load -force HHNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHNR.sm xs ys;
+printf("HHNR.sm\n");
+load -force HHNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHNS.sm xs ys;
+printf("HHNS.sm\n");
+load -force HHOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHOA.sm xs ys;
+printf("HHOA.sm\n");
+load -force HHOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHOC.sm xs ys;
+printf("HHOC.sm\n");
+load -force HHOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHOD.sm xs ys;
+printf("HHOD.sm\n");
+load -force HHOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHOE.sm xs ys;
+printf("HHOE.sm\n");
+load -force HHOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHOR.sm xs ys;
+printf("HHOR.sm\n");
+load -force HHOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHOS.sm xs ys;
+printf("HHOS.sm\n");
+load -force HHP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHP_.sm xs ys;
+printf("HHP_.sm\n");
+load -force HHSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHSA.sm xs ys;
+printf("HHSA.sm\n");
+load -force HHSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHSD.sm xs ys;
+printf("HHSD.sm\n");
+load -force HHV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHV_.sm xs ys;
+printf("HHV_.sm\n");
+load -force HHZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHZn.sm xs ys;
+printf("HHZn.sm\n");
+load -force HHcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHcF.sm xs ys;
+printf("HHcF.sm\n");
+load -force HHcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii HHcP.sm xs ys;
+printf("HHcP.sm\n");
+load -force NCBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCBr.sm xs ys;
+printf("NCBr.sm\n");
+load -force NCC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCC0.sm xs ys;
+printf("NCC0.sm\n");
+load -force NCC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCC3.sm xs ys;
+printf("NCC3.sm\n");
+load -force NCCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCCF.sm xs ys;
+printf("NCCF.sm\n");
+load -force NCCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCCL.sm xs ys;
+printf("NCCL.sm\n");
+load -force NCCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCCN.sm xs ys;
+printf("NCCN.sm\n");
+load -force NCCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCCO.sm xs ys;
+printf("NCCO.sm\n");
+load -force NCCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCCP.sm xs ys;
+printf("NCCP.sm\n");
+load -force NCCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCCW.sm xs ys;
+printf("NCCW.sm\n");
+load -force NCF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCF_.sm xs ys;
+printf("NCF_.sm\n");
+load -force NCFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCFe.sm xs ys;
+printf("NCFe.sm\n");
+load -force NCHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCHL.sm xs ys;
+printf("NCHL.sm\n");
+load -force NCMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCMg.sm xs ys;
+printf("NCMg.sm\n");
+load -force NCMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCMn.sm xs ys;
+printf("NCMn.sm\n");
+load -force NCN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCN0.sm xs ys;
+printf("NCN0.sm\n");
+load -force NCNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCNA.sm xs ys;
+printf("NCNA.sm\n");
+load -force NCNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCNC.sm xs ys;
+printf("NCNC.sm\n");
+load -force NCND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCND.sm xs ys;
+printf("NCND.sm\n");
+load -force NCNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCNP.sm xs ys;
+printf("NCNP.sm\n");
+load -force NCNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCNR.sm xs ys;
+printf("NCNR.sm\n");
+load -force NCNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCNS.sm xs ys;
+printf("NCNS.sm\n");
+load -force NCOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCOA.sm xs ys;
+printf("NCOA.sm\n");
+load -force NCOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCOC.sm xs ys;
+printf("NCOC.sm\n");
+load -force NCOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCOD.sm xs ys;
+printf("NCOD.sm\n");
+load -force NCOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCOE.sm xs ys;
+printf("NCOE.sm\n");
+load -force NCOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCOR.sm xs ys;
+printf("NCOR.sm\n");
+load -force NCOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCOS.sm xs ys;
+printf("NCOS.sm\n");
+load -force NCP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCP_.sm xs ys;
+printf("NCP_.sm\n");
+load -force NCSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCSA.sm xs ys;
+printf("NCSA.sm\n");
+load -force NCSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCSD.sm xs ys;
+printf("NCSD.sm\n");
+load -force NCV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCV_.sm xs ys;
+printf("NCV_.sm\n");
+load -force NCZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCZn.sm xs ys;
+printf("NCZn.sm\n");
+load -force NCcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCcF.sm xs ys;
+printf("NCcF.sm\n");
+load -force NCcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NCcP.sm xs ys;
+printf("NCcP.sm\n");
+load -force NDBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDBr.sm xs ys;
+printf("NDBr.sm\n");
+load -force NDC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDC0.sm xs ys;
+printf("NDC0.sm\n");
+load -force NDC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDC3.sm xs ys;
+printf("NDC3.sm\n");
+load -force NDCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDCF.sm xs ys;
+printf("NDCF.sm\n");
+load -force NDCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDCL.sm xs ys;
+printf("NDCL.sm\n");
+load -force NDCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDCN.sm xs ys;
+printf("NDCN.sm\n");
+load -force NDCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDCO.sm xs ys;
+printf("NDCO.sm\n");
+load -force NDCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDCP.sm xs ys;
+printf("NDCP.sm\n");
+load -force NDCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDCW.sm xs ys;
+printf("NDCW.sm\n");
+load -force NDF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDF_.sm xs ys;
+printf("NDF_.sm\n");
+load -force NDFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDFe.sm xs ys;
+printf("NDFe.sm\n");
+load -force NDHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDHL.sm xs ys;
+printf("NDHL.sm\n");
+load -force NDMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDMg.sm xs ys;
+printf("NDMg.sm\n");
+load -force NDMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDMn.sm xs ys;
+printf("NDMn.sm\n");
+load -force NDN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDN0.sm xs ys;
+printf("NDN0.sm\n");
+load -force NDNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDNA.sm xs ys;
+printf("NDNA.sm\n");
+load -force NDNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDNC.sm xs ys;
+printf("NDNC.sm\n");
+load -force NDND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDND.sm xs ys;
+printf("NDND.sm\n");
+load -force NDNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDNP.sm xs ys;
+printf("NDNP.sm\n");
+load -force NDNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDNR.sm xs ys;
+printf("NDNR.sm\n");
+load -force NDNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDNS.sm xs ys;
+printf("NDNS.sm\n");
+load -force NDOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDOA.sm xs ys;
+printf("NDOA.sm\n");
+load -force NDOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDOC.sm xs ys;
+printf("NDOC.sm\n");
+load -force NDOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDOD.sm xs ys;
+printf("NDOD.sm\n");
+load -force NDOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDOE.sm xs ys;
+printf("NDOE.sm\n");
+load -force NDOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDOR.sm xs ys;
+printf("NDOR.sm\n");
+load -force NDOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDOS.sm xs ys;
+printf("NDOS.sm\n");
+load -force NDP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDP_.sm xs ys;
+printf("NDP_.sm\n");
+load -force NDSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDSA.sm xs ys;
+printf("NDSA.sm\n");
+load -force NDSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDSD.sm xs ys;
+printf("NDSD.sm\n");
+load -force NDV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDV_.sm xs ys;
+printf("NDV_.sm\n");
+load -force NDZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDZn.sm xs ys;
+printf("NDZn.sm\n");
+load -force NDcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDcF.sm xs ys;
+printf("NDcF.sm\n");
+load -force NDcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NDcP.sm xs ys;
+printf("NDcP.sm\n");
+load -force NRBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRBr.sm xs ys;
+printf("NRBr.sm\n");
+load -force NRC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRC0.sm xs ys;
+printf("NRC0.sm\n");
+load -force NRC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRC3.sm xs ys;
+printf("NRC3.sm\n");
+load -force NRCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRCF.sm xs ys;
+printf("NRCF.sm\n");
+load -force NRCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRCL.sm xs ys;
+printf("NRCL.sm\n");
+load -force NRCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRCN.sm xs ys;
+printf("NRCN.sm\n");
+load -force NRCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRCO.sm xs ys;
+printf("NRCO.sm\n");
+load -force NRCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRCP.sm xs ys;
+printf("NRCP.sm\n");
+load -force NRCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRCW.sm xs ys;
+printf("NRCW.sm\n");
+load -force NRF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRF_.sm xs ys;
+printf("NRF_.sm\n");
+load -force NRFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRFe.sm xs ys;
+printf("NRFe.sm\n");
+load -force NRHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRHL.sm xs ys;
+printf("NRHL.sm\n");
+load -force NRMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRMg.sm xs ys;
+printf("NRMg.sm\n");
+load -force NRMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRMn.sm xs ys;
+printf("NRMn.sm\n");
+load -force NRN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRN0.sm xs ys;
+printf("NRN0.sm\n");
+load -force NRNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRNA.sm xs ys;
+printf("NRNA.sm\n");
+load -force NRNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRNC.sm xs ys;
+printf("NRNC.sm\n");
+load -force NRND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRND.sm xs ys;
+printf("NRND.sm\n");
+load -force NRNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRNP.sm xs ys;
+printf("NRNP.sm\n");
+load -force NRNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRNR.sm xs ys;
+printf("NRNR.sm\n");
+load -force NRNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRNS.sm xs ys;
+printf("NRNS.sm\n");
+load -force NROA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NROA.sm xs ys;
+printf("NROA.sm\n");
+load -force NROC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NROC.sm xs ys;
+printf("NROC.sm\n");
+load -force NROD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NROD.sm xs ys;
+printf("NROD.sm\n");
+load -force NROE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NROE.sm xs ys;
+printf("NROE.sm\n");
+load -force NROR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NROR.sm xs ys;
+printf("NROR.sm\n");
+load -force NROS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NROS.sm xs ys;
+printf("NROS.sm\n");
+load -force NRP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRP_.sm xs ys;
+printf("NRP_.sm\n");
+load -force NRSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRSA.sm xs ys;
+printf("NRSA.sm\n");
+load -force NRSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRSD.sm xs ys;
+printf("NRSD.sm\n");
+load -force NRV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRV_.sm xs ys;
+printf("NRV_.sm\n");
+load -force NRZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRZn.sm xs ys;
+printf("NRZn.sm\n");
+load -force NRcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRcF.sm xs ys;
+printf("NRcF.sm\n");
+load -force NRcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii NRcP.sm xs ys;
+printf("NRcP.sm\n");
+load -force OABr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OABr.sm xs ys;
+printf("OABr.sm\n");
+load -force OAC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAC0.sm xs ys;
+printf("OAC0.sm\n");
+load -force OAC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAC3.sm xs ys;
+printf("OAC3.sm\n");
+load -force OACF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OACF.sm xs ys;
+printf("OACF.sm\n");
+load -force OACL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OACL.sm xs ys;
+printf("OACL.sm\n");
+load -force OACN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OACN.sm xs ys;
+printf("OACN.sm\n");
+load -force OACO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OACO.sm xs ys;
+printf("OACO.sm\n");
+load -force OACP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OACP.sm xs ys;
+printf("OACP.sm\n");
+load -force OACW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OACW.sm xs ys;
+printf("OACW.sm\n");
+load -force OAF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAF_.sm xs ys;
+printf("OAF_.sm\n");
+load -force OAFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAFe.sm xs ys;
+printf("OAFe.sm\n");
+load -force OAHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAHL.sm xs ys;
+printf("OAHL.sm\n");
+load -force OAMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAMg.sm xs ys;
+printf("OAMg.sm\n");
+load -force OAMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAMn.sm xs ys;
+printf("OAMn.sm\n");
+load -force OAN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAN0.sm xs ys;
+printf("OAN0.sm\n");
+load -force OANA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OANA.sm xs ys;
+printf("OANA.sm\n");
+load -force OANC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OANC.sm xs ys;
+printf("OANC.sm\n");
+load -force OAND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAND.sm xs ys;
+printf("OAND.sm\n");
+load -force OANP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OANP.sm xs ys;
+printf("OANP.sm\n");
+load -force OANR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OANR.sm xs ys;
+printf("OANR.sm\n");
+load -force OANS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OANS.sm xs ys;
+printf("OANS.sm\n");
+load -force OAOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAOA.sm xs ys;
+printf("OAOA.sm\n");
+load -force OAOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAOC.sm xs ys;
+printf("OAOC.sm\n");
+load -force OAOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAOD.sm xs ys;
+printf("OAOD.sm\n");
+load -force OAOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAOE.sm xs ys;
+printf("OAOE.sm\n");
+load -force OAOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAOR.sm xs ys;
+printf("OAOR.sm\n");
+load -force OAOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAOS.sm xs ys;
+printf("OAOS.sm\n");
+load -force OAP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAP_.sm xs ys;
+printf("OAP_.sm\n");
+load -force OASA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OASA.sm xs ys;
+printf("OASA.sm\n");
+load -force OASD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OASD.sm xs ys;
+printf("OASD.sm\n");
+load -force OAV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAV_.sm xs ys;
+printf("OAV_.sm\n");
+load -force OAZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAZn.sm xs ys;
+printf("OAZn.sm\n");
+load -force OAcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAcF.sm xs ys;
+printf("OAcF.sm\n");
+load -force OAcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OAcP.sm xs ys;
+printf("OAcP.sm\n");
+load -force OCBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCBr.sm xs ys;
+printf("OCBr.sm\n");
+load -force OCC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCC0.sm xs ys;
+printf("OCC0.sm\n");
+load -force OCC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCC3.sm xs ys;
+printf("OCC3.sm\n");
+load -force OCCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCCF.sm xs ys;
+printf("OCCF.sm\n");
+load -force OCCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCCL.sm xs ys;
+printf("OCCL.sm\n");
+load -force OCCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCCN.sm xs ys;
+printf("OCCN.sm\n");
+load -force OCCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCCO.sm xs ys;
+printf("OCCO.sm\n");
+load -force OCCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCCP.sm xs ys;
+printf("OCCP.sm\n");
+load -force OCCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCCW.sm xs ys;
+printf("OCCW.sm\n");
+load -force OCF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCF_.sm xs ys;
+printf("OCF_.sm\n");
+load -force OCFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCFe.sm xs ys;
+printf("OCFe.sm\n");
+load -force OCHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCHL.sm xs ys;
+printf("OCHL.sm\n");
+load -force OCMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCMg.sm xs ys;
+printf("OCMg.sm\n");
+load -force OCMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCMn.sm xs ys;
+printf("OCMn.sm\n");
+load -force OCN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCN0.sm xs ys;
+printf("OCN0.sm\n");
+load -force OCNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCNA.sm xs ys;
+printf("OCNA.sm\n");
+load -force OCNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCNC.sm xs ys;
+printf("OCNC.sm\n");
+load -force OCND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCND.sm xs ys;
+printf("OCND.sm\n");
+load -force OCNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCNP.sm xs ys;
+printf("OCNP.sm\n");
+load -force OCNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCNR.sm xs ys;
+printf("OCNR.sm\n");
+load -force OCNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCNS.sm xs ys;
+printf("OCNS.sm\n");
+load -force OCOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCOA.sm xs ys;
+printf("OCOA.sm\n");
+load -force OCOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCOC.sm xs ys;
+printf("OCOC.sm\n");
+load -force OCOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCOD.sm xs ys;
+printf("OCOD.sm\n");
+load -force OCOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCOE.sm xs ys;
+printf("OCOE.sm\n");
+load -force OCOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCOR.sm xs ys;
+printf("OCOR.sm\n");
+load -force OCOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCOS.sm xs ys;
+printf("OCOS.sm\n");
+load -force OCP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCP_.sm xs ys;
+printf("OCP_.sm\n");
+load -force OCSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCSA.sm xs ys;
+printf("OCSA.sm\n");
+load -force OCSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCSD.sm xs ys;
+printf("OCSD.sm\n");
+load -force OCV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCV_.sm xs ys;
+printf("OCV_.sm\n");
+load -force OCZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCZn.sm xs ys;
+printf("OCZn.sm\n");
+load -force OCcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCcF.sm xs ys;
+printf("OCcF.sm\n");
+load -force OCcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OCcP.sm xs ys;
+printf("OCcP.sm\n");
+load -force ODBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODBr.sm xs ys;
+printf("ODBr.sm\n");
+load -force ODC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODC0.sm xs ys;
+printf("ODC0.sm\n");
+load -force ODC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODC3.sm xs ys;
+printf("ODC3.sm\n");
+load -force ODCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODCF.sm xs ys;
+printf("ODCF.sm\n");
+load -force ODCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODCL.sm xs ys;
+printf("ODCL.sm\n");
+load -force ODCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODCN.sm xs ys;
+printf("ODCN.sm\n");
+load -force ODCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODCO.sm xs ys;
+printf("ODCO.sm\n");
+load -force ODCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODCP.sm xs ys;
+printf("ODCP.sm\n");
+load -force ODCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODCW.sm xs ys;
+printf("ODCW.sm\n");
+load -force ODF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODF_.sm xs ys;
+printf("ODF_.sm\n");
+load -force ODFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODFe.sm xs ys;
+printf("ODFe.sm\n");
+load -force ODHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODHL.sm xs ys;
+printf("ODHL.sm\n");
+load -force ODMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODMg.sm xs ys;
+printf("ODMg.sm\n");
+load -force ODMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODMn.sm xs ys;
+printf("ODMn.sm\n");
+load -force ODN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODN0.sm xs ys;
+printf("ODN0.sm\n");
+load -force ODNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODNA.sm xs ys;
+printf("ODNA.sm\n");
+load -force ODNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODNC.sm xs ys;
+printf("ODNC.sm\n");
+load -force ODND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODND.sm xs ys;
+printf("ODND.sm\n");
+load -force ODNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODNP.sm xs ys;
+printf("ODNP.sm\n");
+load -force ODNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODNR.sm xs ys;
+printf("ODNR.sm\n");
+load -force ODNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODNS.sm xs ys;
+printf("ODNS.sm\n");
+load -force ODOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODOA.sm xs ys;
+printf("ODOA.sm\n");
+load -force ODOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODOC.sm xs ys;
+printf("ODOC.sm\n");
+load -force ODOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODOD.sm xs ys;
+printf("ODOD.sm\n");
+load -force ODOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODOE.sm xs ys;
+printf("ODOE.sm\n");
+load -force ODOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODOR.sm xs ys;
+printf("ODOR.sm\n");
+load -force ODOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODOS.sm xs ys;
+printf("ODOS.sm\n");
+load -force ODP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODP_.sm xs ys;
+printf("ODP_.sm\n");
+load -force ODSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODSA.sm xs ys;
+printf("ODSA.sm\n");
+load -force ODSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODSD.sm xs ys;
+printf("ODSD.sm\n");
+load -force ODV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODV_.sm xs ys;
+printf("ODV_.sm\n");
+load -force ODZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODZn.sm xs ys;
+printf("ODZn.sm\n");
+load -force ODcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODcF.sm xs ys;
+printf("ODcF.sm\n");
+load -force ODcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii ODcP.sm xs ys;
+printf("ODcP.sm\n");
+load -force OWBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWBr.sm xs ys;
+printf("OWBr.sm\n");
+load -force OWC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWC0.sm xs ys;
+printf("OWC0.sm\n");
+load -force OWC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWC3.sm xs ys;
+printf("OWC3.sm\n");
+load -force OWCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWCF.sm xs ys;
+printf("OWCF.sm\n");
+load -force OWCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWCL.sm xs ys;
+printf("OWCL.sm\n");
+load -force OWCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWCN.sm xs ys;
+printf("OWCN.sm\n");
+load -force OWCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWCO.sm xs ys;
+printf("OWCO.sm\n");
+load -force OWCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWCP.sm xs ys;
+printf("OWCP.sm\n");
+load -force OWCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWCW.sm xs ys;
+printf("OWCW.sm\n");
+load -force OWF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWF_.sm xs ys;
+printf("OWF_.sm\n");
+load -force OWFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWFe.sm xs ys;
+printf("OWFe.sm\n");
+load -force OWHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWHL.sm xs ys;
+printf("OWHL.sm\n");
+load -force OWMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWMg.sm xs ys;
+printf("OWMg.sm\n");
+load -force OWMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWMn.sm xs ys;
+printf("OWMn.sm\n");
+load -force OWN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWN0.sm xs ys;
+printf("OWN0.sm\n");
+load -force OWNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWNA.sm xs ys;
+printf("OWNA.sm\n");
+load -force OWNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWNC.sm xs ys;
+printf("OWNC.sm\n");
+load -force OWND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWND.sm xs ys;
+printf("OWND.sm\n");
+load -force OWNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWNP.sm xs ys;
+printf("OWNP.sm\n");
+load -force OWNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWNR.sm xs ys;
+printf("OWNR.sm\n");
+load -force OWNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWNS.sm xs ys;
+printf("OWNS.sm\n");
+load -force OWOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWOA.sm xs ys;
+printf("OWOA.sm\n");
+load -force OWOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWOC.sm xs ys;
+printf("OWOC.sm\n");
+load -force OWOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWOD.sm xs ys;
+printf("OWOD.sm\n");
+load -force OWOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWOE.sm xs ys;
+printf("OWOE.sm\n");
+load -force OWOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWOR.sm xs ys;
+printf("OWOR.sm\n");
+load -force OWOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWOS.sm xs ys;
+printf("OWOS.sm\n");
+load -force OWP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWP_.sm xs ys;
+printf("OWP_.sm\n");
+load -force OWSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWSA.sm xs ys;
+printf("OWSA.sm\n");
+load -force OWSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWSD.sm xs ys;
+printf("OWSD.sm\n");
+load -force OWV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWV_.sm xs ys;
+printf("OWV_.sm\n");
+load -force OWZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWZn.sm xs ys;
+printf("OWZn.sm\n");
+load -force OWcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWcF.sm xs ys;
+printf("OWcF.sm\n");
+load -force OWcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii OWcP.sm xs ys;
+printf("OWcP.sm\n");
+load -force SABr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SABr.sm xs ys;
+printf("SABr.sm\n");
+load -force SAC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAC0.sm xs ys;
+printf("SAC0.sm\n");
+load -force SAC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAC3.sm xs ys;
+printf("SAC3.sm\n");
+load -force SACF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SACF.sm xs ys;
+printf("SACF.sm\n");
+load -force SACL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SACL.sm xs ys;
+printf("SACL.sm\n");
+load -force SACN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SACN.sm xs ys;
+printf("SACN.sm\n");
+load -force SACO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SACO.sm xs ys;
+printf("SACO.sm\n");
+load -force SACP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SACP.sm xs ys;
+printf("SACP.sm\n");
+load -force SACW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SACW.sm xs ys;
+printf("SACW.sm\n");
+load -force SAF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAF_.sm xs ys;
+printf("SAF_.sm\n");
+load -force SAFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAFe.sm xs ys;
+printf("SAFe.sm\n");
+load -force SAHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAHL.sm xs ys;
+printf("SAHL.sm\n");
+load -force SAMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAMg.sm xs ys;
+printf("SAMg.sm\n");
+load -force SAMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAMn.sm xs ys;
+printf("SAMn.sm\n");
+load -force SAN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAN0.sm xs ys;
+printf("SAN0.sm\n");
+load -force SANA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SANA.sm xs ys;
+printf("SANA.sm\n");
+load -force SANC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SANC.sm xs ys;
+printf("SANC.sm\n");
+load -force SAND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAND.sm xs ys;
+printf("SAND.sm\n");
+load -force SANP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SANP.sm xs ys;
+printf("SANP.sm\n");
+load -force SANR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SANR.sm xs ys;
+printf("SANR.sm\n");
+load -force SANS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SANS.sm xs ys;
+printf("SANS.sm\n");
+load -force SAOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAOA.sm xs ys;
+printf("SAOA.sm\n");
+load -force SAOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAOC.sm xs ys;
+printf("SAOC.sm\n");
+load -force SAOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAOD.sm xs ys;
+printf("SAOD.sm\n");
+load -force SAOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAOE.sm xs ys;
+printf("SAOE.sm\n");
+load -force SAOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAOR.sm xs ys;
+printf("SAOR.sm\n");
+load -force SAOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAOS.sm xs ys;
+printf("SAOS.sm\n");
+load -force SAP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAP_.sm xs ys;
+printf("SAP_.sm\n");
+load -force SASA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SASA.sm xs ys;
+printf("SASA.sm\n");
+load -force SASD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SASD.sm xs ys;
+printf("SASD.sm\n");
+load -force SAV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAV_.sm xs ys;
+printf("SAV_.sm\n");
+load -force SAZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAZn.sm xs ys;
+printf("SAZn.sm\n");
+load -force SAcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAcF.sm xs ys;
+printf("SAcF.sm\n");
+load -force SAcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SAcP.sm xs ys;
+printf("SAcP.sm\n");
+load -force SDBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDBr.sm xs ys;
+printf("SDBr.sm\n");
+load -force SDC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDC0.sm xs ys;
+printf("SDC0.sm\n");
+load -force SDC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDC3.sm xs ys;
+printf("SDC3.sm\n");
+load -force SDCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDCF.sm xs ys;
+printf("SDCF.sm\n");
+load -force SDCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDCL.sm xs ys;
+printf("SDCL.sm\n");
+load -force SDCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDCN.sm xs ys;
+printf("SDCN.sm\n");
+load -force SDCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDCO.sm xs ys;
+printf("SDCO.sm\n");
+load -force SDCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDCP.sm xs ys;
+printf("SDCP.sm\n");
+load -force SDCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDCW.sm xs ys;
+printf("SDCW.sm\n");
+load -force SDF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDF_.sm xs ys;
+printf("SDF_.sm\n");
+load -force SDFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDFe.sm xs ys;
+printf("SDFe.sm\n");
+load -force SDHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDHL.sm xs ys;
+printf("SDHL.sm\n");
+load -force SDMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDMg.sm xs ys;
+printf("SDMg.sm\n");
+load -force SDMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDMn.sm xs ys;
+printf("SDMn.sm\n");
+load -force SDN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDN0.sm xs ys;
+printf("SDN0.sm\n");
+load -force SDNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDNA.sm xs ys;
+printf("SDNA.sm\n");
+load -force SDNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDNC.sm xs ys;
+printf("SDNC.sm\n");
+load -force SDND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDND.sm xs ys;
+printf("SDND.sm\n");
+load -force SDNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDNP.sm xs ys;
+printf("SDNP.sm\n");
+load -force SDNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDNR.sm xs ys;
+printf("SDNR.sm\n");
+load -force SDNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDNS.sm xs ys;
+printf("SDNS.sm\n");
+load -force SDOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDOA.sm xs ys;
+printf("SDOA.sm\n");
+load -force SDOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDOC.sm xs ys;
+printf("SDOC.sm\n");
+load -force SDOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDOD.sm xs ys;
+printf("SDOD.sm\n");
+load -force SDOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDOE.sm xs ys;
+printf("SDOE.sm\n");
+load -force SDOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDOR.sm xs ys;
+printf("SDOR.sm\n");
+load -force SDOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDOS.sm xs ys;
+printf("SDOS.sm\n");
+load -force SDP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDP_.sm xs ys;
+printf("SDP_.sm\n");
+load -force SDSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDSA.sm xs ys;
+printf("SDSA.sm\n");
+load -force SDSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDSD.sm xs ys;
+printf("SDSD.sm\n");
+load -force SDV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDV_.sm xs ys;
+printf("SDV_.sm\n");
+load -force SDZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDZn.sm xs ys;
+printf("SDZn.sm\n");
+load -force SDcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDcF.sm xs ys;
+printf("SDcF.sm\n");
+load -force SDcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii SDcP.sm xs ys;
+printf("SDcP.sm\n");
+load -force cFBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFBr.sm xs ys;
+printf("cFBr.sm\n");
+load -force cFC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFC0.sm xs ys;
+printf("cFC0.sm\n");
+load -force cFC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFC3.sm xs ys;
+printf("cFC3.sm\n");
+load -force cFCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFCF.sm xs ys;
+printf("cFCF.sm\n");
+load -force cFCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFCL.sm xs ys;
+printf("cFCL.sm\n");
+load -force cFCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFCN.sm xs ys;
+printf("cFCN.sm\n");
+load -force cFCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFCO.sm xs ys;
+printf("cFCO.sm\n");
+load -force cFCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFCP.sm xs ys;
+printf("cFCP.sm\n");
+load -force cFCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFCW.sm xs ys;
+printf("cFCW.sm\n");
+load -force cFF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFF_.sm xs ys;
+printf("cFF_.sm\n");
+load -force cFFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFFe.sm xs ys;
+printf("cFFe.sm\n");
+load -force cFHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFHL.sm xs ys;
+printf("cFHL.sm\n");
+load -force cFMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFMg.sm xs ys;
+printf("cFMg.sm\n");
+load -force cFMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFMn.sm xs ys;
+printf("cFMn.sm\n");
+load -force cFN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFN0.sm xs ys;
+printf("cFN0.sm\n");
+load -force cFNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFNA.sm xs ys;
+printf("cFNA.sm\n");
+load -force cFNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFNC.sm xs ys;
+printf("cFNC.sm\n");
+load -force cFND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFND.sm xs ys;
+printf("cFND.sm\n");
+load -force cFNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFNP.sm xs ys;
+printf("cFNP.sm\n");
+load -force cFNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFNR.sm xs ys;
+printf("cFNR.sm\n");
+load -force cFNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFNS.sm xs ys;
+printf("cFNS.sm\n");
+load -force cFOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFOA.sm xs ys;
+printf("cFOA.sm\n");
+load -force cFOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFOC.sm xs ys;
+printf("cFOC.sm\n");
+load -force cFOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFOD.sm xs ys;
+printf("cFOD.sm\n");
+load -force cFOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFOE.sm xs ys;
+printf("cFOE.sm\n");
+load -force cFOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFOR.sm xs ys;
+printf("cFOR.sm\n");
+load -force cFOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFOS.sm xs ys;
+printf("cFOS.sm\n");
+load -force cFP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFP_.sm xs ys;
+printf("cFP_.sm\n");
+load -force cFSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFSA.sm xs ys;
+printf("cFSA.sm\n");
+load -force cFSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFSD.sm xs ys;
+printf("cFSD.sm\n");
+load -force cFV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFV_.sm xs ys;
+printf("cFV_.sm\n");
+load -force cFZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFZn.sm xs ys;
+printf("cFZn.sm\n");
+load -force cFcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFcF.sm xs ys;
+printf("cFcF.sm\n");
+load -force cFcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cFcP.sm xs ys;
+printf("cFcP.sm\n");
+load -force cPBr.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPBr.sm xs ys;
+printf("cPBr.sm\n");
+load -force cPC0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPC0.sm xs ys;
+printf("cPC0.sm\n");
+load -force cPC3.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPC3.sm xs ys;
+printf("cPC3.sm\n");
+load -force cPCF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPCF.sm xs ys;
+printf("cPCF.sm\n");
+load -force cPCL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPCL.sm xs ys;
+printf("cPCL.sm\n");
+load -force cPCN.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPCN.sm xs ys;
+printf("cPCN.sm\n");
+load -force cPCO.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPCO.sm xs ys;
+printf("cPCO.sm\n");
+load -force cPCP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPCP.sm xs ys;
+printf("cPCP.sm\n");
+load -force cPCW.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPCW.sm xs ys;
+printf("cPCW.sm\n");
+load -force cPF_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPF_.sm xs ys;
+printf("cPF_.sm\n");
+load -force cPFe.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPFe.sm xs ys;
+printf("cPFe.sm\n");
+load -force cPHL.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPHL.sm xs ys;
+printf("cPHL.sm\n");
+load -force cPMg.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPMg.sm xs ys;
+printf("cPMg.sm\n");
+load -force cPMn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPMn.sm xs ys;
+printf("cPMn.sm\n");
+load -force cPN0.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPN0.sm xs ys;
+printf("cPN0.sm\n");
+load -force cPNA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPNA.sm xs ys;
+printf("cPNA.sm\n");
+load -force cPNC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPNC.sm xs ys;
+printf("cPNC.sm\n");
+load -force cPND.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPND.sm xs ys;
+printf("cPND.sm\n");
+load -force cPNP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPNP.sm xs ys;
+printf("cPNP.sm\n");
+load -force cPNR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPNR.sm xs ys;
+printf("cPNR.sm\n");
+load -force cPNS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPNS.sm xs ys;
+printf("cPNS.sm\n");
+load -force cPOA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPOA.sm xs ys;
+printf("cPOA.sm\n");
+load -force cPOC.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPOC.sm xs ys;
+printf("cPOC.sm\n");
+load -force cPOD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPOD.sm xs ys;
+printf("cPOD.sm\n");
+load -force cPOE.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPOE.sm xs ys;
+printf("cPOE.sm\n");
+load -force cPOR.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPOR.sm xs ys;
+printf("cPOR.sm\n");
+load -force cPOS.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPOS.sm xs ys;
+printf("cPOS.sm\n");
+load -force cPP_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPP_.sm xs ys;
+printf("cPP_.sm\n");
+load -force cPSA.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPSA.sm xs ys;
+printf("cPSA.sm\n");
+load -force cPSD.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPSD.sm xs ys;
+printf("cPSD.sm\n");
+load -force cPV_.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPV_.sm xs ys;
+printf("cPV_.sm\n");
+load -force cPZn.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPZn.sm xs ys;
+printf("cPZn.sm\n");
+load -force cPcF.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPcF.sm xs ys;
+printf("cPcF.sm\n");
+load -force cPcP.dat;
+[xs,ys]=savgol(x,y,5,5,4);
+save -ascii cPcP.sm xs ys;
+printf("cPcP.sm\n");