SWAN使用说明(By杨武).docx

SWAN使用说明SWAN可以用矩形网格和三角形网格进行计算，矩形网格计算与SWASH类似，使用方法可参考SWASH使用说明，故此处仅三角形网格使用方法进行简单阐述。SWAN计算波浪主要需要两个文件，一个是水深地形文件，即“.f14”文件，另一个是INPUT文件，即“.swn”文件，本说明主要从水深地形文件的获取，INPUT文件的输入，数据后处理，以及使用过程中的一些细节进行阐述。水深地形文件“.f14”的获取：adcirc用于潮流计算，首先需要得到计算网格，它采用三角形网格，用sms里面的adcirc模块绘制网格。绘制网格需要有岸线文件 .cst，和水深点即水下地形文件 .pts，水深用正值表示。其中岸线文件有一个固定题头如下： COAST必须大写 1表示只有一条岸线 393表示岸线上的点数 0表示岸线不封闭。有了这两个文件后，分别拖入sms中，选择adcirc模块，即可绘制网格，步骤如下：1. 首先手绘一条代表外海边界的线，并将边界上的点细分，根据需要分段。2. 分别双击代表岸线和外海边界的线，分别选择mainland和ocean，这时外海边界会变成蓝色。3. 选择中的 build polygon 这时会变成可选模式，单击它后双击封闭区域，出现该对话框分别设置如下：分别点ok后，选择中的map—2d mesh，等待后即可生成adcirc所需的三角形网格。4. 点击mesh data 进入mesh模式下，选择file中的save as 将保存类型改为，即可得到计算所需的fort.14文件INPUT文件：INPUT 示例：***********MODEL INPUT*********************************SET LEVEL 0.30MODE STATIONARYCGRID UNSTRUCTURED CIRCLE 36 0.0521 1. 31（定义网格类型，此处定义为三角形）READGRID UNSTRUC ADCIRC（定义水深地形文件类型）INPGRID BOTTOM UNSTRUCTURED（读入水深地形文件类型）READINP BOTTOM 1.0 fort.14 \（读入水深地形文件f14文件）BOUN SEGMENT XY 383138. 4146497. 383660. 4191660. CONSTANT PAR 2. 10. 0. 30（定义造波边界，需特别注意，三角形网格造波与矩形网格造波不同，不是用SIDE N语句，而是通过SEGMENT 加一条线上起点和终点的X和Y坐标来定的，如上383138 4146497. 383660. 4191660指点（383138 4146497）和点（383660. 4191660）构成的线段，CONSTANT后面是波浪的一些参数）GEN3（表面在三代模型中运行）FRIC 0.001（定义底摩阻）BREA 1 0.73 （定义波浪破碎条件）DIFFRAC（考虑绕射）OFF QUAD（OFF语句可以把忽略一些你不需要的波浪作用，比如此处是忽略波浪之间的相互作用）NUM ACCUR STAT MXITST=500 ALFA=0.01（设置迭代的精度，为了精确起见，可以设置的相对大些）************ OUTPUT REQUESTS *************************BLOCK COMPGRID NOHEAD S.mat LAY 3 XP YP TSEC HS RTPBLOCK COMPGRID NOHEAD test.mat LAY 3 XP YP WATL （输出语句）TEST 0,0COMPUTE （计算语句）STOP（结束语句）数据后处理：用SWAN波浪模型计算得到的数据有很多种后处理方式，此处仅介绍我所使用的处理方式。“.mat”文件是MATLAB文件，用MATLAB打开“.mat”后，选中HSIG，会显示一行矩阵，其所对应的个数正对应你所用的三角形网格的节点数，为了方便起见，将其转置一下，得到一列数据，这一列数据就是有效波高，把“.f14”地形水深文件中的水深数据替换掉，此时再用SMS打开就可以显示有效波高的分布了。一些注意事项： 1.三角形网格造波的方法与矩形网格不同，你所选取的起点和终点必须在计算范围内； 2.用有效波高替换水深数据时要保证替换前后数据格式一致，否则SMS打不开替换后的数据；