获取网格位置信息

getcontour_data.m

写入输入文件

fort.33的存储规则

• 1. 第一行总三角网格的vertices数目（nnode）
• 2. 每一行四个位置点，共（nnode）/2+1行（最后一行）
• 3. 读取规则如下

nnodes = fscanf(fid,'%d',1);
%% 4998*2 
nodes  = zeros(nnodes,2);
nodes(:,1) = fscanf(fid,'%f',nnodes);
nodes(:,2) = fscanf(fid,'%f',nnodes);
% Convert from cm to m
%% 单位转换
nodes = nodes*1e-2;

即，先存放r坐标，然后存放z坐标

写入fort.33文件

fort.33的vertices排列顺序

• 1. 最外层边界（index：-1），初始点，然后加密

     
     image.png

• 2. 从内靶板开始的经过b2 mesh boundary 到达外靶板（index：2）

     
     image.png

• 3. structure (index:3)

     
     image.png

• 4. pump & dome

     
     image.png

• 5. lower PFR

     
     image.png

• 6. 边界内填充

     
     image.png

• 7. plasma region

     
     image.png

排列顺序为：极向优先排列（98*38）

限制区域三角网格的产生

遵循，先初始node，随后逐渐加密。
以为