6       subroutine itrSetup ( y,  acold )   argument
10       real*8     y(nshg,ndof),  acold(nshg,ndof)  local
32      &           y=y-(one-alfi)*Delt(1)*CFLfl(1)*acold
56       subroutine itrPredict ( yold,          acold,  argument
61       real*8        yold(nshg,ndof),            acold(nshg,ndof),  local
72         ac = acold*(one-almi/gami)
79         y  = yold+alfi/Dtgl*acold*(one-gami)
85         ac = acold*(one-almi)
92         y  = yold+alfi/Dtgl*acold
98         ac = acold
114         ac = acold*fct2+(y-yold)*fct3
128       subroutine itrCorrect ( y, ac,yold, acold,Dy)  argument
133       real*8   yold(nshg,ndof), acold(nshg,ndof)  local
146       ac= acold*fct1+(y-yold)*fct2
155       subroutine itrCorrectSclr (y, ac, yold, acold, Dyt)  argument
160       real*8   yold(nshg,ndof), acold(nshg,ndof)  local
167       ac(:,is)= acold(:,is)*fct1+(y(:,is)-yold(:,is))*fct2 
176       subroutine itrUpdate( yold,          acold,  argument
181       real*8        yold(nshg,ndof),            acold(nshg,ndof),  local
193          acold(:,6)= acold(:,6)*(1-1/gami) 
197          acold(:,1:nflow) = acold(:,1:nflow) 
198      &                    + (ac(:,1:nflow)-acold(:,1:nflow))*fct2
202          acold(:,7)= zero  
206          acold = acold + (ac-acold)*fct2