Make generic zeros function work on GPU

[glwagner]

Resolves #719.

[codecov]

Codecov Report

Merging #720 into master will decrease coverage by 3.16%.
The diff coverage is 100.00%.

@@            Coverage Diff             @@
##           master     #720      +/-   ##
==========================================
- Coverage   79.75%   76.58%   -3.17%     
==========================================
  Files         122      122              
  Lines        2489     2422      -67     
==========================================
- Hits         1985     1855     -130     
- Misses        504      567      +63     

Impacted Files	Coverage Δ
src/Fields/field.jl	71.42% <100.00%> (-8.34%)	⬇️
src/Coriolis/no_rotation.jl	0.00% <0.00%> (-100.00%)	⬇️
src/Forcing/Forcing.jl	50.00% <0.00%> (-50.00%)	⬇️
src/SurfaceWaves.jl	0.00% <0.00%> (-47.06%)	⬇️
src/AbstractOperations/AbstractOperations.jl	33.33% <0.00%> (-33.34%)	⬇️
src/AbstractOperations/interpolation_utils.jl	63.33% <0.00%> (-33.34%)	⬇️
src/Forcing/model_forcing.jl	66.66% <0.00%> (-33.34%)	⬇️
src/Buoyancy/Buoyancy.jl	61.11% <0.00%> (-23.10%)	⬇️
src/Architectures.jl	66.66% <0.00%> (-22.23%)	⬇️
src/TimeSteppers/TimeSteppers.jl	80.00% <0.00%> (-20.00%)	⬇️
... and 37 more

---

Continue to review full report at Codecov.

Legend - Click here to learn more
Δ = absolute <relative> (impact), ø = not affected, ? = missing data
Powered by Codecov. Last update 6502f07...e778d2d. Read the comment docs.