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

美丽心情

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

 
 
 

日志

 
 

lats4d---->GRIB  

2010-10-18 09:40:24|  分类: WRF有关 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
因为我用的是HadISST.nc月平均的数据,而我需要的是daily的GRIB格式的SST数据。下面是我的解决思路,暂时来看没有什么问题。(到./ungrib.exe一步)
1,nc数据用GrADS写成dat  ----->nc_dat.gs
2,dat用FORTRAN将月平均资料进行线性插值,插值成每天的,并分daily输出n个dat文件。----->sst_daily.f90
3,then,写好ctl描述文件,用gradsnc中的lats4d将所有的dat文件转换成GRIB数据。此期间,有很多重要的问题
(1) 命令:lats4d -v -i sst_daily.ctl -o sst_daily_YYYYMMDD -format grads_grib -ftype ctl -table my.table
(2) my.table写法,起初不清楚怎么写,让lats4自然形成了一个.grads.table之后按照其中的格式写了一个,table主要包括三项:#!variable,#!vert,#!center,根据我搜索的所有资料以及在运行lats4d中显示的错误中改正,最终------->my.table
#!variable
sst     | 11  | ** Daily 1 degree resolution SST [K]     |      | float | sfc | -999 | 16 |   |
#!vert
plev     | Pressure level         | hPa    | multi  | down |  100 | 0 |  0 | 0
sfc      | Earth surface             |          | single |   up |    1 | 0 |  0 | 0
#!center
gsfc      | 13  | 100 | 2
ugamp     | 27  | 100 | 2
#!qc
至少目前我觉得没有错误。
生成的sst_daily_YYYYMMDD.grb在./ungrib.exe之后显示:(只举一个时次例子)
 *** Starting program ungrib.exe ***
Start_date =  2007-12-04_00:00:00 ,      End_date = 2007-12-05_00:00:00
output format is WPS
Path to intermediate files is ./
 ungrib - grib edition num            1
###############################################################################
Inventory for date = 2007-12-04 00:00:00
PRES   SST    
-------------------------------------------------------------------------------
2001.0  X     
-------------------------------------------------------------------------------
 Begin rrpr
##############################################################################
Inventory for date = 2007-12-04 00:00:00
PRES   SST    
-------------------------------------------------------------------------------
2001.0  X     
-------------------------------------------------------------------------------
Subroutine DATINT: Interpolating 3-d files to fill in any missing data...
Looking for data at time 2007-12-04_00
Found file:      SST:2007-12-04_00
Looking for data at time 2007-12-04_06
ERROR: Data not found: 2007-12-04_06:00:00.0000
FORTRAN STOP
(STOP原因:由于只有一个时次,所以它没法插值)
但愿以后不要出错。可以多弄几个时次看看结果。OK!哈哈!
另外,运行完./metgrid.exe之后还要对东海海区进行增加一度的试验,又是一个问题,区域改值!!!!!!!!!
继续研究,现在相当有信心啊,哈哈!
  评论这张
 
阅读(659)| 评论(1)
推荐 转载

历史上的今天

评论

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

页脚

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