SUM (array [,dim] [,mask])

  Class:  Transformational function - Generic

  Returns the sum of all the elements in an entire array or in a
  specified dimension of an array.

  The "array" can be of integer or real type.  The "dim" is optional
  and must be a scalar integer with a value in the range 1 to n,
  where "n" is the rank of "array".  The "mask" is optional and must
  be a logical array that is conformable with "array".

  The result is the same data type as "array".  The result is a
  scalar if "dim" is absent or "array" has rank one.  Otherwise, the
  result is an array with rank that is one less than "array", and
  shape (d1, d2,..., d"dim"-1, d"dim"+1,..., dn), where (d1, d2,...,
  dn) is the shape of "array".

  If only "array" appears, the result is the sum of all elements of
  "array".  If "array" has size zero, the result is zero.

  If "array" and "mask" both appear, the result is the sum of all
  elements of "array" corresponding to true elements of "mask".  If
  "array" has size zero, or every element of "mask" has the value
  .FALSE., the result is zero.

  If "dim" also appears and "array" has rank one, the value is the
  same as SUM (array [,mask=mask]).  Otherwise, the value of element
  (s1, s2,..., s"dim"-1, s"dim"+1,..., sn) of SUM (array, dim,
  [,mask]) is equal to SUM (array (s1, s2,..., s"dim"-1, :, s"dim"+1,
  ..., sn)) [mask=mask (s1, s2, ..., s"dim"-1, :, s "dim"+1 ...,
  sn)].

  Examples:

  SUM ((/2, 3, 4/)) and SUM ((/2, 3, 4/), DIM=1) returns the value 9.

  SUM (B, MASK=B .LT.  0.0) returns the arithmetic sum of the
  negative elements of B.

  Consider that C is the array:

   |1 2 3|
   |4 5 6|.

  SUM (C, DIM=1) returns the value (5, 7, 9).

  SUM (C, DIM=2) returns the value (6, 15).