EFDC操作手冊(cè)資料

EFDC操作手冊(cè)

?綜述

通過(guò)了解EFDC-EPA的全局及對(duì)EFDC模型的全過(guò)程有一個(gè)大概的了解。

?1介紹

?起源

VirhiniaInstituteofMarineScience(Hamrick)

?特性

?矩形或正交曲線(xiàn)網(wǎng)格

水平方向采用

?Sigma垂向坐標(biāo)

垂向

?歐拉運(yùn)輸方程

?拉格朗日粒子

可以模擬拉格朗日的粒子

?干濕網(wǎng)格

模擬淺灘的情況

?閘壩控制

?水生植被

?與其他模型耦合

?數(shù)值計(jì)算

?有限差分方法

?二階精度的有限差分空間度散

?二階精度、三層時(shí)間步長(zhǎng)的時(shí)間差分

?前處理

模型有前處理的網(wǎng)格程序，能夠生成網(wǎng)格及深度的插值計(jì)算。

?網(wǎng)格工具

?深度插值工具

?初始條件

?輸入格式采用文本格式

?運(yùn)行

?實(shí)時(shí)的演示

?診斷文件

?重啟文件

用于熱啟動(dòng)，就是不用從初始計(jì)算重新開(kāi)始計(jì)算，相當(dāng)于半中間的直接啟動(dòng)

?后處理

?時(shí)間序列

?空間分布圖

?矢量圖

?等值線(xiàn)

?剖面圖

?2全局概覽

?運(yùn)行程序

翁要介紹了efdc代碼的執(zhí)行程序的組成部分，這部分內(nèi)容后面會(huì)詳細(xì)敘述.

?efdc.for

源代碼文件，主要代碼文件

?

通用聲明變量文件，一般模型使用

?efdc.par

特別的用于最小化內(nèi)存占用管理

?模型配置文件（輸入文件）

模型通過(guò)后綴名為inp的文本文件來(lái)進(jìn)行輸入文件的配置,按照功能分為5中類(lèi)型。

?空間網(wǎng)格地形組

用于構(gòu)成水休模型的空間網(wǎng)格及地形條件。里面有一些概念可能會(huì)難以理解，簡(jiǎn)單的衡量

空間為什么會(huì)有這么多文件，比如網(wǎng)格方向的概念，對(duì)于圖形界面來(lái)說(shuō)，很難理解，這些

內(nèi)容我自己覺(jué)得是使用字符來(lái)表達(dá)空間所必須的。我們看到的空間實(shí)體，采用文本文件來(lái)

抽象表示，本身就很復(fù)雜。為什么一圖勝前言，信息傳遞給計(jì)算機(jī)能夠讀懂的語(yǔ)言本來(lái)就

比較麻煩。后續(xù)章節(jié)到網(wǎng)格時(shí)會(huì)詳細(xì)描述。

?cell

將水體外部邊界線(xiàn)數(shù)字化，在字符時(shí)代，用于不同的數(shù)字來(lái)代表水體和陸地界面，從而

刻畫(huà)出水體的形狀。

?celllt

cell文件的補(bǔ)充文件，通常為cell文件中不包括入流和出流口的文件。

?dxdy

網(wǎng)格的幾何特性

?空間網(wǎng)格的分辨率

水平方向上每個(gè)網(wǎng)格的大小間距

?深度

這里的深度是初始水深。水深：即水面到底部地形的距離，是水位的表征量。在計(jì)

算過(guò)程中隨著水位變動(dòng)。而初始水深代表了初始水位。

?底部高程

垂直方向的上的網(wǎng)格的高程（一般是網(wǎng)格中心點(diǎn)計(jì)），這里代表實(shí)際的地形。和上

述概念不同的是,這里代表實(shí)際的地形，一般而言，其在計(jì)算過(guò)程中不會(huì)改變。

?底部粗糙率

?水生植被

考慮水底的植被作用，影響水動(dòng)力作用。需要在主控文件中開(kāi)啟植被阻隔。

?Ixly

網(wǎng)格的地理特性

?網(wǎng)格的中心坐標(biāo)

網(wǎng)格中心點(diǎn)的坐標(biāo),EFDC一般采用的是UTM坐標(biāo)系。

?網(wǎng)格的方向

網(wǎng)格的方向，后續(xù)計(jì)算流向時(shí)的轉(zhuǎn)換。

?depth

僅僅用于笛卡爾網(wǎng)格來(lái)指定水深、粗糙度等。為笛卡爾網(wǎng)格單獨(dú)開(kāi)了一道門(mén)。

?mask

阻擋物，相當(dāng)于固定邊界，只不過(guò)其可以設(shè)定在水體中心，比如湖中的島可以采用這個(gè)

來(lái)描述。但需要說(shuō)明的是為了計(jì)算方便，島嶼通常采用挖空網(wǎng)格來(lái)表示。

?主控及運(yùn)行控制組

模型的總指揮文件，控制一些模塊的開(kāi)關(guān)及參數(shù)的設(shè)定。

?efdc

該文件是所有控制文件的核心文件，管控了幾乎所有的參數(shù)及模塊的開(kāi)關(guān)。

?模塊的開(kāi)關(guān)

如是否開(kāi)啟水動(dòng)力模塊，是否模擬染料。

?邊界條件位置及其序列數(shù)量

對(duì)于河道來(lái)說(shuō)，支流入干流位置點(diǎn)，支流數(shù)量等。

?控制輸出格式

是否生成流速等指標(biāo)

