2018-06-19_10-59
/work2/am/2017.LLC.Pacific/HAGIWARA_MASK/MAP.HISTO.CLOUD.V.BIN
am@localhost
$ srcdump.sh MAP.NUM.DAT.V.BIN.PL.BOX.sh
#------------------------------
# List of the following files:
#------------------------------
MAP.NUM.DAT.V.BIN.PL.BOX.sh
#------------------------------
# Machine info
#------------------------------
localhost
/work2/am/2017.LLC.Pacific/HAGIWARA_MASK/MAP.HISTO.CLOUD.V.BIN
Tue Jun 19 11:03:28 JST 2018
#======================
# MAP.NUM.DAT.V.BIN.PL.BOX.sh
#======================
#!/bin/bash
# Description:
#
# Author: am
#
# Host: localhost
# Directory: /work2/am/2017.LLC.Pacific/HAGIWARA_MASK/MAP.HISTO.CLOUD
#
# Revision history:
# This file is created by /usr/local/mybin/ngmt.sh at 11:51 on 06-07-2018.
. ./gmtpar.sh
echo "Bash script $0 starts."
range=140/180/30/50
size=Q130/5
xanot=a10f5g5
yanot=a5f5g5
indir_root=.
indir=${indir_root}/OUTPUT_MAP.HISTO.CLOUD.V.BIN
sdate=20160606
edate=20160625
wlon=140
elon=180
slat=030
nlat=050
inlist=$(ls -1 ${indir}/BOX_DAY_${sdate}-${edate}_${wlon}-${elon}_${slat}-${nlat}*)
figdir=Fig_MAP.HISTO.CLOUD.V.BIN/${sdate}-${edate}_${wlon}-${elon}_${slat}-${nlat}
mkdir -vp $figdir
cpt=$(basename $0 .sh).cpt
makecpt -Cseis -I -T1000/3000/500 > $cpt
dirsstave=/work2/am/2017.LLC.Pacific/HAGIWARA_MASK/OISST.MONTH.AVE.STD
sdate=20160605
insstave=OISST.MONTH.AVE.STD_${sdate}-${edate}_${wlon}-${elon}_${slat}-${nlat}.txt
insst=${dirsstave}/${insstave}
if [ ! -f $insst ]; then
echo Error in $0 : No such file, $insstave
exit 1
fi
grdsst=$(basename $insst .txt).AVE.grd
reso=0.25
awk '{if($1!="#"&&$3>0)print $1,$2,$3}' $insst |\
surface -R$range -I$reso/$reso -T1 -G$grdsst
# ▒▒▒n▒̃}▒X▒N
landmask=landmask.grd
maskedgrd=masked.grd
grdlandmask -R$range -Dl -I$reso -N1/NaN -G$landmask #-V
grdmath $grdsst $landmask OR = $maskedgrd
sdate=20160606
# SOUNDINGS
indir3=/work2/am/2017.LLC.Pacific/ERA-I.SYNOPTIC
infle3=ERA-I.SYNOPTIC.siteinfo.txt
in3=${indir3}/${infle3}
if [ ! -f $in3 ]; then
echo ERROR in $0
echo NO SUCH FILE, $infle3
echo IN $indir3
exit 1
fi
for in in $inlist; do
if [ ! -f $in ]; then
echo Error in $0 : No such file, $in
exit 1
fi
in2=$(echo $in | sed -e "s/DAY/NIGHT/")
out=${figdir}/NUM_$(echo $(basename $in .txt) | sed -e "s/DAY_//").ps
title="DAY"
anot=${xanot}/${yanot}:.$title:WSne
awk '{ if($1==">") print $1,$2,$5; else if($1!="#" && $1!=">") print $0 }' $in|\
psxy -R$range -J$size -C$cpt -L -M -A -K -X1.5 -Y6 -P >$out
grdcontour $maskedgrd -R$range -J$size -C1 -W10/255/255/255 -O -K>>$out
grdcontour $maskedgrd -R$range -J$size -W1 -C1 -O -K>>$out
grdcontour $maskedgrd -R$range -J$size -W3 -A5 -O -K>>$out
# SOUNDING
awk '{if($1!="#") print $8+$9/60., $6+$7/60.,\
" 10 0 0 LT ", substr($3,9,2)}' $in3|\
pstext -D0.02/0.02 -R$range -J$size -O -K >>$out
awk '{if($1!="#") print $8+$9/60.,$6+$7/60.,substr($3,9,2)}' $in3|\
psxy -R$range -J$size -Sa0.1 -W10/255/255/255 -K -O >>$out
awk '{if($1!="#") print $8+$9/60.,$6+$7/60.,substr($3,9,2)}' $in3|\
psxy -R$range -J$size -Sa0.1 -W3 -K -O >>$out
pscoast -R$range -J$size -W10/255/255/255 -O -K >>$out
pscoast -R$range -J$size -B$anot -W3 -O -K >>$out
#psscale -D2.5/-0.5/5/0.1h -C$cpt -E -Ba5f5 -O -K >>$out
title="NIGHT"
anot=${xanot}/${yanot}:.$title:WSne
awk '{ if($1==">") print $1,$2,$5; else if($1!="#" && $1!=">") print $0 }' $in2|\
psxy -R$range -J$size -C$cpt -L -M -A -K -Y-3.5 -O >>$out
grdcontour $maskedgrd -R$range -J$size -C1 -W10/255/255/255 -O -K>>$out
grdcontour $maskedgrd -R$range -J$size -W1 -C1 -O -K>>$out
grdcontour $maskedgrd -R$range -J$size -W3 -A5 -O -K>>$out
# SOUNDING
awk '{if($1!="#") print $8+$9/60., $6+$7/60.,\
" 10 0 0 LT ", substr($3,9,2)}' $in3|\
pstext -D0.02/0.02 -R$range -J$size -O -K >>$out
awk '{if($1!="#") print $8+$9/60.,$6+$7/60.,substr($3,9,2)}' $in3|\
psxy -R$range -J$size -Sa0.1 -W10/255/255/255 -K -O >>$out
awk '{if($1!="#") print $8+$9/60.,$6+$7/60.,substr($3,9,2)}' $in3|\
psxy -R$range -J$size -Sa0.1 -W3 -K -O >>$out
pscoast -R$range -J$size -W10/255/255/255 -O -K >>$out
pscoast -R$range -J$size -B$anot -W3 -O -K >>$out
psscale -D2.5/-0.5/5/0.1h -C$cpt -E -Ba500f500 -O -K >>$out
xoffset=
yoffset=6.25
export LANG=C
curdir1=$(pwd)
now=$(date)
host=$(hostname)
time=$(ls -l ${in} | awk '{print $6, $7, $8}')
time2=$(ls -l ${in2} | awk '{print $6, $7, $8}')
infle=$(basename $in)
infle2=$(basename $in2)
pstext <<EOF -JX6/1.5 -R0/1/0/1.5 -N -X${xoffset:-0} -Y${yoffset:-0} -O >> $out
0 1.50 9 0 1 LM $0 $@
0 1.35 9 0 1 LM ${now}
0 1.20 9 0 1 LM ${host}
0 1.05 9 0 1 LM ${curdir1}
0 0.90 9 0 1 LM Indir: ${indir}
0 0.75 9 0 1 LM Input: ${infle} (${time})
0 0.60 9 0 1 LM Input: ${infle2} (${time2})
0 0.45 9 0 1 LM Input:
EOF
echo "OUTPUT: $out"
done
#rm -fv $cpt
echo "Done $0"
#----------------------
# End of MAP.NUM.DAT.V.BIN.PL.BOX.sh
#----------------------