EFDC操作手册资料

EFDC操作手册

•综述

通过了解EFDC-EPA的全局及对EFDC模型的全过程有一个大概的了解。

•1介绍

•起源

VirhiniaInstituteofMarineScience(Hamrick)

•特性

•矩形或正交曲线网格

水平方向采用

•Sigma垂向坐标

垂向

•欧拉运输方程

•拉格朗日粒子

可以模拟拉格朗日的粒子

•干湿网格

模拟浅滩的情况

•闸坝控制

•水生植被

•与其他模型耦合

•数值计算

•有限差分方法

•二阶精度的有限差分空间度散

•二阶精度、三层时间步长的时间差分

•前处理

模型有前处理的网格程序，能够生成网格及深度的插值计算。

•网格工具

•深度插值工具

•初始条件

•输入格式采用文本格式

•运行

•实时的演示

•诊断文件

•重启文件

用于热启动，就是不用从初始计算重新开始计算，相当于半中间的直接启动

•后处理

•时间序列

•空间分布图

•矢量图

•等值线

•剖面图

•2全局概览

•运行程序

翁要介绍了efdc代码的执行程序的组成部分，这部分内容后面会详细叙述.

•efdc.for

源代码文件，主要代码文件

•

通用声明变量文件，一般模型使用

•efdc.par

特别的用于最小化内存占用管理

•模型配置文件（输入文件）

模型通过后缀名为inp的文本文件来进行输入文件的配置,按照功能分为5中类型。

•空间网格地形组

用于构成水休模型的空间网格及地形条件。里面有一些概念可能会难以理解，简单的衡量

空间为什么会有这么多文件，比如网格方向的概念，对于图形界面来说，很难理解，这些

内容我自己觉得是使用字符来表达空间所必须的。我们看到的空间实体，采用文本文件来

抽象表示，本身就很复杂。为什么一图胜前言，信息传递给计算机能够读懂的语言本来就

比较麻烦。后续章节到网格时会详细描述。

•cell

将水体外部边界线数字化，在字符时代，用于不同的数字来代表水体和陆地界面，从而

刻画出水体的形状。

•celllt

cell文件的补充文件，通常为cell文件中不包括入流和出流口的文件。

•dxdy

网格的几何特性

•空间网格的分辨率

水平方向上每个网格的大小间距

•深度

这里的深度是初始水深。水深：即水面到底部地形的距离，是水位的表征量。在计

算过程中随着水位变动。而初始水深代表了初始水位。

•底部高程

垂直方向的上的网格的高程（一般是网格中心点计），这里代表实际的地形。和上

述概念不同的是,这里代表实际的地形，一般而言，其在计算过程中不会改变。

•底部粗糙率

•水生植被

考虑水底的植被作用，影响水动力作用。需要在主控文件中开启植被阻隔。

•Ixly

网格的地理特性

•网格的中心坐标

网格中心点的坐标,EFDC一般采用的是UTM坐标系。

•网格的方向

网格的方向，后续计算流向时的转换。

•depth

仅仅用于笛卡尔网格来指定水深、粗糙度等。为笛卡尔网格单独开了一道门。

•mask

阻挡物，相当于固定边界，只不过其可以设定在水体中心，比如湖中的岛可以采用这个

来描述。但需要说明的是为了计算方便，岛屿通常采用挖空网格来表示。

•主控及运行控制组

模型的总指挥文件，控制一些模块的开关及参数的设定。

•efdc

该文件是所有控制文件的核心文件，管控了几乎所有的参数及模块的开关。

•模块的开关

如是否开启水动力模块，是否模拟染料。

•边界条件位置及其序列数量

对于河道来说，支流入干流位置点，支流数量等。

•控制输出格式

是否生成流速等指标

•show

在运行过程中显示的指标及信息（提供计算时的实时演示）

・时间序列边界条件组

•aser

气象条件的时间序列文件，包括气压、温度、湿度等

•qser

流量时间序列文件

•dser

染料序列文件

•pser

开边界水位序列文件

•sser

盐度时间序列

•物理过程组

•qctl

水工构筑物控制文件

•vege

