fort. 15 常用参数定义

RUNDES =字母数字运行说明1（<= 32个字符）

 

RUNID =字母数字运行说明2（<= 24个字符）

 

NFOVER =非致命错误覆盖选项;
= 0，输入参数不一致将导致程序终止。
= 1，不一致的输入参数（如果可能）将被自动更正为默认值或一致的值，并继续执行。确保阅读非致命警告消息，以查看是否已修改任何参数。请注意，并非所有不一致的参数都可以自动更正，因此仍然可能导致致命错误消息和程序终止。
NFOVER注意事项：
有时，由于数值参数设置不正确，违反时间步长稳定性标准，高程解会变得物理上过大。 ADCIRC基于对水位高度的用户特定限制终止可能很有用。要启用此功能，必须使用DEBUG_WARN_ELEV编译器指令来编译ADCIRC。这在cmplrflags.mk文件中设置。例如，对于串行模型，
DA：= -DREAL8 -DLINUX -DCSCA –DDEBUG_WARN_ELEV
对于并行模型：
DP：= -DREAL8 -DLINUX -DCSCA -DCMPI -DDEBUG_WARN_ELEV
这将在NFOVER行中指定的Fort.15文件中启用其他参数：
NFOVER，WarnElev，iWarnElevDump，WarnElevDumpLimit，ErrorElev
ADCIRC然后监视最大水位，其行为如下：
如果高程超过WarnElev，则会发出警告。
如果超出了WarnElev并且iWarnElevDump == 1，则将写入全局高程文件（写入fort.69）。
如果由于上述警告限制而写入了WarnElevDumpLimit全局高程文件，则执行将终止。
如果海拔高度超过ErrorElev，则执行将终止。
默认值为：
WarnElev = 20.0！在20米处警告
iWarnElevDump = 0！不要写全局高程文件
WarnElevDumpLimit = 50！如果达到警告级别50次，则终止执行
ErrorElev = 1000.0！终止执行海拔超过1000米的水
示例：要覆盖默认值，请按上述方法重新编译ADCIRC，然后修改fort.15文件中的NFOVER行：
1 10. 0 100 30。
当水位高程超过10米时，这将导致ADCIRC发出警告，不会写入全局高程文件，并且如果生成100条警告或水位超过30米，则ADCIRC将终止。

 

NABOUT =从ADCIRC输出到屏幕或控制台的日志记录级别以及ADCIRC日志文件（fort.16）。 ADCIRC以5个严重级别写入日志消息：DEBUG，ECHO，INFO，WARNING和ERROR（从最低或最低重要到最高或最重要）。选择日志记录级别表示应该记录该级别及更高级别的消息。将日志记录级别设置为WARNING或ERROR将减少日志文件的大小（和混乱），但是可能会丢失重要的信息。错误消息通常是由导致运行停止的问题引起的。
NABOUT = -1 :: DEBUG级别的日志消息和更高版本。这可能会占用大量磁盘空间，并且可能会使ADCIRC变慢。通常仅对ADCIRC开发人员有用。
NABOUT = 0 :: ECHO级日志消息和更高版本。 ECHO级日志消息包括大多数输入文件的回显打印，包括Fort.13，Fort.14和Fort.22文件。
NABOUT = 1 ::: INFO级日志消息和更高版本。这些消息通知用户ADCIRC已完成的重要操作，但不是问题或结果的结果。
NABOUT = 2 ::警告级别的日志消息和更高版本。这些消息表示通常不会对运行造成致命影响的潜在问题。
NABOUT = 3 ::仅错误级别的日志消息。这些消息表明严重的问题，通常会导致运行停止。
在ADCIRC版本49之前，0（ECHO）和1（INFO）是唯一可用的选项。 NABOUT最初代表“缩写输出”（用于记录文件）。

 

NSCREEN =控制日志消息输出到屏幕（即标准输出）。 时间步记录将在每个Abs（NSCREEN）时间步中写入。 屏幕输出主要包括时间步记录。
NSCREEN <0 ::通常将写在屏幕上的日志消息写到名为adcirc.log的文件中。
NSCREEN = 0 ：：日志消息将不会写入屏幕。
NSCREEN> 0 ::将日志消息写入屏幕（也称为标准输出）。

 

IHOT =控制模型是否热启动的参数。 热启动功能可用于2D和3D运行。 如果正在进行谐波分析，则热启动文件还将包含谐波分析数据，因此也可以热启动谐波分析。
= 0冷启动模型
= 17从ascii文件fort.17热启动
= 67使用热启动文件fort.67中的输入信息进行热启动
= 68热启动模型，使用热启动文件fort.68中的输入信息
= 367热启动模型，使用netCDF热启动文件fort.67.nc中的输入信息
= 368个热启动模型，使用netCDF热启动文件fort.68.nc中的输入信息
= 567个热启动模型，使用netCDF4热启动文件fort.67.nc中的输入信息
= 568热启动模型，使用netCDF4热启动文件fort.68.nc中的输入信息

 