?show

在運(yùn)行過(guò)程中顯示的指標(biāo)及信息（提供計(jì)算時(shí)的實(shí)時(shí)演示）

?時(shí)間序列邊界條件組

?aser

氣象條件的時(shí)間序列文件，包括氣壓、溫度、濕度等

?qser

流量時(shí)間序列文件

?dser

染料序列文件

?pser

開(kāi)邊界水位序列文件

?sser

鹽度時(shí)間序列

?物理過(guò)程組

?qctl

水工構(gòu)筑物控制文件

?vege

植被阻礙流動(dòng)特征文件

?wave

表面重力波文件

?初始條件及啟動(dòng)組

模型啟動(dòng)時(shí)的初始條件的設(shè)定

?dye

初始染料空間分布

?restart

熱啟動(dòng)指定文件

?salt

初始鹽度空間分布

?模型結(jié)果文件（輸出文件）

一般以后綴名為.out的文件輸巴

?診斷輸出文件

?重啟文件

?網(wǎng)格點(diǎn)的時(shí)間序列文件

指標(biāo)隨時(shí)間變化的值

?最小平方諧波分析文件

?二維及三維圖形可視化文件

?前處理

前處理,及構(gòu)建模型及參數(shù)的設(shè)置

?3網(wǎng)格

模型構(gòu)建的第一步是生成研究區(qū)的網(wǎng)格文件。

?基本概念

?四邊形網(wǎng)格（三角網(wǎng)格）

EFDC采用四邊形及三角形網(wǎng)格（為四邊形網(wǎng)格的對(duì)角線(xiàn)分割后得到）

?坐標(biāo)系xy->ij

網(wǎng)格的編號(hào)為ij值，對(duì)應(yīng)直角坐標(biāo)系中的ij°x的正方向?yàn)闁|，y方向?yàn)楸狈较颉?/p>

?EFDC采用相對(duì)坐標(biāo)系

地理中的實(shí)際的正東方向用x"正北方向用y\

?國(guó)家坐標(biāo)系

?UTM

需要注意的是，EFDC采用的是國(guó)際標(biāo)準(zhǔn)單位，MKS,所以推薦使用UTM坐標(biāo)系。

?用戶(hù)自定義

?網(wǎng)格文件

?CELL文件

使用數(shù)字陣列表來(lái)代表網(wǎng)格分布情況。

?頭文件

4行頭文件，為固定格式。

?坐標(biāo)軸區(qū)域

水體單元格上下部分為坐標(biāo)軸的標(biāo)識(shí)，i為x方向，j為y方向。

?水髀元格

中間整塊數(shù)字部分就是區(qū)域網(wǎng)格了，采用不同的數(shù)字代表不同的網(wǎng)格類(lèi)型，這里為了明

確表示水體邊界，是加入了陸地類(lèi)型的。邊界的概念為，水體和陸地的交界處。舉個(gè)例

子，windows的掃雷游戲玩過(guò)吧，那么多格子就代表網(wǎng)格，有雷的區(qū)域就用1表示，

那么一整塊雷就全部是L

?CELLT文件

可指定為CELL文件的子集,用于特殊的目的，如表示除入流口以外的水體。則原來(lái)的邊界

可能變?yōu)殛懙?，原?lái)的水體變?yōu)檫吔缁蜿懙亍?/p>

?DXDY

指定網(wǎng)格的幾何特定，單位為m,具體前面說(shuō)明前面已經(jīng)有了不再重復(fù)。

?空間網(wǎng)格的分辨率

水平方向上每個(gè)網(wǎng)格的大小間距

?深度

這里的深度是初始水深，水深：即水面到底部高程的距離，是水位的表征量。在計(jì)算過(guò)

程中隨著水位變動(dòng)。而初始水深代表了初始水位。水位;深度+底部高程

?底部高程

垂直方向的上的網(wǎng)格的高程（一般是網(wǎng)格中心點(diǎn)計(jì)），這里代表實(shí)際的地形。和上述概

念不同的是，這里代表實(shí)際的地形，其在計(jì)算過(guò)程中不會(huì)改變，是地理量。

?底部粗糙率

?水生植被

考慮水底的植被作用，影響水動(dòng)力作用。需要在主控文件中開(kāi)啟植被阻隔。

?LXDY

?網(wǎng)格的中心坐標(biāo)（位置）

網(wǎng)格中心點(diǎn)的坐標(biāo)，EFDC一般采用的是UTM坐標(biāo)系，用于后處理時(shí)的輸出結(jié)果的可

視化，依照方便可設(shè)定為m,KM等。

?旋轉(zhuǎn)矩陣（方向）

網(wǎng)格的方向，用于將網(wǎng)格的偽方向轉(zhuǎn)換為實(shí)際地理中的方向。

?速度方向

結(jié)果文件中的流向轉(zhuǎn)換為真正的流向。（地理坐標(biāo)）

?風(fēng)場(chǎng)方向

實(shí)際中的風(fēng)場(chǎng)利用逆矩陣轉(zhuǎn)換為模型中的風(fēng)向。

?網(wǎng)格生成gefdc.f

實(shí)際操作中網(wǎng)格的生成是采用gefdc.f來(lái)生成的。

?準(zhǔn)備文件

cell

前面已經(jīng)敘述，數(shù)字輪彝

?depdat

指定深度及底部高程。xyz格式，xy為坐標(biāo)，z為值。

?gcell

曲線(xiàn)網(wǎng)格時(shí)進(jìn)行輔助四邊形網(wǎng)格文件同cell類(lèi)似格式

