Allocate IOB arrays for nitrate and algae fluxes when -DCSIRO_BGC

Author: dougiesquire

src/access/shared/gtracer_flux.F90

@@ -89,6 +89,11 @@ subroutine flux_exchange_init (Time, Ocean, Ocean_state, Ice_ocean_boundary, atm
 #if defined(ACCESS_CM)
     allocate(Ice_ocean_boundary%co2(isc:iec,jsc:jec))
 #endif
+#if defined(ACCESS_OM) && defined(CSIRO_BGC)
+    allocate( &
+      Ice_ocean_boundary%iof_nit(isc:iec,jsc:jec), &
+      Ice_ocean_boundary%iof_alg(isc:iec,jsc:jec))
+#endif
 
     Ice_ocean_boundary%u_flux = 0.0
     Ice_ocean_boundary%v_flux = 0.0
@@ -115,6 +120,10 @@ subroutine flux_exchange_init (Time, Ocean, Ocean_state, Ice_ocean_boundary, atm
 #if defined(ACCESS_CM)
     Ice_ocean_boundary%co2 = 0.0
 #endif
+#if defined(ACCESS_OM) && defined(CSIRO_BGC)
+    Ice_ocean_boundary%iof_nit = 0.0
+    Ice_ocean_boundary%iof_alg = 0.0
+#endif
 
     ! Spawn 2D Ocean%fields FMS coupler type from 1D gas_fields_ocn
     call coupler_type_spawn(gas_fields_ocn, Ocean%fields, (/isc,isc,iec,iec/), &