#!/bin/csh -f

set here = `pwd`
set site = ""
set ts_dir = "$here/tsdata.xyz"
set psdir = "$here/psfiles"
set jpgdir = "$here/jpegfiles.xyz"
set pngdir = "$here/pngfiles.xyz"
foreach d ( $ts_dir $psdir $jpgdir $pngdir )
  if ( ! -d $ts_dir ) mkdir -p $d
end
if ( $1 != "" ) set site = `echo $1| tr '[A-Z]' '[a-z]'`
if ( $2 != "" ) set ts_dir = $2
set xrange = "-xr 2004.92:2005"
set yrange = "-yr -20:20"
set xrange = ""
set yrange = ""
#set options = "-term ps -jpg -yu mm -eb   -event slow,step,eqk -dim 500x650 -fontsize 10 $xrange $yrange "
#set options = "-term ps -jpg -yu mm -eb   -event slow,step,eqk -dim 600x800 -fontsize 10 $xrange $yrange "
set options = "-term sp -png -yu mm -eb -event slow,step,eqk -dim 600x800 -fontsize 10 $xrange $yrange "
\rm $psdir/$site* $jpgdir/$site* $pngdir/$site* >& /dev/null


foreach s ( `ls $ts_dir/????.data.x|sort -r |grep "$site"` )

	set stn = `echo $s:t | awk -F. '{print $1}'| tr '[A-Z]' '[a-z]'`
	set STN = `echo $stn | tr '[a-z]' '[A-Z]'`
	echoc "plotting $s:t ..."

#	fit_nev $stn l 0 0 >& $fitdir/$stn.lfit.all.log
	
	#gplot -i $ts_dir/$stn.data.lon  -pf line,offsets -rmf jumps $options
	gplot -i $ts_dir/$stn.data.x  -pf line,offsets,jumps  $options
	gplot -i $ts_dir/$stn.resd.x    $options 
	gplot -i $ts_dir/$stn.ldet.x   -pf annual,jumps $options -psfile $ts_dir/$stn.ldet.png
	
	#mv $ts_dir/$stn.data.jp* $jpgdir
	#mv $ts_dir/$stn.data.ps  $psdir
        mv $ts_dir/$stn.data.png $pngdir
	#mv $ts_dir/$stn.resd.jp* $jpgdir
	#mv $ts_dir/$stn.resd.ps  $psdir
        mv $ts_dir/$stn.resd.png $pngdir
	#mv $ts_dir/$stn.ldet.jp* $jpgdir
	#mv $ts_dir/$stn.ldet.ps  $psdir
        mv $ts_dir/$stn.ldet.png $pngdir
    \rm $ts_dir/*.ps >& /dev/null
	echo " Done"
end