?gridext

水體網(wǎng)格角坐標(biāo)文件（左下角），三角形網(wǎng)格也必須等同四邊形設(shè)定。在設(shè)定參數(shù)

NTYPE=O時(shí)采用。

?vege

同depdat設(shè)定植被

?zrouth

同上

?gefdc設(shè)定文件

網(wǎng)格生成工具的主控文件，總共13個(gè)卡片組，起始efdc.inp文件也是efdc的主控文件，

在這里有異曲同工之妙，所以簡(jiǎn)單學(xué)習(xí)下為后面鋪墊。文件的格式是文本格式，每組以C*

開(kāi)頭，第一行為標(biāo)題，第二行為參數(shù)標(biāo)題，第三行為參數(shù)值。

?C2咔2

生成方式，必選

?NTYPE

生成網(wǎng)格的方式

?0

采用gridext和cell生成網(wǎng)格

?1-5

正交曲線(xiàn)網(wǎng)格，一般采用5

?6

不保證正交性，需要指定邊界點(diǎn)

?7

需要用戶(hù)指定函數(shù)，矩形區(qū)域網(wǎng)格

?8

采用cell和卡4的設(shè)定深度插值法生成網(wǎng)格

?9

深度插值法，但坐標(biāo)系和8不同

?NBPP(NTYPE=l-6)

邊界點(diǎn)數(shù)目

?IJ

IJ的基本信息，最大最小值等

?C3

在正交曲線(xiàn)網(wǎng)格時(shí)，保留一個(gè)規(guī)整的矩形網(wǎng)格備用，用于后期演示

?C4

網(wǎng)格坐標(biāo)信息

?C5

網(wǎng)格迭代選項(xiàng)

?C6

數(shù)值方面的參數(shù)

?C7

單位變換

?C8

插值方式

?C9-10

設(shè)定矩形網(wǎng)格的參數(shù)

?C11

d叩data的插值方式

?C12-13

邊界的設(shè)定

?網(wǎng)格生成VOGG

visualorthogonalgridgenerationsystem,正交曲線(xiàn)網(wǎng)格生成工具

?正交曲線(xiàn)網(wǎng)格

正交曲線(xiàn)網(wǎng)格能夠貼合邊界，彌補(bǔ)了直角網(wǎng)格在復(fù)雜邊界處的鋸齒狀網(wǎng)格，從而減少了網(wǎng)

格數(shù)量，

?VOGG使用

?生成思路

?概化模型

首先確定研究區(qū)域的邊界，也就是水體和陸地的邊界區(qū)域。這部分需要將其概化，

依據(jù)研究目標(biāo)和實(shí)際區(qū)域形狀進(jìn)行概化。

?邊界線(xiàn)應(yīng)該相對(duì)平滑（整體接近矩形塊）

?邊界區(qū)域的角盡量為直角

?生成計(jì)算區(qū)域-矩形網(wǎng)格

?節(jié)點(diǎn)（邊角點(diǎn)）-角度

?線(xiàn)（點(diǎn)和點(diǎn)之間的線(xiàn)）-長(zhǎng)度

?網(wǎng)格數(shù)量梳理

手動(dòng)或者自動(dòng)定義網(wǎng)格數(shù)量，生成U值

?網(wǎng)格生成-由矩形轉(zhuǎn)換為正交曲線(xiàn)網(wǎng)格

?源代碼介紹

?cell.inp

同efdc所需文件類(lèi)似，但其僅僅只有0,9,5,后期生成后再進(jìn)行修改為三角網(wǎng)

格

?gridside

網(wǎng)格點(diǎn)或邊數(shù)量控制，坐標(biāo)位置輸入

?Nside

節(jié)點(diǎn)數(shù)量或者邊數(shù)量

?Ijndex

左下角i的標(biāo)號(hào)

?jjndex

左下角J的標(biāo)號(hào)

?inter

網(wǎng)格插值分辨率

?vogg

控制網(wǎng)格生成參數(shù)及地形插值

?可視化操作

?界面簡(jiǎn)介

?左側(cè)工具欄

工具排列的順序，基本就是網(wǎng)格制作的流程。

?線(xiàn)控制

用于草圖繪制，繪制水體邊界的草圖

?點(diǎn)控制

精確修改草圖的位置點(diǎn)，通過(guò)節(jié)點(diǎn)調(diào)整草圖吻合實(shí)體

?邊界編輯

修改節(jié)點(diǎn)的方向，實(shí)質(zhì)是調(diào)整網(wǎng)格線(xiàn)的走勢(shì)，全部必須是單一方向，一般

按順時(shí)針。

?U排序工具

這里起始就是網(wǎng)格生成前一步，將節(jié)點(diǎn)轉(zhuǎn)換為網(wǎng)格節(jié)點(diǎn)

?網(wǎng)格工具

調(diào)整網(wǎng)格，刪除

?上側(cè)工具欄

?視圖控制

移動(dòng)操作區(qū)域的位置方向

?生成網(wǎng)格

也用于導(dǎo)入生成好的網(wǎng)格

?菜單欄

全局控制及導(dǎo)入文件的菜單

?導(dǎo)入水體邊界（shp）

?導(dǎo)入草圖邊界（shp）

?導(dǎo)入網(wǎng)格文件（shp）

?導(dǎo)入坐標(biāo)文件

?網(wǎng)格生成步驟

?底圖準(zhǔn)備

