mike21 操作實(shí)例-矩形網(wǎng)格

1、一數(shù)據(jù)文件格式l.land.xyz（陸地邊界文件）格式說明：緯度tab經(jīng)度tabtab陸地高程值TOC o 1-5 h z12.3030756.1276911012.3089556.1290911012.3089556.1278411012.3030756.12769010OOOOOO說明：第三列中。表示陸地邊界的結(jié)束點(diǎn)，該點(diǎn)與陸地起始點(diǎn)重疊；每兩列之間用TAB鍵隔開，不要用空格。.water.xyz（水深文件）格式說明：TOC o 1-5 h z12.1680556.17699-20.00312.1666556.17715-20.00312.1660956.17762-20.00312.16

2、63756.17808-20.003說明：第一、二列分別為緯度、經(jīng)度，第三列為該點(diǎn)處水深，第四列暫時(shí)不詳；其他說明同陸地文件；對(duì)于上述文件格式也可用直角坐標(biāo)（大地坐標(biāo)）表示，但是在DefineworkingArea對(duì)話框中要選定DecimalDeg項(xiàng)。如下圖所示：.WL1（水位邊界條件文件）格式說明:WaterlevelrecordingsfromSt1timeElevation1993-12-0200:00:00-0.331993-12-0200:30:00-0.3621993-12-0201:00:00-0.3965說明：第一列日期，第二列時(shí)間，第三列水位；第一列和第二列之間用空格隔開，第

3、二列和第三列之間用TAB隔開；前兩行為文件說明語句應(yīng)保留。.windcastrup（wind）（風(fēng)邊界條件文件）格式說明：WinddataTimeSpeedDirection10000324010184.2610000324010184.26186.689189.1671993-12-0200:00:009.2941993-12-0200:30:0010.0661993-12-0201:00:0010.655說明：第一列為日期，第二列為時(shí)間，第三列為風(fēng)速，第四列為風(fēng)向（與北向順時(shí)針夾角）;第一、二行為文件說明語句應(yīng)保留；第三行暫時(shí)不詳；第一列和第二列之間用空格隔開，第二行和第三行及第四行之間用

4、TAB鍵隔開。.WaterlevelNdrRoese（水位驗(yàn)證文件）格式說明：NdrRoese:WaterLevelmTimeNrdRoese:WaterLevelmUnit100000100001993-12-0200:00:00-0.3443751993-12-0200:15:00-0.3490631993-12-0200:30:00-0.3485881993-12-0200:45:00-0.397647說明：該文件為“用出0。5?！闭緦?shí)測水位；第一、二列之間用空格隔開，第二、三列之間用TAB隔開；第三列為水位。.WindKastrup（風(fēng)過程文件）格式說明二以實(shí)例演示MIKE21的水動(dòng)力

5、模型及輸移擴(kuò)散模型的操作步驟1.網(wǎng)格數(shù)據(jù)（1）命令操作流程水深網(wǎng)格文件設(shè)置菜單流程打開mike21-filenewBathymetries打開“DefineWorkingArea”對(duì)話框后設(shè)置選擇菜單“WorkArea”一選擇“BackgroundManagement”一點(diǎn)擊“import”按鈕一選擇文件類型為“*.xyz”找到陸地邊界文件“加”一再選擇convertform中的“LONG/LAT”后打開陸地邊界文件一在點(diǎn)擊“import”按鈕一按照陸地邊界文件的操作打開水深數(shù)據(jù)文件一后ok在選擇“WorkArea”中的“BathymetryManagement”后出一個(gè)對(duì)話框一點(diǎn)擊“New”

6、出現(xiàn)計(jì)算區(qū)域設(shè)置對(duì)話框，設(shè)置后一OK此時(shí)圖形已經(jīng)框出計(jì)算區(qū)域，在點(diǎn)擊工具欄中工按鈕選擇要插值的水深數(shù)據(jù)范圍一選定后在點(diǎn)擊一次一再打開“WorkArea”中的“BathymetryManagement”一點(diǎn)擊插值按鈕“interpolate”一設(shè)置后OK保存已經(jīng)設(shè)置好的計(jì)算區(qū)域網(wǎng)格文件（例如文件名配前2”文件）“File”一“open”打開已經(jīng)生成的“汩212-0”文件一可以應(yīng)用工具欄中的工具及“Tools”菜單中的“5。萬比”編輯水深及陸地值一存盤，至此已經(jīng)完成網(wǎng)格數(shù)據(jù)文件的編輯。嵌套網(wǎng)格數(shù)據(jù)文件的設(shè)置流程精心整理第一步：建立子網(wǎng)格打開mike21ffilenew選定mikezero中的“Mi

