source:
LMDZ5/branches/LF-private/tools/Max_diff_nc_with_lib/NR_util/unit_matrix.f90
@
3069
Last change on this file since 3069 was 1765, checked in by , 11 years ago | |
---|---|
File size: 326 bytes |
Line | |
---|---|
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.