植被阻碍流动特征文件

•wave

表面重力波文件

•初始条件及启动组

模型启动时的初始条件的设定

•dye

初始染料空间分布

•restart

热启动指定文件

•salt

初始盐度空间分布

•模型结果文件（输出文件）

一般以后缀名为.out的文件输巴

•诊断输出文件

•重启文件

•网格点的时间序列文件

指标随时间变化的值

•最小平方谐波分析文件

•二维及三维图形可视化文件

•前处理

前处理,及构建模型及参数的设置

•3网格

模型构建的第一步是生成研究区的网格文件。

•基本概念

•四边形网格（三角网格）

EFDC采用四边形及三角形网格（为四边形网格的对角线分割后得到）

•坐标系xy->ij

网格的编号为ij值，对应直角坐标系中的ij°x的正方向为东，y方向为北方向。

•EFDC采用相对坐标系

地理中的实际的正东方向用x"正北方向用y\

•国家坐标系

•UTM

需要注意的是，EFDC采用的是国际标准单位，MKS,所以推荐使用UTM坐标系。

•用户自定义

•网格文件

•CELL文件

使用数字阵列表来代表网格分布情况。

•头文件

4行头文件，为固定格式。

•坐标轴区域

水体单元格上下部分为坐标轴的标识，i为x方向，j为y方向。

•水髀元格

中间整块数字部分就是区域网格了，采用不同的数字代表不同的网格类型，这里为了明

确表示水体边界，是加入了陆地类型的。边界的概念为，水体和陆地的交界处。举个例

子，windows的扫雷游戏玩过吧，那么多格子就代表网格，有雷的区域就用1表示，

那么一整块雷就全部是L

•CELLT文件

可指定为CELL文件的子集,用于特殊的目的，如表示除入流口以外的水体。则原来的边界

可能变为陆地，原来的水体变为边界或陆地。

•DXDY

指定网格的几何特定，单位为m,具体前面说明前面已经有了不再重复。

•空间网格的分辨率

水平方向上每个网格的大小间距

•深度

这里的深度是初始水深，水深：即水面到底部高程的距离，是水位的表征量。在计算过

程中随着水位变动。而初始水深代表了初始水位。水位;深度+底部高程

•底部高程

垂直方向的上的网格的高程（一般是网格中心点计），这里代表实际的地形。和上述概

念不同的是，这里代表实际的地形，其在计算过程中不会改变，是地理量。

•底部粗糙率

•水生植被

考虑水底的植被作用，影响水动力作用。需要在主控文件中开启植被阻隔。

•LXDY

•网格的中心坐标（位置）

网格中心点的坐标，EFDC一般采用的是UTM坐标系，用于后处理时的输出结果的可

视化，依照方便可设定为m,KM等。

•旋转矩阵（方向）

网格的方向，用于将网格的伪方向转换为实际地理中的方向。

•速度方向

结果文件中的流向转换为真正的流向。（地理坐标）

•风场方向

实际中的风场利用逆矩阵转换为模型中的风向。

•网格生成gefdc.f

实际操作中网格的生成是采用gefdc.f来生成的。

•准备文件

cell

前面已经叙述，数字轮彝

•depdat

指定深度及底部高程。xyz格式，xy为坐标，z为值。

•gcell

曲线网格时进行辅助四边形网格文件同cell类似格式

•gridext

水体网格角坐标文件（左下角），三角形网格也必须等同四边形设定。在设定参数

NTYPE=O时采用。

•vege

同depdat设定植被

•zrouth

同上

•gefdc设定文件

网格生成工具的主控文件，总共13个卡片组，起始efdc.inp文件也是efdc的主控文件，

在这里有异曲同工之妙，所以简单学习下为后面铺垫。文件的格式是文本格式，每组以C*

开头，第一行为标题，第二行为参数标题，第三行为参数值。

•C2咔2

生成方式，必选

•NTYPE

生成网格的方式

•0

采用gridext和cell生成网格

・1-5

正交曲线网格，一般采用5

•6

不保证正交性，需要指定边界点

•7

需要用户指定函数，矩形区域网格

•8

