# Gridgen Journal File to generate mesh for modified NACA4415
# Developed by Omar Lopez
# Thesis title: "Computational study of a NACA4415 airfoil using synthetic jet control"
# Supervisor: Robert D. Moser
# Department of Mechanical Engineering.
# University of Texas at Austin. 2009
package require PWI_Glyph 1.6.9

# Delete any existing grids and database entities.  Reset AS/W, defaults, and
# tolerances.
gg::memClear
gg::aswDeleteBC -glob "*"
gg::aswDeleteVC -glob "*"
gg::aswSet GENERIC -dim 3
gg::defReset
gg::tolReset
# Delay screen updates and checking for user input until script is finished.
gg::updatePolicy DELAYED

#Read the grid file that contains the modified NACA4415 profile. 
#Make sure the file is in the correct path.
set _ggTemp_(1) [gg::conImport "./modifiedNACA4415.grd"]
set _CN(1) [lindex $_ggTemp_(1) 0]
set _CN(2) [lindex $_ggTemp_(1) 1]
set _CN(3) [lindex $_ggTemp_(1) 2]
set _CN(4) [lindex $_ggTemp_(1) 3]
set _CN(5) [lindex $_ggTemp_(1) 4]
set _CN(6) [lindex $_ggTemp_(1) 5]
set _CN(7) [lindex $_ggTemp_(1) 6]
set _CN(8) [lindex $_ggTemp_(1) 7]
set _CN(9) [lindex $_ggTemp_(1) 8]
set _CN(10) [lindex $_ggTemp_(1) 9]
set _CN(11) [lindex $_ggTemp_(1) 10]
set _CN(12) [lindex $_ggTemp_(1) 11]
set _CN(13) [lindex $_ggTemp_(1) 12]
set _CN(14) [lindex $_ggTemp_(1) 13]
set _CN(15) [lindex $_ggTemp_(1) 14]
set _CN(16) [lindex $_ggTemp_(1) 15]
set _CN(17) [lindex $_ggTemp_(1) 16]
set _CN(18) [lindex $_ggTemp_(1) 17]
set _CN(19) [lindex $_ggTemp_(1) 18]
set _CN(20) [lindex $_ggTemp_(1) 19]
set _CN(21) [lindex $_ggTemp_(1) 20]
set _CN(22) [lindex $_ggTemp_(1) 21]
set _CN(23) [lindex $_ggTemp_(1) 22]
set _CN(24) [lindex $_ggTemp_(1) 23]
unset _ggTemp_(1)
#########################################################

#Create other connectors needed for the RANS region
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(6) -arc 0]
    gg::segAddControlPt [gg::conGetPt $_CN(5) -arc 0]
  gg::segEnd
set _CN(25) [gg::conEnd]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(6) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(5) -arc 1]
  gg::segEnd
set _CN(26) [gg::conEnd]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(8) -arc 0]
    gg::segAddControlPt [gg::conGetPt $_CN(7) -arc 0]
  gg::segEnd
set _CN(27) [gg::conEnd]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(8) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(7) -arc 1]
  gg::segEnd
set _CN(28) [gg::conEnd]
gg::dispSmallText FALSE
set _ggTemp_(1) [list 1.107 -7.52e-4 0]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(23) -arc 1]
    gg::segAddControlPt $_ggTemp_(1)
  gg::segEnd
set _CN(29) [gg::conEnd]
unset _ggTemp_(1)
set _ggTemp_(2) [list -0.1 0 0]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(1) -arc 1]
    gg::segAddControlPt $_ggTemp_(2)
  gg::segEnd
set _CN(30) [gg::conEnd]
unset _ggTemp_(2)
set _ggTemp_(1) [list 0.5 0.2638 0]
gg::conBegin
  gg::segBegin -type CONIC -rho 0.01
    gg::segAddControlPt [gg::conGetPt $_CN(30) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(29) -arc 1]
    gg::segAddControlPt $_ggTemp_(1)
  gg::segEnd
set _CN(31) [gg::conEnd]
unset _ggTemp_(1)
set _ggTemp_(1) [list 0.5 -0.2003 0]
gg::conBegin
  gg::segBegin -type CONIC -rho 0.01
    gg::segAddControlPt [gg::conGetPt $_CN(29) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(30) -arc 1]
    gg::segAddControlPt $_ggTemp_(1)
  gg::segEnd
