SWAN 使用说明(By 杨武)

《SWAN 使用说明(By 杨武)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、SWAN使用说明SWAN可以用矩形网格和三角形网格进行计算，矩形网格计算与SWASH类似，使用方法可参考SWASH使用说明，故此处仅三角形网格使用方法进行简单阐述。SWAN计算波浪主要需要两个文件，一个是水深地形文件，即“.f14”文件，另一个是INPUT文件，即“.swn”文件，本说明主要从水深地形文件的获取，INPUT文件的输入，数据后处理，以及使用过程中的一些细节进行阐述。一．水深地形文件“.f14”的获取：adcirc用于潮流计算，首先需要得到计算网格，它采用三角形网格，用sms里面的adcirc模块绘制网格。绘制网格需要有岸线文件.cst，和水深点即

2、水下地形文件.pts，水深用正值表示。其中岸线文件有一个固定题头如下：COAST必须大写1表示只有一条岸线393表示岸线上的点数0表示岸线不封闭。有了这两个文件后，分别拖入sms中，选择adcirc模块，即可绘制网格，步骤如下：1.首先手绘一条代表外海边界的线，并将边界上的点细分，根据需要分段。2.分别双击代表岸线和外海边界的线，分别选择mainland和ocean，这时外海边界会变成蓝色。3.选择中的buildpolygon这时会变成可选模式，单击它后双击封闭区域，出现该对话框分别设置如下：分别点ok后，选择中的map—2dmesh，等待后即可生成adcir

3、c所需的三角形网格。4.点击meshdata进入mesh模式下，选择file中的saveas将保存类型改为，即可得到计算所需的fort.14文件一．INPUT文件：INPUT示例：***********MODELINPUT*********************************SETLEVEL0.30MODESTATIONARYCGRIDUNSTRUCTUREDCIRCLE360.05211.31（定义网格类型，此处定义为三角形）READGRIDUNSTRUCADCIRC（定义水深地形文件类型）INPGRIDBOTTOMUNSTRUCTURED（读入

4、水深地形文件类型）READINPBOTTOM1.0'fort.14'（读入水深地形文件f14文件）BOUNSEGMENTXY383138.4146497.383660.4191660.CONSTANTPAR2.10.0.30（定义造波边界，需特别注意，三角形网格造波与矩形网格造波不同，不是用SIDEN语句，而是通过SEGMENT加一条线上起点和终点的X和Y坐标来定的，如上3831384146497.383660.4191660指点（3831384146497）和点（383660.4191660）构成的线段，CONSTANT后面是波浪的一些参数）GEN3（表面

5、在三代模型中运行）FRIC0.001（定义底摩阻）BREA10.73（定义波浪破碎条件）DIFFRAC（考虑绕射）OFFQUAD（OFF语句可以把忽略一些你不需要的波浪作用，比如此处是忽略波浪之间的相互作用）NUMACCURSTATMXITST=500ALFA=0.01（设置迭代的精度，为了精确起见，可以设置的相对大些）************OUTPUTREQUESTS*************************BLOCK'COMPGRID'NOHEAD'S.mat'LAY3XPYPTSECHSRTPBLOCK'COMPGRID'NOHEAD'test

6、.mat'LAY3XPYPWATL（输出语句）TEST0,0COMPUTE（计算语句）STOP（结束语句）一．数据后处理：用SWAN波浪模型计算得到的数据有很多种后处理方式，此处仅介绍我所使用的处理方式。“.mat”文件是MATLAB文件，用MATLAB打开“.mat”后，选中HSIG，会显示一行矩阵，其所对应的个数正对应你所用的三角形网格的节点数，为了方便起见，将其转置一下，得到一列数据，这一列数据就是有效波高，把“.f14”地形水深文件中的水深数据替换掉，此时再用SMS打开就可以显示有效波高的分布了。一．一些注意事项：1.三角形网格造波的方法与矩形网格不同

7、，你所选取的起点和终点必须在计算范围内；2.用有效波高替换水深数据时要保证替换前后数据格式一致，否则SMS打不开替换后的数据；