采用cell和卡4的设定深度插值法生成网格

•9

深度插值法，但坐标系和8不同

•NBPP(NTYPE=l-6)

边界点数目

•IJ

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

・C3

在正交曲线网格时，保留一个规整的矩形网格备用，用于后期演示

・C4

网格坐标信息

•C5

网格迭代选项

•C6

数值方面的参数

・C7

单位变换

・C8

插值方式

•C9-10

设定矩形网格的参数

•C11

d叩data的插值方式

•C12-13

边界的设定

•网格生成VOGG

visualorthogonalgridgenerationsystem,正交曲线网格生成工具

•正交曲线网格

正交曲线网格能够贴合边界，弥补了直角网格在复杂边界处的锯齿状网格，从而减少了网

格数量，

•VOGG使用

•生成思路

•概化模型

首先确定研究区域的边界，也就是水体和陆地的边界区域。这部分需要将其概化，

依据研究目标和实际区域形状进行概化。

•边界线应该相对平滑（整体接近矩形块）

•边界区域的角尽量为直角

•生成计算区域-矩形网格

•节点（边角点）-角度

•线（点和点之间的线）-长度

•网格数量梳理

手动或者自动定义网格数量，生成U值

•网格生成-由矩形转换为正交曲线网格

•源代码介绍

•cell.inp

同efdc所需文件类似，但其仅仅只有0,9,5,后期生成后再进行修改为三角网

格

•gridside

网格点或边数量控制，坐标位置输入

•Nside

节点数量或者边数量

•Ijndex

左下角i的标号

•jjndex

左下角J的标号

•inter

网格插值分辨率

•vogg

控制网格生成参数及地形插值

•可视化操作

•界面简介

•左侧工具栏

工具排列的顺序，基本就是网格制作的流程。

•线控制

用于草图绘制，绘制水体边界的草图

•点控制

精确修改草图的位置点，通过节点调整草图吻合实体

・边界编辑

修改节点的方向，实质是调整网格线的走势，全部必须是单一方向，一般

按顺时针。

•U排序工具

这里起始就是网格生成前一步，将节点转换为网格节点

•网格工具

调整网格，删除

•上侧工具栏

•视图控制

移动操作区域的位置方向

•生成网格

也用于导入生成好的网格

•菜单栏

全局控制及导入文件的菜单

•导入水体边界（shp）

•导入草图边界（shp）

•导入网格文件（shp）

•导入坐标文件

•网格生成步骤

•底图准备

依照程序需要将底图准备好

•导入边界

水体边界的底图导入，一般我们要将水体作为底图放入网格生成器中

•草图绘制

简略的框架绘制，绘制为一个矩形框架，控制网格的边界，不同的网格生成程

序中的操作不同，但均要保证网格的边界平滑，能够在宏观上贴合水体边界

•修改节点属性信息

•生成网格

网格在边界基础上生成，通过控制参数来控制网格的分辨率

•调整网格

网格不可能一次就能够完美，需要进行不断微调，因为网格后期基本不会改动，

所以一般在网格生成阶段就需要花费较多精力来调整到足够好。

•地形插值

一般网格的工具都会附带地形工具，也就是给网格垂直方向深度赋值。

•4主控文件

efdc.inp文件

•概述

参数一般是通过整型变量设定,0即为关闭状态，非0为打开状态。

•运行控制参数

控制模型运行时的开关，比如是否开启水动力模块。

•输出控制参数

输出结果的内容控制，控制模型输出什么样的结果文件。

•物理过程参数

对于水动力、水质模型本身的一些参数如扩散系数等。

•卡组说明

对于卡组中的很多值并非都需要很熟悉，所以仅仅先说明一些我比较熟悉的值，后续对其他进

行补充。

•卡1

模型运行的标题文件

•卡1A

设定有关空间和时间计算的参数

•卡2

重启文件，通用设置，诊断切换

•重新启动文件参数hiiD：///China3S/D/3192807.hlml

重新启动文件的使用是用于热启动，也就是说在运行后，下次运行时则从某一时刻开始，

而不重新从最初的状态开始,这里是为了跳过冷启动的时间。在EFDC的源代码中，

