GMT version! Requires GMT (netcdf/gmt/ghostscript) 
  
 Calculate clock parameters from   k-file(s) 
  
    -- REQUIRES  x-files 
  
 Usage      :  sh_plotk -res -obs -model -jump -njump -polyj <order> -polynj <order> -files <kfile(s)>
 Where      :  -res   = plot residual clock [DEFAULT] 
            :  -obs   = plot observed clock  
            :  -model = plot model clock  
            :  Both -obs and -model can be plotted on the same page 
  
            :  -jump  = plot with jumps removed from estimates [DEFAULT] 
            :  -njump = plot without jumps removed from estimates 
            :  Both -jump and -njump can be plotted on the same page 
  
            :  -polyj  <order> of polynomial to be estimated for clock with jumps [DEFAULT c] 
            :  -polynj <order> of polynomial to be estimated for clock without jumps [DEFAULT l]
            :  Specify the order of polynomials... linear, quadratic or cubic (l,q,c or 1,2,3)
  
            :  -files  <kfile(s) list of kfiles to be plotted 
  
 Examples   : sh_plotk -res -jump -njump -polyj c -polynj l -files kkata2.140 
            : sh_plotk -model -obs -jump -njump -polyj 3 -polynj 1 -files kkata2.140 
            : sh_plotk -obs -files kkata2.140   
  
 Statistics :   cat CLOCKS 
  
 Max resid. :   cat CLOCKS | grep max  
  
 creates       plotk.site.ps postcript files