#!/bin/csh

# add this directory to path
set path = ( . $path )

# mapflag global setting; mapflag=0 creates only region-level plots
# which is sufficient for most purposes, but the VERTCON3 archive was
# built using mapflag=2 for all zones except Alaska, so you might want to
# change that value if you need more detailed plots:
#set mapflag = 0 
set mapflag = 2

echo -------------------------------------------------------------------
echo finalRuns - BEGIN
echo -------------------------------------------------------------------

# possible combinations (list taken from 'doit.bat')
#   For CONUS:
#                 ngvd29 navd88 conus #
#   For ALASKA:
#                 ngvd29 navd88 alaska #
#   For HAWAII:
#                 No combinations exist
#   For PR:
#                 lt prvd02 pr #
#   For VI:
#                 lt vivd09 vi #
#   For GUAM:
#                 guvd63 guvd04 guam #
#   For CNMI:
#                 lt nmvd03 cnmi #
#   For AS:
#                 lt asvd02 as #


# NOTE: for CNMI and VI areas, we decided on a manual approach
# where, for each individual island with data, 
#  transformation = mean
#  error = std
# Therefore these areas will not use the doit2/doit3 for transformation generation
#
# The grid generation process for these areas is as follows
# 1) doit/doit2 to get coverage plots for report, and run filters on data
# 2) finalize.* for that region, to build flat-plane transformations
#     and splat the transformations into a grid

echo -------------------------------------------------------------------
echo finalRuns - VI
echo -------------------------------------------------------------------

set frame1 = lt
set frame2 = vivd09
set region = vi
set filter1 = 0

doit.bat ${frame1} ${frame2} ${region} ${filter1}
doit2.bat ${frame1} ${frame2} ${region} ${filter1} ${mapflag}
finalize.vi


echo -------------------------------------------------------------------
echo finalRuns - CNMI
echo -------------------------------------------------------------------

set frame1 = lt
set frame2 = nmvd03
set region = cnmi
set filter1 = 0

doit.bat ${frame1} ${frame2} ${region} ${filter1}
doit2.bat ${frame1} ${frame2} ${region} ${filter1} ${mapflag}
finalize.cnmi



# ====================================================================

# the following areas use gridded transformations (in part or all of the region),
# so need to run all the doit's

# for all of the island regions (pr, guam, as), the "finalize" step
# will take care of splatting the transformation into a grid with large
# "no transformation available" areas


echo -------------------------------------------------------------------
echo finalRuns - GUAM
echo -------------------------------------------------------------------

# GUAM
set frame1 = guvd63
set frame2 = guvd04
set region = guam
set filter1 = 0
set rz = 60
set cod = $rz

doit.bat ${frame1} ${frame2} ${region} ${filter1}
doit2.bat ${frame1} ${frame2} ${region} ${filter1} ${mapflag}
doit3.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}
doit4.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}
finalize.guam


echo -------------------------------------------------------------------
echo finalRuns - AS
echo -------------------------------------------------------------------

set frame1 = lt
set frame2 = asvd02
set region = as
set filter1 = 3
set rz = 120
set cod = $rz

doit.bat ${frame1} ${frame2} ${region} ${filter1}
doit2.bat ${frame1} ${frame2} ${region} ${filter1} ${mapflag}
doit3.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}
doit4.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}
finalize.as


echo -------------------------------------------------------------------
echo finalRuns - PR
echo -------------------------------------------------------------------

# PR PUERTO RICO
set frame1 = lt
set frame2 = prvd02
set region = pr
set filter1 = 0
set rz = 60
set cod = $rz

doit.bat ${frame1} ${frame2} ${region} ${filter1}
doit2.bat ${frame1} ${frame2} ${region} ${filter1} ${mapflag}
doit3.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}
doit4.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}
finalize.pr



# alaska and conus have complete grids, so there is not a "finalize" step

echo -------------------------------------------------------------------
echo finalRuns - CONUS
echo -------------------------------------------------------------------

# CONUS
set frame1 = ngvd29
set frame2 = navd88
set region = conus
set filter1 = 3
set rz = 300
set cod = 180

doit.bat ${frame1} ${frame2} ${region} ${filter1}
doit2.bat ${frame1} ${frame2} ${region} ${filter1} ${mapflag}
doit3.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}
doit4.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}


echo -------------------------------------------------------------------
echo finalRuns - ALASKA
echo -------------------------------------------------------------------

# ALASKA
set frame1 = ngvd29
set frame2 = navd88
set region = alaska
set filter1 = 3
set mapflag = 1
set rz = 60
set cod = $rz

doit.bat ${frame1} ${frame2} ${region} ${filter1}
doit2.bat ${frame1} ${frame2} ${region} ${filter1} ${mapflag}
doit3.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}
doit4.bat ${frame1} ${frame2} ${region} ${filter1} $rz $cod ${mapflag}


echo -------------------------------------------------------------------
echo finalRuns - END
echo -------------------------------------------------------------------


exit
exit