Блог

Я пытаюсь извлечь географический субрегион из глобальных файлов NetCDF, доступных через opendap. Требуемый регион простирается от 13 градусов западной долготы до 10 градусов восточной долготы: долготы от -13 до 10.

Обычно NCO может сделать это следующим образом:

 ncks -d lat,40.,70.-d lon,-10.,13. https://dataserver.nccs.nasa.gov/thredds/dodsC/CMIP5/NASA/GISS/rcp85/E2-H_rcp85_r2i1p1_day/tos_day_GISS-E2-H_rcp85_r2i1p1_20510101-20751231.nc out.nc
 

Однако долготы в результирующем файле NetCDF изменяются от 0 до 13 градусов, а не от -10 до 13. Это связано с тем, что долготы записываются как «градусы востока», поэтому NCO, похоже, не обрабатывает отрицательные значения. Кто-нибудь знает способ настройки вышеуказанного для принятия отрицательных результатов?

Примечание: CDO может обрабатывать отрицательные долготы. Однако в этом случае CDO не является подходящим решением, поскольку он а) должен сначала загрузить глобальный набор данных перед обрезкой, и б) не способен обрабатывать этот конкретный файл opendap.

NCO обрабатывает эти «обернутые координаты», как описано здесь . Пожалуйста, повторите попытку с вашими ограничивающими рамками в [0,360], например,

 ncks -d lat,40.,70. -d lon,347.,10. https://dataserver.nccs.nasa.gov/thredds/dodsC/CMIP5/NASA/GISS/rcp85/E2-H_rcp85_r2i1p1_day/tos_day_GISS-E2-H_rcp85_r2i1p1_20510101-20751231.nc out.nc