2014-10-18 9 views
0

Я загрузил данные из погодного модуля american gfs (от ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.2014101812/master, дата редактирования, если вы читаете более 1 месяца спустя). Это grib-файлы.Градусы, отображающие данные из разных файлов (разные моменты времени) в одном 2D-графике (высота - время)

Я использую grads для построения данных, и он отлично работает, если я накладываю данные из 1 timestep (= 1 файл) на карту. Итак, моя первая попытка состояла в том, чтобы открыть первые 2 файла и сделать затененный участок, показывающий области с одинаковой скоростью ветра в том же цвете. Я получаю пустой сюжет с правильными осями.

'open gfs.t12z.mastergrb2f00.ctl' 

'open gfs.t12z.mastergrb2f03.ctl' 

'set dfile 1' 

'set time 12Z17OCT2014' 
'set gxout shaded' 
'set lat 50' 
'set lon 0' 
'set lev 1000 200' 
'set t 1 2' 
'd sqrt(ugrdprs*ugrdprs+vgrdprs*vgrdprs)*1.9438' 

'set dfile 2' 
'set time 15Z17OCT2014' 
'set t 1 2' 
'd sqrt(ugrdprs*ugrdprs+vgrdprs*vgrdprs)*1.9438' 

Я использую opengrads - Сетка анализа и системы отображения (GrADS) Версия 2.0.2.oga.2.

Любой, кто может помочь мне распечатать данные из разных временных меток (= разных файлов) в один сюжет?

ответ

1

Вы можете использовать один .ctl-файл для построения нескольких файлов данных. Я сделал это с .gdat файлами, поэтому, я думаю, вы сможете сделать это с помощью GRIB. Вы должны создать новый .ctl-файл, используя шаблоны grads.

http://www.iges.org/grads/gadoc/templates.html

В принципе, в DSET линии вы поставите имя файла, похожий на

dset gfs.t12z.mastergrb2f%h2.grb 

где% h2 представляет 2 цифры час.