mike21-操作实例---矩形网格

MIKE21HD操作简要说明1一数据文件格式1.land.xyz（陆地边界文件）格式说明：纬度tab经度tabtab陆地高程值12.3030756.1276911012.3089556.1290911012.3089556.1278411012.3030756.12769010。。。。。。说明：第三列中0表示陆地边界的结束点，该点与陆地起始点重叠；每两列之间用TAB键隔开，不要用空格。2.water.xyz（水深文件）格式说明：12.1680556.17699-20.00312.1666556.17715-20.00312.1660956.17762-20.00312.1663756.17808-20.003。。。。。。说明：第一、二列分别为纬度、经度，第三列为该点处水深，第四列暂时不详；其他说明同陆地文件；对于上述文件格式也可用直角坐标（大地坐标）表示，但是在DefineworkingArea对话框中要选定DecimalDeg项。如下图所示：3.WL1（水位边界条件文件）格式说明：WaterlevelrecordingsfromSt1timeElevation1993-12-0200:00:00-0.33选择该项21993-12-0200:30:00-0.3621993-12-0201:00:00-0.3965。。。。。。说明：第一列日期，第二列时间，第三列水位；第一列和第二列之间用空格隔开，第二列和第三列之间用TAB隔开；前两行为文件说明语句应保留。4.windcastrup(wind)（风边界条件文件）格式说明：WinddataTimeSpeedDirectionUnit10000220000100003240101993-12-0200:00:009.294184.261993-12-0200:30:0010.066186.6891993-12-0201:00:0010.655189.167。。。。。。说明：第一列为日期，第二列为时间，第三列为风速，第四列为风向（与北向顺时针夹角）；第一、二行为文件说明语句应保留；第三行暂时不详；第一列和第二列之间用空格隔开，第二行和第三行及第四行之间用TAB键隔开。5.WaterlevelNdrRoese（水位验证文件）格式说明：NdrRoese:WaterLevel[m]TimeNrdRoese:WaterLevel[m]Unit100000100001993-12-0200:00:00-0.3443751993-12-0200:15:00-0.3490631993-12-0200:30:00-0.3485881993-12-0200:45:00-0.397647。。。。。。说明：该文件为“NdrRoese”站实测水位；第一、二列之间用空格隔开，第二、三列之间用TAB隔开；第三列为水位。6.WindKastrup（风过程文件）格式说明二以实例演示MIKE21的水动力模型及输移扩散模型的操作步骤1.网格数据（1）命令操作流程①水深网格文件设置菜单流程打开mike21→file→new→Bathymetries→打开“DefineWorkingArea”对话框后设置→选择菜单“WorkArea”→选择“BackgroundManagement…”→点击“import…”按钮→选择文件类型为“*.xyz”→找到陆地边界文件“land”→再选择convertform中的“LONG/LAT”后打开陆地边界文件→在点击“import…”按钮→按照陆地边界文件的操作打开水深数据文件→后ok→在选择“WorkArea”中的“BathymetryManagement…”后出一个对话框→点击“New”出现计算区域设置对话框，设置后→OK→此时图形已经框出计算区域，在点击工具栏中按钮选择要插值的水深数据范围→选定后在点击一次→再打开3“WorkArea”中的“BathymetryManagement…”→点击插值按钮“interpolate…”→设置后OK→保存已经设置好的计算区域网格文件（例如文件名“Bat2”文件）→“File”→“open”打开已经生成的“Bat2-0”文件→可以应用工具栏中的工具及“Tools”菜单中的“SetValue…”编辑水深及陆地值→存盘，至此已经完成网格数据文件的编辑。②嵌套网格数据文件的设置流程第一步：建立子网格打开mike21→file→new→选定mikezero中的“MikezeroToolbox”→展开“Transformation”→选定“RotateGrid”→此时弹出一个对话框“setupName”设置子网格区域的名称→下一步→出现“selectinputfile”选择要嵌入子网格的网格文件（例如Bat2-0文件）→下一步→出现“subseriesselection”，在此可以设置子网格需要计算的时间段，也可跳过→下一步→“selectitem”该项暂时不详．．．．．．，跳过→下一步→出现“GridSize”，该项用于设置子网格的尺寸（第一项：设置子网格的结点数，结点数=3×大网格数+1，且子区域边界要和大区域中的网格线重叠；第二项：子网格的网格宽度，取上一级网格宽度的1/3；第三项为子区域的左下角在大网格区域中的位置；最后一项为子网格区域在大网格中的转角，与直角坐标X向的夹角）→下一步→因子网格区域是否包括陆地和水域而提供的选项，本例中包含陆地和水域故选择“dateisBathymetry…”项→下一步→“outputspecification”输出子网格的文件名称（以Bat2-0zi为例）→下一步→点击“Execute”按钮→完成。第二步：子网格边界处理打开mike21→file→new→选定mike21中的“Mike21Toolbox”→展开“Hydrodynamics”→选定“BoderAdjustment”→设置名称→下一步→出现“InputFile”界面，（第一项：输入大网格文件；第二项：输入小网格文件；最后两项“X-coordina1和Y-coordina1”为子区域左下角在大网格中的位置）→下一步→出现“SelectOutputFile”界面（第一项：输入处理后的大网格文件名称；第二项：输出处理后的小网格文件名称）→下一步→点击“Execute”按钮→完成。经过上述两个步骤就完成了对嵌套网格的设置。对于嵌套两个同级网格（如下图所示）的操作说明：将“大”网格和“小1”网格按上述两步处理后，得到“大-处理”和“小1-处理”，后在将“大-处理”和“小2”网格按上述两步处理后，得到“大-处理1”和“小2-处理”，则在水动力模型计算时，输入“大-处理1、小1-处理、小2-处理”三个网格文件。如果三重嵌套，例如“小2”有嵌套在了“小1”中，则要两、两发生关系，既“大”和“小1”处理后的“小1-处理”在和“小2”处理，“小2-处理”在和“大-处理”进行处理。多重嵌套以此类推。大小1小24（2）图形流程示意第一步：见图1图1第二步：见图2图2第三步：加入背景文件，见图3、图4图3选定UTM带号：Utmzone=（经度÷6+31）的整数部分定义区域左下角在大地坐标系下的位置显示区域的宽显示区域的高在显示的区域内添加，陆地及水深背景5图4将陆地边界和水深数据加入之后图形显示如图5所示：图5第四步：选定计算区域，见图6、图7图62.文件类型为：（*.xyz）1.输入3.陆地和水域文件坐标点均为经纬度，故要选择将坐标由经纬度转化为直角坐标选项！选择该项，编辑计算区域6图7经过上述设定后既圈定了计算区域，见图8图8第五步：在圈定计算区域之后，要选定要插值的数据范围。应用工具栏中的按忸，在图形中框取，选定后再点击一次．．．．．按忸．．，之后在调出如图9的对话框。图91.新建一个计算区域2.指定计算区域左下角在显示区域中的坐标5.指定计算区域的转角（与北向夹角）3.指定计算网格的宽度4.指定计算区域在X、Y方向的结点数。6.陆地区域高程值（要和陆地边界文件一致）该框圈定区域既为计算区域1.选择该项，对计算区域的数据进行插值！2.指定插值时寻找结点周围数据的个数！3.选择插值方法！7插值结束之后会出现消息框，该框显示插值所用的时间信息。待全部ok之后，存盘！！（文件名为bat1），此时系统自动生成一个文件名为bat1_0的文件，该文件既为计算用的网格文件。然后打开该网格文件，如图10所示。图10第六步：网格数据调整。可用工具栏中如图11中所示的按钮选定要编辑的网格，然后选择菜单“Tool”下的“Setvalue…”（见图12、图13），调整后如图14所示，保存。图11图12图13图14处理后的网格数据图需要设陆地高程值的区域错误区域在此设置选定网格的值。8至此计算区域的网格数据文件已经整理完毕。如果要在本计算区域内内嵌网格，操作步骤如下：第一步：建立子网格图15图16图18图17图19图20（该项暂不详．．．．．）选择该项选择该项设置子网格名称选择被嵌入的大网格文件设置嵌套的小网格中计算的起始和结束时间，默认时间同上一级网格的计算时间9图21图22图24图23至此已经完成了小网格的编辑，小网格的文件名为Bat1_0(qiantao)，图25所示为嵌套后的小网格。图25嵌套好的计算区域指定X、Y向的小网格结点数：结点数=3×大网格数+1小网格宽度（网格宽度必须为上一级网格宽度的1/3）指定小网格左下角在大网格中的位置嵌套网格在大网格中的转角（与X方向）选择该项，指明小网格中包括陆地和水深输出编辑好的小网格文件结束之前先运行，点击该键10第二步：子网格的边界处理图26图27图29图28图30图31至此已经完成了嵌套网格的准备工作，也产生了两个用于模型计算的网格数据文件，在本例中分别为：大网格文件“Bat1_0(chuli)”，小网格文件“Bat1_0(qiantao)chuli”。选择该项选择该项输入大网格文件名输入小网格文件名输入小网格左下角在大网格中的位置输出调整后的大网格文件名输出调整后的小网格文件名设定都结束后，运行完成112.边界条件的设置（1）命令操作流程第一步：设置时间序列文件（*.dfs0）打开mike21→file→new→选择mikezero中的“Timeseries”（时间序列）→选择“FromAsciiFile”→出现“TimeseriesEidetor:importfromAscii”对话框→在此对话框中选择“FileName”（本例为WL1.txt）→将“Timediscription”项设置成“EquidistantCalendarAxis”（等时间序列轴）→ok→出现水位边界条件“WL1.TXT”的过程线，此时鼠标右键单击此图→选择“properties…”，出现“Fileproperties”对话框→将AxisType设置成“EquidistantcalendarAxis”；将“Iteminformation”中“Type”项中的“undefined”设置成“waterlevel”（水位）→ok→保存文件（本例为Ts1.dfs0）。重复上述的操作，直至将上下两个边界的4个点时间序列文件处理完毕。验证文件也按上述操作过程设置。第二步：建立边界条件打开mike21→file→new→选择“profileseries”（断面序列）→选择“blankTimeseries”（空白时间序列）→此时弹出“Fileproperties”对话框→设置“title”（上边界条件名称，本例为shangbianjie）；设置AxisType为“EquidistantcalendarAxis”；设置Timestep为1800s（因为边界条件中的时间间隔为30分钟）；设置No.ofpoints为13（本例中，上边界的网格数为13个）；设置Gridstep（网格步长）为900；将iteminformation中的Type项由undefined该为waterlevel→ok→此时出现一个空白的断面边界条件编辑器→打开北边界的左侧的点时间文件TS1.dfs