依照程序需要將底圖準(zhǔn)備好

?導(dǎo)入邊界

水體邊界的底圖導(dǎo)入，一般我們要將水體作為底圖放入網(wǎng)格生成器中

?草圖繪制

簡(jiǎn)略的框架繪制，繪制為一個(gè)矩形框架，控制網(wǎng)格的邊界，不同的網(wǎng)格生成程

序中的操作不同，但均要保證網(wǎng)格的邊界平滑，能夠在宏觀上貼合水體邊界

?修改節(jié)點(diǎn)屬性信息

?生成網(wǎng)格

網(wǎng)格在邊界基礎(chǔ)上生成，通過(guò)控制參數(shù)來(lái)控制網(wǎng)格的分辨率

?調(diào)整網(wǎng)格

網(wǎng)格不可能一次就能夠完美，需要進(jìn)行不斷微調(diào)，因?yàn)榫W(wǎng)格后期基本不會(huì)改動(dòng)，

所以一般在網(wǎng)格生成階段就需要花費(fèi)較多精力來(lái)調(diào)整到足夠好。

?地形插值

一般網(wǎng)格的工具都會(huì)附帶地形工具，也就是給網(wǎng)格垂直方向深度賦值。

?4主控文件

efdc.inp文件

?概述

參數(shù)一般是通過(guò)整型變量設(shè)定,0即為關(guān)閉狀態(tài)，非0為打開(kāi)狀態(tài)。

?運(yùn)行控制參數(shù)

控制模型運(yùn)行時(shí)的開(kāi)關(guān)，比如是否開(kāi)啟水動(dòng)力模塊。

?輸出控制參數(shù)

輸出結(jié)果的內(nèi)容控制，控制模型輸出什么樣的結(jié)果文件。

?物理過(guò)程參數(shù)

對(duì)于水動(dòng)力、水質(zhì)模型本身的一些參數(shù)如擴(kuò)散系數(shù)等。

?卡組說(shuō)明

對(duì)于卡組中的很多值并非都需要很熟悉，所以?xún)H僅先說(shuō)明一些我比較熟悉的值，后續(xù)對(duì)其他進(jìn)

行補(bǔ)充。

?卡1

模型運(yùn)行的標(biāo)題文件

?卡1A

設(shè)定有關(guān)空間和時(shí)間計(jì)算的參數(shù)

?卡2

重啟文件，通用設(shè)置，診斷切換

?重新啟動(dòng)文件參數(shù)hiiD：///China3S/D/3192807.hlml

重新啟動(dòng)文件的使用是用于熱啟動(dòng)，也就是說(shuō)在運(yùn)行后，下次運(yùn)行時(shí)則從某一時(shí)刻開(kāi)始，

而不重新從最初的狀態(tài)開(kāi)始,這里是為了跳過(guò)冷啟動(dòng)的時(shí)間。在EFDC的源代碼中，

有關(guān)熱啟動(dòng)的相關(guān)代碼包括：RESTOUT.for.RESTINlO.for、RESTIN1和

RESTMOD.for,其中RESTOUT.for程序處理寫(xiě)熱啟動(dòng)文件RESTART.OUT和

RSTWD.OUT文件，而RESTINl.for分別處理熱啟動(dòng)文件RESTART.inp和

RSTWD.inp的讀入,RESTINlO.for只讀取RESTART.inp,不讀入RSTWD.inp文件，

且一般的水動(dòng)力模擬不調(diào)用RESTINlO.for,只調(diào)用RESTINl.for。RESTMOD.for只

有HDMT和HDMT2T.for兩個(gè)子程序才會(huì)調(diào)用，一般只模擬水動(dòng)力不會(huì)調(diào)用該子程序。

RESTIN2.for也會(huì)涉及RESTART.inp分類(lèi)：EFDC模型EFDC熱啟動(dòng)是在運(yùn)行EFDC經(jīng)

常用到的方法，例如，在運(yùn)行EFDC出現(xiàn)意外中斷時(shí)，熱啟動(dòng)可以幫助用戶(hù)繼續(xù)從中斷

處運(yùn)行EFDC,而不需要從頭開(kāi)始運(yùn)行。但熱啟動(dòng)經(jīng)常會(huì)碰到負(fù)值問(wèn)題而使EFDC無(wú)法

運(yùn)行下去，這種問(wèn)題的解決辦法是：用Uedit32對(duì)RestartJnp文件進(jìn)行修改，將-0.

替換為0.即可。

?ISRESTI

重啟文件的讀取

?1從inp文件中讀取初始條件

?-1和上述一致，但根據(jù)底部高程進(jìn)行調(diào)整

?2對(duì)于多層的處理，當(dāng)新KC大于4時(shí)，可以從KC/2層中初始化變?yōu)镵C層

?10舊版本中讀??？

?ISRESTO

重啟文件的生成

?-1運(yùn)行結(jié)束后生成重啟文件

?N每N個(gè)基準(zhǔn)參考時(shí)間段生成重啟文件

?多線(xiàn)程運(yùn)算

?輸出到屏幕的內(nèi)容

?寫(xiě)出一些平衡值到.out輸出文件中

?卡3

一些數(shù)值計(jì)算的方法及迭代次數(shù)的設(shè)定

?卡4

對(duì)對(duì)流擴(kuò)散方程進(jìn)行自定義

?卡5

?水平動(dòng)量擴(kuò)散設(shè)定

?WASP輸出設(shè)定

?干濕模式設(shè)定

?關(guān)閉干濕

