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

Make figure for example 1 #7

Closed
2 of 4 tasks
NoraLoose opened this issue Jan 15, 2021 · 7 comments
Closed
2 of 4 tasks

Make figure for example 1 #7

NoraLoose opened this issue Jan 15, 2021 · 7 comments
Assignees

Comments

@NoraLoose
Copy link
Member

NoraLoose commented Jan 15, 2021

Example 1 in #1:

POP 0.1 degree data, 5-day average, fixed-scale and fixed-factor filters (with different shapes: Gaussian, sharp, boxcar?) applied to salinity. Salinity does go near 0 in some places, so maybe we can use this to give an example of how the sharp filter can generate negative values from positive data. If we implement boxcar it will have to be ad-hoc because the Laplacian-based ones can't do it.

To summarize, highlight the following aspects with POP 0.1 degree data:

  • fixed-scale vs. fixed-factor filters
  • Gaussian vs. sharp filter shape
  • Sharp filter can generate negative values from positive data
  • filter which takes into account boundaries vs. filter which ignores boundaries
@NoraLoose
Copy link
Member Author

NoraLoose commented Jan 20, 2021

A simple figure that shows differences in Gaussian vs. sharp filter:

SST_gulfstream

We might want to substitute this figure by the model-obs figure (#16).

@NoraLoose
Copy link
Member Author

NoraLoose commented Jan 20, 2021

Boundary-aware vs. boundary-unaware filter. For a change, we are now looking at the Agulhas Current:

SST_Agulhas_boundaries_latlon

Notice the "cold" band/shadow hugging the coast, which arises from filling land values with 0. This figure shows how fill values on land artificially affect the filtered fields in oceanic coastal regions. Apart from the differences near the coast, the filtered fields are identical.

@iangrooms
Copy link
Member

"Sharp filter can generate negative values from positive data" we can actually demonstrate this using the alternate definition of EKE \bar{(u')^2}. I think Jake's along-track example might do this already; in any case it might be easier to demonstrate this using a "variance" than using salinity

@NoraLoose
Copy link
Member Author

"Sharp filter can generate negative values from positive data" we can actually demonstrate this using the alternate definition of EKE \bar{(u')^2}. I think Jake's along-track example might do this already; in any case it might be easier to demonstrate this using a "variance" than using salinity

Agreed. I was also planning to show the alternative definition \bar{(u')^2} for the POP data, and maybe compute a global integral, as you suggest here. Unless we don't want to do the same thing twice, because Jake has kind of done this for the obs data...

There is something odd about salinity, see cells 9 and 10 in this notebook.

@iangrooms
Copy link
Member

Yes, negative values of SSS is very suspicious. Might as well just do the alternative version of EKE instead. Jake's notebook shows local negatives using the alternative definition of EKE, but we could still use the alternative definition on POP data to demonstrate the global integrals so it's not redundant.

@NoraLoose
Copy link
Member Author

NoraLoose commented Mar 5, 2021

I updated the figure that illustrates the effect of boundary-unaware filters. I switched to a warmer region with maximal shoreline length to make the contamination by zero fill values on land even more clear:

SST_boundary_unaware

In the right panel, SST values close to the coast are as low as 6 degree C.

@iangrooms
Copy link
Member

This looks great! Thanks @NoraLoose. I think we can close this issue, and I'll update my post on #1 to reflect that this figure is done.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants