- 积分
- 3287
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-8-2
- 最后登录
- 1970-1-1
![[舒美佳人] 粉丝数:6911 微博数:10772 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
绘图出错信息:gs文件如下:请帮忙纠错!
*light white to dark gray
'set rgb 80 255 255 255'
'set rgb 81 216 216 216'
'set rgb 82 191 191 191'
'set rgb 83 165 165 165'
'set rgb 84 145 145 145'
'set rgb 85 128 128 128'
'set rgb 86 100 100 100'
'set rgb 87 88 88 88'
'set rgb 88 77 77 77'
'set rgb 89 65 65 65'
'reinit'
'open D:\2012tshan\newpic\20130617\20130617.ctl'
'enable print D:\2012tshan\newpic\20130617\11\pm41N.gmf'
t=1
while(t<=17)
'c'
'set t 't
'open D:\2012tshan\newpic\20130617\11\topo10.ctl'
'query time'
res=subwrd(result,3)
year=substr(res,11,2)
month=substr(res,6,3)
day=substr(res,4,2)
hour=substr(res,1,2)
if(month='JAN'); month='01'; mtype='big'; endif
if(month='FEB'); month='02'; mtype='feb'; endif
if(month='MAR'); month='03'; mtype='big'; endif
if(month='APR'); month='04'; mtype='sma'; endif
if(month='MAY'); month='05'; mtype='big'; endif
if(month='JUN'); month='06'; mtype='sma'; endif
if(month='JUL'); month='07'; mtype='big'; endif
if(month='AUG'); month='08'; mtype='big'; endif
if(month='SEP'); month='09'; mtype='sma'; endif
if(month='OCT'); month='10'; mtype='big'; endif
if(month='NOV'); month='11'; mtype='sma'; endif
if(month='DEC'); month='12'; mtype='big'; endif
endmonth=FALSE
if(hour<16)
hour=hour+8
else
hour=hour+8-24
if((mtype='big' & day=31) | (mtype='sma' & day=30))
month=month+1
if(month<10); month='0'%month; endif
if(month=13); month='01'; endif
endmonth=TRUE
endif
if(mtype='feb' & day=28)
month='03'
endmonth=TRUE
endif
if(endmonth=FALSE)
day=day+1
if(day<10); day='0'%day; endif
else
day='01'
endif
endif
if(hour<10); hour='0'%hour; endif
'set ccolor 4 '
'set strsiz 0.3 0.8 '
'draw title 'year'-'month'-'day' 'hour':00'
'set parea 1.0 10.0 1.2 8.0'
'set lev 1000 100'
'set lat 30 50'
'set lon 60 100'
'define vor=hcurl(ugrdprs,vgrdprs)*1e5'
************************thetase***********************
'define tt=TMPprs'
'define es=6.11*exp(17.26*(tt-273.16)/(tt-35.86))'
'define qs=0.622*es/(lev-0.378*es)'
'define qq=RHprs/100.*qs'
'define prs=lev*100'
'define el=log((qq/0.62197)*(prs/100.)/(1.+(qq/0.62197)))'
'define td=273.15+(243.5*el-440.8)/(19.48-el)'
'define tc=tt-1.21*(tt-td)'
'define e=qq/(qq+0.622)*prs'
'define cpm=1004*(1+0.86*qq)'
'define rm=287.04*(1+0.608*qq)'
'define pc=prs*pow(tc/tt,cpm/rm)'
'define st=tc*pow(100000/(pc-e),287.04/1004)'
'define lvc=2500.8*1000-2274.4*(tt-273.15)'
'define qse=st*exp(qq*lvc/1004/tc)'
'set ylab on'
'set xlab on'
'set zlog on'
'set lat 41'
'set lon 75 85'
'set xlopts 1 6 0.16'
'set ylopts 1 6 0.16'
'set mpdset hires'
'set grads off'
'set grid off'
'set lev 1000 100'
'set gxout shaded'
'd qse'
'set gxout contour'
'set cthick 6'
*'set cstyle 4'
'set ccolor 1'
*'set cmax 0'
'd vvelprs*100'
'set cstyle 4'
'set cthick 9'
'set ccolor 1'
'd qse'
'set gxout vector'
'set ccolor 1'
'set arrscl 0.2 20'
'd ugrdprs;vgrdprs'
'set dfile 2'
'set ylab off'
'set xlab off'
'set x 1'
'set y 1'
'set vrange 0 17000'
'set rgb 16 144 144 144'
'set gxout linefill'
'set lfcols 16 0'
*'set lat 37 45'
'set lon 75 85'
'set z 1'
'set t 1'
*'define cc=ave(tp,lon=70,lon=80)'
*'set lon 77'
'define cc=ave(tp,lat=40.5,lat=41.5)'
'set lat 41'
'define a0=tp-tp'
'd cc;a0'
'close 2'
'print'
'c'
t=t+1
endwhile
'disable print'
return
'reinit'
|
|