Science/clag-agn
1
0
Fork 0
mirror of https://asciireactor.com/otho/clag-agn.git synced 2024-11-21 22:45:06 +00:00
Code Issues Packages Projects Releases Wiki Activity

found convergence on 7 models with 4 bins

Browse Source
This commit is contained in:
caes 2017-07-05 09:57:28 -04:00
parent 834e5f0192
commit f4b3f44f90
4 changed files with 29 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
notes
View File

@ -6,3 +6,11 @@ Ok, a quick summary of what we talked about doing for next week:
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
4368A is suspect
4392A is suspect

17
scripts/process_output.pl
View File

@ -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 "

9
scripts/templates/lag_atlas.gp
View File

@ -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'

13
scripts/templates/psd_atlas.gp
View File

@ -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