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

美丽心情

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

 
 
 

日志

 
 

引用 Grads---lats4d.gs(转动力论坛funny)  

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

  下载LOFTER 我的照片书  |

 

引用

ahualian引用 Grads转成Netcdf(转动力论坛funny)

 

引用

zhangljGrads转成Netcdf(转动力论坛funny)
用GrADS写NetCDF、 HDF-SDS或GRIB等格式文件,请使用lats4d.gs:

说明:
可对5种格式的数据进行转换:GrADS IEEE、GSFC Phoenix、GRIB、NetCDF、HDF-SDS。
输出格式有4种:GRIB、NetCDF、HDF-SDS、flat binary。

选项:
-i filename:指定输入文件,
-o filename:指定输出文件(缺省grads.lats),
-cal calender: standard(缺省)、noleap、clim、climleap,
-center ctr:指定资料处理中心(如 PCMDI、GSFC、NCEP),
-de filename:指定维数设置文件,
-format fmt:指定输出文件的格式(输出GRIB推荐使用grads-grib选项),
-ftype ctl|sdf|xdf:指定输入文件类型,
-freq [n] unit:指定输入数据文件的时间步长,unit不能为分钟,
-func expr:在输出前做资料处理的运算,
-grid type: 指定数据格类型(linear缺省、gaussian、generic),
-h:显示帮助信息,
-lat y1 y2:指定纬度范围(缺省用维数环境),
-levs lev1…levN:指定层次(缺省:所有层),
-lon x1 x2:指定经度范围(缺省用维数环境),
-mean:做时间平均,做平均的时间步长用-time选项设置,
-model mod:生成数据资料的模式名称,
-precision nbits:指定grib格式精度(bit数,缺省 nbits=16),与-table一起用,
-table tab: 指定要素列表文件名,
-time t1 t2 [tincr]:指定数据的时间采样,tincr指定跳过的时间间隔,
-title text:输出grib文件时的数据title,
-v:详细方式,
-vars var1…varN:指定要素(缺省为所有要素),
-xsfc:剔除所有地面要素,
-xupper:剔除所有高空要素,
-xvars var1…varN:指定要剔除的要素,
-zrev:垂直方向逆序,
-q:执行完毕退出grads。

扩展:
可以与用户自定义函数(UDF) regrid2() 配合使用,实现差值场输出成特定格式。

范例:
1、在"gradsnc" or "gradshdf" 的grads环境下
    lats4d -h                                                      
    lats4d -v -q -i model -o ex1                                    
    lats4d -v -q -i model -o ex2a -format grads_grib                 
    lats4d -v -q -i model -o ex2b -format grads_grib -precision 8   
    lats4d -v -q -i model -o ex2c -format grads_grib -precision 32   
    lats4d -v -q -i model -o ex3 -levs 700 500 -vars ua va           
    lats4d -v -q -i model -o ex4 -time 1jan1987 3jan1987            
    lats4d -v -q -i model -o ex5 -time = = 2                        
    lats4d -v -q -i model -o ex6 -mean                              
    lats4d -v -q -i model -o ex7 -mean -time = = 2                  
    lats4d -v -q -i model -o ex8 -lat 20 70 -lon -140 -60   
        
2、将 model 文件中的4X5(latxlon)场输出为 1X1数据场:
    lats4d -i model -de geos1x1 -func regrid2(var,1,1,bs_p1,-180,-90)

BUGS:
1、  Sometimes lats4d will only work if you exit and                 
    restart GrADS.                                                   
                                                                    
2、  The option "-precision 32" does not quite work. This            
    appears to be a LATS bug.                                       
                                                                    
3、  Because of a limitation in the GRIB format, "grib" or           
    "grads_grib" output cannot have levels where p<1.               
    To circumvent this problem, a hybrid level number is            
    is used in such cases.                                        
  评论这张
 
阅读(326)| 评论(0)
推荐 转载

历史上的今天

评论

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

页脚

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