有关热启动的相关代码包括：RESTOUT.for.RESTINlO.for、RESTIN1和

RESTMOD.for,其中RESTOUT.for程序处理写热启动文件RESTART.OUT和

RSTWD.OUT文件，而RESTINl.for分别处理热启动文件RESTART.inp和

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

且一般的水动力模拟不调用RESTINlO.for,只调用RESTINl.for。RESTMOD.for只

有HDMT和HDMT2T.for两个子程序才会调用，一般只模拟水动力不会调用该子程序。

RESTIN2.for也会涉及RESTART.inp分类：EFDC模型EFDC热启动是在运行EFDC经

常用到的方法，例如，在运行EFDC出现意外中断时，热启动可以帮助用户继续从中断

处运行EFDC,而不需要从头开始运行。但热启动经常会碰到负值问题而使EFDC无法

运行下去，这种问题的解决办法是：用Uedit32对RestartJnp文件进行修改，将-0.

替换为0.即可。

•ISRESTI

重启文件的读取

•1从inp文件中读取初始条件

•-1和上述一致，但根据底部高程进行调整

•2对于多层的处理，当新KC大于4时，可以从KC/2层中初始化变为KC层

•10旧版本中读取？

•ISRESTO

重启文件的生成

•-1运行结束后生成重启文件

・N每N个基准参考时间段生成重启文件

•多线程运算

•输出到屏幕的内容

•写出一些平衡值到.out输出文件中

•卡3

一些数值计算的方法及迭代次数的设定

•卡4

对对流扩散方程进行自定义

•卡5

•水平动量扩散设定

•WASP输出设定

•干湿模式设定

•关闭干湿

•常量湿深度非线性迭代

•变量湿深度非线性迭代

•常量深度

•常量迭代存在mask

•使用网格面的变量干湿模式

•使用mask的常量干湿模式

•植被阻流作用

•涡旋强度公式

•波浪动力学

・C6

物质的垂向对流和扩散以及物质的源汇数量域。

•每组变量用一个数值代替，且有各自的参数

•ISTRAN是扩散和源汇项的开关

•ISTOPT

选择指标的初始条件的设定方式，以及其的参与方式

•ISADAC

数值计算方式

•C7#时间参数

时间设定参数整数变量，通过和C8配合使用来设定时间参数。

•NTC

efdc采用设定迭代倍数来设定运行时常，也就是说其设定一个运行的时间单位，这里

设定次数。

•NLTC

时间步长，同上面，是通过一个参照值来设定时间步长的。

•变化的步长和固定步长

•输出的频率

•C8#时间参数

实际参数

•TCON

开始时间的转换系数，当相对标准时间不为秒时，要进行转换。

•TBEGIN

模型启动时刻，通过TCON来转换为秒。

•TREF

模型的相对标准时间,相当于一个基准时间，卡7中的变量都是参照此进行计算的。

•科氏力设定

•变步长的设定

•C9#空间参数

水平空间参数设定

•水平网格参数

设定水平网格的数量及结构

•IC\JC

设定I\J的方向的网格数量

•LC\LVC

由于计算方便，EFDC采用了一个数来代替（i,j）数组，将其转换为L,但为了循

环方便，其增加了L=1和L=LC,完全是程序需要。

•ISCLO

正交曲线网格

•并行计算的网格区域#并行计算

有关并行计算的设定

•ISMASK

设定不透水网格面的开关，不透水网格就是可以指定网格某一面是不透水的。

•ISPGNS#开边界设定

设定北部和南部开边界中的周期性信息

•初始条件的空间插值

光滑参数，初始条件的插值是efdc执行的，可以设定其的插值参数。

•C9a

垂向空间参数设定

•垂向层数设定KC

•CIO

设定具体的sigma层数比例分配

•K层

K是从下向上计算,多层时,K=1为底层

•DZC层比例，总和必须为1

•C11

网格、糙率、深度值参数

•矩形网格的固定参数

•糙率的调整转换

•底部深度的调整转换

•干湿深度设定

•C11A

动量修正系数设定

•C11B

角网格修正参数

•C12