ICS =参数，用于控制模型是在球面坐标系还是在笛卡尔坐标系中运行。
= 1 ADCIRC控制方程式在标准笛卡尔坐标中。假定网格文件（fort.14）中的坐标的长度单位与重力单位（如下所述）一致。对于这种类型的运行，在不太可能的情况下需要潮汐强迫（NTIP = 1或2）和/或空间可变的科里奥利系数（NCOR = 1）的情况下，可以使用反向地图投影（Carte Parallelo-grammatique）获得网格的经度和纬度值。但是，我们强烈建议，如果模型域足够大，以至于空间变量科里奥利或潮汐强迫都被认为很重要，则应使用经度，纬度网格使用球形控制方程（ICS = 2）运行模型。
= 2个ADCIRC控制方程在球坐标中转换为笛卡尔坐标，然后使用地图投影离散化（Carte Parallelo-grammatique – CPP）。网格文件（fort.14）中的坐标以经度和纬度为十进制度。

 

IM =型号类型
= 0使用新的GWCE和动量方程公式运行正压2DDI
使用新的GWCE和基于速度的动量方程式进行1次正压3D运行
使用新的GWCE和基于速度的动量方程式进行21次斜压3D运行
= 111112正压2DDI使用集总GWCE（而不是默认的完全一致GWCE）运行。 需要使用此选项以集总+显式模式运行ADCIRC，从而绕过迭代求解器。 显式模式还需要在此文件中指定系数A00，B00，C00（= 0、1、0）。 （见下文）
= 611112使用集总GWCE（而不是默认的完全一致GWCE）运行正压3D。 需要使用此选项以集总+显式模式运行ADCIRC，从而绕过迭代求解器。 显式模式还需要在此文件中指定系数A00，B00，C00（= 0、1、0）。 （见下文）

 

IDEN = 3D运行中的密度强迫形式（对于所有斜压模型运行，从UNIT 11读取初始密度，温度和/或盐度字段）; ADCIRC不支持2D斜压模型运行。
在盐度和温度强迫下运行的= -4诊断斜压ADCIRC
= -3 诊断斜压ADCIRC在温度强制下运行
= -2 在盐度强迫下运行诊断斜压ADCIRC
= -1带有Sigma T强制的诊断斜压ADCIRC
= 0 正压模型运行
= 1 使用Sigma T强制运行的预后斜压ADCIRC
= 2 用盐度强迫运行预后的斜压ADCIRC
= 3 预热斜压ADCIRC在温度强制下运行
= 4 在盐度和温度强迫下运行的预后斜压ADCIRC

NOLIBF =用于控制2DDI ADCIRC运行中使用的底部应力参数化类型的参数。必须指定此参数，但在3D运行中将忽略该参数。注意：在``NWP''部分中，如果用户选择quadratic_friction_coefficient_at_sea_floor，mannings_n_at_sea_floor或chezy_friction_coefficient_at_sea_floor，则NOLIBF必须为1（非线性摩擦公式），因为所有这些公式都是非线性的。如果NOLIBF的值不是1，则是错误会导致ADCIRC停止。
= 0 线性底部摩擦定律。摩擦系数（FFACTOR）在下面指定。
= 1 二次底摩擦定律。摩擦系数（FFACTOR）在下面指定。
= 2 混合非线性底部摩擦定律。在深水中，摩擦系数是恒定的，并产生二次底摩擦定律。在浅水中，摩擦系数随着深度的减小而增加（例如曼宁型摩擦定律）。摩擦系数确定为：FFACTOR ＝ FFACTORMIN * [1+（HBREAK / H）** FTHETA] **（FGAMMA / FTHETA）。必需的参数（FFACTORMIN，HBREAK，FTHETA，FGAMMA）
在下面指定。
对于3D ADCIRC运行，应使用bottom_roughness_length节点属性指定空间变化的底部摩擦。

 

NOLIFA =控制ADCIRC中的有限幅度项的参数。 NOLIFA的值影响最小水深参数（H0）的含义，并要求指定其他参数以及H0（请参见下文）。打开有限幅项时，对流项的时间导数部分也应打开，以保持适当的质量守恒和一致性（即，当NOLIFA> 0时，则NOLICAT = 1）。
= 0有限振幅项不包括在模型运行中（即，除连续性方程中的瞬变项外，所有深度均使用测深而非总深度对深度进行了线性化处理）
禁用元素的润湿和干燥。假定初始水深等于网格文件中指定的测深水深（fort.14）。
模型运行中包含= 1个有限振幅项，并且禁用了元素的润湿和干燥。假定初始水深等于网格文件中指定的测深水深（fort.14）。
= 2个有限幅度项包含在模型运行中，并且启用了元素的润湿和干燥。假定初始水深等于网格文件中指定的测深水深（fort.14）。

 

