1- include ( " steady_state_velocity.jl " )
1+ import ClimaCore : Grids
22
33struct AtmosSimulation{TT, S1 <: AbstractString , S2 <: AbstractString , OW, OD}
44 job_id:: S1
@@ -53,10 +53,10 @@ future work:
5353
5454function AtmosSimulation {FT} (;
5555 model = AtmosModel (),
56- domain:: AtmosDomain = SphereDomain {FT} (),
5756 params:: Parameters.ClimaAtmosParameters = ClimaAtmosParameters (FT),
58- initial_condition:: ICs.InitialCondition = InitialConditions. DecayingProfile (),
5957 comms_ctx:: ClimaComms.AbstractCommsContext = ClimaComms. context (),
58+ grid:: Grids.AbstractGrid = SphereGrid (FT, params, comms_ctx),
59+ initial_condition:: ICs.InitialCondition = InitialConditions. DecayingProfile (),
6060 dt = 600 ,
6161 start_date = DateTime (2010 , 1 , 1 ),
6262 t_start = 0 ,
@@ -119,7 +119,7 @@ function AtmosSimulation{FT}(;
119119 )
120120 spaces = get_spaces_restart (Y)
121121 else
122- spaces = get_spaces (domain , params, comms_ctx)
122+ spaces = get_spaces (grid , params, comms_ctx)
123123 Y = ICs. atmos_state (
124124 initial_condition (params), model,
125125 spaces. center_space,
@@ -130,9 +130,11 @@ function AtmosSimulation{FT}(;
130130 )
131131 end
132132
133- steady_state_velocity = get_steady_state_velocity (
134- domain, initial_condition, params, Y; check_steady_state,
135- )
133+ # TODO : Deal with not having topography
134+ steady_state_velocity = nothing
135+ # get_steady_state_velocity(
136+ # params, Y, topography, initial_condition, mesh_warp_type,
137+ # )
136138
137139 p = build_cache (
138140 Y, model, params, surface_setup, dt, start_date, tracers,
0 commit comments