?常量濕深度非線(xiàn)性迭代

?變量濕深度非線(xiàn)性迭代

?常量深度

?常量迭代存在mask

?使用網(wǎng)格面的變量干濕模式

?使用mask的常量干濕模式

?植被阻流作用

?渦旋強(qiáng)度公式

?波浪動(dòng)力學(xué)

?C6

物質(zhì)的垂向?qū)α骱蛿U(kuò)散以及物質(zhì)的源匯數(shù)量域。

?每組變量用一個(gè)數(shù)值代替，且有各自的參數(shù)

?ISTRAN是擴(kuò)散和源匯項(xiàng)的開(kāi)關(guān)

?ISTOPT

選擇指標(biāo)的初始條件的設(shè)定方式，以及其的參與方式

?ISADAC

數(shù)值計(jì)算方式

?C7#時(shí)間參數(shù)

時(shí)間設(shè)定參數(shù)整數(shù)變量，通過(guò)和C8配合使用來(lái)設(shè)定時(shí)間參數(shù)。

?NTC

efdc采用設(shè)定迭代倍數(shù)來(lái)設(shè)定運(yùn)行時(shí)常，也就是說(shuō)其設(shè)定一個(gè)運(yùn)行的時(shí)間單位，這里

設(shè)定次數(shù)。

?NLTC

時(shí)間步長(zhǎng)，同上面，是通過(guò)一個(gè)參照值來(lái)設(shè)定時(shí)間步長(zhǎng)的。

?變化的步長(zhǎng)和固定步長(zhǎng)

?輸出的頻率

?C8#時(shí)間參數(shù)

實(shí)際參數(shù)

?TCON

開(kāi)始時(shí)間的轉(zhuǎn)換系數(shù)，當(dāng)相對(duì)標(biāo)準(zhǔn)時(shí)間不為秒時(shí)，要進(jìn)行轉(zhuǎn)換。

?TBEGIN

模型啟動(dòng)時(shí)刻，通過(guò)TCON來(lái)轉(zhuǎn)換為秒。

?TREF

模型的相對(duì)標(biāo)準(zhǔn)時(shí)間,相當(dāng)于一個(gè)基準(zhǔn)時(shí)間，卡7中的變量都是參照此進(jìn)行計(jì)算的。

?科氏力設(shè)定

?變步長(zhǎng)的設(shè)定

?C9#空間參數(shù)

水平空間參數(shù)設(shè)定

?水平網(wǎng)格參數(shù)

設(shè)定水平網(wǎng)格的數(shù)量及結(jié)構(gòu)

?IC\JC

設(shè)定I\J的方向的網(wǎng)格數(shù)量

?LC\LVC

由于計(jì)算方便，EFDC采用了一個(gè)數(shù)來(lái)代替（i,j）數(shù)組，將其轉(zhuǎn)換為L(zhǎng),但為了循

環(huán)方便，其增加了L=1和L=LC,完全是程序需要。

?ISCLO

正交曲線(xiàn)網(wǎng)格

?并行計(jì)算的網(wǎng)格區(qū)域#并行計(jì)算

有關(guān)并行計(jì)算的設(shè)定

?ISMASK

設(shè)定不透水網(wǎng)格面的開(kāi)關(guān)，不透水網(wǎng)格就是可以指定網(wǎng)格某一面是不透水的。

?ISPGNS#開(kāi)邊界設(shè)定

設(shè)定北部和南部開(kāi)邊界中的周期性信息

?初始條件的空間插值

光滑參數(shù)，初始條件的插值是efdc執(zhí)行的，可以設(shè)定其的插值參數(shù)。

?C9a

垂向空間參數(shù)設(shè)定

?垂向?qū)訑?shù)設(shè)定KC

?CIO

設(shè)定具體的sigma層數(shù)比例分配

?K層

K是從下向上計(jì)算,多層時(shí),K=1為底層

?DZC層比例，總和必須為1

?C11

網(wǎng)格、糙率、深度值參數(shù)

?矩形網(wǎng)格的固定參數(shù)

?糙率的調(diào)整轉(zhuǎn)換

?底部深度的調(diào)整轉(zhuǎn)換

?干濕深度設(shè)定

?C11A

動(dòng)量修正系數(shù)設(shè)定

?C11B

角網(wǎng)格修正參數(shù)

?C12

渦流及湍流的相關(guān)擴(kuò)散系數(shù)

?AHO

水平動(dòng)量和質(zhì)量擴(kuò)散系數(shù)，常量值

?AHD

無(wú)量綱水平動(dòng)量擴(kuò)散系數(shù)，變量的附加在常量AHO上,Smagorinsky亞尺度水平擴(kuò)

散公式。

?AVO

背景，常量或分子動(dòng)力粘度

?ABO

背景，常量或分子擴(kuò)散系數(shù)

?AVMN

最小渦流粘度

?ABMN

最小渦流擴(kuò)散系數(shù)

?VISMUD

恒定浮泥粘度

?AVBCON

設(shè)定為0時(shí)，渦流未被激活，采用ABO和AVO設(shè)定為垂向擴(kuò)散系數(shù)和粘度系數(shù),設(shè)

定為1時(shí)采用模型計(jì)算值

?ZBRWALL

邊壁對(duì)數(shù)粗糙度，在水平動(dòng)量擴(kuò)散系數(shù)激活后，并且AHD和AHO不為0.

?C12A

?C13

兩個(gè)應(yīng)該是高級(jí)應(yīng)用時(shí)采用的，沒(méi)有詳細(xì)說(shuō)明。略

