/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  8
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.01;  // centimetres.  According to https://cfd.direct/openfoam/user-guide/v8-blockmesh/

#include "../scale";  // gives global N
rx 3.5355339059327373;
arcNx #calc "3 * $N";
arcNy #calc " 2 * $N";

// = // rx #calc "5.0 / Foam::sqrt(2.0)";  // length of short side of rectangular triangle with long side =5

vertices
(
    // F = front (z = 0), B = back   (z = 1)
    // bottom = b, t = top, n = middle (left/right of center of cylinder)
    // l = left, r = right, m = middle (under center of cylinder), s = other square around cylinder corner
    // c = cylinder
    // (0, 0, 0) = center of cylinder in front
    name Fbl (-20 -20 0)
    name Bbl (-20 -20 1)
    name Fbs (20 -20 0)
    name Bbs (20 -20 1)
    name Fbr (200 -20 0)
    name Bbr (200 -20 1)
    name Ftl (-20 21 0)
    name Btl (-20 21 1)
    name Fts (20 21 0)
    name Bts (20 21 1)
    name Ftr (200 21 0)
    name Btr (200 21 1)
    name Fcl (-5 0 0)
    name Bcl (-5 0 1)
    name Fcb (0 -5 0)
    name Bcb (0 -5 1)
    name Fcr (5 0 0)
    name Bcr (5 0 1)
    name Fct (0 5 0)
    name Bct (0 5 1)
    name Fcbl (#calc "-$rx" #calc "-$rx" 0) 
    name Bcbl (#calc "-$rx" #calc "-$rx" 1) 
    name Fcbr ($rx #calc "-$rx" 0)
    name Bcbr ($rx #calc "-$rx" 1)
    name Fctl (#calc "-$rx" $rx 0) 
    name Bctl (#calc "-$rx" $rx 1) 
    name Fctr ($rx $rx 0)
    name Bctr ($rx $rx 1)
    name Fnl (-20 0 0)
    name Bnl (-20 0 1)
    name Fnr (200 0 0)
    name Bnr (200 0 1)
    name Fbm (0 -20 0)
    name Bbm (0 -20 1)
    name Ftm (0 21 0)
    name Btm (0 21 1)
    name Fns (20 0 0)
    name Bns (20 0 1)
    name Bnr (200 0 0)
    name Fnr (200 0 1)
    name Fnm (0 0 0)  
    name Bnm (0 0 0)
);

blocks
(
//    name total hex (Fbl Fbr Ftr Ftl Bbl Bbr Btr Btl) (50 5 1) simpleGrading (1 1 1)
      name blbarc hex (Fbl Fbm Fcb Fcbl Bbl Bbm Bcb Bcbl) ($arcNx $arcNy 1) simpleGrading (1 1 1)
      name bltarc hex (Fnl Fbl Fcbl Fcl Bnl Bbl Bcbl Bcl) ($arcNx $arcNy 1) simpleGrading (1 1 1)
      name tlbarc hex (Ftl Fnl Fcl Fctl Btl Bnl Bcl Bctl) ($arcNx $arcNy 1) simpleGrading (1 1 1)
      name tltarc hex (Ftm Ftl Fctl Fct Btm Btl Bctl Bct) ($arcNx $arcNy 1) simpleGrading (1 1 1)
      name trtarc hex (Fts Ftm Fct Fctr Bts Btm Bct Bctr) ($arcNx $arcNy 1) simpleGrading (1 1 1)
      name trbarc hex (Fns Fts Fctr Fcr Bns Bts Bctr Bcr) ($arcNx $arcNy 1) simpleGrading (1 1 1)
      name brtarc hex (Fbs Fns Fcr Fcbr Bbs Bns Bcr Bcbr) ($arcNx $arcNy 1) simpleGrading (1 1 1)
      name brbarc hex (Fbm Fbs Fcbr Fcb Bbm Bbs Bcbr Bcb) ($arcNx $arcNy 1) simpleGrading (1 1 1)
      name pipetr hex (Fns Fnr Ftr Fts Bns Bnr Btr Bts) (#calc "12 * $arcNy" #calc "$arcNx" 1) simpleGrading (1 1 1)
      name pipebr hex (Fbs Fbr Fnr Fns Bbs Bbr Bnr Bns) (#calc "12 * $arcNy" #calc "$arcNx" 1) simpleGrading (1 1 1)
);

edges
(
    name Fblbarc arc Fcb Fcbl origin Fnm 
    name Fbltarc arc Fcbl Fcl origin Fnm 
    name Ftlbarc arc Fcl Fctl origin Fnm 
    name Ftltarc arc Fctl Fct origin Fnm 
    name Ftrtarc arc Fct Fctr origin Fnm 
    name Ftrbarc arc Fctr Fcr origin Fnm 
    name Fbrtarc arc Fcr Fcbr origin Fnm 
    name Fbrbarc arc Fcbr Fcb origin Fnm 

    name Bblbarc arc Bcb Bcbl origin Bnm 
    name Bbltarc arc Bcbl Bcl origin Bnm 
    name Btlbarc arc Bcl Bctl origin Bnm 
    name Btltarc arc Bctl Bct origin Bnm 
    name Btrtarc arc Bct Bctr origin Bnm 
    name Btrbarc arc Bctr Bcr origin Bnm 
    name Bbrtarc arc Bcr Bcbr origin Bnm 
    name Bbrbarc arc Bcbr Bcb origin Bnm 

);

boundary
(
/*
    movingWall
    {
        type wall;
        faces
        (
            (3 7 6 2)
        );
    }
    fixedWalls
    {
        type wall;
        faces
        (
            (0 4 7 3)
            (2 6 5 1)
            (1 5 4 0)
        );
    }
    frontAndBack
    {
        type empty;
        faces
        (
            (0 3 2 1)
            (4 5 6 7)
        );
    }
*/
);

mergePatchPairs
(
);

// ************************************************************************* //