- 积分
- 669
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-3-24
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
求助 之前没打点的时候是可以每页显示6张图的 打点之后就只显示第一张图了TAT
图在下面
gs文件:
(只看到d smath9(a2)就行 后面的都是一样的)
'reinit '
'open C:\Users\pc.pc-PC\Desktop\xg\haha.ctl'
'open C:\Users\pc.pc-PC\Desktop\xg\Had_60-11-3-8.ctl'
'enable print C:\Users\pc.pc-PC\Desktop\tcorrjscs3-8.gmf'
'define_colors'
*****Define
'set x 1'
'set y 1'
'set z 1'
'set t 1 52'
'define js1=p1'
'define js2=p2'
'define js3=p3'
'define js4=p4'
'define js5=p5'
'define js6=p6'
'set dfile 2'
'set lon -177.5 177.5'
'set lat -87.5 87.5'
'set z 1'
'set t 1'
'define a1=tcorr(js1,t3,t=1,t=19)'
'define a2=tcorr(js2,t4,t=1,t=19)'
'define a3=tcorr(js3,t5,t=1,t=19)'
'define a4=tcorr(js4,t6,t=1,t=19)'
'define a5=tcorr(js5,t7,t=1,t=19)'
'define a6=tcorr(js6,t8,t=1,t=19)'
***************
'set grads off'
'set grid off'
***************************************************
'run page 12 32'
'set map 1 1 4'
'set mproj scaled'
*'set grads off'
'set xlopts 1'
'set gxout grfill'
'set clevs -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5'
'set ccols 47 46 45 44 43 42 22 23 24 25 26 27 '
'd smth9(a1)'
****************************
*'set gxout contour'
*'set ccolor 1'
*'d smth9(a1)'
*****************************
ilat=-87.5
while(ilat<=87.5)
'set lat 'ilat
ilon=-177.5
while(ilon<=177.5)
'set lon 'ilon
'd smth9(a1)'
b=subwrd(result,4)
*say result
*say b
if( (b >= 0.273 & b <= 1) | (b>=-1 &b<=-0.273) )
'q w2xy 'ilon' 'ilat''
x=subwrd(result,3)
y=subwrd(result,6)
* say x
* say y
'set line 1 1 4'
'draw mark 1 'x' 'y' 0.08'
*'q xy2w 'x' 'y''
endif
*endif
ilon=ilon+5
endwhile
ilat=ilat+5
endwhile
************
'run page 22 32'
'set map 1 1 4'
'set mproj scaled'
'set grads off'
'set xlopts 1'
'set gxout grfill'
'set clevs -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5'
'set ccols 47 46 45 44 43 42 22 23 24 25 26 27 '
'd smth9(a2)'
************
'set gxout contour'
'set ccolor 1'
'd smth9(a2)'
*'cbarn'
************
ilat=-87.5
while(ilat<=87.5)
'set lat 'ilat
ilon=-177.5
while(ilon<=177.5)
'set lon 'ilon
'd smth9(a1)'
b=subwrd(result,4)
*ay result
*say b
if( (b >= 0.273 & b <= 1) | (b>=-1 &b<=-0.273) )
'q w2xy 'ilon' 'ilat''
x=subwrd(result,3)
y=subwrd(result,6)
*ay x
*ay y
'set line 1 1 4'
'draw mark 1 'x' 'y' 0.08'
*'q xy2w 'ilon' 'ilat''
endif
*endif
ilon=ilon+5
endwhile
ilat=ilat+5
endwhile
******************
'run page 32 32'
'set map 1 1 4'
'set mproj scaled'
'set xlopts 1'
'set gxout grfill'
'set clevs -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5'
'set ccols 47 46 45 44 43 42 22 23 24 25 26 27 '
'd smth9(a3)'
************
'set gxout contour'
'set ccolor 1'
'd smth9(a3)'
************
'run page 11 32'
'set map 1 1 4'
'set mproj scaled'
'set grads off'
'set xlopts 1'
'set gxout grfill'
'set clevs -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5'
'set ccols 47 46 45 44 43 42 22 23 24 25 26 27 '
'd smth9(a4)'
************
'set gxout contour'
'set ccolor 1'
'd smth9(a4)'
*********************
************
'run page 21 32'
'set map 1 1 4'
'set mproj scaled'
'set grads off'
'set xlopts 1'
'set gxout grfill'
'set clevs -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5'
'set ccols 47 46 45 44 43 42 22 23 24 25 26 27 '
'd smth9(a5)'
************
'set gxout contour'
'set ccolor 1'
'd smth9(a5)'
************
'run page 31 32'
'set map 1 1 4'
'set mproj scaled'
'set grads off'
'set xlopts 1'
'set gxout grfill'
'set clevs -0.5 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5'
'set ccols 47 46 45 44 43 42 22 23 24 25 26 27 '
'd smth9(a6)'
************
'set gxout contour'
'set ccolor 1'
'd smth9(a6)'
********************
*'draw title the relationship between extreme precipitaton and temperature field'
'printim C:\Users\pc.pc-PC\Desktop\jlq.jpg white'
'print'
'disable print'
'reinit'
|
|