?C14

潮汐、大氣強(qiáng)迫條件、地下水和渠道子網(wǎng)格的參數(shù)設(shè)定

?潮汐強(qiáng)迫條件數(shù)目

?風(fēng)序列數(shù)目

序列代表了時(shí)間和值構(gòu)成了一個(gè)時(shí)間序列，如流量序列。

?大氣序列數(shù)目

?地下水設(shè)定

?土壤濕度平衡

?地下水和底床交互

?激活河道子網(wǎng)格模式

?波浪文件選擇

?潮汐高程同化#疑問(wèn)

?對(duì)于超出干網(wǎng)格的水開(kāi)啟滲透或消除作用

?C15

潮汐調(diào)和設(shè)定符號(hào)設(shè)定

?C16

水位邊界條件和壓力邊界條件參數(shù)設(shè)置

?水位邊界條件或開(kāi)邊界數(shù)目（東南西北）

?潮汐調(diào)和邊界數(shù)目

?強(qiáng)迫類(lèi)型

?邊界時(shí)間序列數(shù)目

?C17

周期水位邊界及壓力邊界具體設(shè)定

?編號(hào)

?振幅

?相位

?C18

南部開(kāi)邊界設(shè)定

?SB

?方式

?第一類(lèi)型開(kāi)邊界

?指定高程

?調(diào)和常數(shù)

?第二類(lèi)型

?RADIATION-SEPARATION輻射分離條件0切向速度

?RADIAHON-SEPARATION條件自由切向速度

?時(shí)間序列方式

?應(yīng)用潮汐調(diào)和常數(shù)

?C19

西部開(kāi)邊界設(shè)定

?C20

東部開(kāi)邊界設(shè)定

?C21

北部開(kāi)邊界設(shè)定

?C21A

水位和速度數(shù)據(jù)

?C22

模型讀取的時(shí)間序列的數(shù)量和污染物的種類(lèi)設(shè)定

?C23

速度、體積源匯項(xiàng)、流量控制、取水和退水邊界的設(shè)定

?體積源匯項(xiàng)的位置數(shù)目

常量或者變化的體積源匯項(xiàng)（如流量），對(duì)于正向入流是考慮污染物隨流的輸入，而出

流則考慮了物質(zhì)的損失，

?體積源匯項(xiàng)的時(shí)間序列數(shù)目

?壓力控制措施（閘壩）的位置數(shù)目

控制水位的水工構(gòu)筑物

?壓力控制措施（閘壩）的時(shí)間序列數(shù)目

?取水'回流邊界的位置數(shù)目

常量或者隨時(shí)間變化的取水和回水邊界，為電站和冷卻系統(tǒng)設(shè)計(jì)。

?取水'回流邊界的時(shí)間序列數(shù)目

?C24

體積源匯項(xiàng)的位置指定，濃度序列鏈接。

?位置設(shè)定

通過(guò)設(shè)定網(wǎng)格ij來(lái)確定邊界條件的位置

?為位置指定序列號(hào)

為該位置指定相應(yīng)的流量序號(hào)、污染物濃度序號(hào)，對(duì)應(yīng)dser中的組號(hào)（從上到下）

?比例因子QSFACTOR

設(shè)定比例

?恒定入流和出流設(shè)定

實(shí)際當(dāng)中并不直接使用，而采用時(shí)間序列的恒定格式。若設(shè)定，則是此值加上時(shí)間序列

值。

?側(cè)向入流，體積量由m2/s轉(zhuǎn)換到m3/s

?動(dòng)量通量的朝向指定

?C25-C26

恒定序列的恒量設(shè)定。

?C27-C31

射流和羽狀流的設(shè)置。

?C32

水位和壓力流的參數(shù)設(shè)定信息（包括閘壩水位設(shè)置）

?位置（上游和下游設(shè)定）

支持單向流和雙向流，同時(shí)水位控制類(lèi)則只需要設(shè)定上游，如溢流壩等，但是如果是有

控制的閘，則需要上下游同時(shí)設(shè)定。

?水流控制類(lèi)型

?水位-流量曲線(xiàn)設(shè)定

常用，一般用于河流的下游邊界，通過(guò)W-Q的關(guān)系來(lái)進(jìn)行下游邊界條件的控制。

如果有水文部門(mén)的公式則使用構(gòu)造水位和流量值，如果沒(méi)有，這可以用曼寧公式推

的大概的關(guān)系。

?上游水深壓差控制型

?潮汐入口加速流

?動(dòng)量通量設(shè)定

?動(dòng)量通量寬度

?C33

流量取水、回水邊界

?位置（上游和下游設(shè)定）

支持單向流和雙向流，司時(shí)水位控制類(lèi)則只需要設(shè)定上游，如溢流壩等

?取水回水序列編號(hào)

?恒定取水回水量值

?動(dòng)量通量設(shè)定

?動(dòng)量通量寬度

?C34-C35

取水回水水質(zhì)指標(biāo)及泥沙恒定濃度設(shè)定

?C36-C42

沉積物模塊設(shè)置

?C43

有毒污染物的初始條件和參數(shù)設(shè)定

?初始條件

?空間均勻濃度分布

?水柱中的空間變化的濃度分布

?在沉積物中空間變化的初始條件

?在水柱中和沉積物中空間變化濃度分布

?在沉積物中的單位表征

?水中有毒污染物的降解率及參考溫度

?底泥中的有毒污染物降解率及參考溫度

