@@ -725,8 +725,8 @@ contains
725725
726726 if (q == 0 ) then
727727 jac_rhs(j, k, l) = alf_igr* (2._wp * (dvel(1 , 2 )* dvel(2 , 1 )) &
728- + dvel(1 , 1 )** 2_wp + dvel(2 , 2 )** 2_wp &
729- + (dvel(1 , 1 ) + dvel(2 , 2 ))** 2_wp )
728+ + dvel(1 , 1 )** 2._wp + dvel(2 , 2 )** 2._wp &
729+ + (dvel(1 , 1 ) + dvel(2 , 2 ))** 2._wp )
730730 end if
731731 end do
732732
@@ -838,8 +838,8 @@ contains
838838 dvel(:, 2 ) = dvel_small
839839
840840 jac_rhs(j, k, l) = alf_igr* (2._wp * (dvel(1 , 2 )* dvel(2 , 1 )) &
841- + dvel(1 , 1 )** 2_wp + dvel(2 , 2 )** 2_wp &
842- + (dvel(1 , 1 ) + dvel(2 , 2 ))** 2_wp )
841+ + dvel(1 , 1 )** 2._wp + dvel(2 , 2 )** 2._wp &
842+ + (dvel(1 , 1 ) + dvel(2 , 2 ))** 2._wp )
843843
844844 end if
845845 end do
@@ -1257,9 +1257,9 @@ contains
12571257 jac_rhs(j, k, l) = alf_igr* (2._wp * (dvel(1 , 2 )* dvel(2 , 1 ) &
12581258 + dvel(1 , 3 )* dvel(3 , 1 ) &
12591259 + dvel(2 , 3 )* dvel(3 , 2 )) &
1260- + dvel(1 , 1 )** 2_wp + dvel(2 , 2 )** 2_wp &
1261- + dvel(3 , 3 )** 2_wp &
1262- + (dvel(1 , 1 ) + dvel(2 , 2 ) + dvel(3 , 3 ))** 2_wp )
1260+ + dvel(1 , 1 )** 2._wp + dvel(2 , 2 )** 2._wp &
1261+ + dvel(3 , 3 )** 2._wp &
1262+ + (dvel(1 , 1 ) + dvel(2 , 2 ) + dvel(3 , 3 ))** 2._wp )
12631263 end if
12641264 end do
12651265
@@ -1429,9 +1429,9 @@ contains
14291429 jac_rhs(j, k, l) = alf_igr* (2._wp * (dvel(1 , 2 )* dvel(2 , 1 ) &
14301430 + dvel(1 , 3 )* dvel(3 , 1 ) &
14311431 + dvel(2 , 3 )* dvel(3 , 2 )) &
1432- + dvel(1 , 1 )** 2_wp + dvel(2 , 2 )** 2_wp &
1433- + dvel(3 , 3 )** 2_wp &
1434- + (dvel(1 , 1 ) + dvel(2 , 2 ) + dvel(3 , 3 ))** 2_wp )
1432+ + dvel(1 , 1 )** 2._wp + dvel(2 , 2 )** 2._wp &
1433+ + dvel(3 , 3 )** 2._wp &
1434+ + (dvel(1 , 1 ) + dvel(2 , 2 ) + dvel(3 , 3 ))** 2._wp )
14351435 end if
14361436 end do
14371437 end do
0 commit comments