Line data Source code
1 : module mo_lu_factor 2 : private 3 : public :: lu_fac 4 : contains 5 262148400 : 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 262148400 : lu(1) = 1._r8 / lu(1) 13 262148400 : lu(2) = 1._r8 / lu(2) 14 262148400 : lu(3) = 1._r8 / lu(3) 15 262148400 : lu(4) = 1._r8 / lu(4) 16 262148400 : lu(5) = 1._r8 / lu(5) 17 262148400 : lu(6) = lu(6) * lu(5) 18 262148400 : lu(7) = 1._r8 / lu(7) 19 262148400 : lu(8) = lu(8) * lu(7) 20 262148400 : lu(9) = 1._r8 / lu(9) 21 262148400 : lu(10) = 1._r8 / lu(10) 22 262148400 : lu(11) = 1._r8 / lu(11) 23 262148400 : lu(12) = 1._r8 / lu(12) 24 262148400 : lu(13) = lu(13) * lu(12) 25 262148400 : lu(14) = 1._r8 / lu(14) 26 262148400 : lu(15) = 1._r8 / lu(15) 27 262148400 : lu(16) = 1._r8 / lu(16) 28 262148400 : lu(17) = 1._r8 / lu(17) 29 262148400 : lu(18) = 1._r8 / lu(18) 30 262148400 : lu(19) = 1._r8 / lu(19) 31 262148400 : lu(20) = 1._r8 / lu(20) 32 262148400 : lu(21) = 1._r8 / lu(21) 33 262148400 : lu(22) = 1._r8 / lu(22) 34 262148400 : lu(23) = 1._r8 / lu(23) 35 262148400 : lu(24) = 1._r8 / lu(24) 36 262148400 : lu(26) = 1._r8 / lu(26) 37 262148400 : lu(27) = 1._r8 / lu(27) 38 262148400 : lu(28) = 1._r8 / lu(28) 39 262148400 : lu(29) = 1._r8 / lu(29) 40 262148400 : lu(30) = 1._r8 / lu(30) 41 262148400 : lu(31) = 1._r8 / lu(31) 42 262148400 : lu(32) = 1._r8 / lu(32) 43 262148400 : lu(33) = lu(33) * lu(32) 44 262148400 : lu(34) = 1._r8 / lu(34) 45 262148400 : lu(35) = 1._r8 / lu(35) 46 262148400 : end subroutine lu_fac01 47 262148400 : 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 262148400 : call lu_fac01( lu ) 55 262148400 : end subroutine lu_fac 56 : end module mo_lu_factor