Skip to content

Commit

Permalink
address comments from @acordonez
Browse files Browse the repository at this point in the history
  • Loading branch information
lee1043 committed Jan 29, 2024
1 parent 4fa5bcd commit c16b3fb
Showing 1 changed file with 4 additions and 4 deletions.
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

0 comments on commit c16b3fb

Please sign in to comment.