# VERTCON 3 build script for region:  CNMI
# See DRU-12, p. 143-145
#
# This batch file is to be run AFTER the following scripts have already been run:
#      doit.bat
#      doit2.bat
#      (You do NOT run doit3.bat nor doit4.bat for the CNMI region)
# It calls various standard "grid manipulation" programs to 
# create "flat plane" transformations, and "flat plane" error grids, 
# surrounding the three main CNMI islands (the northern islands get none)
# Anything outside of the three main islands gets flag "-999.0"
#
# All code used is listed below:
#  /home/dru/Linux/Mylib/gfill    (Source:  /home/dru/Linux/Mylib/Source/gfill.f)
#  /home/dru/Linux/Mylib/gpull    (Source:  /home/dru/Linux/Mylib/Source/gpull.f)
#  /home/dru/Linux/Mylib/gsplat   (Source:  /home/dru/Linux/Mylib/Source/gsplat.f)
#                          
##########################
# Transformation:
##########################
#
# Step 1:  Fill in the full grid area with flag -999.d0 at 1' x 1' spacing
#
gfill << !
1
143  
147
241
13.8
22
493
cnmiflag.b
-999.d0
!
#
# Step 2: Create "flat plane" (constant) transformation for Rota
#
gfill << !
1
145.1
145.3
13
14.1
14.21666666666666666666666666666666666666666666666666666666
8
cnmi01.b
-0.0325d0
!
#
# Step 3: Create "flat plane" (constant) transformation for Tinian
#
gfill << !
1
145.5666666666666666666666666666666666666666666
145.6833333333333333333333333333333333333333333
8 
14.9
15.11666666666666666666666666666666666666666666
14
cnmi02.b
0.0360d0
!
#
# Step 4: Create "flat plane" (constant) transformation for Saipan
#
gfill << !
1
145.6833333333333333333333333333333333333333333
145.85
11
15.08333333333333333333333333333333333333333333
15.3
14
cnmi03.b
-0.3050d0
!
#
# Step 5: Drop (splat) the Rota part into the full size flag-filled grid
#
gsplat << !
cnmiflag.b
cnmi01.b
cnmi04.b
!
#
# Step 6: Drop (splat) the Tinian part into the grid from step 5
#
gsplat << !
cnmi04.b
cnmi02.b
cnmi05.b
!
#
# Step 7: Drop (splat) the Saipan part into the grid from step 6
#
gsplat << !
cnmi05.b
cnmi03.b
lt.nmvd03.cnmi.trn.b
!
##########################
# Error Grid (uses same "flag" grid as transformation):
##########################
#
# Step 1: Create "flat plane" (constant) error for Rota
#
gfill << !
1
145.1
145.3
13
14.1
14.21666666666666666666666666666666666666666666666666666666
8
cnmi01e.b
0.0025d0
!
#
# Step 2: Create "flat plane" (constant) error for Tinian
#
gfill << !
1
145.5666666666666666666666666666666666666666666
145.6833333333333333333333333333333333333333333
8 
14.9
15.11666666666666666666666666666666666666666666
14
cnmi02e.b
0.0000d0
!
#
# Step 3: Create "flat plane" (constant) error for Saipan
#
gfill << !
1
145.6833333333333333333333333333333333333333333
145.85
11
15.08333333333333333333333333333333333333333333
15.3
14
cnmi03e.b
0.0000d0
!
#
# Step 4: Drop (splat) the Rota part into the full size flag-filled grid
#
gsplat << !
cnmiflag.b
cnmi01e.b
cnmi04e.b
!
#
# Step 5: Drop (splat) the Tinian part into the grid from step 4
#
gsplat << !
cnmi04e.b
cnmi02e.b
cnmi05e.b
!
#
# Step 6: Drop (splat) the Saipan part into the grid from step 5
#
gsplat << !
cnmi05e.b
cnmi03e.b
lt.nmvd03.cnmi.err.b
!
##########################
# Clean Up
##########################
#
# Step 1: Remove preliminary files
#
rm -f cnmi01.b
rm -f cnmi02.b
rm -f cnmi03.b
rm -f cnmi04.b
rm -f cnmi05.b
rm -f cnmi01e.b
rm -f cnmi02e.b
rm -f cnmi03e.b
rm -f cnmi04e.b
rm -f cnmi05e.b
rm -f cnmiflag.b