Use GMT to make a multiple panel GAMIT atmospheric parameters  plots  with uncertainty scaling
  
  requirements  : GMT(netcdf/gmt/ghostscript)
                : gamit/utils/poly01
  
  USAGE :    sh_gamit_atmos   <options>. 
  
  EXAMPLE:   sh_gamit_atmos  -o o* -sites PIN2 VNDP PIN2 YAM2
  
=+=+=+=+=+=+=+=+=+=+=+=+=+ OPTIONS +=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+
 
 -o[files]     list    : o-file names. MUST have this.
   		     
 -f[ile]       file    : Site names. 
   		     
 -g[rad]               : command line argument for plotting graidents instead of atmospheric zenith delays.    
   		     
 -sites   list    : command line argument for selected sites.    
   		     
   		      IF -f and -b are both omitted, all sites are plotted.
   		     
 -free  		    : Biases free  solution Default.
 -fixed  		: Biases fixed solution.
 -both 		    : Biases free with fixed solution. Fixed is dashed. stats refere to free sol.
   		     
 -d[elete]    <sites etc> :  names of unwanted sites.
   		     
 -t[otal]        : total atmospheric zenith delay (Default is adjustment to a priori model).
  
 -df              file   : file contains the names of unwanted sites. M
 -[delete_file]          : Makes life easier for -d
 
 -u[nc_scale]   value   : scale all uncertainties with this number. Default is 1.
 
 -y[scale]   min max   : vertical   scale.  If not issued it will be calculated.
  
 -x[scale]   min max   : horizontal scale.  If not issued it will be calculated.
  
 -frame      value     : gmt border day-axis frame  ticks. Default is  1. 
  
 -anot       value     : gmt border day-axis label intervals. Default is  1. 
  
 -a         ext        : Add more descriptors to postscript file name. Default is psgamit.#