NOLICA =控制ADCIRC中对流项的参数（时间导数部分以连续性方程的GWCE形式出现，并由NOLICAT控制）。
当包括对流项的这些（空间导数）部分时，GWCE中对流项的时间导数部分也应包括在内（即，NOLICA = 1，NOLICAT = 1时）。
= 0对流术语不包括在计算中
= 1个对流词包含在计算中
NOLICA和NOLICAT参数本身将在整个域上激活或取消对流条款; 但是，这些术语也可以在逐个元素的基础上激活或停用。 有关详细信息，请参阅节点属性文件（fort.13）的文档中advection_state节点属性的定义。

 

NOLICAT =参数，用于控制对流项的时间导数部分，该对数项以ADCIRC中的连续性方程的GWCE形式出现。 GWCE中的其余对流项和动量方程中的整个对流项均由NOLICA控制。如果包括有限振幅或对流项的其余部分以保持质量守恒和解的一致性，则应包括这些项。
= 0在GWCE连续性方程式中出现的对流项的时间导数部分不包括在计算中。
= 1在计算中包括在GWCE连续性方程中出现的对流项的时间导数部分。
NOLICA和NOLICAT参数本身将在整个域上激活或取消对流条款;但是，这些术语也可以在逐个元素的基础上激活或停用。有关详细信息，请参阅节点属性文件（fort.13）的文档中advection_state节点属性的定义。

 

NWP =运行中使用的节点属性的数量。 节点属性是网格中每个节点的属性，在空间上变化，但时间上恒定。 有关示例，请参见AttrName。 用户必须在节点属性文件（fort.13）中提供节点属性数据本身。 如果NWP不为零，则节点属性的名称出现在以下各行中，每行一个：

FOR j=1 to NWP

AttrName(j)

end j loop

 

NCOR =参数，用于控制科里奥利参数在空间中是恒定的，还是从网格中的节点的y坐标（假定为纬度）中计算得出的，是向下读取还是在空间上变化。 栅格坐标系由ICS参数指定（请参见上文）。
= 0，以读取空间常数的科里奥利参数
= 1，以计算空间可变的科里奥利参数

 

NTIP =参数，用于控制是否使用潮汐势和自吸/负载潮汐强迫
强制ADCIRC。
= 0，未使用潮汐势和自吸力/潮汐力
= 1，使用了潮汐强迫作用
= 2，使用了潮汐势和自吸力/潮汐力。 在这种情况下，会从``自我吸引力/地球负荷潮汐强迫文件''中读取网格中每个节点上每个组成部分的自我吸引力/负荷潮汐信息。

 

NWS

 

NRAMP =斜坡选项参数，用于控制是否将斜坡应用于ADCIRC强制功能。
= 0强制功能不使用斜坡功能；冷启动后立即施加全力强迫。
= 1指定持续时间（DRAMP，相对于冷启动时间的天数）的单个双曲正切斜率函数将应用于所有强制。有关斜坡函数的更多信息，请参见DRAMP行的描述。
= 2与NRAMP = 1相同，除了指定的持续时间（DRAMPExtFlux，相对于冷启动时间加上FluxSettlingTime的天数）的第二个单独的双曲正切斜率，还将专门读取外部通量强迫（例如，河流边界条件） DRAMP生产线。此外，还将在DRAMP线上指定IBTYPE = 52河沿的FluxSettlingTime参数。如果网格（fort.14）文件中没有IBTYPE = 52边界，则将读取但忽略FluxSettlingTime。有关更多信息，请参见DRAMP的描述。
= 3与NRAMP = 2相同，除了指定持续时间（DRAMPIntFlux，相对于冷启动时间加上FluxSettlingTime的天数）的第三个单独的双曲正切斜率，还将专门用于内部通量强迫（例如，在堤坝和涵洞下的流量）在DRAMP线上读取。有关更多信息，请参见DRAMP行的描述。
= 4与NRAMP = 3相同，不同之处在于，还将读取指定持续时间的第四个单独的双曲正切斜率（DRAMPElev，相对于冷启动时间加上FluxSettlingTime的天数），专门用于标高指定的边界强迫（例如潮汐边界）。 DRAMP生产线。有关更多信息，请参见DRAMP行的描述。
= 5与NRAMP = 4相同，不同之处在于，还将在DRAMP线上读取专门用于潮汐强迫的指定持续时间（DRAMPTip，相对于冷启动时间加上FluxSettlingTime的天数）的第五个单独的双曲正切斜率。有关更多信息，请参见DRAMP行的描述。
= 6与NRAMP = 5相同，除了指定的持续时间（DRAMPMete，相对于冷启动时间加上FluxSettlingTime的天数）的第六个单独的双曲正切斜率也将被读取（特别是气象强迫（即风和大气压）） DRAMP生产线。有关更多信息，请参见DRAMP行的描述。
= 7与NRAMP = 6相同，不同之处在于，还将在DRAMP线上读取指定持续时间的第七个单独的双曲线正切斜率（DRAMPWRad，相对于冷启动时间加上FluxSettlingTime的天数），专门用于波浪辐射应力强迫。有关更多信息，请参见DRAMP行的描述。
= 8与NRAMP = 7相同，除了也会从DRAMP行读取延迟参数（相对于冷启动时间加上FluxSettlingTime的天数，DUnRampMete）。气象坡道延迟参数DUnRampMete在以下情况下很有用：从冷启动开始，从长期无气象的潮汐加速中启动气象强制运行。相对于ADCIRC冷启动时间，气象斜坡延迟参数将气象斜坡的应用启动延迟了指定的时间长度。有关更多信息，请参见DRAMP行的描述。

 

