------------------------------
List of the following files:
------------------------------
time-height_obs.run.sh
sonde.km.ctl
------------------------------
Machine info
------------------------------
aofd165.bio.mie-u.ac.jp
/work2/am/WRF/WRF.Post/GrADS/PL1101/Time-Height_OBS
Tue Sep 19 14:09:14 JST 2017
======================
time-height_obs.run.sh
======================
#!/bin/bash
gs=$(basename $0 .run.sh).gs
if [ ! -f $gs ]; then
echo ERROR in $0: No such file, $gs
exit 1
fi
opt=""
cat<<EOF>$gs
function timehgt( args )
#
# Default values
#
ctl='./sonde.km.ctl'
outdir='Fig'
lev1=0 ;#100
lev2=12 ;#5000
lev1f2=lev1/1000 ;* m -> km
lev2f2=lev2/1000 ;* m -> km
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
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,3)
i=i+1
endwhile
### SET-UP ###
'set time 'stime' 'etime
'set lev 'lev1' 'lev2
'cc'
'set grads off'
'set grid off'
'set parea 1 7.8 7 10'
'set xlopts 1 0 0'
### T ####
'set gxout grfill'
#'set cmin -35'
#'set cmax 5'
#'set cint 5'
levs='-70 -65 -60 -55 -50 -45 -40 -35 -30 -25 -20 -15 -10 -5 0'
kind='white->grainbow'
xcbar='7.9 8.0 7.0 9.8 -fs 2 -line on -edge circle'
'color -levs ' levs ' -kind ' kind ' -xcbar ' xcbar
'set xlab off'
'set ylab off'
'd t'
### UV ####
'set arrlab off'
'set arrscl 0.5 50'
'set ccolor 0'
'set cthick 10'
'd skip(u,5,1);v'
'set arrscl 0.5 50'
'set ccolor 1'
'set cthick 2'
'set xlab on'
'set ylab on'
'vec.gs skip(u,5,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
x=subwrd(result,3)
y=subwrd(result,6)
say x' 'y
'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.8 3.8 6.8'
'colorcb 7 b'
'set xlopts 1 3 0.12'
'set xlab on'
'set ylab on'
'set gxout grfill'
levs='0 10 20 40 60 80 100'
kind='white->white->grainbow'
xcbar=' 7.9 8.0 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 Height [km]'
'draw string 0.3 5.0 Height [km]'
'set string 1 c 4 0'
'draw string 8.2 9.9 [\`ao\`nC]'
'draw string 8.2 6.7 [%]'
'set string 1 c 4 0'
'set strsiz 0.18 0.23'
'draw string 4.5 10.7 'title
'allclose'
out='PL1101_Time-Height_OBS.km.eps'
'gxprint 'out
say
say 'Output: 'out
say
if (quitopt="yes")then
quit
endif
EOF
grads -bcp "$gs $opt"
exit 0
----------------------
End of time-height_obs.run.sh
----------------------
======================
sonde.km.ctl
======================
dset /work2/am/2016.PolarLow/Sonde.Profile.03.Height/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 125 linear 0.1 0.1
tdef 34 linear 00z18Jan2011 6hr
vars 18
Height 125 99 height (m)
Pscl 125 99 Air Temprature (K)
t 125 99 Air Temerature (degC)
rh 125 99 Relative Humidity (%)
v 125 99 y-wind component (m s-1)
u 125 99 x-wind component (m s-1)
p 125 99 Pressure (hPa)
TD 125 99 K
MR 125 99 g/kg
DD 125 99 dgr
FF 125 99 m/s
AZ 125 99 dgr
Range 125 99 m
Lon 125 99 dgr
Lat 125 99 dgr
SpuKey 125 99 bitfield
UsrKey 125 99 bitfield
RadarH 125 99 m
endvars
----------------------
End of sonde.km.ctl
----------------------