注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

美丽心情

相信自己!&&Happy everyday! &LOVE

 
 
 

日志

 
 

lats4d写GRIB,新问题_done  

2010-10-17 11:05:59|  分类: WRF有关 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
GRIB (GRIdded Binary) data
这几天,因为有HadISST的nc数据,要转成GRIB格式导入WPS中,所以钻研了几天的lats4d。
学了最基本的转换命令,lats4d -i hadisst -o ss_tgrib -format grads_grib.这是第一步,写成后发现WPS中运行./ungrib.exe时出错,找不到该时次的数据……通过上网查原因,原来是monthly数据的问题,于是又费劲周折的把nc写成dat再用FORTRAN插值成daily的数据(因为我不会用GrADS插值),最后再把写成的某几个daily.dat用gradsnc(lats4d)转成GRIB。还是,到了./ungrib.exe这一步又出错了,同样的错误……我无语
./ungrib.exe之后的错误提示:
Subroutine DATINT: Interpolating 3-d files to fill in any missing data...
Looking for data at time 2007-12-01_00
ERROR: Data not found: 2007-12-01_00:00:00.0000
FORTRAN STOP
个人觉得还是写的GRIB数据有问题,因为同样的我用rtg_sst_grb_0.5.20050625的数据就可以实现。到底问题出在哪里呢?
===========================================================================
[wxj@localhost WPS]$ ./util/g1print.exe ../DATA/SST/rtg_sst_grb_0.5.20050625
Copen: File = ../DATA/SST/rtg_sst_grb_0.5.20050625                                                                                   
Fortran Unit = 0
UNIX File descriptor: 3
----------------------------------------------------
 rec GRIB GRIB  Lvl  Lvl  Lvl         Time      Fcst
 Num Code name  Code one  two                   hour
----------------------------------------------------
   1  11          1    0    0  2005-06-25_00:00 +00
***** End-Of-File on C unit   3
==========================================================================以上是查看rtg_st。
[wxj@localhost WPS]$ ./util/g1print.exe GRIBFILE.AAA
Copen: File = GRIBFILE.AAA                                                                                                           
Fortran Unit = 0
UNIX File descriptor: 3
----------------------------------------------------
 rec GRIB GRIB  Lvl  Lvl  Lvl         Time      Fcst
 Num Code name  Code one  two                   hour
----------------------------------------------------
   1   1          1    0    0  2007-12-01_00:00 +00
***** End-Of-File on C unit   3
===============================================这是查看我写的 sst_daily_grib_20071201.grb
只是GRIB Code不一样。
另外:Vtable
GRIB1| Level| From |  To     | metgrid  | metgrid | metgrid                              |GRIB2|GRIB2|GRIB2|GRIB2|
Param| Type |Level1|Level2| Name     | Units    | Description                         |Discp|Catgy|Param|Level|
-------+--------+-------+---------+-----------+----------+-----------------------------------+-------+------+--------+------+
  11 |          1   |   0   |            | SST      |     K        | Sea Surface Temperature   |  0     |     0 |      0  |       1 |
-----+---------+--------+---------+-----------+----------+--------------------------------------------+------------------------------------+
写GRIB数据的时候自动生成的vtable则是:
   name | id | title | units   | datatype | surface | decimal_scale_factor | precision | comments_1 | comments_2
   sst     | 1  | Daily 1 degree resolution SST [K]  |    | float | sfc | -999 | 16  |   |
 其中:#name = variable name (no blanks)
          # id = GRIB parameter number (>127 => AMIP-2 specific)
         # title = long name (description)
         # units = variable units
        # datatype = float or int
         # level_type = level_type in vertical dimension table, or blank if values must be defined via lats_vert_dim
         # decimal_scale_factor = GRIB decimal scale factor, or -999 if no decimal scaling
        # precision = number of bits of precision if stored in GRIB,
        #             or -999 for level-dependent bit length (ignored if decimal_scale_factor is set)
# comments_1 = comments, ignored by LATS
# comments_2 = comments, ignored by LATS
是不是GRIB parameter number 和GRIB Code (Gridded Information in Binary)的原因啊?他们到底什么意思?

GRIB codes:   http://badc.nerc.ac.uk/help/formats/grib/

Some of the information in the GRIB record information sections are represented by integer codes. For instance, the name of the parameter in the GRIB record is indicated using an integer known as the GRIB parameter code, or parameter identifier. Code tables, which are not part of the GRIB file but are held separately for each dataset using data in GRIB format.

我又用wgrib看了看两个数据:

[wxj@localhost SST]$ wgrib sst_daily_grib_20071201.grb
Undefined parameter table (center 100-2 table 128), using NCEP-opn
1:0:d=07120100:PRES:kpds5=1:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0
[wxj@localhost SST]$ wgrib rtg_sst_grb_0.5.20050511
1:0:d=05051100:TMP:kpds5=11:kpds6=1:kpds7=0:TR=1:P1=0:P2=0:TimeU=2:sfc:anl:NAve=0

还是对这些东西不了解,对比下来我写的数据就是与可用的不同。关键是问题出在哪里,怎么改啊,要在lats4d中添加什么参数,或者还是自己写个vtable?我仿照WPS中的Vtable写了个,gradsnc中不支持,老是出现错误提示……他不认识我写的,唉!
我再想想~~~~
Who can help me?

http://forum.wrfforum.com/viewtopic.php?f=22&t=233

http://bbs.lasg.ac.cn/bbs/thread-34060-2-37.html
写table.




  评论这张
 
阅读(696)| 评论(1)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017