Line data Source code
1 : module mo_lu_factor 2 : private 3 : public :: lu_fac 4 : contains 5 4348047600 : subroutine lu_fac01( lu ) 6 : use shr_kind_mod, only : r8 => shr_kind_r8 7 : implicit none 8 : !----------------------------------------------------------------------- 9 : ! ... dummy args 10 : !----------------------------------------------------------------------- 11 : real(r8), intent(inout) :: lu(:) 12 4348047600 : lu(1) = 1._r8 / lu(1) 13 4348047600 : lu(2) = 1._r8 / lu(2) 14 4348047600 : lu(3) = 1._r8 / lu(3) 15 4348047600 : lu(4) = 1._r8 / lu(4) 16 4348047600 : lu(5) = 1._r8 / lu(5) 17 4348047600 : lu(6) = lu(6) * lu(5) 18 4348047600 : lu(7) = 1._r8 / lu(7) 19 4348047600 : lu(8) = lu(8) * lu(7) 20 4348047600 : lu(9) = 1._r8 / lu(9) 21 4348047600 : lu(10) = 1._r8 / lu(10) 22 4348047600 : lu(11) = 1._r8 / lu(11) 23 4348047600 : lu(12) = 1._r8 / lu(12) 24 4348047600 : lu(13) = lu(13) * lu(12) 25 4348047600 : lu(14) = 1._r8 / lu(14) 26 4348047600 : lu(15) = 1._r8 / lu(15) 27 4348047600 : lu(16) = 1._r8 / lu(16) 28 4348047600 : lu(17) = 1._r8 / lu(17) 29 4348047600 : lu(18) = 1._r8 / lu(18) 30 4348047600 : lu(19) = 1._r8 / lu(19) 31 4348047600 : lu(20) = 1._r8 / lu(20) 32 4348047600 : lu(21) = 1._r8 / lu(21) 33 4348047600 : lu(22) = 1._r8 / lu(22) 34 4348047600 : lu(23) = 1._r8 / lu(23) 35 4348047600 : lu(24) = 1._r8 / lu(24) 36 4348047600 : lu(26) = 1._r8 / lu(26) 37 4348047600 : lu(27) = 1._r8 / lu(27) 38 4348047600 : lu(28) = 1._r8 / lu(28) 39 4348047600 : lu(29) = 1._r8 / lu(29) 40 4348047600 : lu(30) = 1._r8 / lu(30) 41 4348047600 : lu(31) = 1._r8 / lu(31) 42 4348047600 : lu(32) = 1._r8 / lu(32) 43 4348047600 : lu(33) = lu(33) * lu(32) 44 4348047600 : lu(34) = 1._r8 / lu(34) 45 4348047600 : lu(35) = 1._r8 / lu(35) 46 4348047600 : end subroutine lu_fac01 47 4348047600 : subroutine lu_fac( lu ) 48 : use shr_kind_mod, only : r8 => shr_kind_r8 49 : implicit none 50 : !----------------------------------------------------------------------- 51 : ! ... dummy args 52 : !----------------------------------------------------------------------- 53 : real(r8), intent(inout) :: lu(:) 54 4348047600 : call lu_fac01( lu ) 55 4348047600 : end subroutine lu_fac 56 : end module mo_lu_factor