SetDirectory["add the location of mathematica module directory of pencil code"];
Needs["pc`"]
dir="add the run directory address";
outdir="add the output dir address";

Figure 5

runs={"P1024_k01_kf10e","P1024_k1_kf30d","P1024_k1_kf10c","V1024_k01_kf10i","V1024_k1_kf30d","V1024_k1_kf10c"};
Do[
{t[sim,1],specgw[sim]}=read1D[sim,"power_GWs.dat"];
{t[sim,1],speckin[sim]}=read1D[sim,"power_kin.dat"];
{t[sim,1],specstr[sim]}=read1D[sim,"power_Str.dat"],{sim,runs}]
Do[n1=readParamNml[runs[[j]],"start.in","WAV1"];n2=readParamNml[runs[[j+3]],"start.in","WAV1"];plotgw[j]=ListLogLogPlot[{Table[{(i-1)*n1,Mean[specgw[runs[[j]]][[1;;Dimensions[t[runs[[j]],1]][[1]],i]]]},{i,1,512}],Table[{(i-1)*n2,Mean[specgw[runs[[j+3]]][[1;;Dimensions[t[runs[[j+3]],1]][[1]],i]]]},{i,1,512}]},Joined->True,PlotStyle-> {Black,{Dashed,Black},Red,{Dashed,Red},Blue,{Dashed,Blue}},Frame->True,FrameStyle->Directive[Black,Bold,10],FrameLabel->{Style["k",Black,Bold,12],Style["Sp((OverDot[h]))",Black,Bold,12]},PlotLegends->Placed[LineLegend[{StringTake[Style[runs[[j]],10]//ToString,-15],StringTake[Style[runs[[j+3]],10]//ToString,-15]},LegendFunction->"Frame",LegendMarkerSize->15],{0.32,0.28}]];plotkin[j]=ListLogLogPlot[{Table[{(i-1)*n1,speckin[runs[[j]]][[1,i]]},{i,1,512}],Table[{(i-1)*n2,speckin[runs[[j+3]]][[1,i]]},{i,1,512}],Table[{(i-1)*n1,speckin[runs[[j]]][[Dimensions[t[runs[[j]],1]][[1]],i]]},{i,1,512}],Table[{(i-1)*n2,speckin[runs[[j+3]]][[Dimensions[t[runs[[j+3]],1]][[1]],i]]},{i,1,512}]},Joined->True,PlotStyle-> {Black,{Dashed,Black},Red,{Dashed,Red},Blue,{Dashed,Blue}},Frame->True,FrameStyle->Directive[Black,Bold,10],FrameLabel->{Style["k",Black,Bold,12],Style["(Subscript[E, k])",Black,Bold,12]},PlotLegends->Placed[LineLegend[{StringTake[Style[runs[[j]],10]//ToString,-15],StringTake[Style[runs[[j+3]],10]//ToString,-15]},LegendFunction->"Frame",LegendMarkerSize->15],{0.32,0.28}]];plotstr[j]=ListLogLogPlot[{Table[{(i-1)*n1,specstr[runs[[j]]][[1,i]]},{i,1,512}],Table[{(i-1)*n2,specstr[runs[[j+3]]][[1,i]]},{i,1,512}],Table[{(i-1)*n1,specstr[runs[[j]]][[Dimensions[t[runs[[j]],1]][[1]],i]]},{i,1,512}],Table[{(i-1)*n2,specstr[runs[[j+3]]][[Dimensions[t[runs[[j+3]],1]][[1]],i]]},{i,1,512}]},Joined->True,PlotStyle-> {Black,{Dashed,Black},Red,{Dashed,Red},Blue,{Dashed,Blue}},Frame->True,FrameStyle->Directive[Black,Bold,10],FrameLabel->{Style["k",Black,Bold,12],Style["Sp(T)",Black,Bold,12]},PlotLegends->Placed[LineLegend[{StringTake[Style[runs[[j]],10]//ToString,-15],StringTake[Style[runs[[j+3]],10]//ToString,-15]},LegendFunction->"Frame",LegendMarkerSize->15],{0.32,0.28}]],{j,1,3}]

Clear[omega0];
tt={80,80,802,82,80,802};(*These values are chosen such that the spectrum is at the same time for all runs*)
kpeak={10,30,10,10,30,10};
Do[omega0[j]=((((Total[speckin[runs[[j]]][[1,All]]]*2))^2)/((kpeak[[j]])^2))*(54/(π^2))((((4/3)))^2),{j,1,6}];

gwnorm=1;
tti={1,1,11,3,1,11};
plotRangesgw={{{0.1,10^2},{10^-9*(0.4/(6*omega0[1])),10^-5*(((*1*)0.4)/(6*omega0[1]))}},{{1,10^3},{10^-12*(1/(6*omega0[5])),10^-5*(1/(6*omega0[5]))}},{{1,10^3},{10^-12*(1/(6*omega0[3])),10^-5*(1/(6*omega0[3]))}}};
plotRangeskin={{{0.1,10^2},{10^-21,10^-2}},{{1,10^3},{10^-21,10^-2}},{{1,10^3},{10^-21,10^-2}}};
plotRangesstr={{{0.1,10^2},{10^-23,10^-2}},{{1,10^3},{10^-23,10^-2}},{{1,10^3},{10^-23,10^-2}}};
Do[n1=readParamNml[runs[[j]],"start.in","WAV1"];n2=readParamNml[runs[[j+3]],"start.in","WAV1"];plotgw[j]=ListLogLogPlot[{Table[{(i-1)*n1,(1/(6*omega0[j]))*gwnorm*specgw[runs[[j]]][[tt[[j]],i]]},{i,1,512}],Table[{(i-1)*n2,(1/(6*omega0[j+3]))*gwnorm*specgw[runs[[j+3]]][[tt[[j+3]],i]]},{i,1,512}]},Joined->True,PlotStyle-> {{Black,Thickness[0.005]},{Dotted,Black,Thickness[0.005]},{Red,Thickness[0.005]},{Dotted,Red,Thickness[0.005]}},PlotRange->plotRangesgw[[j]],Frame->True,FrameStyle->Directive[Black,4],FrameLabel->{Style["k",Black,6,FontSlant->Italic],If[j==1,Style["(Subscript[Ω, GW])/k (Subscript[Ω, 0])",Black,6]]},PlotLegends->Placed[LineLegend[{Style["A"<>ToString[j],4],Style["V"<>ToString[j],4]},Spacings->0,LegendMarkerSize->{10,0.1}],{0.52,0.28}],Epilog->{If[j==1,Text[Style["(g)",5],{4,-5}]],If[j==2,Text[Style["(h)",5],{6,-4}]],If[j==3,Text[Style["(i)",5],{6.3,-5}]]}];plotkin[j]=ListLogLogPlot[{Table[{(i-1)*n1,speckin[runs[[j]]][[tt[[j]],i]]},{i,1,512}],Table[{(i-1)*n2,speckin[runs[[j+3]]][[tt[[j+3]],i]]},{i,1,512}],Table[{(i-1)*n1,speckin[runs[[j]]][[tti[[j]],i]]},{i,1,512}],Table[{(i-1)*n2,speckin[runs[[j+3]]][[tti[[j+3]],i]]},{i,1,512}]},Joined->True,PlotRange->plotRangeskin[[j]],PlotStyle->{{Black,Thickness[0.005]},{Dotted,Black,Thickness[0.005]},{Red,Thickness[0.005]},{Dotted,Red,Thickness[0.005]}},Frame->True,FrameStyle->Directive[Black,4],FrameLabel->{If[j==4,Style["k",Black,6,FontSlant->Italic]],If[j==1,Style["(Subscript[E, K])",Black,6]]},PlotLegends->Placed[LineLegend[{Style["A"<>ToString[j],4],Style["V"<>ToString[j],4]},Spacings->0,LegendMarkerSize->{10,0.1}],{0.32,0.28}],Epilog->{If[j==1,Text[Style["(a)",5],{4,-23.5}]],If[j==2,Text[Style["(b)",5],{6,-9.5}]],If[j==3,Text[Style["(c)",5],{6.3,-9.5}]]}];plotstr[j]=ListLogLogPlot[{Table[{(i-1)*n1,specstr[runs[[j]]][[tt[[j]],i]]},{i,1,512}],Table[{(i-1)*n2,specstr[runs[[j+3]]][[tt[[j+3]],i]]},{i,1,512}],Table[{(i-1)*n1,specstr[runs[[j]]][[tti[[j]],i]]},{i,1,512}],Table[{(i-1)*n2,specstr[runs[[j+3]]][[tti[[j+3]],i]]},{i,1,512}]},Joined->True,PlotRange->plotRangesstr[[j]],PlotStyle-> {{Black,Thickness[0.005]},{Dotted,Black,Thickness[0.005]},{Red,Thickness[0.005]},{Dotted,Red,Thickness[0.005]}},Frame->True,FrameStyle->Directive[Black,4],FrameLabel->{If[j==4,Style["k",Black,6]],If[j==1,Style["Sp(T)",Black,6]]},PlotLegends->Placed[LineLegend[{Style["A"<>ToString[j],4],Style["V"<>ToString[j],4]},Spacings->0,LegendMarkerSize->{10,0.1}],{0.32,0.28}],Epilog->{If[j==1,Text[Style["(d)",5],{4,-32}]],If[j==2,Text[Style["(e)",5],{6,-12}]],If[j==3,Text[Style["(f)",5],{6.3,-12}]]}],{j,1,3}]

plott=GraphicsGrid[{{plotkin[1],plotkin[2],plotkin[3]},{plotstr[1],plotstr[2],plotstr[3]},{plotgw[1],plotgw[2],plotgw[3]}},Spacings->-50,AspectRatio->1/GoldenRatio]

Export[outdir<>"acoustic_vs_vortical.eps",plott]

Figure 10

runs={"P1024_k01_kf10e_rho_nonuni","P1024_k1_kf10c_rho_nonuni","P1024_k01_kf10e","P1024_k1_kf10c"};
Do[
{t[sim,1],specgw[sim]}=read1D[sim,"power_GWs.dat"];
{t[sim,1],speckin[sim]}=read1D[sim,"power_kin.dat"];
{t[sim,1],specstr[sim]}=read1D[sim,"power_Str.dat"],{sim,runs}]
Do[n1=readParamNml[runs[[j]],"start.in","WAV1"];n2=readParamNml[runs[[j+2]],"start.in","WAV1"];plotgw[j]=ListLogLogPlot[{Table[{(i-1)*n1,Mean[specgw[runs[[j]]][[1;;Dimensions[t[runs[[j]],1]][[1]],i]]]},{i,1,512}],Table[{(i-1)*n2,Mean[specgw[runs[[j+2]]][[1;;Dimensions[t[runs[[j+2]],1]][[1]],i]]]},{i,1,512}]},Joined->True,PlotStyle-> {Black,{Dashed,Black},Red,{Dashed,Red},Blue,{Dashed,Blue}},Frame->True,FrameStyle->Directive[Black,Bold,10],FrameLabel->{Style["k",Black,Bold,12],Style["Sp((OverDot[h]))",Black,Bold,12]},PlotLegends->Placed[LineLegend[{StringTake[Style[runs[[j]],10]//ToString,-15],StringTake[Style[runs[[j+2]],10]//ToString,-15]},LegendFunction->"Frame",LegendMarkerSize->15],{0.32,0.28}]];plotkin[j]=ListLogLogPlot[{Table[{(i-1)*n1,speckin[runs[[j]]][[1,i]]},{i,1,512}],Table[{(i-1)*n2,speckin[runs[[j+2]]][[1,i]]},{i,1,512}],Table[{(i-1)*n1,speckin[runs[[j]]][[Dimensions[t[runs[[j]],1]][[1]],i]]},{i,1,512}],Table[{(i-1)*n2,speckin[runs[[j+2]]][[Dimensions[t[runs[[j+2]],1]][[1]],i]]},{i,1,512}]},Joined->True,PlotStyle-> {Black,{Dashed,Black},Red,{Dashed,Red},Blue,{Dashed,Blue}},Frame->True,FrameStyle->Directive[Black,Bold,10],FrameLabel->{Style["k",Black,Bold,12],Style["(Subscript[E, k])",Black,Bold,12]},PlotLegends->Placed[LineLegend[{StringTake[Style[runs[[j]],10]//ToString,-15],StringTake[Style[runs[[j+2]],10]//ToString,-15]},LegendFunction->"Frame",LegendMarkerSize->15],{0.32,0.28}]];plotstr[j]=ListLogLogPlot[{Table[{(i-1)*n1,specstr[runs[[j]]][[1,i]]},{i,1,512}],Table[{(i-1)*n2,specstr[runs[[j+2]]][[1,i]]},{i,1,512}],Table[{(i-1)*n1,specstr[runs[[j]]][[Dimensions[t[runs[[j]],1]][[1]],i]]},{i,1,512}],Table[{(i-1)*n2,specstr[runs[[j+2]]][[Dimensions[t[runs[[j+2]],1]][[1]],i]]},{i,1,512}]},Joined->True,PlotStyle-> {Black,{Dashed,Black},Red,{Dashed,Red},Blue,{Dashed,Blue}},Frame->True,FrameStyle->Directive[Black,Bold,10],FrameLabel->{Style["k",Black,Bold,12],Style["Sp(T)",Black,Bold,12]},PlotLegends->Placed[LineLegend[{StringTake[Style[runs[[j]],10]//ToString,-15],StringTake[Style[runs[[j+2]],10]//ToString,-15]},LegendFunction->"Frame",LegendMarkerSize->15],{0.32,0.28}]],{j,1,2}]

Clear[omega0]
tt={80,605,80,605};
kpeak={10,10,10,10};
Do[omega0[j]=((((Total[speckin[runs[[j]]][[1,All]]]*2))^2)/((kpeak[[j]])^2))*(54/(π^2))((((4/3)))^2),{j,1,4}];

gwnorm=1;
tti={1,11,1,11};
plotRangesgw={{{0.1,10^2},{10^-9*(0.4/(6*omega0[1])),10^-5*(((*1*)0.4)/(6*omega0[1]))}},{{1,10^3},{10^-12*(1/(6*omega0[2])),10^-5*(1/(6*omega0[2]))}}};
plotRangeskin={{{0.1,10^2},{10^-21,10^-2}},{{1,10^3},{10^-21,10^-2}}};
plotRangesstr={{{0.1,10^2},{10^-23,10^-2}},{{1,10^3},{10^-23,10^-2}}};
Do[n1=readParamNml[runs[[j]],"start.in","WAV1"];n2=readParamNml[runs[[j+2]],"start.in","WAV1"];plotgw[j]=ListLogLogPlot[{Table[{(i-1)*n1,(1/(6*omega0[j]))*gwnorm*specgw[runs[[j]]][[tt[[j]],i]]},{i,1,512}],Table[{(i-1)*n2,(1/(6*omega0[j+2]))*gwnorm*specgw[runs[[j+2]]][[tt[[j+2]],i]]},{i,1,512}]},Joined->True,PlotStyle-> {{Green,Thickness[0.005]},{Black,Thickness[0.005]},{Red,Thickness[0.005]},{Blue,Thickness[0.005]},{Dotted,Red,Thickness[0.005]}},PlotRange->plotRangesgw[[j]],Frame->True,FrameStyle->Directive[Black,4],FrameLabel->{Style["k",Black,6,FontSlant->Italic],If[j==1,Style["(Subscript[Ω, GW])/k (Subscript[Ω, 0])",Black,6]]},PlotLegends->Placed[LineLegend[{Style["A"<>ToString[If[j==1,j,j+1]]<>"'",4],Style["A"<>ToString[If[j==1,j,j+1]],4]},Spacings->0,LegendMarkerSize->{10,0.1}],{0.52,0.28}],Epilog->{If[j==1,Text[Style["(c)",5],{4,-5}]],If[j==2,Text[Style["(d)",5],{6.3,-5}]]}];plotkin[j]=ListLogLogPlot[{Table[{(i-1)*n1,speckin[runs[[j]]][[tt[[j]],i]]},{i,1,512}],Table[{(i-1)*n2,speckin[runs[[j+2]]][[tt[[j+2]],i]]},{i,1,512}],Table[{(i-1)*n1,speckin[runs[[j]]][[tti[[j]],i]]},{i,1,512}],Table[{(i-1)*n2,speckin[runs[[j+2]]][[tti[[j+2]],i]]},{i,1,512}]},Joined->True,PlotRange->plotRangeskin[[j]],PlotStyle->{{Green,Thickness[0.005]},{Black,Thickness[0.005]},{Blue,Thickness[0.005]},{Red,Thickness[0.005]}},Frame->True,FrameStyle->Directive[Black,4],FrameLabel->{If[j==3,Style["k",Black,6,FontSlant->Italic]],If[j==1,Style["(Subscript[E, K])",Black,6]]},PlotLegends->Placed[LineLegend[{Style["A"<>ToString[If[j==1,j,j+1]]<>"'",4],Style["A"<>ToString[If[j==1,j,j+1]],4]},Spacings->0,LegendMarkerSize->{10,0.1}],{0.32,0.28}],Epilog->{If[j==1,Text[Style["(a)",5],{4,-23.5}]],If[j==2,Text[Style["(b)",5],{6,-9.5}]]}],{j,1,2}]

plott=GraphicsGrid[{{plotkin[1],plotkin[2]},{plotgw[1],plotgw[2]}},Spacings->-30,AspectRatio->1/GoldenRatio]

Export[outdir<>"acoustic_vs_acoustic.eps",plott]

Figure 3

sim1=SetDirectory["/home/ramkishor/cvs_axel/jani/GW/kinflow/"];
spec={"kin","GWs"};
{th1,spec1}=read1D["512a_2nd_dt5c_cos_dt","power_"<>spec[[1]]<>".dat"];
{th2,spec2}=read1D["512a_2nd_dt5c_cos_dt","power_"<>spec[[2]]<>".dat"];

cs=(1/(Sqrt[3]));
kp=30;
omegabar=1;
nn=0;
power[k_]:=(1/(k^2))(((((k/kp)))^(4+nn))/(1+((((k/kp)))^(6+nn))))
func3[k_,t_]:=((4 (omegabar^2))/(4(π^2)k))*(1/2)NIntegrate[((q^3)/q1)(((1-((((((q^2)+(k^2)-(q1^2))/(2 k q))))^2)))^2)power[q]power[q1](-((2 (-(k^6)+(cs^2) (k^4) (q^2)+(cs^4) (k^2) (q^4)-(cs^6) (q^6)+(cs^2) (k^4) (q1^2)-10 (cs^4) (k^2) (q^2) (q1^2)+(cs^6) (q^4) (q1^2)+(cs^4) (k^2) (q1^4)+(cs^6) (q^2) (q1^4)-(cs^6) (q1^6)+2 cs k q ((k^4)+2 (cs^2) (k^2) (-(q^2)+(q1^2))+(cs^4) ((q^4)+2 (q^2) (q1^2)-3 (q1^4))) Cos[cs q1 t] Sin[k t] Sin[cs q t]+2 cs (k^5) q1 Cos[cs q t] Sin[k t] Sin[cs q1 t]+4 (cs^3) (k^3) (q^2) q1 Cos[cs q t] Sin[k t] Sin[cs q1 t]-6 (cs^5) k (q^4) q1 Cos[cs q t] Sin[k t] Sin[cs q1 t]-4 (cs^3) (k^3) (q1^3) Cos[cs q t] Sin[k t] Sin[cs q1 t]+4 (cs^5) k (q^2) (q1^3) Cos[cs q t] Sin[k t] Sin[cs q1 t]+2 (cs^5) k (q1^5) Cos[cs q t] Sin[k t] Sin[cs q1 t]+Cos[k t] (((k^6)-(cs^2) (k^4) ((q^2)+(q1^2))+(cs^6) ((((q^2)-(q1^2)))^2) ((q^2)+(q1^2))-(cs^4) (k^2) ((q^4)-10 (q^2) (q1^2)+(q1^4))) Cos[cs q t] Cos[cs q1 t]+2 (cs^2) q q1 (-3 (k^4)+(cs^4) ((((q^2)-(q1^2)))^2)+2 (cs^2) (k^2) ((q^2)+(q1^2))) Sin[cs q t] Sin[cs q1 t])))/((((k+cs (q-q1)))^2) (((k+cs (-q+q1)))^2) (((k-cs (q+q1)))^2) (((k+cs (q+q1)))^2)))),{q,0,10 *kp},{q1,Abs[q-k],q+k},Method->"AdaptiveMonteCarlo"]
dataa[t_]:=ParallelTable[{(10^n),((((10^n)))^2)*(π^2)*func3[(10^n),t]},{n,0,3,0.01}];
n={1,1400};
n={18,2000};
gwnorm=1;
kp=30;
omega0={((((Total[spec1[[n[[1]],All]]]*2))^2)/(kp^2))*(54/(π^2))((((4/3)))^2),((((Total[spec1[[n[[2]],All]]]*2))^2)/(kp^2))*(54/(π^2))((((4/3)))^2)}

plot1=GraphicsGrid[{{Show[{ListLogLogPlot[Table[{(i-1),spec1[[n[[1]],i]]},{i,1,256}],Joined->True,PlotRange->{{1,300},{(10^(-7)),(10^(-1))}},PlotStyle->{{Blue,Dotted,Thickness[0.01]}},Frame->True,FrameStyle->Directive[Black,8],FrameLabel->{Style["k",Black,8],Style["(Subscript[E, K])",Black,8]},PlotLegends->Placed[LineLegend[{Style["η="<>ToString[N[th2[[n[[1]]]],1]],8]},LegendMarkerSize->8],{0.48,0.2}]],LogLogPlot[5.5*(10^(-2))(k^2)power[k],{k,1,300},PlotStyle->{Red,Thickness[0.01]}],LogLogPlot[(10^(-6.5))(k^4),{k,1,15},PlotStyle->Gray],Epilog->{Text["(k^4)",{1.1,-7.9}],Text[Style["(a)",7],{0.5,-3}]}}],Show[{ListLogLogPlot[Table[{(i-1),(1/(6omega0[[1]]))*spec2[[n[[1]],i]]},{i,1,256}],Joined->True,PlotRange->{{1,300},{(10^(-8)),0.5}},PlotStyle->{{Blue,Dotted,Thickness[0.01]}},Frame->True,FrameStyle->Directive[Black,8],FrameLabel->{Style["k",Black,8,FontSlant->Italic],Style["(Subscript[Ω, GW])/k (Subscript[Ω, 0])",Black,8]},PlotLegends->Placed[LineLegend[{Style["η="<>ToString[N[th2[[n[[1]]]],1]],8]},LegendMarkerSize->8],{0.25,0.85}]],ListLogLogPlot[Transpose[{dataa[th2[[n[[1]]]]][[All,1]],gwnorm*dataa[th2[[n[[1]]]]][[All,2]]}],PlotStyle->{Red,Thickness[0.01]}]},Epilog->{Text[Style["(b)",7],{5,-2}]}]},{Show[{ListLogLogPlot[Table[{(i-1),spec1[[n[[2]],i]]},{i,1,256}],Joined->True,PlotRange->{{1,3*(10^2)},{(10^(-7)),(10^(-1))}},PlotStyle->{{Blue,Dotted,Thickness[0.01]}},Frame->True,FrameStyle->Directive[Black,8],FrameLabel->{Style["k",Black,8],Style["(Subscript[E, K])",Black,8]},PlotLegends->Placed[LineLegend[{Style["η="<>ToString[N[th2[[n[[2]]]],4]],8]},LegendMarkerSize->8],{0.48,0.2}]],LogLogPlot[5.5*(10^(-2))(k^2)power[k],{k,1,300},PlotStyle->{Red,Thickness[0.01]}],LogLogPlot[(10^(-6.5))(k^4),{k,1,15},PlotStyle->Gray],Epilog->{Text["(k^4)",{1.1,-7.9}],Text[Style["(c)",7],{0.5,-3}]}}],Show[{ListLogLogPlot[Table[{(i-1),(1/(6omega0[[2]]))*spec2[[n[[2]],i]]},{i,1,256}],Joined->True,PlotRange->{{1,300},{(10^(-6)),2*10}},PlotStyle->{{Blue,Dotted,Thickness[0.01]}},Frame->True,FrameStyle->Directive[Black,8],FrameLabel->{Style["k",Black,8,FontSlant->Italic],Style["(Subscript[Ω, GW])/k (Subscript[Ω, 0])",Black,8]},PlotLegends->Placed[LineLegend[{Style["η="<>ToString[N[th2[[n[[2]]]],4]],8]},LegendMarkerSize->8],{0.48,0.12}]],ListLogLogPlot[Transpose[{dataa[th2[[n[[2]]]]][[All,1]],gwnorm*dataa[th1[[n[[2]]]]][[All,2]]}],PlotStyle->{Red,Thickness[0.01]}]},Epilog->{Text[Style["(d)",7],{5,2}]}]}},Spacings->0,AspectRatio->1/GoldenRatio]

Export[outdir<>"spectrum2_512a_2nd_dt5c_cos_dt.eps",plot1]