search

willaguiar / ASC_and_heat_transport

Github repository for Analysis of ASC speed and cross slope heat transport on Panan simulation
0 stars 0 forks source link

New visualisation of cross-slope properties in Antarctic margins? #38

Open taimoorsohail opened 2 months ago

taimoorsohail commented 2 months ago

I have written a piece of example code that plots ocean properties binned by isobath, as shown below. I have remapped the isobath bins (which are the same as st_ocean) by the cumulative area held by each bin, normalised. This produces a 'normalised distance from Antarctica', creating a more realistic slope representation. We can also multiply the normalised cumulative area by latitude to create a "pseudo-Latitude". Anyway, I was wondering what we all think about this? We can mask the ocean properties by regime before binning, which would give us a clear representation of the average T, S and rho in each regime, without having to resort to using individual sections. Let me know!

Avg_thermohaline_props

adele-morrison commented 1 month ago

Love it! Would be great to apply this to u_along, temperature and density for temporal composites of months with strong or weak ASC.

taimoorsohail commented 1 month ago

Yep, it is easy enough to do with scalar tracer properties like temperature, and density, but u_along is harder. We would have to rotate the velocity vectors to be along- and across-slope for each isobath (the same isobaths as there are st_ocean values). I'm guessing this is an expensive code, considering we only have some selected isobaths for which u_along is available? The same applies for CSHT, by the way, with the added complexity of the zonal convergence term. I'm not sure how to handle u_along and CSHT, if there is any code that does this for multiple isobaths that would be great to see and incorporate.

adele-morrison commented 1 month ago

I think for ACC speed we could just use u in these zonally averaged transect plots. It won't give values that match u_along, but it should still show when / where the ACC is strong/weak.

And maybe we can use CSHT to define the temporal composites (i.e. all the times when CSHT is southward), and then just show a schematic arrow for CSHT, but not actual data.

taimoorsohail commented 1 month ago

Cool! I'll try this! By the way, the xwmb package seems promising in quickly assessing transport across and along any given arbitrary section - might be faster than the methods we are currently using to calculate u_along and CSHT+ZC.

adele-morrison commented 1 month ago

Oh, good to know about xwmb! I'm looking forward to Henri's upcoming COSIMA talk about it.

adele-morrison commented 1 month ago

@taimoorsohail did you upload the code for this? I was going to make a new Figure 4 version using this method now, but am struggling to find the code.

taimoorsohail commented 4 weeks ago

Sorry @adele-morrison, I neglected to upload the relevant file before going on leave. I've uploaded it to my branch (taimoor) under the name "Circumpolar_average_example.ipynb". Let me know if you need any more help.