home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: Science / Science.zip / gmt_os2.zip / src / surface.out < prev    next >
Encoding:
Text File  |  1996-09-25  |  3.5 KB  |  63 lines
  1. GMT Warning: Could not determine home directory!
  2. surface 3.0a:OS/2 - Adjustable tension continuous curvature surface gridding
  3.  
  4. usage: surface [xyz-file] -G<output_grdfile_name> -I<xinc>[m|c][/<yinc>[m|c]]
  5.     -R<west>/<east>/<south>/<north> [-A<aspect_ratio>] [-C<convergence_limit>]
  6.     [-H] [-Ll<limit>] [-Lu<limit>] [-N<n_iterations>] ] [-S<search_radius>[m]]
  7.     [-T<tension>[i][b] ] [-Q] [-V[l]] [-Z<over_relaxation_parameter>] [-:] [-b[d]]
  8.  
  9.     surface will read from standard input or [xyz-file].
  10.  
  11.     Required arguments to surface:
  12.     -G sets output grd File name
  13.     -I sets the Increment of the grid; enter xinc, optionally xinc/yinc.
  14.         Default is yinc = xinc.  Append an m [or c] to xinc or yinc to indicate minutes [or seconds]
  15.         e.g.  -I10m/5m grids longitude every 10 minutes, latitude every 5 minutes.
  16.  
  17.     -R specifies the min/max coordinates of data region in user units.
  18.        Use dd:mm[:ss] format for regions given in degrees and minutes [and seconds].
  19.        Append r if -R specifies the longitudes/latitudes of the lower left
  20.        and upper right corners of a rectangular area
  21.  
  22.     OPTIONS:
  23.     -A<aspect_ratio>  = 1.0  by default which gives an isotropic solution.
  24.         i.e. xinc and yinc assumed to give derivatives of equal weight; if not, specify
  25.         <aspect_ratio> such that yinc = xinc / <aspect_ratio>.
  26.         e.g. if gridding lon,lat use <aspect_ratio> = cosine(middle of lat range).
  27.     -C<convergence_limit> iteration stops when max abs change is less than <c.l.>
  28.         default will choose 0.001 of the range of your z data (1 ppt precision).
  29.         Enter your own convergence limit in same units as z data.
  30.     -H means input/output file has 0 Header record(s) [OFF]
  31.        Optionally, append number of header records
  32.     -L constrain the range of output values:
  33.         -Ll<limit> specifies lower limit; forces solution to be >= <limit>.
  34.         -Lu<limit> specifies upper limit; forces solution to be <= <limit>.
  35.         <limit> can be any number, or the letter d for min (or max) input data value,
  36.         or the filename of a grdfile with bounding values.  [Default solution unconstrained].
  37.         Example:  -Ll0 gives a non-negative solution.
  38.     -N sets max <n_iterations> in each cycle; default = 250.
  39.     -S sets <search_radius> to initialize grid; default = 0 will skip this step.
  40.         This step is slow and not needed unless grid dimensions are pathological;
  41.         i.e., have few or no common factors.
  42.         Append m to give <search_radius> in minutes.
  43.     -T adds Tension to the gridding equation; use a value between 0 and 1.
  44.         default = 0 gives minimum curvature (smoothest; bicubic) solution.
  45.         1 gives a harmonic spline solution (local max/min occur only at data points).
  46.         typically 0.25 or more is good for potential field (smooth) data;
  47.         0.75 or so for topography.  Experiment.
  48.         Append B or b to set tension in boundary conditions only;
  49.         Append I or i to set tension in interior equations only;
  50.         No appended letter sets tension for both to same value.
  51.     -Q Query for grid sizes that might run faster than your -R -I give.
  52.     -V Run in verbose mode [OFF].
  53.         Append l for long verbose
  54.     -Z sets <over_relaxation parameter>.  Default = 1.4
  55.         Use a value between 1 and 2.  Larger number accelerates convergence but can be unstable.
  56.         Use 1 if you want to be sure to have (slow) stable convergence.
  57.  
  58.     -: Expect lat/lon input rather than lon/lat [OFF].
  59.     -b means input triplets are binary.  Default is ascii.
  60.        Append d for double precision [Default is single precision].
  61.     (See man page for gmtdefaults to set these and other defaults to ON or OFF)
  62.     (For additional details, see Smith & Wessel, Geophysics, 55, 293-305, 1990.)
  63.