From f4b3f44f90697a1af26c5677b53ef7f2227ea45a Mon Sep 17 00:00:00 2001
From: caes <caes>
Date: Wed, 5 Jul 2017 09:57:28 -0400
Subject: [PATCH] found convergence on 7 models with 4 bins

---
 notes                          | 10 +++++++++-
 scripts/process_output.pl      | 17 ++++++++++-------
 scripts/templates/lag_atlas.gp |  9 +++++----
 scripts/templates/psd_atlas.gp | 13 +++++--------
 4 files changed, 29 insertions(+), 20 deletions(-)

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