grads处理JRA逐日资料问题

carpenter3

     
      我想利用JRA-25逐日资料（每隔6h）做三年夏季每个月平均的比湿垂直分布图，就是画三张垂直剖面图，分别是6、7、8月份，我之前参考http://blog.sina.com.cn/s/blog_4d5c3ba901000bm1.html 写出来的gs文件一直有错误，而且利用cmd显示能生成idx文件，但是找不着在哪。纠结了好久还是各种，请教高人指教下gs文件要怎么写还有那个问题是怎么回事！谢谢了！
  ctl文件已上传附件。

carpenter3

我自己顶下

kongfeng0824

请问JRA-25是关于什么的资料呢？有降水资料吗？可以给一个链接地址码？谢谢

善人/jw

LZ把ctl直接贴出来，还有你用gs时的错误提示也贴出来吧

carpenter3

善人/jw 发表于 2013-5-28 00:11

                               
登录/注册后可看大图

LZ把ctl直接贴出来，还有你用gs时的错误提示也贴出来吧

我的gs文件错误很多就不帖了，直接做个1996年6月份的也行，主要是我的批处理时循环结构还不是很熟悉，麻烦指教下
ctl文件：

dset  I:\anl_p.%y4%m2%d2%h2
options template
index I:\RCM\REGCM_GLOBALDAT\JRA_25\anl_p_daily\anl_p.idx
undef 9.999E+20
title I:\RCM\REGCM_GLOBALDAT\JRA_25\anl_p_daily\anl_p.1996
* produced by grib2ctl v0.9.12.5p32d
dtype grib 255
ydef  145  linear  -90.0  1.25
xdef  288  linear    0.0  1.25
tdef 1104 linear 00Z01jun1996 6hr
zdef 23  levels
1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10 7 5 3 2 1 0.4
vars 13
HGTprs   23 7,100,0   ** Geopotential height            [gpm]
TMPprs   23 11,100,0  ** Air temperature                  [K]
SPFHprs  12 51,100,0  ** Specific humidity            [kg/kg]
DEPRprs   8 18,100,0  ** Dew point depression             [K]
UGRDprs  23 33,100,0  ** Zonal wind                     [m/s]
VGRDprs  23 34,100,0  ** Meridional wind                [m/s]
CWATprs  12 76,100,0  ** Cloud water content          [kg/kg]
PRMSLmsl  0 2,102,0   ** Pressure reduced to MSL         [Pa]
TMPsfc    0 11,105,2  ** Surface 2m air temperature       [K]
SPFHsfc   0 51,105,2  ** Surface 2m specific humidity [kg/kg]
DEPRsfc   0 18,105,2  ** Surface 2m dew point depression  [K]
UGRDsfc   0 33,105,10 ** Surface 10m zonal wind         [m/s]
VGRDsfc   0 34,105,10 ** Surface 10m meridional wind    [m/s]
ENDVARS

carpenter3

kongfeng0824 发表于 2013-5-27 22:01

                               
登录/注册后可看大图

请问JRA-25是关于什么的资料呢？有降水资料吗？可以给一个链接地址码？谢谢

有降雨，直接搜JRA-25就能找着

carpenter3

求指教啊
这是我写的求一个月的
'reinit'
'open   I:\anl_p.ctl'
tm=1
while(tm<121)
'set t' tm
'set gxout contour'
'set grid on'
'set grads off'
'set xlopts 1 5 0.15'
'set ylopts 1 5 0.15'
'set lon 70 106'
'set lat 30 '
'set xlopts 1 5 0.15'
'set ylopts 1 5 0.15'
'set cthick 6'
'set ccolor 1'
'set z 1000 100'
'set string  1 tr 6 0'
'set strsiz 0.15'
'define tt=ave(SPFHprs,t='tm',t='tm+3',1)'
‘define ttt=ave(tt,lon=70,lat=106)’
'draw string 5.5 4.1  JRA '
tm=tm+4
endwhile
'd ttt'
'print'
'printim I:\199606.png x1000y800 white '
'disable print'

;

善人/jw

把经度也设置成一个值
set lon 70

guoxiaoguo

有没有大神处理过JRA-55的资料的啊！！！

笨蛋学徒

guoxiaoguo 发表于 2014-8-12 17:31
有没有大神处理过JRA-55的资料的啊！！！

请问JRA-55资料怎么处理啊？？？