1173       subroutine CalcRCRConvCoef (stepn, numSrfs)  argument
1179       integer numSrfs, stepn      local
1182       if (stepn .eq. 0) then
1190       if (stepn .ge. 1) then
1191         RCRConvCoef(1,:) =-ValueListRCR(3,:)*exp(-dtRCR(:)*(stepn+alfi))
1193         RCRConvCoef(stepn+1,:) = ValueListRCR(1,:)*(1-alfi) 
1196         RCRConvCoef(stepn+2,:) = ValueListRCR(1,:)*alfi 
1200       if (stepn .ge. 2) then
1201         do j=2,stepn
1203      &        exp(-dtRCR(:)*(stepn + alfi + 2 - j))*
1209       RCRConvCoef(stepn+1,:)= RCRConvCoef(stepn+1,:)
1210      &                  - RCRConvCoef(stepn+2,:)*(1.0-alfi)/alfi 
1211       RCRConvCoef(stepn+2,:)= RCRConvCoef(stepn+2,:)/alfi 
1219       subroutine CalcHopRCR (timestepRCR, stepn, numSrfs)  argument
1225       integer numSrfs, stepn        local
1229       call RCRint(timestepRCR*(stepn + alfi),PdistCur)
1231      &     *exp(-dtRCR(1:numSrfs)*(stepn + alfi)) + PdistCur(1:numSrfs)