7、kezeroToolbox”展開aTransformation”選定“RotateGrid”一此時(shí)彈出一個(gè)對(duì)話框“setupName設(shè)置子網(wǎng)格區(qū)域的名稱一下一步一出現(xiàn)selectinputfile選擇要嵌入子網(wǎng)格的網(wǎng)格文件（例如Bat2-0文件）一下一步一出現(xiàn)“subseriesselection”，在此可以設(shè)置子網(wǎng)格需要計(jì)算的時(shí)間段，也可跳過一下一步一“selectitem”闔陶寸不詳，跳過一下一步-出現(xiàn)“GridSize”，該項(xiàng)用于設(shè)置子網(wǎng)格的尺寸（第一項(xiàng)：設(shè)置子網(wǎng)格的結(jié)點(diǎn)數(shù)，結(jié)點(diǎn)數(shù)=3X大網(wǎng)格數(shù)+1,且子區(qū)域邊界要和大區(qū)域中的網(wǎng)格線重疊；第二項(xiàng)：子網(wǎng)格的網(wǎng)格寬度，取上一級(jí)網(wǎng)格寬度的1/3

8、;第三項(xiàng)為子區(qū)域的左下角在大網(wǎng)格區(qū)域中的位置；最后一項(xiàng)為子網(wǎng)格區(qū)域在大網(wǎng)格中的轉(zhuǎn)角，與直角坐標(biāo)X向的夾角）一下一步一因子網(wǎng)格區(qū)域是否包括陸地和水域而提供的選項(xiàng)，本例中包含陸地和水域故選擇“dateisBathymetry”項(xiàng)一下一步aoutputspecification”輸出子網(wǎng)格的文件名稱（以Bat2-0zi為例）一下一步一點(diǎn)擊“Execute”按鈕一完成。第二步：子網(wǎng)格邊界處理打開mike21-filenew選定mike21中的“Mike21Toolbox”f展開“Hydrodynamics”選定“BoderAdjustment”f設(shè)置名稱一下一步一出現(xiàn)“彼西江”界面，（第一項(xiàng)：輸入大網(wǎng)

9、格文件;第二項(xiàng)：輸入小網(wǎng)格文件；最后兩項(xiàng)“X-coordina1和Y-coordina1”為子區(qū)域左下角在大網(wǎng)格中的位置）一下一步一出現(xiàn)“SelectOutputFile”界面（第一項(xiàng)：輸入處理后的大網(wǎng)格文件名稱；第二項(xiàng):輸出處理后的小網(wǎng)格文件名稱）一下一步一點(diǎn)擊“Execute”按鈕一完成。經(jīng)過上述兩個(gè)步驟就完成了對(duì)嵌套網(wǎng)格的設(shè)置。對(duì)于嵌套兩個(gè)同級(jí)網(wǎng)格（如下圖所示）的操作說明：將“小大”網(wǎng)格和“小”網(wǎng)格按上述兩步處理后，得到“大處理”和“小1-處理”，后在將“大-處理”和“小2”網(wǎng)格按上述兩步處理后，得到“大-處理1”和“小2-處理”，則在水動(dòng)力模型計(jì)算時(shí)，輸入“人-處理1、小1-處理、小2

10、-處理”三個(gè)網(wǎng)格文件。如果三重嵌套，例如“小2”有嵌套在了“小1”中，則要兩、兩發(fā)生關(guān)系，既“大”和“小1”處理后的“小1-處理”在和“小2”處理，“小2-處理”在和“大-處理”進(jìn)行處理。多重嵌套以此類推。（2）圖形流程示意第一步：見圖1圖1第二步：見圖2精心整理精心整理精心整理精心整理第四步：選定計(jì)算區(qū)域，見圖第四步：選定計(jì)算區(qū)域，見圖6、圖7圖2第三步：加入背景文件，見圖3、圖4廣1輸入在顯示的區(qū)域內(nèi)添加,廣1輸入2.文件類型為:（*.xyz2.文件類型為:（*.xyz）圖4將陸地邊界和水深數(shù)據(jù)加入之后圖形顯示如圖5所示：圖5鎏MIPEZero-但讓1-MndifiedWorkArcoWi

