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