PFEflux_l = (PFEflux_u+(Fe_spm(i,j,k)+Fe_ads_inorg(i,j,k)
Fe_burial(i,j) = PFEflux_l
PFEflux_l = 0. _d 0
& *hFacC(i,j,k,bi,bj)-PFEflux_l)*recip_drF(k)
PFEflux_u = PFEflux_l
_RL PFEflux_l