96 #define L (0.5 * (1. - G))  macro
101   {H, 0, 0, 0, L, 0, 0, 0},
102   {0, H, 0, 0, 0, L, 0, 0},
103   {0, 0, H, 0, 0, 0, L, 0},
104   {0, 0, 0, H, 0, 0, 0, L},
105   {L, 0, 0, 0, H, 0, 0, 0},
106   {0, L, 0, 0, 0, H, 0, 0},
107   {0, 0, L, 0, 0, 0, H, 0},
108   {0, 0, 0, L, 0, 0, 0, H}
111 …0},    {0, 0, 0},        {0, P *H, 0},     {M * L, M *L, P}, {P * L, 0, 0},    {0, 0, 0},        {…
112 … * H, M *H, M}, {0, P *H, 0},     {0, 0, 0},        {M * L, 0, 0},    {P * L, M *L, P}, {0, P *L, …
113 … * H, P *H, M}, {M * H, 0, 0},    {0, 0, 0},        {0, M *L, 0},     {P * L, P *L, P}, {M * L, 0,…
114 …    {P * H, 0, 0},    {M * H, P *H, M}, {0, M *L, 0},     {0, 0, 0},        {P * L, 0, 0},    {M *…
115 …{{M * L, M *L, M}, {P * L, 0, 0},    {0, 0, 0},        {0, P *L, 0},     {M * H, M *H, P}, {P * H,…
116 …{{M * L, 0, 0},    {P * L, M *L, M}, {0, P *L, 0},     {0, 0, 0},        {M * H, 0, 0},    {P * H,…
117 …{{0, 0, 0},        {0, M *L, 0},     {P * L, P *L, M}, {M * L, 0, 0},    {0, 0, 0},        {0, M *…
118 …{{0, M *L, 0},     {0, 0, 0},        {P * L, 0, 0},    {M * L, P *L, M}, {0, M *H, 0},     {0, 0, …
122   {H * H * H, L * H * H, L * L * H, H * L * H, H * H * L, L * H * L, L * L * L, H * L * L},
123   {L * H * H, H * H * H, H * L * H, L * L * H, L * H * L, H * H * L, H * L * L, L * L * L},
124   {L * L * H, H * L * H, H * H * H, L * H * H, L * L * L, H * L * L, H * H * L, L * H * L},
125   {H * L * H, L * L * H, L * H * H, H * H * H, H * L * L, L * L * L, L * H * L, H * H * L},
126   {H * H * L, L * H * L, L * L * L, H * L * L, H * H * H, L * H * H, L * L * H, H * L * H},
127   {L * H * L, H * H * L, H * L * L, L * L * L, L * H * H, H * H * H, H * L * H, L * L * H},
128   {L * L * L, H * L * L, H * H * L, L * H * L, L * L * H, H * L * H, H * H * H, L * H * H},
129   {H * L * L, L * L * L, L * H * L, H * H * L, H * L * H, L * L * H, L * H * H, H * H * H}
132 …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 * …
133 …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 * …
134 …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}…
135 …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}…
136 …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…
137 …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…
138 …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}…
139 …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}…
144   {H * H, L *H, L *L, H *L},
145   {L * H, H *H, H *L, L *L},
146   {L * L, H *L, H *H, L *H},
147   {H * L, L *L, L *H, H *H}
150   {{M * H, M *H}, {P * H, M *L}, {P * L, P *L}, {M * L, P *H}},
151   {{M * H, M *L}, {P * H, M *H}, {P * L, P *H}, {M * L, P *L}},
152   {{M * L, M *L}, {P * L, M *H}, {P * H, P *H}, {M * H, P *L}},
153   {{M * L, M *H}, {P * L, M *L}, {P * H, P *L}, {M * H, P *H}}
157 #undef L
471     PetscReal      L, m = 1.0;  in THICreate()  local
473     L = thi->Lx;  in THICreate()
474     PetscCall(PetscOptionsReal("-thi_L", "Domain size (m)", "", L, &L, &flg));  in THICreate()
475     if (flg) thi->Lx = thi->Ly = L;  in THICreate()