source:
LMDZ5/trunk/tools/Max_diff_nc_with_lib/NR_util/unit_matrix.f90
@
1874
Last change on this file since 1874 was 1765, checked in by , 11 years ago | |
---|---|
File size: 326 bytes |
Rev | Line | |
---|---|---|
[1765] | 1 | module unit_matrix_m |
2 | ||
3 | implicit none | |
4 | ||
5 | contains | |
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 | ||
18 | end module unit_matrix_m |
Note: See TracBrowser
for help on using the repository browser.