時間ー高度断面図（縦軸は気圧）

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

List of the following files:

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

time-p_obs.run.sh

sonde.ctl

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

Machine info

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

aofd165.bio.mie-u.ac.jp

/work2/am/WRF/WRF.Post/GrADS/PL1101/Time-Height_OBS_P

Tue Sep 19 13:50:18 JST 2017

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

time-p_obs.run.sh

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

#!/bin/bash

gs=$(basename $0 .run.sh).gs

export LANG=C

host=$(hostname)

cwd=$(pwd)

timestamp=$(date)

runname="Observation"

out='PL1101_Time-P_OBS.eps'

cat<<EOF>$gs

function timehgt( args )

#

# Default values

#

ctl='sonde.ctl'

outdir='Fig'

lev1=95

lev2=1010

title='PL1101 (In-situ observation)'

#

# Decode options

#

i = 1

while( 1 )

arg = subwrd( args, i )

i = i + 1;

if( arg = '' ); break; endif

while( 1 )

if( arg = '-dm' ); dirmodel= subwrd(args,i);i=i+1; break;endif

if( arg = '-rn' ); runname= subwrd(args,i) ;i=i+1; break;endif

if( arg = '-do' ); dirobs= subwrd(args,i) ;i=i+1; break;endif

if( arg = '-co' ); ctlobs= subwrd(args,i) ;i=i+1; break;endif

if( arg = '-dtime' ); datetime= subwrd(args,i);i=i+1; break;endif

if( arg = '-lon' ); lon= subwrd(args,i) ;i=i+1; break;endif

if( arg = '-lat' ); lat= subwrd(args,i) ;i=i+1; break;endif

if( arg = '-q' ); quitopt=yes ;i=i+1; break;endif

say 'Syntax error : arg= 'arg

return

endwhile

endwhile

'!mkdir -vp 'outdir

say

say 'Open 'ctl

say

'open 'ctl

stime='00z18jan2011'

etime='12z26jan2011'

'set time 'stime' 'etime

### SET TROPOPAUSE HEIGHT ###

indir='/work2/am/2016.PolarLow/Sonde.Profile.04.Lapse.rate.smooth/output'

'q dims'

say result

line=sublin(result,5)

its=subwrd(line,11)

ite=subwrd(line,13)

i=its

#say 'its= 'its

#say 'ite= 'ite

while (i<=ite)

'set t 'i

'q dims'

line=sublin(result,5)

datetime=subwrd(line,6)

hh=substr(datetime,1,2)

dd=substr(datetime,4,2)

mmm=substr(datetime,6,3)

yy=substr(datetime,11,2)

if (mmm='JAN') ;mm=01 ;endif

if (mmm='FEB') ;mm=02 ;endif

if (mmm='MAR') ;mm=03 ;endif

if (mmm='APR') ;mm=04 ;endif

if (mmm='MAY') ;mm=05 ;endif

if (mmm='JUN') ;mm=06 ;endif

if (mmm='JUL') ;mm=07 ;endif

if (mmm='AUG') ;mm=08 ;endif

if (mmm='SEP') ;mm=09 ;endif

if (mmm='OCT') ;mm=10 ;endif

if (mmm='NOV') ;mm=11 ;endif

if (mmm='DEC') ;mm=12 ;endif

#say hh' 'dd' 'mm' 'yy

input=indir'/'yy''mm''dd''hh'_TPH.txt'

res = read(input) ;* SKIP HEADER

res = read(input)

line=sublin( res, 2)

pdatetime.i=subwrd(line,1)

tph.i=subwrd(line,4)

#say 'pdatetime.i='pdatetime.i

#say 'tph.i='tph.i

i=i+1

endwhile

'cc'

'set grads off'

'set grid off'

'set parea 1 7.5 7 10'

'set xlopts 1 0 0'

'set ylopts 1 3 0.15'

'set ylevs 1000 800 700 600 500 400 300 200 100'

