Grads转成Netcdf

zhanglj 的 Grads转成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.