Конвертировать календарь с помощью ириса - PullRequest
Новая
       25

Конвертировать календарь с помощью ириса

0 голосов
/ 08 октября 2018

Как я могу конвертировать календарь с Iris?Например, преобразование «григорианского» календаря в «365_day»?

В данный момент я делаю 

original_tcoord = cube_in.coords()[tcoord_idx]

tmp_time = cf_units.num2date(
    original_tcoord.points, 
    original_tcoord.units.name, 
    target_calendar,
)

new_time = cf_units.date2num(
    tmp_time, 
    original_tcoord.units.name, 
    target_calendar,
)
new_unit = cf_units.Unit(
    original_tcoord.units.name,
    calendar=target_calendar
)
new_tcoord = iris.coords.DimCoord(
    new_time,
    standard_name=original_tcoord.standard_name, 
    long_name=original_tcoord.long_name, 
    var_name=original_tcoord.var_name,
    units=new_unit,
)

cube_adjusted = cube_in.copy()
cube_adjusted.remove_coord('time')
cube_adjusted.add_dim_coord(
    new_tcoord, 
    tcoord_idx
)

Конечно, это не самый быстрый способ?

1 Ответ

0 голосов
/ 19 ноября 2018

Вероятно, это можно сделать более просто, просто заменив единицы измерения на вашей существующей координате: 

tcoord = cube.coord('time')
tcoord.units = cf_units.Unit(tcoord.units.origin, calendar='gregorian')

Это предложение исходит от Рут Комер, а не от меня, но я публикую здесь, как и предлагалось.

...