Real*8 avgField
field(I) = field(I) + avgField
O avgField,
avgField = 0.
avgField = avgField+field(I)
avgField = avgField/FLOAT(nValid)
field(I) = field(I)-avgField