diff --git a/src/atchem2.f90 b/src/atchem2.f90
index 22af31b4b..594fb125c 100644
--- a/src/atchem2.f90
+++ b/src/atchem2.f90
@@ -351,6 +351,7 @@ end subroutine FCVFUN
   else if ( solverType == 3 ) then
     call FSUNDENSEMATINIT( 1, neq,  neq, ier )
     ! call FCVDENSE( neq, ier )
+    call FSUNDENSELINSOLINIT( 1, ier )
     call FCVDLSINIT( ier )
     ! UNEXPECTED SOLVER TYPE
   else