Skip to content

Commit

Permalink
Add debugging info
Browse files Browse the repository at this point in the history
  • Loading branch information
dennisYatunin committed Jul 9, 2024
1 parent 3e30fd8 commit 23ca08e
Show file tree
Hide file tree
Showing 2 changed files with 8 additions and 3 deletions.
3 changes: 0 additions & 3 deletions src/analytic_solutions/analytic_solutions.jl
Original file line number Diff line number Diff line change
Expand Up @@ -43,9 +43,6 @@ function initial_thermo_state(params, coord)
T_min = FT(100)
z_iso = log((T_min / T_sfc - a) / (1 - a)) / β
p_iso = p_sfc * ((1 - a) / (1 - a * T_sfc / T_min))^(cp_d / R_d)
z > z_iso &&
@warn "Extending constant-N profile with isothermal profile above T = \
$T_min (z = $z_iso)" maxlog = 1

# The perturbation in approximate_FΔuvw is computed around a background
# state that is not hydrostatically balanced and is unphysical above 30 km,
Expand Down
8 changes: 8 additions & 0 deletions src/utils/common_spaces.jl
Original file line number Diff line number Diff line change
Expand Up @@ -120,6 +120,14 @@ function make_hybrid_spaces(
hypsography =
Hypsography.LinearAdaption(Geometry.ZPoint.(z_surface))
end
@show device
@show ClimaComms.device(Fields.coordinate_field(h_space))
@show ClimaComms.device(z_grid)
@show ClimaComms.device(z_surface)
@show ClimaComms.device(Geometry.ZPoint.(z_surface))
@show typeof(h_space)
@show typeof(z_grid)
@show typeof(z_surface)
end
grid = Grids.ExtrudedFiniteDifferenceGrid(h_grid, z_grid, hypsography; deep)
# TODO: return the grid
Expand Down

0 comments on commit 23ca08e

Please sign in to comment.