G =重力常数。 该常数的单位决定了ADCIRC运算的距离单位。 （ADCIRC始终以秒为单位运行，因此G的时间单位必须为秒。）ICS = 2时，要求G = 9.81 m / sec2。 无论ICS是什么，当NTIP = 1或NCOR = 1时，都要求G = 9.81 m / sec2。

 

TAU0 = =广义波连续性方程（GWCE）加权因子，用于加权GWCE的原始部分和波动部分的相对贡献。如果在上面的Fort.15文件中将“ primitive_weighting_in_continuity_equation”指定为节点属性，则将读入此行，但将忽略该行。如果未使用节点属性文件，或者在节点属性（fort.13）文件中包含“ primitive_weighting_in_continuity_equation”，但在Fort.15文件中未指定该参数，则将使用此TAU0参数。

= 0，GWCE是纯波动方程。
<1，GWCE的行为就像一个纯原始连续性方程。设置TAU0的一个好的经验法则是将其设置为等于等效线性摩擦系数的最大值（例如，对于线性摩擦TAU0 = TAU；对于二次摩擦TAU0 =最大值（速度* CF /深度）。典型值TAU0在0.005 – 0.1的范围内。
= -1，TAU0在空间上变化但时间恒定;根据深度计算如下：如果深度> = 10，则TAU0设置为0.005，如果深度<10，则TAU0设置为0.020。
= -2，TAU0在空间上变化但时间恒定;根据深度计算如下：如果深度> = 200，则TAU0设置为0.005;如果深度<200但> 1，则TAU0设置为1 /深度;如果深度<1，则设置TAU0。到1.0。
= -3，则TAU0随时间变化。 TAU0是从节点属性文件中读取的TAU0Base计算得出的。
如果TAU0Base <0.025; TAU0 = TAU0Base（时间常数）
如果TAU0Base> = 0.025; TAU0 = TAU0Base + 1.5 TK（i）其中TK（i）= Cd | U | / H
可以使用ADCIRC实用程序tau0_gen.f生成TAU0Base值。该程序基于以下逻辑进行生成：
如果平均。节点与其邻居之间的距离<1750 m TAU0Base = 0.03
如果平均。节点与其邻居之间的距离> 1750 m且深度<10m； TAU0Base = 0.02
如果平均。节点与其邻居之间的距离> 1750 m，深度> 10m； TAU0Base = 0.005
= -5，FullDomainTimeVaryingTau0 =。 TAU0在空间和时间上变化，并取决于局部摩擦;它限制在Tau0FullDomainMin和Tau0FullDomainMax指定的范围内。
Tau0 = Tau0Min + 1.5 * TK（i）
对于随时间变化的tau0公式，ADCIRC能够写出内部计算的tau0值。这些值将写入到Fort.90文件中，该文件的格式和输出频率与水面高程输出文件（fort.63）相同。通过在Fort.15文件中tau0输入值的十分之一处放置1，来指定fort.90文件的产生。例如，如果tau0 = -3.1，仍然根据上面对tau0 = -3的描述进行tau0的计算，并且还会生成fort.90输出文件。
 

Tau0FullDomainMin，Tau0FullDomainMax –仅当TAU0 = -5 –时空变化的TAU0方案必须在其之间的指定值时才包括此行。 建议值为Tau0FullDomainMin = 0.005和Tau0FullDomainMax = 0.2。

 

DTDP = ADCIRC时间步长（以秒为单位）。 注意：模型中的时间计算为：TIME = STATIM * 86400. + DTDP * IT。
> 0 =不使用预测器-校正器算法。
<0 =使用预测器-校正器算法。