------------------------------
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
----------------------