add CRS example in proj documentation

docs/source/apps/proj.rst

@@ -174,7 +174,7 @@ also used for supporting files like datum shift files.
 
 .. versionadded:: 9.3.0
 
-    *{crs}* is one of the possibilities accepted by proj_create(), provided it
+    *{crs}* is one of the possibilities accepted by :c::func:`proj_create()`, provided it
     expresses a projected CRS, like a WKT string, an object code (like "EPSG:32632")
     a PROJJSON string, etc.
     The projection computed will be those of the map projection implied by
@@ -212,6 +212,41 @@ data will appear as three lines of::
 
     460770.43     5011865.86
 
+This other example
+
+.. code-block:: console
+
+    proj EPSG:6421 -V <<EOF
+    -120 35.8
+    EOF
+
+Will perform the projection of the coordinates in "NAD83(2011) / California zone 4"
+(`EPSG:6421`) into its geographic system, "NAD83(2011)", showing the expanded annotated listing.
+The output will appear as:
+
+.. code-block:: console
+
+    #Lambert Conformal Conic
+    #	Conic, Sph&Ell
+    #	lat_1= and lat_2= or lat_0, k_0=
+    # +proj=lcc +lat_0=35.3333333333333 +lon_0=-119 +lat_1=37.25 +lat_2=36
+    # +x_0=2000000 +y_0=500000 +ellps=GRS80
+    #Final Earth figure: ellipsoid
+    #  Major axis (a): 6378137.000
+    #  1/flattening: 298.257222
+    #  squared eccentricity: 0.006694380023
+    Longitude: 120dW [ -120 ]
+    Latitude:  35d48'N [ 35.8 ]
+    Easting (x):   1909606.87
+    Northing (y):  552253.58
+    Meridian scale (h) : 1.00004382  ( 0.004382 % error )
+    Parallel scale (k) : 1.00004382  ( 0.004382 % error )
+    Areal scale (s):     1.00008765  ( 0.008765 % error )
+    Angular distortion (w): 0.000
+    Meridian/Parallel angle: 90.00000
+    Convergence : -0d35'47.714" [ -0.59658715 ]
+    Max-min (Tissot axis a-b) scale error: 1.00004 1.00004
+
 .. only:: man
 
     Other programs