wilcrofter, Simon Byrne, CliMA lead programmer, discusses extending Julia's broadcast operator for best performance of Earth system models on GPU clusters. I imagine, or perhaps hope, this is another emerging design pattern for climate modeling.
https://www.youtube.com/watch?v=0c5oHS1865U&list=PLP8iPy9hna6T7PRe2sucSonFsrrH-oEZC&index=63&pp=iAQB