Reading of zero points.

Gnuplotting.
miroslavbroz · Dec 12, 2024 · 7d902af · 7d902af
1 parent 6e96374
commit 7d902af
Show file tree

Hide file tree

Showing 8 changed files with 48 additions and 6 deletions.
diff --git a/EPS/Ecl12.dat b/EPS/Ecl12.dat
@@ -0,0 +1 @@
+../Ecl12.dat
diff --git a/EPS/chi2_ECL.dat b/EPS/chi2_ECL.dat
@@ -0,0 +1 @@
+../chi2_ECL.dat
diff --git a/EPS/chi2_ECL.plt b/EPS/chi2_ECL.plt
@@ -0,0 +1,27 @@
+#!/usr/bin/gnuplot
+
+set colors classic
+set term post eps enh color dashed
+set out "chi2_ECL.eps"
+set size 1.0,0.7
+
+set xl "JD {/Symbol -} 2400000"
+set yl "eclipse duration [day]" offset +0.5,0
+
+y1=0.24
+set yr [y1:]
+set zeroaxis
+set key right font "Helvetica,14"
+
+load "T0.plt"
+set arrow from T0-2400000,graph 0 to T0-2400000,graph 1 nohead lt 0
+
+p \
+  "duration.dat" u ($1-2400000):2   t "synthetic duration" w lp pt 1 lc 'gray',\
+  "Ecl12.dat"    u ($1-2400000):2:3 t "observed duration" w err lc "blue" pt 1 ps 0.5,\
+  "chi2_ECL.dat" u ($1-2400000):2:3 t "residua" w l lc 'red' lw 3,\
+  "Omc12.dat"    u ($1-2400000):(y1) t "minima" w p pt 2 ps 2 lc 'green',\
+
+q
+
+
diff --git a/EPS/chi2_TTV.dat b/EPS/chi2_TTV.dat
@@ -0,0 +1 @@
+../chi2_TTV.dat
diff --git a/EPS/chi2_TTV.plt b/EPS/chi2_TTV.plt
@@ -9,19 +9,24 @@ set size 1.0,0.7
 day = 1.  # day
 
 set xl "JD {/Symbol -} 2400000"
-set yl "ETV [day] wrt. N-body model" offset +0.5,0
+set yl "TTV [day] wrt. N-body model" offset +0.5,0
 
+#set xr [59000:60500]
 tmp = 0.25*day  # day
-set yr [-tmp:tmp]
+#set yr [-tmp:tmp]
 set zeroaxis
+set key left width +5 font "Helvetica,14"
 
 load "T0.plt"
 set arrow from T0-2400000,graph 0 to T0-2400000,graph 1 nohead lt 0
 
-p "minima.dat" u ($1-2400000):(0) t "synthetic minima (primary and secondary)" w p lt 7 pt 1 ps 1.0,\
-  "Omc12.dat" u ($1-2400000):(0) t "observed minima" w p lt 2 pt 2 ps 2,\
-  "chi2_TTV.dat" u ($1-2400000):($3*day):($5*day) t "O-C" w err lt 1 lw 1 pt 1 ps 0.5,\
-  "Omc12.dat" u ($1-2400000):($3*day) w p t "Mayer et al., Fig. 9" lc 'gray'
+p "minima.dat" u ($1-2400000):(0) t "synthetic minima" w p lc 'gray' pt 1 ps 1.0,\
+  "Omc12.dat" u ($1-2400000):(0) t "observed minima" w p lc 'green' pt 2 ps 2,\
+  "chi2_TTV.dat" u ($1-2400000):($3*day):($5*day) t "O-C" w err lc 'blue' lw 1 pt 1 ps 0.5,\
+  "<awk '!/^#/{ print $1,0.0; print $1,$3; print s; }' chi2_TTV.dat" u ($1-2400000):($2*day) t "residua" w l lc 'red' lw 3,\
+  "<awk '($6>10.0)' chi2_TTV.dat" u ($1-2400000):($3*day) t "chi^2 > 10.0" w p pt 6 ps 2 lc 'red',\
 
 q
 
+  "Omc12.dat" u ($1-2400000):($3*day) w p t "Mayer et al., Fig. 9" lc 'gray'
+
diff --git a/EPS/duration.dat b/EPS/duration.dat
@@ -0,0 +1 @@
+../duration.dat
diff --git a/EPS/minima.dat b/EPS/minima.dat
@@ -0,0 +1 @@
+../minima.dat
diff --git a/chi2/chi2_func.f90 b/chi2/chi2_func.f90
@@ -196,6 +196,11 @@ double precision function chi2_func(x)
   psf_param(i) = x_param(j)
 enddo
 
+do i = 1, nband
+  j = j+1
+  zero(i) = x_param(j)
+enddo
+
 j = j+1
 gamma = x_param(j)
 j = j+1