grads计算水汽通量和水汽通量散度的gs

小猴纸

不知道为什么算的不对 来学习对比一下

开门查看

水汽通量散

目美儿

{:eb511:}

dannuo526

好东西 好分享 6666

RickestRick

想看！啊怎么还要十个字要凑够十个字啊

kmzgl

谢谢楼主的分享！！！

hugeal10000

找了很久了！

jinglj2001

*700水汽通量散度（阴影）+等高线+风羽图
'reinit'
*'set mpdset hires cnriver cnworld'
iyear="19880722"
buwei="0"
tt=0
while (tt<13)
'set mpdset d:\grads_read_ncep\cnworld'
'set map 15 1 5'
'set rgb 16 250 250 250'
'set rgb 17 235 235 235'
'set rgb 18 220 220 220'
'set rgb 19 205 205 205'
'set rgb 20 190 190 190'
'set rgb 21 175 175 175'
'set rgb 22 160 160 160'
'set rgb 23 145 145 145'
'set rgb 24 130 130 130'
'set rgb 25 115 115 115'
'set rgb 26 100 100 100'
'set rgb 27  85  85  85'
'set rgb 28  70  70  70'
'set rgb 29  55  55  55'
'set rgb 30  40  40  40'
'set rgb 31  25  25  25'
'set rgb 32  10  10  10'
'set rgb 33   5   5   5'
'set rbcols 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33'
**时间 tt不足2位数时补足2位，否则找不到NC文件
if (tt>=10)  
'sdfopen d:\grads_read_ncep\19880720-0723\19880722\era5_level_'  iyear  ''  tt  '.nc'
endif
if (tt<=9)   
'sdfopen d:\grads_read_ncep\19880720-0723\19880722\era5_level_' iyear ''  buwei  ''  tt  '.nc'
endif
*'set font 3'
'draw mark 3 83.6 45.933 16'
'set parea 1.40 10.47 0.7 7.73'
'SET LON 40 115'
'SET LAT 30 60'
'set t 1'
'set clopts -1 -1 0.12'
'set xlopts 1 4 0.18'
'set ylopts 1 4 0.16'
'set strsiz 0.2 0.22'
'set grads off'
'set grid off'
'set csmooth on'
'set xlpos -10'
'set ylpos -10'
'set clopts -16 0.18 '
'SET LON 70 90'
'SET LAT 40 50'
'set lev 700'
'define uq=u*q/9.8'
'define vq=v*q/9.8'
'define divq=hdivg(uq,vq)'
'SET LON 40 115'
'SET LAT 30 60'
*'set font 1'
'set gxout shaded'
'set cint 0.5'
*'set clevs  5 10 15 20 25 30'
'set clevs  -15 -5 0 5  '
*'set ccols  0 17 20 23 26 29  32'
'set ccols   26 20 0 0'
'set black  -10 0.099'
*'set black  0 100'
'set csmooth on'
'd smth9(divQ*100000000)'
'run d:\opengrads\contents\resources\scripts\cbarn.gs 0.9 0'
'set gxout barb'
'set ccolor 1'
'set cthick 4'
**'d skip(u*2.5,5,5);skip(v*2.5,5,5)'
'd skip(u*2.5,10,10);skip(v*2.5,10,10)'
'set gxout contour'
'set clskip 2'
'set ccolor 1'
'set cthick 6'
'set cint 4'
*'set clevs'
'd smth9(z/98)'
'set clskip 1'
'set cthick 8'
'set ccolor 2'
'set cint 4'
'set cstyle 2e 1'
*'d t-273.15'
*'set mpdset =d:\grads_read_ncep\cnworld'
**shi为北京时间
shi=tt+8
'draw mark 3 83.6 45.933 16'
'run d:\opengrads\contents\resources\scripts\axis.gs -type B  -position o -start 40 -end 115 -label on  -lfont 4 -lsize 0.22 -lthick 3.0 -interval 10 -suffix `3.'
'run d:\opengrads\contents\resources\scripts\axis.gs -type L  -position o -start  30 -end  60 -label on  -lfont 4 -lsize 0.22 -lthick 3.0 -hoffset -0.2 -interval 5 -suffix `3.'
*'gxprint d:\grads_read_ncep\1966052407\700HGTW-1966052407.png white x1024  y768'
**在地图上托里县的位置画一个圆点
'q w2xy 83.36  45.93'
x1 = subwrd(result,3)
y1 = subwrd(result,6)
'set string 1 c 6'
'set strsiz 0.115'
'draw string 'x1' 'y1'  `32'
'gxprint d:\grads_read_ncep\19880720-0723\19880722\700H-WIN-SQTLSD\700htW-' iyear '-' shi '.png  white x1024 y768 '
tt=tt+1
'c'
'reinit'
endwhile
*'c'
*'reinit'
;

999zyy

感谢楼主

dcxswfff1

求真实体会gs自己写才有深刻记忆