涡流及湍流的相关扩散系数

•AHO

水平动量和质量扩散系数，常量值

•AHD

无量纲水平动量扩散系数，变量的附加在常量AHO上,Smagorinsky亚尺度水平扩

散公式。

•AVO

背景，常量或分子动力粘度

•ABO

背景，常量或分子扩散系数

•AVMN

最小涡流粘度

•ABMN

最小涡流扩散系数

•VISMUD

恒定浮泥粘度

•AVBCON

设定为0时，涡流未被激活，采用ABO和AVO设定为垂向扩散系数和粘度系数,设

定为1时采用模型计算值

•ZBRWALL

边壁对数粗糙度，在水平动量扩散系数激活后，并且AHD和AHO不为0.

•C12A

•C13

两个应该是高级应用时采用的，没有详细说明。略

•C14

潮汐、大气强迫条件、地下水和渠道子网格的参数设定

•潮汐强迫条件数目

•风序列数目

序列代表了时间和值构成了一个时间序列，如流量序列。

•大气序列数目

•地下水设定

•土壤湿度平衡

•地下水和底床交互

•激活河道子网格模式

•波浪文件选择

•潮汐高程同化#疑问

•对于超出干网格的水开启渗透或消除作用

•C15

潮汐调和设定符号设定

•C16

水位边界条件和压力边界条件参数设置

•水位边界条件或开边界数目（东南西北）

•潮汐调和边界数目

•强迫类型

•边界时间序列数目

•C17

周期水位边界及压力边界具体设定

•编号

•振幅

•相位

•C18

南部开边界设定

•SB

•方式

•第一类型开边界

•指定高程

•调和常数

•第二类型

•RADIATION-SEPARATION辐射分离条件0切向速度

•RADIAHON-SEPARATION条件自由切向速度

•时间序列方式

•应用潮汐调和常数

•C19

西部开边界设定

•C20

东部开边界设定

•C21

北部开边界设定

•C21A

水位和速度数据

•C22

模型读取的时间序列的数量和污染物的种类设定

•C23

速度、体积源汇项、流量控制、取水和退水边界的设定

•体积源汇项的位置数目

常量或者变化的体积源汇项（如流量），对于正向入流是考虑污染物随流的输入，而出

流则考虑了物质的损失，

•体积源汇项的时间序列数目

•压力控制措施（闸坝）的位置数目

控制水位的水工构筑物

•压力控制措施（闸坝）的时间序列数目

•取水'回流边界的位置数目

常量或者随时间变化的取水和回水边界，为电站和冷却系统设计。

•取水'回流边界的时间序列数目

•C24

体积源汇项的位置指定，浓度序列链接。

•位置设定

通过设定网格ij来确定边界条件的位置

•为位置指定序列号

为该位置指定相应的流量序号、污染物浓度序号，对应dser中的组号（从上到下）

•比例因子QSFACTOR

设定比例

•恒定入流和出流设定

实际当中并不直接使用，而采用时间序列的恒定格式。若设定，则是此值加上时间序列

值。

•侧向入流，体积量由m2/s转换到m3/s

•动量通量的朝向指定

•C25-C26

恒定序列的恒量设定。

•C27-C31

射流和羽状流的设置。

•C32

水位和压力流的参数设定信息（包括闸坝水位设置）

•位置（上游和下游设定）

支持单向流和双向流，同时水位控制类则只需要设定上游，如溢流坝等，但是如果是有

控制的闸，则需要上下游同时设定。

•水流控制类型

•水位-流量曲线设定

常用，一般用于河流的下游边界，通过W-Q的关系来进行下游边界条件的控制。

如果有水文部门的公式则使用构造水位和流量值，如果没有，这可以用曼宁公式推

的大概的关系。

•上游水深压差控制型

•潮汐入口加速流

•动量通量设定

•动量通量宽度

•C33

流量取水、回水边界

•位置（上游和下游设定）

支持单向流和双向流，司时水位控制类则只需要设定上游，如溢流坝等

•取水回水序列编号

•恒定取水回水量值

•动量通量设定

•动量通量宽度

•C34-C35

