[ecco-support] [EXTERNAL] A Question about the ecco_v4_py tutorial

[Wang, Ou (US 329B)]

Hello,

Thank you Hong for the reply.

The discontinuity actually reflects the tiled alignment of the ECCO V4 model grid. The map you posted is the advective flux of potential temperature along the model’s native x-direction, as indicated by the map title. The model’s native x- and y-directions are not necessarily along the zonal and meridional directions. An example of the tile alignment of the ECCO V4 grid is available at https://github.com/ECCO-GROUP/ECCO-v4-Python-Tutorial/blob/master/Tutorials_as_Jupyter_Notebooks/ECCO_v4_Loading_the_ECCOv4_native_model_grid_parameters.ipynb (cell 5). The x-direction is from left to right, while the y-direction is from bottom to top. So from tiles 10-12 to tiles 0-2, the x-direction rotates.

Best,
Ou

From: ecco-support <ecco-support-bounces at mit.edu> on behalf of Li Pan <lipan at udel.edu>
Date: Wednesday, December 20, 2023 at 11:02 PM
To: ecco-support at mit.edu <ecco-support at mit.edu>
Subject: [EXTERNAL] [ecco-support] A Question about the ecco_v4_py tutorial
Dear ECCO team,

Thank you for developing the Python package ecco_v4_py!
I'm using your Python package tutorial to analyze the ECCO estimate Version 4. I find some regions in the global map of the advective flux of the potential temperature are not connected (For example, the eastern and western Atlantic Oceans are disconnected). Could you please provide me with some hints to address the problem?
The figure is as follows:
[cid:ii_lqefyp661]
The code I used is as follows:
import xarray as xr
import glob
from os.path import join
import ecco_v4_py as ecco
import matplotlib.pyplot as plt
## =================================================================
ECCO_dir = ('D:\eccov4r4\python\Downloads')
ecco_grid = xr.open_dataset(glob.glob(join(ECCO_dir,'*GEOMETRY*.nc'))[0])
vol = (ecco_grid.rA*ecco_grid.drF*ecco_grid.hFacC).transpose('tile','k','j','i')
# year_start = 1993
# year_end = 2016
ecco_vars_int = xr.open_mfdataset(join(ECCO_dir,'*tem_flux*','*.nc'),\
                                       data_vars='minimal',coords='minimal',compat='override')
curr_t_ind = 0
tmp = str(ecco_vars_int.time.values[curr_t_ind])
plt.figure(figsize=(15,5));
ecco.plot_proj_to_latlon_grid(ecco_grid.XC, ecco_grid.YC, ecco_vars_int['ADVx_TH'].isel(time=curr_t_ind,k=0), show_colorbar=True,
                              cmin=-1e7, cmax=1e7, cmap='RdBu_r', user_lon_0=-67, dx=0.2, dy=0.2)
plt.title(r'U advection of potential temperature at the sea surface [$^\circ$C s$^{-1}$] during ' +
          str(tmp)[0:4] +'/' + str(tmp)[5:7], fontsize=16)
plt.show()
#=======================================================================
The data file I used has been attached. It is downloaded from: https://search.earthdata.nasa.gov/search/granules?portal=podaac-cloud&p=C1991543740-POCLOUD&pg[0][v]=f&pg[0][gsk]=-start_date&q=ECCO_L4_OCEAN_3D_TEMPERATURE_FLUX_LLC0090GRID_MONTHLY&tl=1703109688!3!!&lat=-0.5625&zoom=0<https://urldefense.us/v3/__https:/search.earthdata.nasa.gov/search/granules?portal=podaac-cloud&p=C1991543740-POCLOUD&pg*0**Av*=f&pg*0**Agsk*=-start_date&q=ECCO_L4_OCEAN_3D_TEMPERATURE_FLUX_LLC0090GRID_MONTHLY&tl=1703109688!3!!&lat=-0.5625&zoom=0__;W11bXVtdW10!!PvBDto6Hs4WbVuu7!Jw9mrk9ZI9veqFEH5atI3IXjnXMp3eq2_mzoeiF0kdIzabsxJMnAt9KZZSmtzUa2pfBQmULN706upQpBm8o$>

Thanks for your time. I'm looking forward to your response.
[Image removed by sender.] OCEAN_3D_TEMPERATURE_FLUX_mon_mean_2017-12_ECCO...<https://urldefense.us/v3/__https:/drive.google.com/file/d/12JTIby5oVpS688QwXCfkmf01ofPWz9yw/view?usp=drive_web__;!!PvBDto6Hs4WbVuu7!Jw9mrk9ZI9veqFEH5atI3IXjnXMp3eq2_mzoeiF0kdIzabsxJMnAt9KZZSmtzUa2pfBQmULN706uuKCxn5Y$>[Image removed by sender.]

Best regards,
Li Pan

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mit.edu/pipermail/ecco-support/attachments/20231221/657b1075/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 249607 bytes
Desc: image.png
URL: <http://mailman.mit.edu/pipermail/ecco-support/attachments/20231221/657b1075/attachment-0001.png>