任意两点间斜剖面做法总结【利用GRADS自带函数】

liuyw

river 发表于 2017-9-7 22:10
不兼容的网格？你地形的资料有插值成0.25的吗？好像弄那个插值不成那么高的分辨率

我是在您的1*1地形资料上又差值成0.25的，如果不行的话，有地方可以下载0.25*0.25的地形数据？

river

liuyw 发表于 2017-9-8 08:18
我是在您的1*1地形资料上又差值成0.25的，如果不行的话，有地方可以下载0.25*0.25的地形数据？

也没有啦，我就是自己猜测的。你的地形都能画出来说明没问题。你是不是其他数据资料不是你说的0.25*0.25的分辨率啊

liuyw

river 发表于 2017-9-8 16:07
也没有啦，我就是自己猜测的。你的地形都能画出来说明没问题。你是不是其他数据资料不是你说的0.25*0.25 ...

在这个过程中，使用了两个资料，一是插值得到的0.25*0.25的地形资料，二是FNL0.25*0.25。都是没有问题的啊，楼主大神，可以留个联系方式，私下交流一下吗？

river

liuyw 发表于 2017-9-8 16:17
在这个过程中，使用了两个资料，一是插值得到的0.25*0.25的地形资料，二是FNL0.25*0.25。都是没有问题的 ...

你把你地形插值的脚本直接贴上来，我试一下。如果你的FNl 资料也不太大，也传上来

liuyw

river 发表于 2017-9-8 18:25
你把你地形插值的脚本直接贴上来，我试一下。如果你的FNl 资料也不太大，也传上来

一、地形脚本如下：
'reinit'
'open E:\gs\dixing\1X1\orog1.ctl'
'set gxout fwrite'
'set fwrite E:\gs\dixing\0.25X0.25\orog0.25.grd'
'set t 1'
zz=1
while(zz<=17)
'set z 'zz''
'set x 1 360'
'set y 1 180'
'd re(orog,1440,'linear',0,0.25,720,'linear',-90,0.25,bl)'
zz=zz+1
endwhile

'disable fwrite'
;
二、FNL资料太大了，我把ctl文件贴上，如下。具体祥见附件1.
dset D:/data/gdas1.fnl0p25.201707%d2%h2.f00.grib2
options template
index D:/data/gdas1.fnl0p25.2017070700.f00.grib2.idx
undef 9.999E+20
title D:/data/gdas1.fnl0p25.2017070700.f00.grib2
*  produced by g2ctl v0.0.3d
* griddef=1:0:(1440 x 721):grid_template=0:winds(N/S): lat-lon grid:(1440 x 721) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 0.250000 lon 0.000000 to 359.750000 by 0.250000 #points=1038240

dtype grib2
ydef 721 linear -90.000000 0.25
xdef 1440 linear 0.000000 0.250000
tdef 8 linear 00Z07jul2017 6hr
*  z has 31 levels -prs
zdef 31 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000 700 500 300 200 100
options pascals
vars 121

三、我在您的基础上，对脚本进行了修改，添加了在平面图上鼠标选择起止点，没有调用地形文件，脚本见附件2.

river

liuyw 发表于 2017-9-8 18:43
一、地形脚本如下：
'reinit'
'open E:\gs\dixing\1X1\orog1.ctl'

你什么都不干，直接打开ctl ，然后画你用的那几个变量，看看能不能出图？

liuyw

river 发表于 2017-9-8 18:48
你什么都不干，直接打开ctl ，然后画你用的那几个变量，看看能不能出图？

可以出图的，就是加上地形文件后，就不出了

river

liuyw 发表于 2017-9-8 18:50
可以出图的，就是加上地形文件后，就不出了

是不是你路径没修改？我看你脚本里地形数据用的还是1X1的，'open E:\gs\dixing\1X1\orog1.ctl'。你自己看看

小其其格

'd maskout(coll2gr(2,-u),-coll2gr(3,-u));maskout(coll2gr(1,-u),-coll2gr(3,-u));mag(coll2gr(2,-u),coll2gr(1,-u))'

这里绘制的为什么一直是-u，而不是具体的u,v或者omega之类的？

river

小其其格 发表于 2017-9-8 19:24
'd maskout(coll2gr(2,-u),-coll2gr(3,-u));maskout(coll2gr(1,-u),-coll2gr(3,-u));mag(coll2gr(2,-u),col ...

自己去看coll2gr用法