取水回水水质指标及泥沙恒定浓度设定

•C36-C42

沉积物模块设置

•C43

有毒污染物的初始条件和参数设定

•初始条件

•空间均匀浓度分布

•水柱中的空间变化的浓度分布

•在沉积物中空间变化的初始条件

•在水柱中和沉积物中空间变化浓度分布

•在沉积物中的单位表征

•水中有毒污染物的降解率及参考温度

•底泥中的有毒污染物降解率及参考温度

•C44

其他有毒污染物的参数

•溶解及颗粒态有机碳吸附作用

•水面蒸发率

•分子量对蒸发率的影响

•光分解速率

•太阳能辐射影响

•孔隙水中的扩散系数

•沉积物上层与水体间的扩散系数（正负值不一样）

•扩散系数

•消散系数

•沉积物中的颗粒混合扩散系数

•沉积物中的颗粒混合深度范围

•C45

有毒污物的蹒设定

•C45A

有毒有机化合物的时间和空间量设定

•时间恒定

•空间均匀分布和空间变化的分布

•C45B

有毒有机化合物的设定

•C45C

有毒物颗粒态在水柱中的分布

•C45D

有毒物颗粒态在沉积层中的分布

•C46

开边界中的浓度序列数目设定。

•C47

南侧开边界入流浓度指定时间序列序号

•C48

南侧浓度边界的底层恒定浓度

•染料

•TOX等

•C49

南侧浓度边界的底层恒定浓度

•黏性泥沙

•非黏性泥沙

•C50

南侧浓度边界的表层恒定浓度

•C51

南侧浓度边界的表层恒定浓度

•C52-C66

其他方向的开边界浓度设定

•C66A-C66B

数据同化，我的理解是将临近的边界条件进行类似设定。

•C67

拉格朗日粒子追踪

•计算方法选择

•点位置释放

•释放粒子数目

•粒子释放间隔

•输出时间步长

•范围内释放，MLRPDRT成组，并可计算全部释放时间的平均值

•边界区域设定

•可选高阶轨迹积分方法

•输出速度矢量控制

•C68

释放点

•C69

笛卡尔网格中心经纬度

•C70

二进制文件及ASCII文件设定

•C71

水平等值线生成

•C71A

沉积层的等值线生成

•C71B

食物链模型控制

•C72

水平面的水位和压力等值线

•C73

速度矢量的输出

•C74

垂向剖面等值线

•C75

剖面单元格的数目

•C76

垂向剖面线的位置

•C77-79

垂向速度矢量

•C80-81

3d文件输出

•C82

调和分析输出

•C83

调和分析的位置切换

•C84

时间序列文件输出

•输出文件的覆盖和追加

•输出点位置数目

•开始输出的时间点

•结束输出的时间点

•输出间隔

•C85-C86

输出模式模式化设定

•C87

输出量位置及指标确定

•检出水位结果时间序列

•输出浓度变量的时间序列

•涡流粘度和扩散速率的时间序列

•外部模式水平速率的时间序列

•输出外部模式水平传输率的时间序列

•每一层的水平速度的时间序列

•输出外部模式净体积源汇的时间序列

•输出净体积源/汇的时间序列

•备注

C88-90

垂向输出设定，垂向采样设定

•C91

沉降速度的范围

•C92

所选网格的最小悬浮物浓度

•C93

空间变化的侵蚀函数

•C94

空间变化的食物链暴露深度

•5水动力模块

水动力文件的附加文件，如流量、浓度时间序列，水位流量曲线等。

•1aser

模型必须，大气及气象条件输入。

•使用参数：卡6的参数控制大气驱动的模式

•转换系数

转换系数是为了将原始数据统一进行转换为模型内部使用，一般而言，为了保持数据统一

性，不推荐修改这些值，而建议的方法是输入之前就转换好，这样防止出错。就是在输入

之前就处理好。

•序列数目

•时间换算为秒，86400代表后续使用时间为天

•数值统一加上一个值

•风速转换

•降雨、蒸发转换

•数据

•风速

•风向

•气压

•相对湿度

•降雨

•蒸发

•短波辐射

•2浓度序列文件（#边界条件）

