When I try to calculate the difference of velocity of ECCOv4 (LLC90 data) using the following python codes: div_uv = (grid_llc.diff(UVELMASS.mean(dim=‘time’) * data_uv.dyG, ‘X’, boundary=‘fill’) + grid_llc.diff(VVELMASS.mean(dim=‘time’) * data_uv.dxG, ‘Y’, boundary=‘fill’))/(data_uv.rA) I get the…

You might get more help if you give us a way to reproduce this? E.g. post a notebook here that has access to the output you use?

[image] navidcy: E.g. post a notebook here that has access to the output you use? I’ve created a topic with information on how to share notebook code [image] How-to: Share your jupyter notebook Technical About How to share your notebook code so that oth…

##load import numpy as np import xarray as xr import pandas as pd import os from xmitgcm import open_mdsdataset import xgcm #ecco_v4_py import sys from os.path import join,expanduser user_ECCO_dir = ('/g/data/p10/rh1363/ECCO/ECCOv4-py/') sys.path.append(join(user_ECCO_dir,'ECCOv4-py')) import ecco_v…

Thank you so much. Now I add the codes sucecessfully.

Thank you for you advice. Now I add the codes here.

The data can be download from here calculate vorticity using diff_2d_vector in LLC90 grid · Issue #599 · xgcm/xgcm · GitHub

Seems like access to p10 is required, right?

No you don’t need to, just because I put the data in p10. You can download from here ( calculate vorticity using diff_2d_vector in LLC90 grid · Issue #599 · xgcm/xgcm · GitHub )

[image] ruhui: You can download from here Ruhui, when you’re asking for help it is best to make it as easy as possible for people to assist you. One way you can do that is to put the data you’re using somewhere others can access it. The /scratch/public/ area on gadi is for just this purpose. …

Thank you so much. Now I put the test data on /scratch/public/rh1363/ Anyone can use the follow code to test it without any change. ##load import numpy as np import xarray as xr import pandas as pd import os from xmitgcm import open_mdsdataset import xgcm #ecco_v4_py import sys from os.path import…

How to avoid anomalous strip when executing the difference in LLC90 data, like ECCOv4

COSIMA Ocean

Aidan (Aidan Heerdegen, ACCESS-NRI Release Team Lead) 24 April 2023 01:01 3

I’ve created a topic with information on how to share notebook code

Topic		Replies	Views
AUS2200 vegetation fraction ancil creation issues AUS2200	23	406	19 November 2024
How to convert LLC (.data and .) output file to regular lat-lon grid data COSIMA	3	281	22 February 2023
Ancils for small test domain fail checking during reconfiguration due to longitude rounding errors Regional Nesting Suite	17	96	17 February 2026
Wind components on different grids AUS2200	4	100	13 December 2024
Issues loading ACCESS-OM2-01 daily data from RYF run Technical	6	334	20 May 2023

How to avoid anomalous strip when executing the difference in LLC90 data, like ECCOv4

Related topics