set _CN(32) [gg::conEnd]
unset _ggTemp_(1)
set _CN(33) [gg::conSplit $_CN(31) [gg::conGetPt $_CN(31) -x 0.8356]]
set _CN(34) [gg::conSplit $_CN(33) [gg::conGetPt $_CN(33) -x 1.004]]
set _CN(35) [gg::conSplit $_CN(32) [gg::conGetPt $_CN(32) -x 0.8574]]
set _CN(36) [gg::conSplit $_CN(32) [gg::conGetPt $_CN(32) -x 0.9721]]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(1) -arc 0]
    gg::segAddControlPt [gg::conGetPt $_CN(31) -arc 1]
  gg::segEnd
set _CN(47) [gg::conEnd]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(4) -arc 0]
    gg::segAddControlPt [gg::conGetPt $_CN(33) -arc 1]
  gg::segEnd
set _CN(48) [gg::conEnd]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(3) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(35) -arc 0]
  gg::segEnd
set _CN(49) [gg::conEnd]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(10) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(32) -arc 1]
  gg::segEnd
set _CN(50) [gg::conEnd]
gg::conBegin
  gg::segBegin -type CONIC -rho 0.5
    gg::segAddControlPt [gg::conGetPt $_CN(10) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(4) -arc 0]
    gg::segAddControlPt [gg::conGetPt $_CN(23) -arc 1]
  gg::segEnd
set _CN(51) [gg::conEnd]
set _CN(52) [gg::conSplit $_CN(51) [gg::conGetPt $_CN(51) -arc 0.49249061616]]

#Create connectors needed for the LES and EULER regions
set _ggTemp_(1) [list 1 2.5 0]
set _ggTemp_(2) [list 6 2.5 0]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt $_ggTemp_(1)
    gg::segAddControlPt $_ggTemp_(2)
  gg::segEnd
set _CN(37) [gg::conEnd]
unset _ggTemp_(1)
unset _ggTemp_(2)
set _ggTemp_(3) [list 6 2 0]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(37) -arc 1]
    gg::segAddControlPt $_ggTemp_(3)
  gg::segEnd
set _CN(38) [gg::conEnd]
unset _ggTemp_(3)
set _ggTemp_(4) [list 6 -2 0]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(38) -arc 1]
    gg::segAddControlPt $_ggTemp_(4)
  gg::segEnd
set _CN(39) [gg::conEnd]
unset _ggTemp_(4)
set _ggTemp_(5) [list 6 -2.5 0]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(39) -arc 1]
    gg::segAddControlPt $_ggTemp_(5)
  gg::segEnd
set _CN(40) [gg::conEnd]
unset _ggTemp_(5)
set _ggTemp_(6) [list 1 -2.5 0]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(40) -arc 1]
    gg::segAddControlPt $_ggTemp_(6)
  gg::segEnd
set _CN(41) [gg::conEnd]
unset _ggTemp_(6)
set _ggTemp_(7) [list -1.5 0 0]
gg::conBegin
  gg::segBegin -type CIRCULAR_ARC
    gg::segAddControlPt [gg::conGetPt $_CN(41) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(37) -arc 0]
    gg::segAddControlPt $_ggTemp_(7)
  gg::segEnd
set _CN(42) [gg::conEnd]
unset _ggTemp_(7)
set _ggTemp_(1) [list 4 2 0]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(38) -arc 1]
    gg::segAddControlPt $_ggTemp_(1)
  gg::segEnd
set _CN(43) [gg::conEnd]
unset _ggTemp_(1)
set _ggTemp_(2) [list 4 -2 0]
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(43) -arc 1]
    gg::segAddControlPt $_ggTemp_(2)
  gg::segEnd
set _CN(44) [gg::conEnd]
unset _ggTemp_(2)
gg::conBegin
  gg::segBegin -type 3D_LINE
    gg::segAddControlPt [gg::conGetPt $_CN(44) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(39) -arc 1]
  gg::segEnd
set _CN(45) [gg::conEnd]
set _ggTemp_(3) [list -0.172 0 0]
gg::conBegin
  gg::segBegin -type CONIC -rho 0.5
    gg::segAddControlPt [gg::conGetPt $_CN(43) -arc 1]
    gg::segAddControlPt [gg::conGetPt $_CN(44) -arc 1]
    gg::segAddControlPt $_ggTemp_(3)
  gg::segEnd
set _CN(46) [gg::conEnd]
unset _ggTemp_(3)
######################################################################

