Hi, I was trying to retrieve the Channel velocity data by using <co

Kernel keep dying while using getCutout about pyjhtdb HOT 1 OPEN

idies commented on September 17, 2024 2

Kernel keep dying while using getCutout

from pyjhtdb.

Comments (1)

shervinsahba commented on September 17, 2024

I'm still getting familiar with the dataset, but the error may be from JHTDB's switch to 1-indexing from 0-indexing in the space and time coordinates. If you think that's the case with your query, try something like

import numpy as np
import pyJHTDB

lJHTDB = pyJHTDB.libJHTDB()
lJHTDB.initialize()
auth_token = "<your auth token>"   #<--- AUTH
lJHTDB.add_token(auth_token)

result = lJHTDB.getbigCutout(
        data_set='isotropic1024coarse', fields='u', 
        t_start=1, t_end=3, t_step=1,
        start=np.array([1, 1, 1], dtype = int),
        end=np.array([1024, 1024, 1], dtype = int),
        step=np.array([16, 16, 1], dtype = int),
        filter_width=4,
        filename='cutout')

lJHTDB.finalize()

If you change t_start or one of the spatial starts to 0, it will fail with error

"System.Web.Services.Protocols.SoapException: Server was unable to process request. ---> System.Exception: JHTDB started using 1-based indexing from Sept 16. Please update your JHTDB library if you have not done so.

If you're using Jupyter, the error may not print but it should be in the console.

from pyjhtdb.

Recommend Projects