RHprs不是26层而是21层，运行不了ff_tlop.gs，进而无法做出Skew-T图。请指点。

kql87919

看了《用GrADS绘制的Skew-T图》的帖子，准备一试。
用ff_tlop.gs准备从fnl资料里fwrite一部分温度，相对湿度，u和v，ff_tlop.gs文件如下：
'reinit'
'open d://201204.ctl'
'set gxout fwrite'
'set fwrite d:/tlop.grd'

'set lat 20 40'
'set lon 110 130'
t=1
while(t<=7)
   'set t ' t
   z=1
   while(z<=26)
     'set z ' z
     'd TMPprs'
     'd RHprs'
     'd UGRDprs'
     'd VGRDprs'
     z=z+1
     endwhile
   t=t+1
endwhile
'disable fwrite'
;

但是由于RHprs不是26层而是21层，运行不了ff_tlop.gs，进而无法做出Skew-T图。请高人指点

传说中的谁

相对湿度只有21层，最高到100hPa，循环的时候你可以只取到第21层就可以啦

kql87919

传说中的谁 发表于 2012-4-25 16:52

                               
登录/注册后可看大图

相对湿度只有21层，最高到100hPa，循环的时候你可以只取到第21层就可以啦

明白啦，谢谢

传说中的谁

kql87919 发表于 2012-4-25 17:12

                               
登录/注册后可看大图

明白啦，谢谢

反正日常的预报中也很少关注100hPa以上的东西

朕小通

传说中的谁 发表于 2012-4-25 17:15

                               
登录/注册后可看大图

反正日常的预报中也很少关注100hPa以上的东西

请教下，grd数据的排放顺序应该是先变量，后高度。他编的gs是先高度，后变量。这样写出的grd能正确么？

易小凯

朕小通 发表于 2012-4-30 21:33

                               
登录/注册后可看大图

请教下，grd数据的排放顺序应该是先变量，后高度。他编的gs是先高度，后变量。这样写出的grd能正确么？

如果只是用在fortran中就没啥问题
要是根据这个写ctl就会出问题
实况哥v587