11、ndowHdp5etCurrentCcntourLevel.,弋trl+L1仆圓逸嚼IBathymetryMariiernent.Ctrl+MBdd：.groiridManageiriStExport.、選擇該過項(xiàng)，編輯計(jì)算Resize.11區(qū)區(qū)域Settings.5etFont.ShowBackgroundIrrages星日LEdtVew意BIT5.指定計(jì)算區(qū)域的轉(zhuǎn)角（與北向夾角）2.指定計(jì)算區(qū)域左下角在顯示區(qū)域中的坐標(biāo)3.指定計(jì)算網(wǎng)格的寬度1.新建一個(gè)計(jì)算區(qū)域4.指定計(jì)算區(qū)域在X、Y方向的結(jié)點(diǎn)數(shù)。6.陸地區(qū)域高程值（要和陸地邊界文件一致）J圖7經(jīng)過上述設(shè)定后既圈定了計(jì)算區(qū)域，見圖8該框圈定

12、區(qū)域既為計(jì)算區(qū)域圖8在圖形中框第五步：在圈定計(jì)算區(qū)域之后，要選定要插值的數(shù)據(jù)范圍。應(yīng)用工具欄中的按忸，取，選定后再點(diǎn)擊一次按忸，之后在調(diào)出如圖9的對(duì)話框。在圖形中框加閉71口口0口。1ZZSOlCOtDRBdb&d笈工（MDg廣X1.選擇該項(xiàng)，對(duì)計(jì)算區(qū)域的數(shù)據(jù)進(jìn)行插值！三圖9三插值結(jié)束之后會(huì)出現(xiàn)消息框，該框顯示插值所用的時(shí)間信息。待全部ok之后，存盤?。ㄎ募麨閎atl）,此時(shí)系統(tǒng)自動(dòng)生成一個(gè)文件名為bat1_0的文件，該文件既為計(jì)算用的網(wǎng)格文件。然后打開該網(wǎng)格文件，如圖10所示。需要設(shè)陸地高程值的區(qū)域圖10第六步:網(wǎng)格數(shù)據(jù)調(diào)整?？捎霉ぞ邫谥腥鐖D11中所示的按鈕選定要編輯的網(wǎng)格，下的“5?！氨?/p>

13、”（見圖12、圖13），調(diào)整后如圖14所示，保存。n-在此設(shè)置選定網(wǎng)格的值。圖11圖12圖13圖14處理后的網(wǎng)格數(shù)據(jù)圖至此計(jì)算區(qū)域的網(wǎng)格數(shù)據(jù)文件已經(jīng)整理完畢。如果要在本計(jì)算區(qū)域內(nèi)內(nèi)嵌網(wǎng)格，操作步驟如下：第一步：建立子網(wǎng)格精心整理精心整理圖圖21圖22設(shè)置子網(wǎng)格名稱圖15圖16設(shè)置子網(wǎng)格名稱選擇被嵌入的大網(wǎng)格文件圖18圖17設(shè)置嵌套的小網(wǎng)格中計(jì)算的起始和結(jié)束時(shí)間，默認(rèn)時(shí)間同上一級(jí)網(wǎng)格的計(jì)算時(shí)間圖19圖20（該項(xiàng)暫不詳）指定X、Y向的小網(wǎng)格結(jié)點(diǎn)數(shù):結(jié)點(diǎn)數(shù)=3義大網(wǎng)格數(shù)+1小網(wǎng)格寬度（網(wǎng)格寬度必須為上一級(jí)網(wǎng)格寬度的1/3）指定小網(wǎng)格一左下角在大一網(wǎng)格中的位置選擇該項(xiàng)，指明小網(wǎng)格中包括陸地和水深嵌套網(wǎng)

14、格在大網(wǎng)格中的轉(zhuǎn)角（與乂方向）精心整理精心整理圖24圖23輸出編輯好的小網(wǎng)格文件至此已經(jīng)完成了小網(wǎng)格的編輯，小網(wǎng)格的文件名為BatLO(qiantao),圖25所示為嵌套后的小網(wǎng)格。圖25嵌套好的計(jì)算區(qū)域第二步：子網(wǎng)格的邊界處理圖26圖27輸入大網(wǎng)格文件名輸入小網(wǎng)格文件名輸入大網(wǎng)格文件名輸入小網(wǎng)格文件名輸入小網(wǎng)格左下角在大網(wǎng)格中的位置圖29圖28=設(shè)定都結(jié)束后，運(yùn)行完成圖30圖31至此已經(jīng)完成了嵌套網(wǎng)格的準(zhǔn)備工作，也產(chǎn)生了兩個(gè)用于模型計(jì)算的網(wǎng)格數(shù)據(jù)文件，在本例中分別為：大網(wǎng)格文件“8211_09點(diǎn)11）”，小網(wǎng)格文件“8211_0口120）g皿1”。2.邊界條件的設(shè)置（1）命令操作流程第一步：

