source: LMDZ5/branches/testing/tools/Max_diff_nc_with_lib/NR_util/unit_matrix.f90 @ 2232

Last change on this file since 2232 was 1910, checked in by Laurent Fairhead, 11 years ago

Merged trunk changes r1860:1909 into testing branch

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 326 bytes
Line 
1module unit_matrix_m
2
3  implicit none
4
5contains
6
7  SUBROUTINE unit_matrix(mat)
8    USE nrtype, only: wp
9    REAL(WP), DIMENSION(:,:), INTENT(OUT) :: mat
10    INTEGER i,n
11    n=min(size(mat,1),size(mat,2))
12    mat(:,:)=0.0_wp
13    do i=1,n
14       mat(i,i)=1.0_wp
15    end do
16  END SUBROUTINE unit_matrix
17
18end module unit_matrix_m
Note: See TracBrowser for help on using the repository browser.