#Set the number of nodes in the connectors
gg::conRedimBegin
  gg::conRedim $_CN(4) 50
  gg::conRedim $_CN(12) 20
  gg::conRedim $_CN(6) 25
  gg::conRedim $_CN(19) 25
  gg::conRedim $_CN(15) 29
  gg::conRedim $_CN(18) 25
  gg::conRedim $_CN(14) 15
  gg::conRedim $_CN(26) 15
  gg::conRedim $_CN(5) 25
  gg::conRedim $_CN(25) 15
  gg::conRedim $_CN(13) 10
  gg::conRedim $_CN(11) 50
  gg::conRedim $_CN(23) 43
  gg::conRedim $_CN(24) 43
  gg::conRedim $_CN(2) 50
  gg::conRedim $_CN(9) 10
  gg::conRedim $_CN(27) 15
  gg::conRedim $_CN(10) 20
  gg::conRedim $_CN(7) 25
  gg::conRedim $_CN(17) 15
  gg::conRedim $_CN(20) 25
  gg::conRedim $_CN(16) 29
  gg::conRedim $_CN(21) 25
  gg::conRedim $_CN(28) 15
  gg::conRedim $_CN(8) 25
  gg::conRedim $_CN(22) 50
  gg::conRedim $_CN(49) 75
  gg::conRedim $_CN(50) 75
  gg::conRedim $_CN(29) 75
  gg::conRedim $_CN(48) 75
  gg::conRedim $_CN(47) 75
  gg::conRedim $_CN(30) 75
  gg::conRedim $_CN(34) 50
  gg::conRedim $_CN(33) 50
  gg::conRedim $_CN(31) 80
  gg::conRedim $_CN(35) 80
  gg::conRedim $_CN(36) 50
  gg::conRedim $_CN(32) 50
  gg::conRedim $_CN(46) 1024
  gg::conRedim $_CN(44) 483
  gg::conRedim $_CN(43) 50
  gg::conRedim $_CN(45) 50
  gg::conRedim $_CN(39) 42
  gg::conRedim $_CN(38) 7
  gg::conRedim $_CN(40) 7
  gg::conRedim $_CN(41) 52
  gg::conRedim $_CN(37) 52
  gg::conRedim $_CN(42) 80
  gg::conRedim $_CN(51) 50
  gg::conRedim $_CN(52) 50
gg::conRedimEnd
####################################################################

#Set the distribution of the nodes in each connector
#5e-5 corresponds to the first cell heigth in the SS of the airfoil.
gg::conBeginSpacing $_CN(4) -sub 1 5e-5
gg::conBeginSpacing $_CN(12) -sub 1 1e-4
gg::conEndSpacing $_CN(12) -sub 1 5e-5
gg::conBeginSpacing $_CN(6) -sub 1 5e-5
gg::conEndSpacing $_CN(6) -sub 1 5e-5
gg::dispSmallText FALSE
gg::conBeginSpacing $_CN(5) -sub 1 5e-5
gg::conEndSpacing $_CN(5) -sub 1 5e-5
gg::conEndSpacing $_CN(19) -sub 1 5e-5
gg::conEndSpacing $_CN(18) -sub 1 5e-5
gg::conBeginSpacing $_CN(15) -sub 1 5e-5
gg::conEndSpacing $_CN(15) -sub 1 5e-5
gg::conBeginSpacing $_CN(14) -sub 1 5e-5
gg::conEndSpacing $_CN(14) -sub 1 5e-5
gg::conBeginSpacing $_CN(26) -sub 1 5e-5
gg::conEndSpacing $_CN(26) -sub 1 5e-5
gg::dispSmallText FALSE
gg::conBeginSpacing $_CN(25) -sub 1 5e-5
gg::conEndSpacing $_CN(25) -sub 1 5e-5
gg::conBeginSpacing $_CN(13) -sub 1 5e-5
gg::conEndSpacing $_CN(13) -sub 1 5e-5
gg::conBeginSpacing $_CN(11) -sub 1 5e-5
gg::conEndSpacing $_CN(11) -sub 1 5e-5
gg::conBeginSpacing $_CN(52) -sub 1 5e-5
gg::conEndSpacing $_CN(52) -sub 1 5e-5
gg::conBeginSpacing $_CN(23) -sub 1 5e-5
gg::conEndSpacing $_CN(23) -sub 1 1e-3
gg::conBeginSpacing $_CN(24) -sub 1 5e-5
gg::conEndSpacing $_CN(24) -sub 1 1e-3
gg::conBeginSpacing $_CN(2) -sub 1 5e-5
gg::conEndSpacing $_CN(2) -sub 1 5e-5
gg::conBeginSpacing $_CN(51) -sub 1 5e-5
gg::conEndSpacing $_CN(51) -sub 1 5e-5
gg::conBeginSpacing $_CN(9) -sub 1 5e-5
gg::conEndSpacing $_CN(9) -sub 1 5e-5
gg::conBeginSpacing $_CN(7) -sub 1 5e-5
gg::conEndSpacing $_CN(7) -sub 1 5e-5
gg::conBeginSpacing $_CN(17) -sub 1 5e-5
gg::conEndSpacing $_CN(17) -sub 1 5e-5
gg::conBeginSpacing $_CN(20) -sub 1 5e-5
gg::conEndSpacing $_CN(21) -sub 1 5e-5
gg::conBeginSpacing $_CN(16) -sub 1 5e-5
gg::conEndSpacing $_CN(16) -sub 1 5e-5
gg::conBeginSpacing $_CN(8) -sub 1 5e-5
gg::conEndSpacing $_CN(8) -sub 1 5e-5
gg::conBeginSpacing $_CN(28) -sub 1 5e-5
gg::conEndSpacing $_CN(28) -sub 1 5e-5
gg::conBeginSpacing $_CN(27) -sub 1 5e-5
gg::conEndSpacing $_CN(27) -sub 1 5e-5
gg::conBeginSpacing $_CN(10) -sub 1 1e-4

