用JRA-55的土壤湿度转WRF模式

金玉姻缘

C:\Users\mayan\Desktop\IMG_20160611_125847.jpg
C:\Users\mayan\Desktop\IMG_20160611_125750.jpg
C:\Users\mayan\Desktop\IMG_20160611_125524.jpg

金玉姻缘

kazami2003 发表于 2016-6-10 14:01
其實, 不是讀不到, 而是....JRA55的grib data裡面的level code有問題. SoilW ---> 225, SoilT ---> 211 (  ...

师兄好 ，这是编写的vtable以及错误
GRIB1 Param    LEVEL1 Type    From level1   To level2   metgrid name   metgrid units    metgrid description
225                       213                0                                 SM                 Proportion         Soil wetness


GRIB2 discp     GRIB2 catgy       GRIB2 param      GRIB2 level
4                          0                    225                     0

错误
#############################
Inventory  for date =
Pres    SM
------------------------------------------------------------
Subroutine DATTINT:Interpolating 3-d files to fill in ang missing ...........
ERROR :Data not found :

kazami2003

金玉姻缘 发表于 2016-6-11 13:09
师兄好 ，这是编写的vtable以及错误
GRIB1 Param    LEVEL1 Type    From level1   To level2   metgr ...

No, 我說的是"必須修改" grib file…改Vtable沒用, 因為依照WMO對grib1的定義, level code 213, 221是"不存在"的, 必須把grib file裡面的level code改成112, 這樣ungrib才讀得到

kazami2003

如果真的非得用JRA55的資料, 而ungrib又讀不到的話, 還有一個方式, 就是把grib1 --> WRF IM file

http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Basics/IM_files/

不過這就要自己coding了.....上面連結自己看一下囉。

金玉姻缘

kazami2003 发表于 2016-6-13 23:49
如果真的非得用JRA55的資料, 而ungrib又讀不到的話, 還有一個方式, 就是把grib1 --> WRF IM file

http:/ ...

再次感谢老师！谢谢你的帮助，我学习学习！

Bozhero

JRA的土壤湿度2015年之后确实分为三层了，但是是交错的，感觉不好插值到NOAH 啊