среда, 7 декабря 2016 г.

Восток - дело тонкое. Реанализ JRA55 download

Задача: скачать данные японского реанализа JRA55.

Краткое описание JRA55 доступно в виде html и pdf файла.

Инструкция

1) Просто так скачать JRA55 нельзя - нужен доступ (application).

На официальном сайте JRA55 нужно пройти регистрацию. ВНИМАНИЕ! Регистрация занимает около трёх дней, так что блицкриг "зарегистрировался и скачал" не пройдёт.

Проходим по ссылке из меню сайта - APPLICATION - ACCEPT.

Дальше заполняем форму и ждём, пока на указанный адрес электронной почты вам не пришлют логин и пароль. Пароль можно (и даже нужно) будет поменять.


2) После получения логин/пароль нужно ... скачать данные через ftp. =\

Подсказка лежит на той же странице сайта, где есть меню с APPLICATION. Она расположена чуть ниже под заголовком "Downloading JRA-55 data".

3) Для работы с ftp вам нужен либо выход к компьютеру с ОС Linux, либо соответствующая вашей ОС программа для работы с ftp. Для Windows это может быть FileZilla.

4) FileZilla

Скачать FileZilla на русском языке

Скачав актуальную версию FileZilla и установив её стандартным способом на ПК, открываем:

а) Файл -> Менеджер сайтов
б)  Хост: ds.data.jma.go.jp (см. пример bash-скрипта с сайта)
     Тип входа: Нормальный
     Логин: ваш логин из email
     Пароль: ваш пароль (из email или ваш, новый)
в) Нажимаем "Соединить"
г) В появившемся списке файлов проходим в JRA-55/Hist
д) Далее пути для суточных данных Daily, для месячных - Monthly
е) Далее нужно выбрать папку с пакетом параметров (surface, isobaric levels и т.д.). Описание папок и вложенных в файлы параметров можно узнать, например, тут.
ж) Суточные данные лежат в папках "годмесяц" (201206). Нужно зайти в такую месячную папку и найти по названию параметра (I), номеру дня (II) и часа (III) соответствующий файл. Например, anl_p125_ugrd.1974032718.
з) Наведя курсор мыши на файл нажимаем правую кнопку и из выпадающего меню выбираем "Скачать".
и) Можно выделить множество файлов и при правом клике выбрать "Добавить файлы в задание". При этом выбранные файлы будут добавлены в список и отображаться внизу. После этого можно нажать правой кнопкой мыши на этом окне и выбрать "Выполнить задание". Все файлы из группы будут загружены.
к) КУДА СКАЧИВАЮТСЯ файлы? Они скачиваются в открытую в левой части окна FileZilla (локальный сайт) папку. По умолчанию это C:/Users/username.

5) На Linux закачка производится либо через консоль, как показано в примере ниже, либо через bash-скрипт для скачивания большого количества данных.

$ ftp ds.data.jma.go.jp # log-on to the JMA Data Dissemination System (JDDS)
????????                # user ID
********                # password
binary                  # changes transfer mode to binary
cd JRA-55/Hist/Monthly/anl_surf125
get anl_surf125.195801
get anl_surf125.195802
...
bye

Вот пример из сети для работы с месячными данными JRA55. Для суточных (6-часовые) данных  предлагаю вам свой рабочий вариант. В нём нужно указать ваши логин/пароль (USER/
PASSWD).


# ===========================================================
HOST='ds.data.jma.go.jp'
USER='your_username'
PASSWD='your_password' 


# Название папки с пакетом параметров. Тут - параметры на изобарических уровнях
STR="anl_p125_hgt."

YEAR=1966   # начальный год


FINALYEAR=2016   # последний год (не включая его)
MONTH=1
DAY=1
HOUR=0

while [ $YEAR -lt $FINALYEAR ]; do

/usr/bin/ftp -inv $HOST<<EOF
user $USER $PASSWD
binary

cd JRA-55/Hist/Daily/anl_p125/$YEAR`printf "%02d" $MONTH`
get $STR$YEAR`printf "%02d" $MONTH``printf "%02d" $DAY``printf "%02d" $HOUR` 


bye
EOF

    let HOUR=HOUR+6
      if [ $HOUR = 24 ]; then
         HOUR=00
     let DAY=DAY+1
       if [ $DAY = 32 ]; then
         DAY=01
             let MONTH=MONTH+1
             if [ $MONTH = 13 ]; then
               MONTH=1
               let YEAR=YEAR+1
             fi
       fi 
      fi
done
exit

# ===========================================================

Сохраните этот код в файл, добавьте расширение ".sh", сделайте его исполняемым с помощью команды

chmod +x script_name.sh

Некоторые замечания:

  •  После скачивания файлы не имеют расширения "*.nc" 
  • Файлы хранят информацию о нескольких или одном параметра за один временной срок
  • Другие японские реанализы (например, JRA25) доступны с общего родительского сайта JRA 

Комментариев нет:

Отправить комментарий

Как перевести UV в направление и скорость ветра? How to convert wind UV-components to direction and velocity?

 Всё просто.  def uv2dir(u, v):     '''     Источник:     https://github.com/blaylockbk/Ute_WRF/blob/master/functions/wind_calc...