Add task for plotting maps of wind stress curl

[xylar]

The wind stress curl is first computed on the native grid (at vertices) and then is remapped to comparison grids.

The reconstruciton of the wind stress at edges is computed with a utility funciton vector_cell_to_edge_isotropic(), analogous to the subroutine with the same name in the MPAS framework.

Checklist

• User's Guide has been updated
• If this PR adds a new analysis task, it has also been added to the user's guide
• Developer's Guide has been updated
• API documentation in the Developer's Guide (api.rst) has any new or modified class, method and/or functions listed
• Documentation has been built locally and changes look as expected
• Testing comment in the PR documents testing used to verify the changes

[xylar]

Testing

I ran this on a SORRMr4 G-Case:
https://web.lcrc.anl.gov/public/e3sm/diagnostic_output/ac.xylar/analysis/FAnSSIE/wind-stress-curl/20250426.GMPAS-JRA1p5-DIB-PISMF-SIS.TL319_SOwISC12to30E3r4_ais4to20.coupled_ismf.chrysalis/clim_yrs5-11/ocean/index.html#antarctic_extended_wsc

It seems to be doing what it should.

[xylar]

@irenavankova, I still need to add the documentation but could you have a look when you can and see if this might be helpful to you, or if you think it will be useful in general even if not?

[irenavankova]

Thank you for implementing this. It works well for me. The only thing is that I think units on the figures should be N/m^3

[irenavankova]

Was your example done on a run with a finer atmospheric mesh? For the V2.1 There are much bigger squares: https://web.lcrc.anl.gov/public/e3sm/diagnostic_output/ac.vankova/analysis/FISMF/hist_curl/clim_2005-2014_ts_1951-2014/html/ocean/index.html#&gid=3&pid=1

[xylar]

Thank you for implementing this. It works well for me. The only thing is that I think units on the figures should be N/m^3

Oh, dear, that's embarrassing. Good call. Should be fixed now.

Was your example done on a run with a finer atmospheric mesh? For the V2.1 There are much bigger squares: https://web.lcrc.anl.gov/public/e3sm/diagnostic_output/ac.vankova/analysis/FISMF/hist_curl/clim_2005-2014_ts_1951-2014/html/ocean/index.html#&gid=3&pid=1

That's probably the difference between a G-Case and a B-Case. I think we're making a really poor choice in how we interpolate wind stress in E3SM in B-Cases but I got the impression we didn't have a choice about it.

[xylar]

Test suite

I ran the developer's test suite:
https://web.lcrc.anl.gov/public/e3sm/diagnostic_output/ac.xylar/analysis_testing/chrysalis/add-wind-stress-curl/
The wind stress curl doesn't look good at 240 km resolution (no surprise):

but it does show up as global, Arctic and Antarctic maps in analysis using --polar_regions and only as a global plot in a water-cycle analysis configuration, as expected.