?C44

其他有毒污染物的參數(shù)

?溶解及顆粒態(tài)有機(jī)碳吸附作用

?水面蒸發(fā)率

?分子量對(duì)蒸發(fā)率的影響

?光分解速率

?太陽(yáng)能輻射影響

?孔隙水中的擴(kuò)散系數(shù)

?沉積物上層與水體間的擴(kuò)散系數(shù)（正負(fù)值不一樣）

?擴(kuò)散系數(shù)

?消散系數(shù)

?沉積物中的顆?；旌蠑U(kuò)散系數(shù)

?沉積物中的顆?；旌仙疃确秶?/p>

?C45

有毒污物的蹣設(shè)定

?C45A

有毒有機(jī)化合物的時(shí)間和空間量設(shè)定

?時(shí)間恒定

?空間均勻分布和空間變化的分布

?C45B

有毒有機(jī)化合物的設(shè)定

?C45C

有毒物顆粒態(tài)在水柱中的分布

?C45D

有毒物顆粒態(tài)在沉積層中的分布

?C46

開(kāi)邊界中的濃度序列數(shù)目設(shè)定。

?C47

南側(cè)開(kāi)邊界入流濃度指定時(shí)間序列序號(hào)

?C48

南側(cè)濃度邊界的底層恒定濃度

?染料

?TOX等

?C49

南側(cè)濃度邊界的底層恒定濃度

?黏性泥沙

?非黏性泥沙

?C50

南側(cè)濃度邊界的表層恒定濃度

?C51

南側(cè)濃度邊界的表層恒定濃度

?C52-C66

其他方向的開(kāi)邊界濃度設(shè)定

?C66A-C66B

數(shù)據(jù)同化，我的理解是將臨近的邊界條件進(jìn)行類(lèi)似設(shè)定。

?C67

拉格朗日粒子追蹤

?計(jì)算方法選擇

?點(diǎn)位置釋放

?釋放粒子數(shù)目

?粒子釋放間隔

?輸出時(shí)間步長(zhǎng)

?范圍內(nèi)釋放，MLRPDRT成組，并可計(jì)算全部釋放時(shí)間的平均值

?邊界區(qū)域設(shè)定

?可選高階軌跡積分方法

?輸出速度矢量控制

?C68

釋放點(diǎn)

?C69

笛卡爾網(wǎng)格中心經(jīng)緯度

?C70

二進(jìn)制文件及ASCII文件設(shè)定

?C71

水平等值線(xiàn)生成

?C71A

沉積層的等值線(xiàn)生成

?C71B

食物鏈模型控制

?C72

水平面的水位和壓力等值線(xiàn)

?C73

速度矢量的輸出

?C74

垂向剖面等值線(xiàn)

?C75

剖面單元格的數(shù)目

?C76

垂向剖面線(xiàn)的位置

?C77-79

垂向速度矢量

?C80-81

3d文件輸出

?C82

調(diào)和分析輸出

?C83

調(diào)和分析的位置切換

?C84

時(shí)間序列文件輸出

?輸出文件的覆蓋和追加

?輸出點(diǎn)位置數(shù)目

?開(kāi)始輸出的時(shí)間點(diǎn)

?結(jié)束輸出的時(shí)間點(diǎn)

?輸出間隔

?C85-C86

輸出模式模式化設(shè)定

?C87

輸出量位置及指標(biāo)確定

?檢出水位結(jié)果時(shí)間序列

?輸出濃度變量的時(shí)間序列

?渦流粘度和擴(kuò)散速率的時(shí)間序列

?外部模式水平速率的時(shí)間序列

?輸出外部模式水平傳輸率的時(shí)間序列

?每一層的水平速度的時(shí)間序列

?輸出外部模式凈體積源匯的時(shí)間序列

?輸出凈體積源/匯的時(shí)間序列

?備注

C88-90

垂向輸出設(shè)定，垂向采樣設(shè)定

?C91

沉降速度的范圍

?C92

所選網(wǎng)格的最小懸浮物濃度

?C93

空間變化的侵蝕函數(shù)

?C94

空間變化的食物鏈暴露深度

?5水動(dòng)力模塊

水動(dòng)力文件的附加文件，如流量、濃度時(shí)間序列，水位流量曲線(xiàn)等。

?1aser

模型必須，大氣及氣象條件輸入。

?使用參數(shù)：卡6的參數(shù)控制大氣驅(qū)動(dòng)的模式

?轉(zhuǎn)換系數(shù)

轉(zhuǎn)換系數(shù)是為了將原始數(shù)據(jù)統(tǒng)一進(jìn)行轉(zhuǎn)換為模型內(nèi)部使用，一般而言，為了保持?jǐn)?shù)據(jù)統(tǒng)一

性，不推薦修改這些值，而建議的方法是輸入之前就轉(zhuǎn)換好，這樣防止出錯(cuò)。就是在輸入

之前就處理好。

?序列數(shù)目

?時(shí)間換算為秒，86400代表后續(xù)使用時(shí)間為天

?數(shù)值統(tǒng)一加上一個(gè)值

?風(fēng)速轉(zhuǎn)換

?降雨、蒸發(fā)轉(zhuǎn)換

?數(shù)據(jù)

?風(fēng)速

?風(fēng)向

?氣壓

?相對(duì)濕度

?降雨

?蒸發(fā)

?短波輻射

?2濃度序列文件（#邊界條件）

輸入點(diǎn)源式的濃度序列文件

?使用參數(shù)

?ISTYP