输入点源式的浓度序列文件

•使用参数

•ISTYP

序列类型

•1（优先采用）

代表了采用权重分配垂向各层的浓度情况。那么后续在转换参数下方多了一行权重

值。权重值，从左到右为垂向的底层到表层，如001则代表分三层，在表层释放。

•0

采用单独的序列组来分配各层的浓度分布情况

•MCSER

序列数目

•TCCSER

时间转换单位为秒，如果后续采用时间为天，则为86400.

•TACSER

时间转换之前，统一增加某值

•RMULADJ

浓度单位换算系数

•ADDADJ

转换单位前，全部浓度值增加某值。

•数据

•时间

•指标值

•3空间分布条1牛（#初始条1牛）

•浓度空间分布的输出条件

通过网格，层，值控制

•4WASP耦合文件

•5贝类浓度模拟

・63D显示输出指定文件

•7土壤湿度文件

•8开边界的周期性测试

•9遮罩设定

设定不过水的网格面，阻挡水流流过网格某一个面

•阻挡方向

•1西方向

•2南侧方向

•3西和南方向

•10渠道网格模型

子网格的渠道模型，允许一维狭窄网格通过二维大的平面网格。

•11子网格变化参数

快速的改变子网格的尺寸

•12开边界水彳滤制pser

开边界的水位控制文件

•参数

同浓度等输入条件一致

•数据

•13水工构筑物

水工构筑物控制，流量水位曲线，这里就是水力学中的水工构筑物控制。

•参数控制

•开启权重

•序列数目

•上下游水位差

•上下游乘积系数

•14流量序列

流量序列文件，同前面浓度序列设置

・15重启文件

设定重新启动模型的初始条件

•16对流扩散文件

EFDC在运移模式中指定运移扩散文件

•17屏幕显示控制

在模型运行过程中，显示某位置的模拟值指标

•18贝类行为信息

•19植被阻力

・20波浪模拟文件

•6水质和富营养化模块

•运行

模型的运行

•7编译和执行代码

这部分涉及到程序编程的内容。为了保证模型的兼容性，模型设置可以在许多平台运行，不同平台

的编译方式有所不同。一般情况我们都使用的是windows系统,所以采用windows系统下的编

译方式。

•文件

•EFDC.ROF

•EFDC.COM

全局common块

•EFDC.PAR

全局声明变量

•自定义参数

为保证模型的最优化，需要将模型设置为较为符合建模区的参数，其设定为越大肯定保证模型

的最大兼容，但是内存就会浪费，小值会导致模型在其他网格下不能运行，如果嫌麻烦就设定

足够大，现在机器的配置应该足够强。

•ICM=MAXIMUMXORICELLINDEXTOSPECIFICGRIDINFILEcell.inp

CELL文件中最大X或者I标号

•IGM=MAXIMUMICELLINDEXFORGRAPHICSGRIDSPECIFIEDBYGCELLMAP.INP

gcell文件中最大I单元标号

•JCM=MAXIMUMYORJCELLINDEXTOSPECIFICGRIDINFILEcell.inp

CELL中最大的V坐标或者J格标号

•JGM=MAXIMUMJCELLINDEXFORGRAPHICSGRIDSPECIFIEDBY

GCELLMAP.INP

Gcell文件中的最大可视化J坐标

•KCM=MAXIMUMNUMBEROFLAYERS,MAXLOOPINDEXKC

垂向最大层数

•KGM=KCM

•KSM=KCM-1

•KPCM

MAXIMUMNUMBEROFCONSTANTELEVATIONLEVELFORTHREEDIMENSION

GRAPHICOUTPUT三维图形输出中高程值的最大值

•LCM=MAXIMUMNUMBEROFWATERCELLS+2OR1+THEMAXLOOPINDEX

LA

网格数目+2或LA值+1

•LCMW=SETTOLCMIFISWAVE.GE.lOTHERWISE=2

波浪文件开启时为1,否则为2

•LCGLM=SETTOLCMIFISLRD.GE.lOTHERWISE=2

•MGM=2*MTM

•MLM=MAXIMUNNUMBEROFH