#2e-5 corresponds to the first cell height in the PS of the airfoil.
gg::conEndSpacing $_CN(10) -sub 1 2e-5
gg::conEndSpacing $_CN(22) -sub 1 5e-5
gg::conBeginSpacing $_CN(29) -sub 1 7e-4
gg::conEndSpacing $_CN(29) -sub 1 7e-3
gg::conBeginSpacing $_CN(48) -sub 1 5e-5
gg::conEndSpacing $_CN(48) -sub 1 7e-3
gg::conBeginSpacing $_CN(47) -sub 1 5e-5
gg::conEndSpacing $_CN(47) -sub 1 7e-3
gg::conBeginSpacing $_CN(30) -sub 1 5e-5
gg::conEndSpacing $_CN(30) -sub 1 7e-3
gg::conBeginSpacing $_CN(49) -sub 1 2e-5
gg::conEndSpacing $_CN(49) -sub 1 7e-3
gg::conBeginSpacing $_CN(50) -sub 1 2e-5
gg::conEndSpacing $_CN(50) -sub 1 7e-3
gg::conBeginSpacing $_CN(31) -sub 1 5e-3
gg::conEndSpacing $_CN(31) -sub 1 1.5e-2
gg::conBeginSpacing $_CN(33) -sub 1 1.5e-2
gg::conEndSpacing $_CN(33) -sub 1 1e-3
gg::conBeginSpacing $_CN(34) -sub 1 1e-3
gg::conEndSpacing $_CN(34) -sub 1 2.6e-3
gg::conEndSpacing $_CN(32) -sub 1 1e-3
gg::conEndSpacing $_CN(32) -sub 1 1.5e-3
gg::conBeginSpacing $_CN(32) -sub 1 3.935e-3
gg::conEndSpacing $_CN(36) -sub 1 7e-3
gg::conBeginSpacing $_CN(36) -sub 1 1e-3
gg::conEndSpacing $_CN(35) -sub 1 5e-3
gg::conBeginSpacing $_CN(35) -sub 1 1e-2

#0.02 is the expected characteristic size of the mesh in the LES region i.e.\Delta_o
gg::conDim $_CN(46) -spacing 0.02
gg::conDim $_CN(44) -spacing 0.02
gg::conBeginSpacing $_CN(45) -sub 1 8e-3
gg::conEndSpacing $_CN(45) -sub 1 1e-1
gg::conEndSpacing $_CN(43) -sub 1 8e-3
gg::conBeginSpacing $_CN(43) -sub 1 1e-1
gg::conDim $_CN(42) -spacing 0.1
gg::conDim $_CN(37) -spacing 0.1
gg::conDim $_CN(38) -spacing 0.1
gg::conDim $_CN(39) -spacing 0.1
gg::conDim $_CN(40) -spacing 0.1
gg::conDim $_CN(41) -spacing 0.1

#########################################################################
#DOMAIN CREATION
#Create RANS region (structured) mesh
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(18)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(14)
    gg::edgeAddCon $_CN(26)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(19)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(15)
  gg::edgeEnd
set _DM(1) [gg::domEnd]
gg::dispConGPS FALSE
gg::dispSmallText FALSE
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(5)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(25)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(6)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(26)
  gg::edgeEnd
set _DM(2) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(11)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(13)
    gg::edgeAddCon $_CN(25)
    gg::edgeAddCon $_CN(12)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(52)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(23)
  gg::edgeEnd
