SWAN-使用说明(By-杨武)(共4页).docx

《SWAN-使用说明(By-杨武)(共4页).docx》由会员分享，可在线阅读，更多相关《SWAN-使用说明(By-杨武)(共4页).docx（4页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

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. 选择中的 build polygon 这时会变成可选模式，单击它后双击封闭区域，出现该对话框分别设置如下：分别点ok后，选择中的map2d mesh，等待后即可生成adcirc所

3、需的三角形网格。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 （读入水

4、深地形文件f14文件）BOUN SEGMENT XY . . . . CONSTANT PAR 2. 10. 0. 30（定义造波边界，需特别注意，三角形网格造波与矩形网格造波不同，不是用SIDE N语句，而是通过SEGMENT 加一条线上起点和终点的X和Y坐标来定的，如上 . . 指点（ ）和点（. ）构成的线段，CONSTANT后面是波浪的一些参数）GEN3（表面在三代模型中运行）FRIC 0.001（定义底摩阻）BREA 1 0.73 （定义波浪破碎条件）DIFFRAC（考虑绕射）OFF QUAD（OFF语句可以把忽略一些你不需要的波浪作用，比如此处是忽略波浪之间的相互作用）NUM AC

5、CUR 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打不开替换后的数据； 专心-专注-专业