#!/usr/bin/gnuplot 
#
# created: 2010May25 /JD
# lastupd: 2010May25 /JD
#          
# purpose: 
# ---------------------------------------------------------
set term post eps color solid "Helvetica" 26
set size 2.3,1.2
`mkdir EPS 2>/dev/null`

#set format x "%2.2f"
#set format y "%2.2f"
set grid y x
set logscale x

set  title "<pair> - ZTD bias dependence on height/latitude" font "HelveticaBold,32"
set ylabel "ZTD bias [mm]" offset 0.5,0
set xlabel "height [log km]     /     latitude [log deg]" offset 0.5,0
set output "EPS/dep_ztd_<pair>.eps"

plot 0 not w lines lt -1,\
     'VAR/ztd_<pair>.dep'  u ($2):6      not w points lt 1 lw 1 pt 9 ps 1.5,\
     'VAR/ztd_<pair>.dep'  u ($2):5      not w points lt 3 lw 3 pt 8 ps 1.5,\
     'VAR/ztd_<pair>.dep'  u ($4/1000):6 not w points lt 1 lw 1 pt 7 ps 1.2,\
     'VAR/ztd_<pair>.dep'  u ($4/1000):5 not w points lt 3 lw 3 pt 6 ps 1.2

`/usr/bin/convert -colorspace RGB  -density 300 EPS/dep_ztd_<pair>.eps EPS/dep_ztd_<pair>.png`
`scp EPS/dep_ztd_<pair>.png moses@pecny.pecny.cz:WWW_FIL/GOP/`
`scp EPS/dep_ztd_<pair>.eps moses@pecny.pecny.cz:WWW_FIL/GOP/`


quit


# GO1-GO0
# -------
set  title "GO1-GO0 - ZTD bias dependence on height/latitude" font "HelveticaBold,32"
set ylabel "ZTD bias [mm]" offset 0.5,0
set xlabel "latitude [log deg] / altitude [log km]" offset 0.5,0
set output "EPS/dep_ztd_R2_GO1-R2_GO0.eps"

plot 0 not w lines lt -1,\
     'VAR/ztd_R2_GO1-R2_GO0.dep'  u ($2):5      not w points lt 3 lw 3 pt 6 ps 1.0,\
     'VAR/ztd_R2_GO1-R2_GO0.dep'  u ($2):6      not w points lt 1 lw 1 pt 7 ps 1.0,\
     'VAR/ztd_R2_GO1-R2_GO0.dep'  u ($4/1000):5 not w points lt 2 lw 3 pt 6 ps 1.0,\
     'VAR/ztd_R2_GO1-R2_GO0.dep'  u ($4/1000):6 not w points lt 4 lw 1 pt 7 ps 1.0

`/usr/bin/convert -colorspace RGB  -density 300 EPS/dep_ztd_R2_GO1-R2_GO0.eps EPS/dep_ztd_R2_GO1-R2_GO0.png`
`scp EPS/dep_ztd_R2_GO1-R2_GO0.png moses@pecny.pecny.cz:WWW_FIL/GOP/`