15、設(shè)置時(shí)間序列文件（*4$0）打開mike21-filenew選擇mikezero中的“Timeseries”（時(shí)間序列）選擇“FromAsciiFile”出現(xiàn)“TimeseriesEidetor:importfromAs對(duì)對(duì)話框在此對(duì)話框中選擇“FileNam”（本例為WL1.txt-將“Timediscription”項(xiàng)設(shè)置成“EquidistantCalendarAx（等時(shí)間序列軸）ok出現(xiàn)水位邊界條件“WL1.TXT”的過程線，此時(shí)鼠標(biāo)右鍵單擊此圖一選擇“properties”，出現(xiàn)“Fileproperties”對(duì)話框?qū)xisType設(shè)置成“EquidistantcalendarAx

16、is；E“Iteminformation”中“Type”項(xiàng)中的“undefined”設(shè)置成“waterlevel”（水位）一ok保存文件（本例為Ts1.dfs0）。重復(fù)上述的操作，直至將上下兩個(gè)邊界的4個(gè)點(diǎn)時(shí)間序列文件處理完畢。驗(yàn)證文件也按上述操作過程設(shè)置。第二步：建立邊界條件打開mike21file-new選擇“profileseries”（斷面序列）選擇“blankTimeseries”（空白時(shí)間序列）一此時(shí)彈出“Fileproperties”對(duì)話框一設(shè)置“title”（上邊界條件名稱，本例為shangbianjie）;設(shè)置AxisType為“EquidistantcalendarAxis

17、”；設(shè)置Timestep為1800s（因?yàn)檫吔鐥l件中的時(shí)間間隔為30分鐘）；設(shè)置No.ofpoints為13（本例中，上邊界的網(wǎng)格數(shù)為13個(gè)）；設(shè)置Gridstep（網(wǎng)格步長）為900;將iteminformation中的Type項(xiàng)由undefined該為waterlevelfok此時(shí)出現(xiàn)一個(gè)空白的斷面邊界條件編輯器一打開北邊界的左側(cè)的點(diǎn)時(shí)間文件TS1.dfs0,并將其復(fù)制拷貝到空白單元格的第一列，打開右側(cè)的點(diǎn)時(shí)間序列文件TS2.dfs0,并將其拷貝復(fù)制到空白單元格的最后一列一選擇“Tools”菜單中的“interpolation”進(jìn)行插值一然后存盤（本例文件名為shangbianjie）。至

18、此已經(jīng)完成了上邊界的準(zhǔn)備，重復(fù)上述操作過程設(shè)置下邊界條件（本例中下邊界條件文件名為xiabianjie）。（2）圖形流程示意第一步：生成*在$0文件精心整理精心整理圖圖34圖32圖33選定要處理的文件選定該項(xiàng)，設(shè)置成等時(shí)間序列軸精心整理精心整理2oomnZoornOutFteWQU江口口mNe)ctZoom將此項(xiàng)設(shè)置成waterlevel將此項(xiàng)設(shè)置成waterlevelGGrid，擊一3Legend，擊一身哪圖的Pr.opertie!皿印加選定此項(xiàng)ForiL.；,圖35經(jīng)過上述設(shè)置后，將文件存盤，重復(fù)設(shè)置直至將所有的邊界處的點(diǎn)文件設(shè)置完畢，設(shè)置后的文件后綴名為.dfs0。第二步：建立邊界條件文

19、件圖36圖37該項(xiàng)設(shè)置成:waterlevel設(shè)置上邊界文件名稱設(shè)置成等時(shí)間序列軸上邊界的網(wǎng)格數(shù)網(wǎng)格步長該項(xiàng)設(shè)置成:waterlevel設(shè)置上邊界文件名稱設(shè)置成等時(shí)間序列軸上邊界的網(wǎng)格數(shù)網(wǎng)格步長圖38圖39建立好空白文檔后，打開左邊界點(diǎn)的.dfs0文件，并將其復(fù)制到空白文檔的第一列，打開右邊界點(diǎn)的.dfs0文件，并將其復(fù)制到空白文檔的最后一列，然后按圖40、圖41所示的操作進(jìn)行插值。圖40選定該項(xiàng)其他各項(xiàng)暫不詳圖41設(shè)置好后存盤，本例的文件名為shangbianjie。然后重復(fù)上述的操作設(shè)置下邊界條件，文件名為xiabianjie。至此，邊界條件已經(jīng)建立完畢。.水動(dòng)力模型的設(shè)置和計(jì)算圖形流程示

