diff --git a/notes b/notes index 3e4f4d7..c0623ca 100644 --- a/notes +++ b/notes @@ -5,4 +5,12 @@ Ok, a quick summary of what we talked about doing for next week: 3) Try 5 bins and 6 bins: how many wavelengths fully coverage for each? Make 18 panel figures for these too, so we can compare everything 4) Write out ascii files with: low frequency, high frequency, PSD and Lag for each wavelength. This will be what I need to try and fit the lags and PSDs….. -Ed \ No newline at end of file +Ed + + + + + + +4368A is suspect +4392A is suspect \ No newline at end of file diff --git a/scripts/process_output.pl b/scripts/process_output.pl index 125b12a..0e6eb89 100755 --- a/scripts/process_output.pl +++ b/scripts/process_output.pl @@ -18,6 +18,9 @@ our $debug=0; if ($debug) { $verbose=1; } +# @ARGV[0] takes the echo wavelength for use in filename +$echo_λ = @ARGV[0]; +$ref_λ = "3465A"; # This section locates the output data of interest in a # psdlab output file. @@ -27,9 +30,9 @@ if (!${^UTF8LOCALE}) { my $charset=$ENV{LANG}; -# @ARGV[0] takes the echo wavelength for use in filename - -$echo_λ = @ARGV[0]; +my $echopsd_tabfile="data/tables/psd_${echo_λ}.tab"; +my $refpsd_tabfile="data/tables/psd_${ref_λ}.tab"; +my $timelag_tabfile="data/tables/lag_${ref_λ}_${echo_λ}.tab"; # This program attempts to fill function_bin with the tabulated # PSDs and time lags. @@ -138,7 +141,7 @@ if($debug) { close($outputfile); -open($datafile,'>',"tmp.refpsd") or die $!; +open($datafile,'>',$refpsd_tabfile) or die $!; while( each %function_bin) { say $datafile $_ . " " . @@ -148,7 +151,7 @@ while( each %function_bin) { } close($datafile); -open($datafile,'>',"tmp.echopsd") or die $!; +open($datafile,'>',$echopsd_tabfile) or die $!; while( each %function_bin) { say $datafile $_ . " " . @@ -169,7 +172,7 @@ close($datafile); # } # close($datafile); -open($datafile,'>',"tmp.lag") or die $!; +open($datafile,'>',$timelag_tabfile) or die $!; while( each %function_bin) { say $datafile $_ . " " . @@ -183,7 +186,7 @@ close($datafile); #open($datafile,'>',"tmp.echoPSD") or die $!; -open($datafile,'>',"cackett_psdlag_" . $echo_λ . ".tab") or die $!; +open($datafile,'>',"data/tables/cackett_psdlag_" . $echo_λ . ".tab") or die $!; say $datafile encode($charset, sprintf("#freqmin freqmax psd ". "psd error lag lag error " diff --git a/scripts/templates/lag_atlas.gp b/scripts/templates/lag_atlas.gp index 266ac69..ca6d739 100644 --- a/scripts/templates/lag_atlas.gp +++ b/scripts/templates/lag_atlas.gp @@ -5,7 +5,7 @@ set termopt enhanced set macros # Placement of the a,b,c,d labels in the graphs - POS = "at graph 0.61,0.85 font 'Times,12'" + POS = "at graph 0.72,0.88 font 'Times,9'" # x- and ytics for each row resp. column NOXNUMS = "unset xlabel;\ @@ -34,10 +34,11 @@ set macros unset key set logscale x set xtics auto font 'Times,9' offset 0,.5 -set ytics (-1,0,1,2,3,4,5) font 'Times,9' -set ytics add ('' -1.5 1,'' -0.5 1,'' 0.5 1,'' 1.5 1,'' 2.5 1,'' 3.5 1,'' 4.5 1,'' 5.5 1) +set ytics (-4,-2,0,2,4) font 'Times,9' +set ytics add ('' -5.5 1,'' -4.5 1,'' -3.5 1,'' -2.5 1,'' -1.5 1,'' -0.5 1,'' 0.5 1,'' 1.5 1,'' 2.5 1,'' 3.5 1,'' 4.5 1,'' 5.5 1) +set ytics add ('' -5 2,'' -3 2,'' -1 2,'' 1 2,'' 3 2,'' 5 2) set xrange [0.005:0.620]; -set yrange [-2:6] +set yrange [-5:5] # Draw line at origin set arrow from 0.005,0 to 0.620,0 nohead lt 3 lc rgb 'black' diff --git a/scripts/templates/psd_atlas.gp b/scripts/templates/psd_atlas.gp index fa7cbc4..086ecd6 100644 --- a/scripts/templates/psd_atlas.gp +++ b/scripts/templates/psd_atlas.gp @@ -5,7 +5,7 @@ set termopt enhanced set macros # Placement of the a,b,c,d labels in the graphs - POS = "at graph 0.61,0.85 font 'Times,12'" + POS = "at graph 0.72,0.88 font 'Times,9'" # x- and ytics for each row resp. column NOXNUMS = "unset xlabel;\ @@ -33,16 +33,13 @@ set macros unset key set logscale x set xtics auto font 'Times,9' offset 0,.5 -set ytics auto font 'Times,9' -set ytics add ('' -0.5 1,'' -1.5 1,'' -2.5 1,'' -3.5 1,'' -4.5 1) - -# xrange from 2016 REU freq bins -#set xrange [0.005:0.620]; -set yrange [-5.5:0.5] +set ytics (-4,-2,0,2,4) font 'Times,9' +set ytics add ('' 2.5 1,'' 1.5 1,'' 0.5 1,'' -0.5 1,'' -1.5 1,'' -2.5 1,'' -3.5 1,'' -4.5 1) +set ytics add ('' -5 2,'' -3 2,'' -1 2,'' 1 2,'' 3 2,'' 5 2) # xrange pulled from example.py from first version of clag set xrange [0.001:5] -set yrange [-5.5:0.5] +set yrange [-5.5:2.5] set multiplot layout 6,3 rowsfirst