序列類(lèi)型

?1（優(yōu)先采用）

代表了采用權(quán)重分配垂向各層的濃度情況。那么后續(xù)在轉(zhuǎn)換參數(shù)下方多了一行權(quán)重

值。權(quán)重值，從左到右為垂向的底層到表層，如001則代表分三層，在表層釋放。

?0

采用單獨(dú)的序列組來(lái)分配各層的濃度分布情況

?MCSER

序列數(shù)目

?TCCSER

時(shí)間轉(zhuǎn)換單位為秒，如果后續(xù)采用時(shí)間為天，則為86400.

?TACSER

時(shí)間轉(zhuǎn)換之前，統(tǒng)一增加某值

?RMULADJ

濃度單位換算系數(shù)

?ADDADJ

轉(zhuǎn)換單位前，全部濃度值增加某值。

?數(shù)據(jù)

?時(shí)間

?指標(biāo)值

?3空間分布條1牛（#初始條1牛）

?濃度空間分布的輸出條件

通過(guò)網(wǎng)格，層，值控制

?4WASP耦合文件

?5貝類(lèi)濃度模擬

?63D顯示輸出指定文件

?7土壤濕度文件

?8開(kāi)邊界的周期性測(cè)試

?9遮罩設(shè)定

設(shè)定不過(guò)水的網(wǎng)格面，阻擋水流流過(guò)網(wǎng)格某一個(gè)面

?阻擋方向

?1西方向

?2南側(cè)方向

?3西和南方向

?10渠道網(wǎng)格模型

子網(wǎng)格的渠道模型，允許一維狹窄網(wǎng)格通過(guò)二維大的平面網(wǎng)格。

?11子網(wǎng)格變化參數(shù)

快速的改變子網(wǎng)格的尺寸

?12開(kāi)邊界水彳濾制pser

開(kāi)邊界的水位控制文件

?參數(shù)

同濃度等輸入條件一致

?數(shù)據(jù)

?13水工構(gòu)筑物

水工構(gòu)筑物控制，流量水位曲線(xiàn)，這里就是水力學(xué)中的水工構(gòu)筑物控制。

?參數(shù)控制

?開(kāi)啟權(quán)重

?序列數(shù)目

?上下游水位差

?上下游乘積系數(shù)

?14流量序列

流量序列文件，同前面濃度序列設(shè)置

?15重啟文件

設(shè)定重新啟動(dòng)模型的初始條件

?16對(duì)流擴(kuò)散文件

EFDC在運(yùn)移模式中指定運(yùn)移擴(kuò)散文件

?17屏幕顯示控制

在模型運(yùn)行過(guò)程中，顯示某位置的模擬值指標(biāo)

?18貝類(lèi)行為信息

?19植被阻力

?20波浪模擬文件

?6水質(zhì)和富營(yíng)養(yǎng)化模塊

?運(yùn)行

模型的運(yùn)行

?7編譯和執(zhí)行代碼

這部分涉及到程序編程的內(nèi)容。為了保證模型的兼容性，模型設(shè)置可以在許多平臺(tái)運(yùn)行，不同平臺(tái)

的編譯方式有所不同。一般情況我們都使用的是windows系統(tǒng),所以采用windows系統(tǒng)下的編

譯方式。

?文件

?EFDC.ROF

?EFDC.COM

全局common塊

?EFDC.PAR

全局聲明變量

?自定義參數(shù)

為保證模型的最優(yōu)化，需要將模型設(shè)置為較為符合建模區(qū)的參數(shù)，其設(shè)定為越大肯定保證模型

的最大兼容，但是內(nèi)存就會(huì)浪費(fèi)，小值會(huì)導(dǎo)致模型在其他網(wǎng)格下不能運(yùn)行，如果嫌麻煩就設(shè)定

足夠大，現(xiàn)在機(jī)器的配置應(yīng)該足夠強(qiáng)。

?ICM=MAXIMUMXORICELLINDEXTOSPECIFICGRIDINFILEcell.inp

CELL文件中最大X或者I標(biāo)號(hào)

?IGM=MAXIMUMICELLINDEXFORGRAPHICSGRIDSPECIFIEDBYGCELLMAP.INP

gcell文件中最大I單元標(biāo)號(hào)

?JCM=MAXIMUMYORJCELLINDEXTOSPECIFICGRIDINFILEcell.inp

CELL中最大的V坐標(biāo)或者J格標(biāo)號(hào)

?JGM=MAXIMUMJCELLINDEXFORGRAPHICSGRIDSPECIFIEDBY

GCELLMAP.INP

Gcell文件中的最大可視化J坐標(biāo)

?KCM=MAXIMUMNUMBEROFLAYERS,MAXLOOPINDEXKC

垂向最大層數(shù)

?KGM=KCM

?KSM=KCM-1

?KPCM

MAXIMUMNUMBEROFCONSTANTELEVATIONLEVELFORTHREEDIMENSION

GRAPHICOUTPUT三維圖形輸出中高程值的最大值

?LCM=MAXIMUMNUMBEROFWATERCELLS+2OR1+THEMAXLOOPINDEX

LA

網(wǎng)格數(shù)目+2或LA值+1

?LCMW=SETTOLCMIFISWAVE.GE.lOTHERWISE=2

波浪文件開(kāi)啟時(shí)為1,否則為2

?LCGLM=SETTOLCMIFISLRD.GE.lOTHERWISE=2

?MGM=2*MTM

?MLM=MAXIMUNNUMBEROFH