'set time 'stime' 'etime

'set yflip on'

'set zlog on'

'set lev 'lev1' 'lev2

### T ####

'set gxout grfill'

#'set cmin -75'

#'set cmax 10'

#'set cint 5'

'set clevs -65 -60 -55 -50 -45 -40 -35 -30 -25 -20 -15 -10 -5 0 5'

'set xlab off'

'set ylab off'

'd t-273'

'xcbar 7.6 7.7 7.0 9.8 -fs 2 -line on -edge circle'

### UV ####

'set arrlab off'

'set arrscl 0.5 50'

'set ccolor 0'

'set cthick 10'

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

'set arrscl 0.5 50'

'set ccolor 1'

'set cthick 2'

'set xlab on'

'set ylab on'

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

### TROPOPAUSE HEIGHT ###

i=its

while (i<=ite)

say pdatetime.i' 'tph.i

#tphkm=tph.i/1000

#say tphkm

#'q w2xy 'pdatetime.i' 'tphkm

'q w2xy 'pdatetime.i' 'tph.i

x=subwrd(result,3)

y=subwrd(result,6)

say x' 'y

say

'set line 0'

'draw mark 3 'x' 'y' '0.12

'set line 1'

'draw mark 3 'x' 'y' '0.08

i=i+1

endwhile

#

# RH

#

'set time 'stime' 'etime

'set parea 1 7.5 3.8 6.8'

'colorcb 7 b'

'set xlopts 1 3 0.15'

'set ylopts 1 3 0.15'

'set ylevs 1000 800 700 600 500 400 300 200 100'

'set xlab on'

'set ylab on'

'set gxout grfill'

levs='0 10 20 40 60 80 100'

kind='white->white->grainbow'

xcbar=' 7.6 7.7 3.8 6.6 -line on -edge circle'

'color -levs ' levs ' -kind ' kind ' -xcbar ' xcbar

'd rh'

'set strsiz 0.13 0.15'

'set string 1 c 4 90'

'draw string 0.3 8.5 P [hPa]'

'draw string 0.3 5.0 P [hPa]'

'set string 1 c 4 0'

'draw string 7.9 9.9 [\`ao\`nC]'

'draw string 7.9 6.7 [%]'

'set string 1 c 4 0'

'set strsiz 0.18 0.2'

'draw string 4.5 10.8 'title

'allclose'

# Footer

'set strsiz 0.08 0.1'

'set string 1 l 2'

'draw string 0.5 2.8 ${timestamp} ${host}'

'draw string 0.5 2.6 ${cwd} ${gs}'

'draw string 0.5 2.4 ${out}'

'gxprint ${out}'

say

say 'Output: ${out}'

say

if (quitopt="yes")then

quit

endif

EOF

opt=""

grads -bcp "$gs $opt"

exit 0

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

End of time-p_obs.run.sh

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

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

sonde.ctl

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

dset /work2/am/WRF/WRF.Post/GrADS/PL1101/Sonde.Profile.02.AM_2/output/%y2%m2%d2%h2.grd

options template

undef -999.9

title MR11 sonde data

xdef 1 linear 0.0 2.5

ydef 1 linear 60.0 2.5

zdef 183 linear 100 5

tdef 35 linear 00z18Jan2011 6hr

vars 18

p 183 99 Pressure (hPa)

Pscl 183 99 ln

t 183 99 Air Temerature (K)

rh 183 99 Relative Humidity (%)

v 183 99 y-wind component (m s-1)

u 183 99 x-wind component (m s-1)

Height 183 99 height (m)

TD 183 99 K

MR 183 99 g/kg

DD 183 99 dgr

FF 183 99 m/s

AZ 183 99 dgr

Range 183 99 m

Lon 183 99 dgr

Lat 183 99 dgr

SpuKey 183 99 bitfield

UsrKey 183 99 bitfield

RadarH 183 99 m

endvars

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

End of sonde.ctl

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