20、意精心整理精心整理圖圖43圖42水動(dòng)力模型水動(dòng)力及輸移擴(kuò)散模型水動(dòng)力及泥沙模型精心整理精心整理圖圖 圖58設(shè)置結(jié)束后，選定該選項(xiàng)運(yùn)行模型精心整理精心整理圖60至此，模型設(shè)置及運(yùn)行全部結(jié)束!.計(jì)算數(shù)據(jù)提取（驗(yàn)證）（1）命令操作流程第一步：提取數(shù)據(jù)打開mike2Jfilenew選擇mikezero中的mikezerotoolbo展開“Extractio”f選擇“Timeseriesform2Dfiles”（從二維模型文件中提取時(shí)間序列）一出現(xiàn)setupname（設(shè)置提取操作的名稱）一下一步一出現(xiàn)inputdate,指定數(shù)據(jù)的來源一下一步一出現(xiàn)subseriesspecification,設(shè)置提取數(shù)

21、據(jù)的時(shí)間范圍f下一步一出現(xiàn)selectitems,選擇要提取的內(nèi)容一下一步一出現(xiàn)pointselect,設(shè)置提取點(diǎn)的數(shù)目和位置f下一步一出現(xiàn)outputfileselection,設(shè)置要輸出的數(shù)據(jù)文件的名稱及標(biāo)題一下一步f運(yùn)行Execute,完成。第二步：對(duì)比計(jì)算結(jié)果及驗(yàn)證資料filenewf選擇mikezero中的“plotcomposer”（圖形設(shè)計(jì)）fok出現(xiàn)一個(gè)空白文檔，選擇菜單中的plot項(xiàng)中“insertnewplotobject”（插入新的圖形文件）一選擇“timeseriesplot”（時(shí)間序列圖形）fokf點(diǎn)擊按鈕（newitem）f出現(xiàn)打開對(duì)話框，選定要添加的文件后，在該對(duì)

22、話框下部的selectitem中設(shè)置要繪制的內(nèi)容fok-確定。對(duì)于鍵可重復(fù)操作，添加多個(gè)圖形文件的輸出，也可按如下操作：鼠標(biāo)右鍵單擊圖形一出現(xiàn)菜單后，選擇“properties”選項(xiàng)f在按上述方法添加。此時(shí)圖形的對(duì)比操作已經(jīng)結(jié)束，至于如何調(diào)試，仁者見仁，智者見智！（2）圖形流程示意圖61圖62設(shè)置提取操作的名稱選擇要被提取的目標(biāo)文件設(shè)置被提取數(shù)據(jù)的時(shí)間段設(shè)置被提取數(shù)據(jù)的時(shí)間段圖63圖64設(shè)置要被提取的內(nèi)容圖65圖66匚設(shè)置要被提取的點(diǎn)數(shù)指定被提取的點(diǎn)的位置設(shè)置輸出文件名和位置圖圖72精心整理圖67圖68GmeriQitedTime?cties:圖69插入新的圖形對(duì)象圖70圖71選擇該項(xiàng)：時(shí)間序

23、列圖形精心整理精心整理SeriesPlotPt1卿旗俄xT-AxiiCurmT-txLAniiE-tLi.aai至蛔困m|律審沮不性靈TLQUD五件名（BLJ1QUUIAtfTiti打開網(wǎng)上部.居F+向CT應(yīng)b111Ml工&1ML115T-ifi-sdInfe-IIL-aniInf.!uaD-Tirssri*i.xd.e-3.勢鎮(zhèn)白stvDirrriiTT：j-dfs67.田生2.選擇已經(jīng)提取的數(shù)據(jù)文件1.添加3.選定圖形中要顯示的內(nèi)容TlS圖73該圖中顯示出所設(shè)置的時(shí)間序列過程線！圖74圖75圖76所示為計(jì)算與實(shí)測值的對(duì)比情況。0.2-0.01-0.2T-口白-0.01-0.8-DO:CO12-04OD;0O12-0600:0012-1000;M12-08P(43D0.33.LiD):SurfaceBevalionmNEReese:施teLeMmm選擇該項(xiàng)，可做圖形輸出?cofla?1993-12-02i-尸.-yi-Help.SavetoBitmap.Font.CopytoCSpbcerdiavetoMetaF