MITgcm / xmitgcm

Read MITgcm mds binary files into xarray
http://xmitgcm.readthedocs.io
MIT License
56 stars 65 forks source link

Assertion Error when .load() or .plot() #291

Closed LizzWebb closed 2 years ago

LizzWebb commented 2 years ago

Hi all,

I am trying to load and/or plot data from the llc4320 model; I have a script that has done this many times in the past but it recently stopped working, and now gives the message: Assertion Error:

The lines I used are below:

model = llcreader.ECCOPortalLLC4320Model()
ds_Ovel  = model.get_dataset(varnames=['U','V'], k_levels=[0])
OceUvel = ds_Ovel.U.isel(time=0,face=5, k=0)
OceUvel.load()

Does anyone know what is causing this, or if there is an update I missed? Thank you in advance!

rabernat commented 2 years ago

Hi @LizzWebb! Sorry for the friction. Could you post the full stack trace from your error? (Like all of the lines that python spits out above AssertionError.) That will help us narrow down the problem.

rabernat commented 2 years ago

FYI, as of right now, it looks like the ECCO data portal is down.

https://data.nas.nasa.gov/ecco/data.php?dir=/eccodata/llc_4320/compressed/0000010368

image

I suspect this is the root of the problem.

LizzWebb commented 2 years ago

Ah, my mistake. The ecco portal previously posted on their homepage when they were down. I didn't think to check the individual data sets.

Thanks, @rabernat. I will wait for the data set to come back, and post the full error if it is not resolved.