MPAS-Analysis icon indicating copy to clipboard operation
MPAS-Analysis copied to clipboard

Add OSNAP Arctic Transects transport time series

Open anirban89 opened this issue 1 year ago • 5 comments

Add Arctic Transects from branch https://github.com/milenaveneziani/geometric_features/tree/flipOSNAP_renameArcticTransectsTag of geometric features. Separate all arctic transects from standard transects using transectGroup flag in config option timeSeriesTransport. When running mpas_analysis without the flag --polar_regions only standard transects are plotted. All Arctic transects including the three new ones "Hudson Bay-Labrador Sea", "OSNAP section East" and "OSNAP section West" are only plotted when --polar_regions is specified. timeSeriesTransport split into transportGroups "Transport Transects" and "Arctic Transport Transects". By default all Arctic transects are removed from standard transects.

User can specify the list of standard transects to plot with transectsToPlot under config option (default.cfg) timeSeriesTransportTransects and leave it as all under timeSeriesArcticTransportTransects. List can not include an Arctic Transect. Similarly user can also specify the list of arctic transects to plot with transectsToPlot under config option timeSeriesArcticTransportTransects but to plot these mpas_analysis needs to be run with the --polar_regions flag.

Files changed: default.cfg, polar_regions.cfg and ocean/time_series_transport.py @xylar @milenaveneziani @alicebarthel

anirban89 avatar Sep 27 '22 20:09 anirban89

@anirban89, I think this still has the default title from the name of the branch. Could you change it to be a little more informative?

xylar avatar Sep 28 '22 21:09 xylar

This will need https://github.com/MPAS-Dev/geometric_features/pull/185 to be merged, a new release of geometric_features, and MPAS-Analysis to be updated to that new version. This isn't as big a deal as it sounds but I'd like to hold off on that process so we can do a bit more testing of https://github.com/MPAS-Dev/geometric_features/pull/185.

xylar avatar Sep 28 '22 21:09 xylar

@xylar yes I agree this PR does not need to be merged right now

anirban89 avatar Sep 28 '22 22:09 anirban89

@alicebarthel and @milenaveneziani, we want to add obs. to this. What is the right upper and lower value we should use? Min and max in time?

xylar avatar Oct 10 '22 18:10 xylar

@xylar: sorry, I never replied to this. I usually compute a standard deviation from all the obs I have (or get it from previous literature) and then plot that as +-sigma around the mean value. @alicebarthel: what do you think?

milenaveneziani avatar Dec 08 '22 15:12 milenaveneziani