79 #define L (0.5 * (1. - G))  macro
84   {H, 0, 0, 0, L, 0, 0, 0},
85   {0, H, 0, 0, 0, L, 0, 0},
86   {0, 0, H, 0, 0, 0, L, 0},
87   {0, 0, 0, H, 0, 0, 0, L},
88   {L, 0, 0, 0, H, 0, 0, 0},
89   {0, L, 0, 0, 0, H, 0, 0},
90   {0, 0, L, 0, 0, 0, H, 0},
91   {0, 0, 0, L, 0, 0, 0, H}
94 …0},    {0, 0, 0},        {0, P *H, 0},     {M * L, M *L, P}, {P * L, 0, 0},    {0, 0, 0},        {…
95 … * H, M *H, M}, {0, P *H, 0},     {0, 0, 0},        {M * L, 0, 0},    {P * L, M *L, P}, {0, P *L, …
96 … * H, P *H, M}, {M * H, 0, 0},    {0, 0, 0},        {0, M *L, 0},     {P * L, P *L, P}, {M * L, 0,…
97 …    {P * H, 0, 0},    {M * H, P *H, M}, {0, M *L, 0},     {0, 0, 0},        {P * L, 0, 0},    {M *…
98 …{{M * L, M *L, M}, {P * L, 0, 0},    {0, 0, 0},        {0, P *L, 0},     {M * H, M *H, P}, {P * H,…
99 …{{M * L, 0, 0},    {P * L, M *L, M}, {0, P *L, 0},     {0, 0, 0},        {M * H, 0, 0},    {P * H,…
100 …{{0, 0, 0},        {0, M *L, 0},     {P * L, P *L, M}, {M * L, 0, 0},    {0, 0, 0},        {0, M *…
101 …{{0, M *L, 0},     {0, 0, 0},        {P * L, 0, 0},    {M * L, P *L, M}, {0, M *H, 0},     {0, 0, …
105   {H * H * H, L * H * H, L * L * H, H * L * H, H * H * L, L * H * L, L * L * L, H * L * L},
106   {L * H * H, H * H * H, H * L * H, L * L * H, L * H * L, H * H * L, H * L * L, L * L * L},
107   {L * L * H, H * L * H, H * H * H, L * H * H, L * L * L, H * L * L, H * H * L, L * H * L},
108   {H * L * H, L * L * H, L * H * H, H * H * H, H * L * L, L * L * L, L * H * L, H * H * L},
109   {H * H * L, L * H * L, L * L * L, H * L * L, H * H * H, L * H * H, L * L * H, H * L * H},
110   {L * H * L, H * H * L, H * L * L, L * L * L, L * H * H, H * H * H, H * L * H, L * L * H},
111   {L * L * L, H * L * L, H * H * L, L * H * L, L * L * H, H * L * H, H * H * H, L * H * H},
112   {H * L * L, L * L * L, L * H * L, H * H * L, H * L * H, L * L * H, L * H * H, H * H * H}
115 …L * M * H, L * H * M}, {P * L * H, L * P * H, L * L * M}, {M * L * H, H * P * H, H * L * M}, {M * …
116 …L * M * H, L * H * M}, {P * H * H, H * M * H, H * H * M}, {P * L * H, H * P * H, H * L * M}, {M * …
117 …L * H, L * M * H, L * L * M}, {P * L * H, H * M * H, H * L * M}, {P * H * H, H * P * H, H * H * M}…
118 …L * H, H * M * H, H * L * M}, {P * L * H, L * M * H, L * L * M}, {P * H * H, L * P * H, L * H * M}…
119 …L, H * M * L, H * H * M}, {P * H * L, L * M * L, L * H * M}, {P * L * L, L * P * L, L * L * M}, {M…
120 …L, L * M * L, L * H * M}, {P * H * L, H * M * L, H * H * M}, {P * L * L, H * P * L, H * L * M}, {M…
121 …L * L, L * M * L, L * L * M}, {P * L * L, H * M * L, H * L * M}, {P * H * L, H * P * L, H * H * M}…
122 …L * L, H * M * L, H * L * M}, {P * L * L, L * M * L, L * L * M}, {P * H * L, L * P * L, L * H * M}…
127   {H * H, L *H, L *L, H *L},
128   {L * H, H *H, H *L, L *L},
129   {L * L, H *L, H *H, L *H},
130   {H * L, L *L, L *H, H *H}
133   {{M * H, M *H}, {P * H, M *L}, {P * L, P *L}, {M * L, P *H}},
134   {{M * H, M *L}, {P * H, M *H}, {P * L, P *H}, {M * L, P *L}},
135   {{M * L, M *L}, {P * L, M *H}, {P * H, P *H}, {M * H, P *L}},
136   {{M * L, M *H}, {P * L, M *L}, {P * H, P *L}, {M * H, P *H}}
140 #undef L
541     PetscReal      L, m = 1.0;  in THICreate()  local
543     L = thi->Lx;  in THICreate()
544     PetscCall(PetscOptionsReal("-thi_L", "Domain size (m)", "", L, &L, &flg));  in THICreate()
545     if (flg) thi->Lx = thi->Ly = L;  in THICreate()