...
elem(ie)%derived%dpdiss_ave(i,j,k)=elem(ie)%derived%dpdiss_ave(i,j,k)+&
rhypervis_subcycle*eta_ave_w*(elem(ie)%state%dp3d(i,j,k,nt)-dp3d_ref(i,j,k,ie))
Update Nov 1, 2022: Please see additional code fix described in CAM github issue #678.