Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Tweak for sea-ice metric demo #1030

Merged
merged 10 commits into from
Jan 29, 2024
Prev Previous commit
address comments from @acordonez
  • Loading branch information
lee1043 committed Jan 29, 2024
commit c16b3fb2663fb7288fd26e469fc15e695d61b7e9
8 changes: 4 additions & 4 deletions doc/jupyter/Demo/sea_ice_line_plots.py
Original file line number Diff line number Diff line change
Expand Up @@ -14,12 +14,12 @@

arctic = (
(ds.siconc.where(ds.lat > 0) * 1e-2 * area.areacello * 1e-6)
.where(ds.siconc > 0.15)
.where(ds.siconc > 15)
.sum(("lat", "lon"))
)
"""
Note for the above line
- where siconc > 0.15: to consider sea ice extent instead of total sea ice area
- where siconc > 15: to consider sea ice extent instead of total sea ice area (criteria: 15%)
- multiply 1e-2: to convert percentage (%) to fraction
"""

Expand All @@ -29,12 +29,12 @@
obs = xc.open_dataset(obs_file)
obs_area = 625
obs_arctic = (
obs.ice_conc.where(obs.lat > 0).where(obs.ice_conc > 0.15) * 1e-2 * obs_area
obs.ice_conc.where(obs.lat > 0).where(obs.ice_conc > 15) * 1e-2 * obs_area
).sum(("xc", "yc"))
"""
Note for the above lines
- obs_area = 625 # area size represented by each grid (625 km^2 = 25 km x 25 km resolution)
- where siconc > 0.15: to consider sea ice extent instead of total sea ice area
- where siconc > 15: to consider sea ice extent instead of total sea ice area (criteria: 15%)
- multiply 1e-2: to convert percentage (%) to fraction
"""

Expand Down
Loading