# VERTCON 3 build script for region:  PR
# 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
#      doit3.bat
#      doit4.bat
# It calls various standard "grid manipulation" programs to cut out the "pulled in tension" transformation
# and error grids around the island of Puerto Rico itself, and fill in the rest of the region with 
# 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
291
294.9
235
17
19
121
prflag.b
-999.d0
!
#
# Step 2: Pull PART of the actual transformation grid out of the full "04.b" grid (western half of island of PR)
#
gpull << !
vmtcdoht.lt.prvd02.pr.000.60.60.04.b
pr01.b
17.9
18.55
292.7
294.4
!
#
# Step 3: Pull PART of the actual transformation grid out of the full "04.b" grid (eastern half of island of PR, avoiding Viequez)
#
gpull << !
vmtcdoht.lt.prvd02.pr.000.60.60.04.b
pr02.b
18.2
18.55
294.4
294.5
!
#
# Step 4: Drop (splat) the western part of island PR into the full size flag-filled grid
#
gsplat << !
prflag.b
pr01.b
pr03.b
!
#
# Step 5: Drop (splat) the eastern part of island PR into the grid from step 4
#
gsplat << !
pr03.b
pr02.b
lt.prvd02.pr.trn.b
!
##########################
# Error Grid (uses same "flag" grid as transformation):
##########################
#
# Step 1: Pull PART of the actual error grid out of the full "ete*b" grid (western half of island of PR)
#
gpull << !
vmeteoht.lt.prvd02.pr.000.60.60.b
pr01e.b
17.9
18.55
292.7
294.4
!
#
# Step 2: Pull PART of the actual error grid out of the full "04.b" grid (eastern half of island of PR, avoiding Viequez)
#
gpull << !
vmeteoht.lt.prvd02.pr.000.60.60.b
pr02e.b
18.2
18.55
294.4
294.5
!
#
# Step 3: Drop (splat) the western part of island PR into the full size flag-filled grid
#
gsplat << !
prflag.b
pr01e.b
pr03e.b
!
#
# Step 4: Drop (splat) the eastern part of island PR into the grid from step 4
#
gsplat << !
pr03e.b
pr02e.b
lt.prvd02.pr.err.b
!
##########################
# Clean Up
##########################
#
# Step 1: Remove preliminary files
#
rm -f pr01.b
rm -f pr02.b
rm -f pr03.b
rm -f pr01e.b
rm -f pr02e.b
rm -f pr03e.b
rm -f prflag.b