As of CHARMM version 32 the COOR HBONd command can also give the distance distribution, from which you should easily be able to compute the average. From corman.doc:
"Distance and lifetime histograms can be computed for all (putative) hydrogen
bonds encountered in the analysis"


Lennart Nilsson
Karolinska Institutet
Stockholm, Sweden