求基于Ncep1.0*1.0资料的流线垂直剖面，要去地形，可批量出图

micolian

求基于Ncep1.0*1.0资料的流线垂直剖面，要求同时加地形，可批量出图。目前可以单独实现平行于经纬度的加地形剖面，可以单独批量出流线图，可以单独出斜线剖面图，但是能够斜线、批量同时加地形的垂直剖面图出不来，主要是本人grads编程不够精通，手头有现有的单独实现程序却无法整合，所以悬赏求大神们实现斜线剖面、批量出流线图、加地形三个功能的gs文件整合，多谢！！！三个独立功能的gs文件如下：
*************加地形**************
'reinit'
*newfile.ctl为多时次合并后的Ncep1.0*1.0资料描述文件
'open e:\1603\newfile.ctl'
*topo10.ctl为地形文件的描述文件，附件中将附带地形文件和描述文件，并且加入合并Ncep1.0*1.0资料的批处理grib.bat
'open e:\1603\topo10.ctl'
'enable print e:\1603\liuxian.gmf'
'set parea 1.0 10.0 1.2 8.0'
'set grads off'
'set grid off'
'set xlopts 1 6 0.18'
'set ylopts 1 6 0.18'
'set xlint 1'

'set lon 75 85'
'set lat 35 47'
'set lev 850 200'
'set zlog on'
'set ylevs 850 700 600 500 400 300 200'

'set lev 850 200'
'set lat 34 50'
'set lon 75'

'set t 9'

'set csmooth on'
'set clab off'
'set clevs 0'
'set gxout stream'
'set clab on'

'set ccolor 1'
'set cthick 8'

'set clopts 1 6 0.18'
'set clab forced'

'set csmooth on'
'd ugrdprs;vgrdprs'

'set dfile 2'
'set ylab off'
'set xlab off'
'set vrange 1500 12000'
'set x 1'
'set y 1'
'set z 1'
'set t 1'

***Y轴数据变化范围（850百帕对应1500米，200百帕对应12000米）

'set rgb 16 144 144 144'
'set gxout linefill'
'set lfcols 15 0'
'set lat 41 45'
'set lon 82.1'
***剖面的经度
'define cc=ave(tp,lon=82.0,lon=82.2)'
***剖面经度相邻的两点
'define a0=tp-tp'
'd cc;a0'
'print'
'disable print'
;

*******************批量出图******************
'reinit'
'set mpdset cnworld cnriver'
'set grads off'
'set csmooth on'
'open E:\1512\newfile.ctl'
'enable print E:\1512\500wind111222.gmf white x600 y450'
'set lon 15 105'
'set lat 0 90'
'set clopts -1 -1 0.12'
'set xlint 15'
'set ylint 15'
'set xlopts 1 4 0.18'
'set ylopts 1 4 0.16'
'set strsiz 0.2 0.22'
'set grads off'

'set lev 500'
n=1
while(n<=23)
'set t 'n
'set grid on'
'set grads off'


'set csmooth on'


'set gxout stream'

'd ugrdprs;vgrdprs'
'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
*'set annot 2'
'draw title 500hPa UV'recs''

'print'
'c'
n=n+1
endwhile
'disable print'

**********斜线剖面**************
'reinit'
'set mpdset cnworld cnriver'
'set grads off'
'set csmooth on'
'set zlog on'
'open E:\1603\newfile.ctl'
'enable print E:\1603\vvvvv1111.gmf white x600 y450'
'set x 1'
'set y 1'
'set lev 1000 100'
lon1 = 60.0
lon2 = 80.0
lat1 = 55.0
lat2 = 35.0
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(vvelprs,'lon','lat')'
  lon = lon + 1
endwhile

'set map 15 1 9'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'

'set x 14 16'
'set xaxis 'lon1' 'lon2
'set clab on'
'set gxout shaded'
'set cmax 0'
'd coll2gr(1,-u)'
'set gxout contour'
*'set cint .5'
'set cthick 11'
'd coll2gr(1,-u)'
*'d vvelprs'
'draw title omega along (55N,95W) to (15N,90W)'
'print'
'disable print'

QQ联系412854357