读取JRA数据

huayazi

各位高手：
本人想处理JRA-55的数据，由于格式是grib的，因此想用ncl_convert2nc先把grib数据转换为nc格式的先。当直接输入命令ncl_convert2nc ./201301/anl_surf.2013013118.grib -v U_GRD_GDS4_HTGL,V_GRD_GDS4_HTGL,g4_lat_0,g4_lon_1所生成的nc数据是正确的，里头含有U和V相关变量及经纬度信息共4个变量。
由于想批量处理，写了个shell
#!/bin/bash
source ~/.bash_profile
cd /home/qxs/dxh/201301
for file in /home/qxs/dxh/201301/anl_surf.20*
do
        ln -s ${file} $file".grib"
        /share/local/ncl_ncarg/bin/ncl_convert2nc $file".grib" -v U_GRD_GDS4_HTGL,V_GRD_GDS4_HTGL,g4_lat_0,g4_lon_1
        ncks -d g4_lon_1,120.0,120.0 -d g4_lat_0,40.0,40.0 -v U_GRD_GDS4_HTGL,V_GRD_GDS4_HTGL $file".nc" $file"_new.nc"
done


结果生成的nc文件里头就缺乏g4_lon_1变量了。


哪位高手帮忙解决一下？

huayazi

藕茭白空心菜

楼主解决了吗