INTEGER expMatrixCounter, impMatrixCounter, expMatrixWriteCount,
COMMON /MATRIX_PARAMS/ expMatrixCounter, impMatrixCounter,
impMatrixCounter=0
impMatrixCounter = impMatrixCounter + 1
& impMatrixWriteCount, impMatrixCounter
recipImpMatrixCounter = (1. _d 0)/dble(impMatrixCounter)