#!/bin/sh

MET_DIR=/d1/johnhg/MET/MET_development/git/MET-develop/met

############################################################################
#
# RSMCs
#
############################################################################
#
# WMO Reference:
#   https://www.wmo.int/pages/prog/www/tcp/images/TCcentres_2.JPG
#
############################################################################
#
# 1: Caribbean Sea, Gulf of Mexico, North Atlantic, Eastern North Pacific
# MIAMI
# 0 -140
# 65 -140
# 65 -20
# 0 -20
# 
# 2: Central North Pacific (Central North Pacific Ocean)
# HONOLULU
# 0 -180
# 60 -180
# 65 -170
# 65 -140
# 0 -140
# 
# 3: Western North Pacific (Western North Pacific Ocean and South China Sea)
# TOKYO
# 0 100
# 60 100
# 60 180
# 0 180
# 
# 4: India (Bay of Bengal and the Arabian Sea)
# NEW_DEHLI
# 0 45
# 60 45
# 60 100
# 0 100
# 
# 5: South-West India (South-West Indian Ocean)
# LA_REUNION
# -40 30
# 0 30
# 0 90
# -40 90
# 
# 6: Australia
# AUSTRALIA
# -36 90
# 0 90
# 0 160
# -32 160
# -32 130
# -36 130
# 
# 7: South Pacific (South-West Pacific Ocean)
# NADI
# -25 160
# 0 160
# 0 -120
# -25 -120
# 
# 8: Tasman Sea
# WELLINGTON
# -40 160
# -25 160
# -25 -120
# -40 -120
# 
############################################################################

# Create a global 0.1 degree field with all 0's
${MET_DIR}/bin/gen_vx_mask \
${MET_DIR}/share/met/tc_data/dland_global_tenth_degree.nc \
${MET_DIR}/share/met/tc_data/dland_global_tenth_degree.nc \
rsmc_global_tenth_degree.nc \
-type data -mask_field 'name="dland"; level="(*,*)";' \
-thresh ne-9999 -value 0 -name RSMC

# Apply poly masks for RSMC regions
for POLY in `ls RSMC/RSMC*.poly`; do
  cp rsmc_global_tenth_degree.nc tmp.nc
  VALUE=`basename $POLY | cut -d'_' -f2 | cut -d'.' -f1`
  ${MET_DIR}/bin/gen_vx_mask \
  tmp.nc $POLY rsmc_global_tenth_degree.nc \
  -type poly -value $VALUE -name RSMC
done

rm tmp.nc

# Plot the result
$MET_DIR/bin/plot_data_plane \
rsmc_global_tenth_degree.nc \
rsmc_global_tenth_degree.ps \
'name="RSMC"; level="(*,*)";'

convert -rotate 90 -background white -flatten \
rsmc_global_tenth_degree.ps rsmc_global_tenth_degree.png

############################################################################
#
# Basins
#
#############################################################################
#
# ATCF Format:
#   https://www.nrlmry.navy.mil/atcf_web/docs/database/new/abrdeck.html
#
# AOML Image:
#   http://www.aoml.noaa.gov/hrd/tcfaq/6RSMCs.jpg
#
############################################################################
#
# 1: Atlantic Basin
# ATLANTIC
# 0.0   -77.5
# 3.4   -76.6
# 6.3   -76.1
# 7.9   -77.3
# 8.9   -78.3
# 9.2   -79.0
# 9.1   -79.7
# 8.6   -80.3
# 8.4   -81.2
# 8.5   -82.2
# 9.0   -82.6
# 9.5   -83.8
# 13.1  -86.4
# 14.1  -87.2
# 15.5  -90.4
# 16.7  -93.5
# 16.8  -95.4
# 18.9  -99.3
# 20.7  -101.4
# 21.6  -103.2
# 26.9  -105.3
# 40.0  -105.3
# 65.0  -105.3
# 65.0   0.0
# 0.0    0.0
# 
# 2: Eastern Pacific
# EASTERN_PACIFIC
# 0.0   -140.0
# 65.0  -140.0
# 65.0  -105.3
# 40.0  -105.3
# 26.9  -105.3
# 21.6  -103.2
# 20.7  -101.4
# 18.9  -99.3
# 16.8  -95.4
# 16.7  -93.5
# 15.5  -90.4
# 14.1  -87.2
# 13.1  -86.4
# 9.5   -83.8
# 9.0   -82.6
# 8.5   -82.2
# 8.4   -81.2
# 8.6   -80.3
# 9.1   -79.7
# 9.2   -79.0
# 8.9   -78.3
# 7.9   -77.3
# 6.3   -76.1
# 3.4   -76.6
# 0.0   -77.5
# 
# 3: Central Pacific
# CENTRAL_PACIFIC
#   0 -180
#  60 -180
#  65 -170
#  65 -140
#   0 -140
#
# 4: Western Pacific
# WESTERN_PACIFIC
#   0  100
#  60  100
#  60  180
#   0  180
# 
# 5: Northern Indian Ocean
# NORTHERN_INDIAN_OCEAN
#   0   45
#  60   45
#  60  100
#   0  100
# 
# 6: Southern Indian Ocean
# SOUTHERN_INDIAN_OCEAN
# -60   20
#   0   20
#   0   90
# -60   90
# 
# 7: Australia
# AUSTRALIA
# -60   90
#   0   90
#   0  160
# -60  160
# 
# 8: South Pacific
# SOUTH_PACIFIC
# -60  160
#   0  160
#   0 -120
# -60 -120
#
############################################################################

# Create a global 0.1 degree field with all 0's
${MET_DIR}/bin/gen_vx_mask \
${MET_DIR}/share/met/tc_data/dland_global_tenth_degree.nc \
${MET_DIR}/share/met/tc_data/dland_global_tenth_degree.nc \
basin_global_tenth_degree.nc \
-type data -mask_field 'name="dland"; level="(*,*)";' \
-thresh ne-9999 -value 0 -name BASIN

# Apply poly masks for BASIN regions
for POLY in `ls BASIN/BASIN*.poly`; do
  cp basin_global_tenth_degree.nc tmp.nc
  VALUE=`basename $POLY | cut -d'_' -f2 | cut -d'.' -f1`
  ${MET_DIR}/bin/gen_vx_mask \
  tmp.nc $POLY basin_global_tenth_degree.nc \
  -type poly -value $VALUE -name BASIN
done

rm tmp.nc

# Plot the result
$MET_DIR/bin/plot_data_plane \
basin_global_tenth_degree.nc \
basin_global_tenth_degree.ps \
'name="BASIN"; level="(*,*)";'

convert -rotate 90 -background white -flatten \
basin_global_tenth_degree.ps basin_global_tenth_degree.png