'
'set t 1'
'set z 1 12'
'set zlog on'
'set gxout shaded'
'set cint 8000'
'set rgb 18 125 125 125'
'set csmooth on'
'set rbcols 18'
'set cmax 0'
'define zz=dz.1'
'd zz'
'set dfile 2'
'set t 'i
'set lat 40'
'set lon 100 120'
'set lev 1000 100'
'define prs=lev'
'define g=9.8'
'define e=6.112*exp(7.5*(air.2-273.15)/(air.2-35.85))'
'define q=0.62197*e/(prs-0.378*e)'
'set gxout contour'
'set cint 0.5'
'set cthick 9'
'define q1=maskout(q*1000,zz)'
'd q1'
'set gxout vector'
'define u1=maskout(uwnd.3,zz)'
'define w1=maskout(omega.4*10,zz)'
'define uw=maskout(mag(uwnd.3,omega.4),zz)'
'd u1;w1;uw'
*rect(95,112.5,6,8)
'printim e:m_dpicuwq0'i'_40N.png white'
'c'
i=i+1
endwhile
'reinit'
*-------------------------------rectangle-----------------------
function rect(lon1,lon2,lev1,lev2)
'q ll2xy 'lon1' 30'
str=sublin(result,1)
xlon1=subwrd(str,1)
'q ll2xy 'lon2' 30'
str=sublin(result,1)
xlon2=subwrd(str,1)
'q gr2xy 1 'lev1
str=sublin(result,1)
ylev1=subwrd(str,6)
'q gr2xy 1 'lev2
str=sublin(result,1)
ylev2=subwrd(str,6)
'q gr2xy 1 1'
str=sublin(result,1)
ybot=subwrd(str,6)
'q gr2xy 1 12'
str=sublin(result,1)
ytop=subwrd(str,6)
'draw line 'xlon1' 'ylev1' 'xlon2' 'ylev1
'draw line 'xlon1' 'ylev2' 'xlon2' 'ylev2
'draw line 'xlon1' 'ybot' 'xlon1' 'ytop
'draw line 'xlon2' 'ybot' 'xlon2' 'ytop
return
===============grads画某点或站点的垂直风廓线.gs===========================
CTL文件
dset G:forummcla-AS-2009-09-16-040000-g2.gra
undef 1.0e30
options big_endian
title RAMS Output
xdef 68 linear 116.2485504 0.0044996
ydef 52 linear 39.8840141 0.0044996
zdef 32 levels 16.889 54.422 97.585 147.223 204.307 269.953 345.446 432.262 532.102 646.917 778.955 930.798 1105.417 1306.230 1537.164 1802.739 2108.150 2459.372 2863.278 3327.770 3861.935 4476.225 5182.659 5995.058 6929.316 8010.843 9196.865 10399.916 11599.916 12799.916 13999.916 15199.916
tdef 25 linear 04:00Z16sep2009 10mn
vars 5
scalar1 32 99 - RAMS : scalar1 concentration [#/m3]
u 32 99 - RAMS : u [m/s]
v 32 99 - RAMS : v [m/s]
topo 1 99 - RAMS : topo [m]
tempc 32 99 - RAMS : temperature [C]
endvars
GS文件
'reinit'
'open g:forummcl1.ctl'
'set lat 39.91'
'set lon 116.39'
'set t 1 25 2'
'set mpdset beijing'
*'set lat 39.89 40'
*'set lon 116.25 116.6'
'set z 1 17'
'set gxout barb'
'set grid off'
'set grads off'
'd u;v'
'printim g:forummcl116.39_39.91.png white'
'c'
;
==============通过点击任意画箭头.gs=======================================
function arrow (args)
say 'Click where you want the arrow-bottom'
'query bpos'
x1 = subwrd(result,3)
y1 = subwrd(result,4)
say 'Click where you want the arrow-top'
'query bpos'
x2 = subwrd(result,3)
y2 = subwrd(result,4)
dx=x2-x1
dy=y2-y1
'difine l=0.2*sqrt(dx*dx+dy*dy)
'define f=atan(