fullsplit.script

Include "param.geo";
Merge "part.msh";
N_DOM=NUM_SLICES;

vol[] = {};         idx_vol[] = {};
sur_scat[] = {};    idx_sur_scat[] = {};
sur_inf[] = {};     idx_sur_inf[] = {};
sigma_left[] = {};  idx_sigma_left[] = {};
sigma_right[] = {}; idx_sigma_right[] = {};
For i In {1:N_DOM}
  idx_vol[] += #vol[];                 vol[] += Physical Volume{4000 + i};
  idx_sur_scat[] += #sur_scat[];       sur_scat[] += Physical Surface{1000 + i};
  idx_sur_inf[] += #sur_inf[];         sur_inf[] += Physical Surface{2000 + i};
  idx_sigma_left[] += #sigma_left[];   sigma_left[] += Physical Surface{3000 + 1000 * (i-2)};
  idx_sigma_right[] += #sigma_right[]; sigma_right[] += Physical Surface{3000 + 1000 * (i-1)};
EndFor
idx_vol[] += #vol[];
idx_sur_scat[] += #sur_scat[];
idx_sur_inf[] += #sur_inf[];
idx_sigma_left[] += #sigma_left[];
idx_sigma_right[] += #sigma_right[];

Mesh.Binary = 1;
Geometry.OrientedPhysicals = 0; // partitions have negative tags!!

For i In {1:N_DOM}
  idom = i-1;
  Delete Physicals;
  Physical Volume(4000+i) = vol[{idx_vol[i-1]:idx_vol[i]-1}];
  Physical Surface(1000+i) = sur_scat[{idx_sur_scat[i-1]:idx_sur_scat[i]-1:1}];
  Physical Surface(2000+i) = sur_inf[{idx_sur_inf[i-1]:idx_sur_inf[i]-1:1}];
  Physical Surface(-(4000 + 1000 * (i-2))) = sigma_left[{idx_sigma_left[i-1]:idx_sigma_left[i]-1:1}];
  Physical Surface(4000 + 1000 * (i-1)) = sigma_right[{idx_sigma_right[i-1]:idx_sigma_right[i]-1:1}];
  Save Sprintf("falcon_mshcut%g.msh", idom);
EndFor