PDAF_generate_rndmat

This page documents the routine PDAF_generate_rndmat of PDAF.

The routine returns a random matrix with special properties.

The routine provides a special functionality and can e.g. be used to generate a random matrix to perturb model forcings.

The interface is the following:

  SUBROUTINE PDAF_generate_rndmat(dim, rndmat, mattype)

with the following arguments:

  INTEGER, INTENT(in) :: dim         ! Size of matrix rndmat
  REAL, INTENT(out)   :: rndmat(dim, dim) ! Matrix
  INTEGER, INTENT(in) :: mattype     ! Select type of random matrix:
                                     !   (1) orthonormal random matrix
                                     !   (2) orthonormal with eigenvector (1,...,1)^T

Note:

• This matrix is e.g. used internally by PDAF_sampleens for second-order exact sampling.