1064       subroutine UpdHistConv(y,nsrfIdList,numSrfs)  argument
1071       integer   nsrfIdList(0:MAXSURF), numSrfs  local
1079       call GetFlowQ(NewQ,y,nsrfIdList,numSrfs) !new flow at time n+1
1085             QHistImp(j,1:numSrfs)=QHistImp(j+1,1:numSrfs)
1087          QHistImp(ntimeptpT+1,1:numSrfs) = NewQ(1:numSrfs)
1115                write(816,*) (QHistImp(j,n),n=1, numSrfs)
1124                write(8166,*) (QHistImp(j,n),n=1, numSrfs)
1134          QHistRCR(lstep+1,1:numSrfs) = NewQ(1:numSrfs)
1146                   write(816,*) j, (QHistRCR(j,n),n=1,numSrfs) ! read from file of previous run
1149             write(816,*) lstep+1, (QHistRCR(lstep+1,n),n=1, numSrfs)
1160                   write(8166,*) (QHistRCR(j,n),n=1, numSrfs)
1173       subroutine CalcRCRConvCoef (stepn, numSrfs)  argument
1179       integer numSrfs, stepn      local
1219       subroutine CalcHopRCR (timestepRCR, stepn, numSrfs)  argument
1225       integer numSrfs, stepn        local
1230       HopRCR(1:numSrfs) = RCRic(1:numSrfs) 
1231      &     *exp(-dtRCR(1:numSrfs)*(stepn + alfi)) + PdistCur(1:numSrfs)
1237       subroutine calcRCRic(y,srfIdList,numSrfs)  argument
1243       integer   srfIdList(0:MAXSURF), numSrfs, irankCoupled  local
1254          call GetFlowQ(Qini,VelOnly,srfIdList,numSrfs) !get initial flow
1255          QHistRCR(1,1:numSrfs)=Qini(1:numSrfs) !initialize QHistRCR
1257          call integrScalar(Pini,POnly,srfIdList,numSrfs) !get initial pressure integral
1259          call integrScalar(CoupleArea,POnly,srfIdList,numSrfs) !get surf area
1260          Pini(1:numSrfs) = Pini(1:numSrfs)/CoupleArea(1:numSrfs)
1262          Qini(1:numSrfs)=QHistRCR(1,1:numSrfs)
1263          Pini(1:numSrfs)=zero    ! hack
1266       RCRic(1:numSrfs) = Pini(1:numSrfs) 
1267      &          - ValueListRCR(1,:)*Qini(1:numSrfs)-PdistIni(1:numSrfs)
1272       subroutine integrScalar(scalInt,scal,srfIdList,numSrfs)  argument
1279       integer   srfIdList(0:MAXSURF), numSrfs, irankCoupled, i, k  local
1284         if(numSrfs.gt.zero) then
1285           do k = 1,numSrfs