Доступ к данным с помощью openDap - PullRequest
0 голосов
/ 17 октября 2019

Я пытаюсь получить доступ к данным с веб-сайта через openDap с python. Для этого нужна библиотека pydap и open_url. Теперь основная проблема заключается в том, что когда я пытаюсь запустить скрипт для получения данных, которые имеют определенное значение глубины, я получаю следующую ошибку:

Error {
    code = 5;
    message = "Range last index less than first index";

Веб-сайт, с которого я пытаюсь получить информацию, являетсяна основе https://tds.hycom.org/thredds/dodsC/datasets/GLBv0.08/expt_93.0/data/forecasts/hycom_glbv_930_2019101512_t024_uv3z.nc.html.

Ниже приведен фрагмент моего кода:

for h in ['1', '4']:
    h=convert_depth(depth)
    depthstr=str(depth)
    hstr=str(h)

    h2=h
    h1=h-1

    dataset = open_url('http://tds.hycom.org/thredds/dodsC/datasets/GLBv0.08/expt_93.0/data/forecasts/hycom_glbv_930_'+sys.argv[1]+'12_t'+sys.argv[2]+'_uv3z.nc?time[0:1:0],depth['+hstr+':1:'+hstr+'],lat[0:1:3000],lon[0:1:4499],water_u[0:1:0]['+hstr+':1:'+hstr+'][0:1:3000][0:1:4499],water_v[0:1:0]['+hstr+':1:'+hstr+'][0:1:3000][0:1:4499]')
    dataset2 = open_url('http://tds.hycom.org/thredds/dodsC/datasets/GLBv0.08/expt_93.0/data/forecasts/hycom_glbv_930_'+sys.argv[1]+'12_t'+sys.argv[2]+'_ts3z.nc?time[0:1:0],depth['+hstr+':1:'+hstr+'],lat[0:1:3000],lon[0:1:4499],water_temp[0:1:0]['+hstr+':1:'+hstr+'][0:1:3000][0:1:4499]')

    u = dataset['water_u']
    v = dataset['water_v']
    tmp = dataset2['water_temp']

    uu=u.water_u[0,h1:h2]/float(1000)
    vv=v.water_v[0,h1:h2]/float(1000)
    tt=tmp.water_temp[0,h1:h2]/float(1000)+20

Как настроить этот код, чтобы я больше не получал это странное сообщение об ошибке и мог получить данные длявсе широты и долготы для конкретных глубин 1 и 4? Я использую Python 2.7 и Pydap == 3.2.2.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...