明星のGPSゾンデデータから時間・高度断面図を作成する その2

2012年5月東シナ海 Leg 1

5/15- 5/17

[2016年 10月 25日 火曜日 09:29:03 JST]

[~/2016.Kawai.Assim/Time-Height.Diagram/grads]

[am@aofd165]

ga-> ronbun.leg1.03.gs

All files closed; all defined objects released;

All GrADS attributes have been reinitialized

Cannot use log scaling when coordinates <= 0

Linear scaling used

Output: Ronbun.leg1.03.eps

$ srcdump.sh ronbun.leg1.03.gs

aofd165.bio.mie-u.ac.jp

/work1/am/2016.Kawai.Assim/Time-Height.Diagram/grads

Tue Oct 25 10:48:20 JST 2016

======================

ronbun.leg1.03.gs

======================

'reinit'

'open ./ctl/Sonde.ctl'

'set display color white'

'c'

'set grads off'

#### RGB set ####

'set rgb 20 250 250 250'

'set rgb 21 230 230 230'

'set rgb 22 210 210 210'

'set rgb 23 190 190 190'

'set rgb 24 160 160 160'

'set rgb 25 130 130 130'

'set rgb 26 100 100 100'

'set rgb 27 70 70 70'

'set rgb 28 40 40 40'

# 'set rgb 29 40 40 40'

#################

'set rgb 31 210 210 210'

'set rgb 32 200 200 200'

'set rgb 33 190 190 190'

'set rgb 34 180 180 180'

'set rgb 35 170 170 170'

'set rgb 36 150 150 150'

'set rgb 37 130 130 130'

'set rgb 38 110 110 110'

'set rgb 39 90 90 90'

'set rgb 40 70 70 70'

'set rgb 41 50 50 50'

'set rgb 42 30 30 30'

###

'set arrlab off'

'set parea 1 7.8 7.4 10.4'

'set xlopts 1 0 0'

'set time 18z15may2012 22z17may2012'

'set lev 0 5000'

'set ylint 500'

### EPT ####

'set gxout grfill'

'set cmin 280'

'set cmax 350'

'set cint 5'

'd ept'

'xcbar 7.9 8.0 7.4 10.2 -line on -edge circle'

############

'set clevs 5 10 15'

'set ccols 1 1 1 2'

'set arrscl 0.5 50'

'set ccolor 0'

'set cthick 10'

'd skip(u,8,1);v'

'set arrscl 0.5 50'

'set ccolor 1'

'set cthick 2'

'vec.gs skip(u,8,1);v -SCL 0.5 50 -P 7 10.75 -SL m/s -HDSZ 0.1'

# 'd skip(u,8,1);v'

'draw title ECS May, 2012'

# 'draw ylab Height [m]'

'close 1'

###

'set parea 1 7.8 4.2 7.2'

'open ./ctl/Xbt.ctl'

# 'set time 00z17may2012 12z17may2012'

'set time 18z15may2012 22z17may2012'

'set lev 0 60'

'set ylint 10'

'set yflip on'

'set zlog on'

'set gxout grfill'

# 'set clevs 20 21 21.5 22 22.5 23 23.5 24 24.5 25 26'

'set clevs 15 16 17 18 19 20 21 22 23 24 25'

# 'set ccols 42 41 40 39 38 37 36 35 34 33 32 31'

'd st'

'xcbar 7.9 8.0 4.2 7.0 -line on -edge circle'

'close 1'

'set yflip off'

'open ./ctl/Gps.ctl'

'set xlopts 1 1 0'

'set parea 1 7.8 2.5 4'

# 'set time 00z17may2012 12z17may2012'

'set time 18z15may2012 22z17may2012'

'set vrange 17 27'

'set ylint 2'

'set ccolor 2'

'set cmark 0'

'set cthick 3'

'set cstyle 1'

'd sst'

'set ccolor 4'

'set cmark 0'

'set cthick 3'

'set cstyle 1'

'd ts'

'set xlopts 1 3 0.15'

'set parea 1 7.8 0.8 2.3 '

# 'set time 00z17may2012 12z17may2012'

'set time 18z15may2012 22z17may2012'

'set vrange 1006 1012'

'set ylint 1'

'set ccolor 8'

'set cmark 0'

'set cthick 3'

'set cstyle 1'

'd p'

'set string 1.2 c 0.15'

'set strsiz 0.15 0.2'

'draw string 0.3 0.3 S'

'draw string 8.0 0.3 N'

'set strsiz 0.1 0.1'

# 'draw string 6.65 6.75 [m/s]'

# 'draw string 1.3 6.93 (b)'

# 'draw string 1.3 4.88 (c)'

# 'draw string 1.3 2.25 (d)'

# 'set strsiz 0.15 0.2'

# 'draw string 3.9 0.4 S1'

# 'draw string 4.4 0.4 S2'

# 'draw string 4.9 0.4 S3'

# 'draw string 5.35 0.4 S4'

# 'draw string 5.85 0.4 S5'

# 'draw string 6.35 0.4 S6'

# 'draw string 6.85 0.4 S7'

# 'draw string 7.3 0.4 S8'

# 'draw string 7.8 0.4 S9'

'set strsiz 0.13 0.15'

'set string 1 c 4 90'

'draw string 0.3 9 Height [m]'

'draw string 0.3 6 Depth [m]'

'draw string 0.3 3.3 [`ao`nC]'

'draw string 0.3 1.5 [hPa]'

'set string 1 c 4 0'

'draw string 8.2 10.3 [K]'

'draw string 8.2 7.1 [`ao`nC]'

'set string 1 c 4 0'

'set strsiz 0.12 0.15'

'set string 2'

'draw string 1.3 3.7 SST'

'set string 4'

'draw string 1.3 3.4 SAT'

out='Ronbun.leg1.03.eps'

'gxprint 'out

say

say 'Output: 'out

say

----------------------

End of ronbun.leg1.03.gs

----------------------