fvcom inlet例子的运行（错误）

来源：互联网发布：js开源特效编辑：程序博客网时间：2024/04/28 14:47

8月9日更新：由于在写此文时所用的estuary的make.inc是错误的，所以，此文有错误，留作纪念算了。

运行inlet例子后显示

!================================================================!
   _______ _     _ _______ _______ _______ ______     _____
(_______)(_)   (_)(_______)(_______)(_______)(_____ \   (_____)
   _____    _     _ _        _     _ _ _ _ _____) ) _ __ _
| ___) | |   | || |      | |   | || ||_|| |(_____ ( | |/ /| |
| |       \ \ / / | |_____ | |___| || |   | | _____) )_|   /_| |
|_|        \___/   \______) \_____/ |_|   |_|(______/(_)\_____/
-- Beta Release
!================================================================!
!                                                                !
!========DOMAIN DECOMPOSITION USING: METIS 4.0.1 ================!
!======Copyright 1998, Regents of University of Minnesota========!
!                                                                !
&NML_NETCDF
NC_ON   = T,
NC_FIRST_OUT    = cycle=0                                                                         ,
NC_OUT_INTERVAL = cycle=72                                                                        ,
NC_OUTPUT_STACK =           0,
NC_SUBDOMAIN_FILES      = FVCOM
               ,
NC_GRID_METRICS = F,
NC_FILE_DATE    = F,
NC_VELOCITY     = T,
NC_SALT_TEMP    = F,
NC_TURBULENCE   = F,
NC_AVERAGE_VEL = T,
NC_VERTICAL_VEL = T,
NC_WIND_VEL     = F,
NC_WIND_STRESS = F,
NC_EVAP_PRECIP = F,
NC_SURFACE_HEAT = F,
NC_GROUNDWATER = F,
NC_BIO = F,
NC_WQM = F,
NC_VORTICITY    = F
/
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
FVCOM Fatal Error!
Can Not Read NameList NML_NETCDF from file: ./inlet_run.nml
Stopping FVCOM
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[rzf@localhost run]$

我知道estuary和river plume的例子是可以运行成功的所以我比较一下这三个namlist中的NETCDF配置的异同，发现estuary 和 river plume的NETCDF配置一样，而inlet的配置多了两个东西，即

NC_WAVE_PARA = T,
NC_WAVE_STRESS = T

删除后，可运行成功。

回头到手册中查找 NC_WAVE_PARA = T, NC_WAVE_STRESS = T的意思，有he一段话“the selection needs to be chosen when the surface wave model is selected for the case with either only waves or current-wave interactions”

所以猜想在生成fvcom可执行文件之前的make.inc 中要取消对于

FLAG_32 = -DWAVE_CURRENT_INTERACTION 或者FLAG_35 = -DWAVE_ONLY的注释

结果如我所料，取消注释之后编译fvcom成功运行成功

0 0