set _DM(3) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(24)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(51)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(10)
    gg::edgeAddCon $_CN(27)
    gg::edgeAddCon $_CN(9)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(2)
  gg::edgeEnd
set _DM(4) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(7)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(28)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(8)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(27)
  gg::edgeEnd
set _DM(5) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(28)
    gg::edgeAddCon $_CN(17)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(20)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(16)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(21)
  gg::edgeEnd
set _DM(6) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(29)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(34)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(48)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(52)
  gg::edgeEnd
set _DM(7) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(4)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(47)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(33)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(48)
  gg::edgeEnd
set _DM(8) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(1)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(30)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(31)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(47)
  gg::edgeEnd
set _DM(9) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(3)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(49)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(35)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(30)
  gg::edgeEnd
set _DM(10) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(49)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(36)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(50)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(22)
  gg::edgeEnd
set _DM(11) [gg::domEnd]
gg::domBegin -type STRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(50)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(32)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(29)
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(51)
  gg::edgeEnd
set _DM(12) [gg::domEnd]

#Create the LES region mesh
gg::domBegin -type UNSTRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(46)
    gg::edgeAddCon $_CN(44)
    gg::edgeReorient
  gg::edgeEnd
  gg::edgeBegin
    gg::edgeAddCon $_CN(35)
    gg::edgeAddCon $_CN(31)
    gg::edgeAddCon $_CN(33)
    gg::edgeAddCon $_CN(34)
    gg::edgeAddCon $_CN(32)
    gg::edgeAddCon $_CN(36)
  gg::edgeEnd
set _DM(13) [gg::domEnd]
gg::domBegin -type UNSTRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(45)
    gg::edgeAddCon $_CN(39)
    gg::edgeAddCon $_CN(43)
    gg::edgeAddCon $_CN(44)
  gg::edgeEnd
set _DM(14) [gg::domEnd]
gg::domBegin -type UNSTRUCTURED
  gg::edgeBegin
    gg::edgeAddCon $_CN(41)
    gg::edgeAddCon $_CN(40)
    gg::edgeAddCon $_CN(45)
    gg::edgeAddCon $_CN(46)
    gg::edgeAddCon $_CN(43)
    gg::edgeAddCon $_CN(38)
    gg::edgeAddCon $_CN(37)
    gg::edgeAddCon $_CN(42)
  gg::edgeEnd
set _DM(15) [gg::domEnd]
set _ggTemp_(1) [list $_DM(14)]
gg::domUnsSolverBegin $_ggTemp_(1)
  gg::domUnsSolverAtt $_DM(14) -boundary_decay 0.92
  gg::domUnsSolverRun REFINE
gg::domUnsSolverEnd
unset _ggTemp_(1)
##################################################################

#Create blocks by extrusion of the domains.
gg::blkExtrusionBegin [list \
    $_DM(1) \
    $_DM(2) \
    $_DM(3) \
    $_DM(4) \
    $_DM(5) \
    $_DM(6) \
    $_DM(7) \
    $_DM(8) \
    $_DM(9) \
    $_DM(10) \
    $_DM(11) \
    $_DM(12) \
    $_DM(13) \
    $_DM(14) \
    $_DM(15) \
  ]  -default HYPERBOLIC
  gg::blkExtrusionMode TRANSLATE
  gg::blkExtrusionAtt -distance 1
  gg::blkExtrusionAtt -direction [list 0 0 1]

  #50 is the number of nodes in the spanwise direction. If 2D the set the number
  #of steps to 1.
  gg::blkExtrusionStep 50 
set _ggTemp_(1) [gg::blkExtrusionEnd]
set _BL(1) [lindex $_ggTemp_(1) 0]
set _BL(2) [lindex $_ggTemp_(1) 1]
set _BL(3) [lindex $_ggTemp_(1) 2]
set _BL(4) [lindex $_ggTemp_(1) 3]
set _BL(5) [lindex $_ggTemp_(1) 4]
set _BL(6) [lindex $_ggTemp_(1) 5]
set _BL(7) [lindex $_ggTemp_(1) 6]
set _BL(8) [lindex $_ggTemp_(1) 7]
set _BL(9) [lindex $_ggTemp_(1) 8]
set _BL(10) [lindex $_ggTemp_(1) 9]
set _BL(11) [lindex $_ggTemp_(1) 10]
set _BL(12) [lindex $_ggTemp_(1) 11]
set _BL(13) [lindex $_ggTemp_(1) 12]
unset _ggTemp_(1)