use default radiation

Author: navidcy

test/test_ocean_sea_ice_model.jl

@@ -69,16 +69,16 @@ using ClimaSeaIce.Rheologies
         #####
 
         # Adding a sea ice model to the coupled model
-        τua = Field{Face, Center, Nothing}(grid) 
+        τua = Field{Face, Center, Nothing}(grid)
         τva = Field{Center, Face, Nothing}(grid)
 
-        dynamics = SeaIceMomentumEquation(grid; 
+        dynamics = SeaIceMomentumEquation(grid;
                                           coriolis = ocean.model.coriolis,
                                           top_momentum_stress = (u=τua, v=τva),
                                           rheology = ElastoViscoPlasticRheology(),
                                           solver = SplitExplicitSolver(120))
 
-        sea_ice  = sea_ice_simulation(grid; dynamics, advection=WENO(order=7)) 
+        sea_ice  = sea_ice_simulation(grid; dynamics, advection=WENO(order=7))
         liquidus = sea_ice.model.ice_thermodynamics.phase_transitions.liquidus
 
         # Set the ocean temperature and salinity
@@ -110,11 +110,9 @@ Return a test-bed coupled simulation with a Checkpointer.
 function testbed_coupled_simulation(grid; stop_iteration=8)
     ocean = ocean_simulation(grid)
 
-    radiation = Radiation(arch)
-
     atmosphere = JRA55PrescribedAtmosphere(arch; backend=JRA55NetCDFBackend(4))
 
-    coupled_model = OceanSeaIceModel(ocean; atmosphere, radiation)
+    coupled_model = OceanSeaIceModel(ocean; atmosphere)
 
     simulation = Simulation(coupled_model; Δt=10, stop_iteration)