Fixing bugs

code/RDownScript.m

@@ -16,7 +16,7 @@
 (*Load package and set root paths*)
 
 
-rootpath="~/git/rdstackingproject/";
+rootpath="/work/francisco.jimenez/sio/git/rdstackingproject/";
 
 
 SetDirectory[rootpath<>"code"];
@@ -24,7 +24,7 @@ SetDirectory[rootpath<>"code"];
 <<mcmc.m
 
 
-mysxscase="/Users/xisco/SXS/BBH_SKS_d14.3_q1.22_sA_0_0_0.330_sB_0_0_-0.440";
+mysxscase=rootpath<>"SXS/BBH_SKS_d14.3_q1.22_sA_0_0_0.330_sB_0_0_-0.440";
 tshift=0;
 npoints=10;
 ntones=1;
@@ -160,7 +160,7 @@ tab=Table[
 ansatz=OvertoneModelV2[ntones,{\[Eta],\[Chi]1,\[Chi]2},t0,"Fit\[Alpha]"->{},"Vary\[Omega]"->True,"\[Omega]val"->{-\[Omega]fact,\[Omega]fact},"Export_\[Omega]val"->True,"ModesData"->Modedata];
 cfit=NonlinearModelFit[data,ansatz[[1]],{x0,x1,x2,x3,x4,x5,x6,x7},t];
 cfitd=Transpose[{data[[All,1]],Normal[cfit]/.t->data[[All,1]]}];
-{ansatz[[2]],cfit["BestFitParameters"],1-EasyMatchT[data,cfitd,t0,t0+90]},{j,npoints}];
+{ansatz[[2]],{x0,x1,x2,x3,x4,x5,x6,x7}/.cfit["BestFitParameters"],1-EasyMatchT[data,cfitd,t0,t0+90]},{j,npoints}];
 
 
-Export[rootpath<>"plots/freq_bfitpars_mmatch.dat",tab]
+Export[rootpath<>"results_data/freq_bfitpars_mmatch.dat",tab]