Permalink
Browse files

Fixed use of undeclared variables. Added check of return code from zg…

…ttrf
  • Loading branch information...
John Donners
John Donners committed Oct 19, 2016
1 parent ecb2df3 commit d5cf6d2c379ebd1c88dd26ecd6c9a1b2973dd280
Showing with 10 additions and 4 deletions.
  1. +4 −4 CreateInitialConditions.F90
  2. +6 −0 SolvePressureCorrection.F90
@@ -45,10 +45,10 @@ subroutine CreateInitialConditions
enddo
!assign the boundary conditions at k=1 and k=nx
do ic=xstart(3),xend(3)
do jc=xstart(2),xend(2)
temp(1 ,jc,ic) = tempbp(jc,ic)
temp(nx,jc,ic) = temptp(jc,ic)
do i=xstart(3),xend(3)
do j=xstart(2),xend(2)
temp(1 ,j,i) = tempbp(j,i)
temp(nx,j,i) = temptp(j,i)
end do
end do
@@ -127,6 +127,12 @@ subroutine SolvePressureCorrection
call zgttrf(nxm, amph(2), acphT, apph(1), appph, phpiv, info)
if (info.gt.0) then
print*,'Singular value found in LAPACK routine zgttrf: info=',info
print*,'Please try to adjust either NX or STR3 in bou.in'
call MPI_Abort(MPI_COMM_WORLD,1,ierr)
endif
call zgttrs('N',nxm,1,amph(2),acphT,apph(1),appph,phpiv, &
dphc(1,j,i), nxm, info)

0 comments on commit d5cf6d2

Please sign in to comment.