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