Com

Project.toml

@@ -20,7 +20,10 @@ JSServe = "824d6782-a2ef-11e9-3a09-e5662e0c26f9"
 NCDatasets = "85f8d34a-cbdd-5861-8df4-14fed0d494ab"
 NodeJS = "2bd173c7-0d6d-553b-b6af-13a54713934c"
 PrettyTables = "08abe8d2-0d0c-5749-adfa-8a2ac140af0d"
-Rasters = "a3a2b9e3-a471-40c9-b274-f788e487c689"
+Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 RasterDataSources = "3cb90ccd-e1b6-4867-9617-4276c8b2ca36"
+Rasters = "a3a2b9e3-a471-40c9-b274-f788e487c689"
 RemoteS = "26020ff4-e92a-488e-a723-09c28e89e04a"
-Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
+SatelliteToolboxPropagators = "c2b69894-ea78-4e2b-9ba6-cedbbc3d14d7"
+SatelliteToolboxTle = "7ff27aeb-5fff-4337-a9ee-a9fe6b7ed35e"
+SatelliteToolboxTransformations = "6b019ec1-7a1e-4f04-96c7-a9db1ca5514d"

documentation/modules/grdgradient.md

@@ -23,7 +23,7 @@ The 2-D gridded data set to be contoured.
 Optional Arguments
 ------------------
 
-- **A** or **azim** : -- *azim=azim* **|** *azim=(azim1, azim2)*\
+- **A** or **azim** or **azimuth** : -- *azim=azim* **|** *azim=(azim1, azim2)*\
  Azimuthal direction for a directional derivative; *azim* is the angle in the x,y plane measured in degrees positive clockwise from north (the +y direction) toward east (the +x direction). The negative of the directional derivative, *-[dz/dx sin(azim) + dz/dy cos(azim)]*, is found; negation yields positive values when the slope of z(x,y) is downhill in the *azim* direction, the correct sense for shading the illumination of an image (see \myreflink{grdimage} and \myreflink{grdview}) by a light source above the x,y plane shining from the *azim* direction. Optionally, supply two azimuths, **azim=(azim1, azim2)**, in which case the gradients in each of these directions are calculated and the one larger in magnitude is retained; this is useful for illuminating data with two directions of lineated structures, e.g., **azim=(0,270)** illuminates from the north (top) and west (left). Finally, if *azim* is a file it must be a grid of the same domain, spacing and registration as *in\_grdfile* and we will update the azimuth at each output node when computing the directional derivatives.
 
 \textinput{common_opts/opt_save_grd}

documentation/modules/grdview.md

@@ -141,3 +141,44 @@ Examples
 --------
 
 See the Example 04 at the Historical Collection gallery.
+
+
+Make a low resolution image plot with texture draping. Use a dpi=50 to make it lighter to load the man page.
+The drape option is used to overlay an image texture (from a cork tree trunk) on the grid data from the
+Gulf of Guinea zone. (Note, the `viz` function calls `grdview` when it detects the `view` option)
+
+\begin{examplefig}{}
+```julia
+using GMT
+
+G = grdcut("@earth_relief_01m", region=(0,10,0,10));
+viz(G, shade="+nt0.5", proj=:merc, view=(145,35), drape="@wood_texture.jpg", zsize=5, surftype=(image=50,))
+```
+\end{examplefig}
+
+This example imagines that a OVNI (UFO) ship is pointing a spot light with a radius of 100 km over Denver.
+
+\begin{examplefig}{}
+```julia
+using GMT
+
+G = grdcut("@earth_relief_30s", region=(-108,-103,35,40));
+
+# Compute the distances to Denver (here we need to use the terse syntax)
+Gr = grdmath("-R? -104.9903 39.7392 SDIST", G);
+
+# Mask distances > 100 km
+Ginside = Gr < 100;
+
+# Compute gradients along two directions to use in shading.
+Gint1 = grdgradient(G, azim=90, norm=(cauchy=true, amp=1));
+Gint2 = grdgradient(G, azim=45, norm=(cauchy=true, amp=1));
+
+# Combine the gradients
+Gint = Gint2 * 0.5 + Gint1 + Ginside - 0.5;
+
+# View the result
+C = makecpt(range=(1000,4000));
+grdview(G, shade=Gint, proj=:merc, zsize=1, view=(155,25), colormap=C, surftype=:image, show=true)
+```
+\end{examplefig}

documentation/modules/histogram.md

@@ -33,6 +33,9 @@ Optional Arguments
  from a predefined function name. Possible names are: *akley, eggs, circle, parabola, rosenbrok, sombrero, x, y, xy, x+y*.
  For further details consult the same option in the plot manual.
 
+- **B** or **axes** or **frame**\
+ Set map boundary frame and axes attributes. Default is to draw and annotate left, bottom and vertical
+ axes and just draw left and top axes. More at \myreflink{frame}
 
 - **C** or **color** or **cmap** : -- *color=cpt*\
  Give a CPT or specify **color="color1,color2 [,color3 ,...]"** or **color=((r1,g1,b1),(r2,g2,b2),...)** to build a linear continuous CPT from those colors automatically. The mid x-value for each bar is used to look-up the bar color. If no argument is given then we select the current CPT.

documentation/modules/legend.md

@@ -14,7 +14,7 @@ file or from a GMTdatset type. Unless otherwise noted, annotations will be made
 annotation font and size in effect (i.e., `FONT_ANNOT_PRIMARY`)
 
 ~~~
-<img src="/assets/docs/seislegend.png" alt="Simple" width="500" class="center"/>
+<img src="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/assets/seislegend.png" alt="Simple" width="500" class="center"/>
 ~~~
 
 Example of a map legend, here for seismicity in the Pacific region.

examples/misc/orbits.md

@@ -19,8 +19,8 @@ orbit descriptor.
 
 \begin{examplefig}{}
 ```julia
-import Pkg; Pkg.add("RemoteS")  # Hide
-using GMT, RemoteS
+import Pkg; Pkg.add("RemoteS"); Pkg.add("SatelliteToolboxTle"); Pkg.add("SatelliteToolboxPropagators"); Pkg.add("SatelliteToolboxTransformations"); # Hide
+using GMT, RemoteS, SatelliteToolboxTle, SatelliteToolboxPropagators, SatelliteToolboxTransformations
 
 # Compute some Aqua orbits for ... 
 orb = sat_tracks(tle=["1 27424U 02022A 23229.43455053 .00001116 00000-0 25175-3 0 9994";