Roms can not write down the .rst file

[zhqianj]

Hello everyone
I am new to ROMS, and I encounter the following issure. It states that: DEF_RST - Illegal output type, io_type = 0. Many many thanks for your help.


The log file states:
--------------------------------------------------------------------------------
Model Input Parameters: ROMS/TOMS version 4.0
Tuesday - July 6, 2021 - 10:22:00 AM
--------------------------------------------------------------------------------

tide_bore

Operating system : Linux
CPU/hardware : x86_64
Compiler system : gfortran
Compiler command : /usr/bin/mpif90
Compiler flags : -frepack-arrays -O3 -ffast-math -
OCN Communicator : 1140850688, PET size = 4

Input Script : ./roms_tidebore.in

SVN Root URL : https://www.myroms.org/svn/src/trunk
SVN Revision : 1077

Local Root : /mnt/e/numerical_model/roms_1053
Header Dir : /mnt/d/tidal_bore/tidalbore_small
Header file : tide_bore.h
Analytical Dir : /mnt/d/tidal_bore/tidalbore_small

Resolution, Grid 01: 698x178x10, Parallel Nodes: 4, Tiling: 2x2


Physical Parameters, Grid: 01
=============================

25000 ntimes Number of timesteps for 3-D equations.
60.000 dt Timestep size (s) for 3-D equations.
5 ndtfast Number of timesteps for 2-D equations between
each 3D timestep.
1 ERstr Starting ensemble/perturbation run number.
1 ERend Ending ensemble/perturbation run number.
0 nrrec Number of restart records to read from disk.
T LcycleRST Switch to recycle time-records in restart file.
30 nRST Number of timesteps between the writing of data
into restart fields.
1 ninfo Number of timesteps between print of information
to standard output.
T ldefout Switch to create a new output NetCDF file(s).
10 nHIS Number of timesteps between the writing fields
into history file.
0 nQCK Number of timesteps between the writing fields
into quicksave file.
F LuvSponge Turning OFF sponge on horizontal momentum.
F LtracerSponge(01) Turning OFF sponge on tracer 01: temp
F LtracerSponge(02) Turning OFF sponge on tracer 02:
5.0000E-06 Akt_bak(01) Background vertical mixing coefficient (m2/s)
for tracer 01: temp
1.0000E-06 Akt_bak(02) Background vertical mixing coefficient (m2/s)
for tracer 02:
5.0000E-05 Akv_bak Background vertical mixing coefficient (m2/s)
for momentum.
5.0000E-06 Akk_bak Background vertical mixing coefficient (m2/s)
for turbulent energy.
5.0000E-06 Akp_bak Background vertical mixing coefficient (m2/s)
for turbulent generic statistical field.
3.000 gls_p GLS stability exponent.
1.500 gls_m GLS turbulent kinetic energy exponent.
-1.000 gls_n GLS turbulent length scale exponent.
7.6000E-06 gls_Kmin GLS minimum value of turbulent kinetic energy.
1.0000E-12 gls_Pmin GLS minimum value of dissipation.
5.4770E-01 gls_cmu0 GLS stability coefficient.
1.4400E+00 gls_c1 GLS shear production coefficient.
1.9200E+00 gls_c2 GLS dissipation coefficient.
-4.0000E-01 gls_c3m GLS stable buoyancy production coefficient.
1.0000E+00 gls_c3p GLS unstable buoyancy production coefficient.
1.0000E+00 gls_sigk GLS constant Schmidt number for TKE.
1.3000E+00 gls_sigp GLS constant Schmidt number for PSI.
1400.000 charnok_alpha Charnok factor for Zos calculation.
0.500 zos_hsig_alpha Factor for Zos calculation using Hsig(Awave).
0.250 sz_alpha Factor for Wave dissipation surface tke flux .
100.000 crgban_cw Factor for Craig/Banner surface tke flux.
3.0000E-04 rdrg Linear bottom drag coefficient (m/s).
1.0000E-03 rdrg2 Quadratic bottom drag coefficient.
2.0000E-02 Zob Bottom roughness (m).
2.0000E-02 Zos Surface roughness (m).
1.0000E-01 Dcrit Minimum depth for wetting and drying (m).
1 Vtransform S-coordinate transformation equation.
1 Vstretching S-coordinate stretching function.
1.0000E+00 theta_s S-coordinate surface control parameter.
1.0000E+00 theta_b S-coordinate bottom control parameter.
1.000 Tcline S-coordinate surface/bottom layer width (m) used
in vertical coordinate stretching.
1025.000 rho0 Mean density (kg/m3) for Boussinesq approximation.
0.000 dstart Time-stamp assigned to model initialization (days).
0.000 tide_start Reference time origin for tidal forcing (days).
0.00 time_ref Reference time for units attribute (yyyymmdd.dd)
0.0000E+00 Tnudg(01) Nudging/relaxation time scale (days)
for tracer 01: temp
0.0000E+00 Tnudg(02) Nudging/relaxation time scale (days)
for tracer 02:
0.0000E+00 Znudg Nudging/relaxation time scale (days)
for free-surface.
0.0000E+00 M2nudg Nudging/relaxation time scale (days)
for 2D momentum.
0.0000E+00 M3nudg Nudging/relaxation time scale (days)
for 3D momentum.
0.0000E+00 obcfac Factor between passive and active
open boundary conditions.
F VolCons(1) NLM western edge boundary volume conservation.
F VolCons(2) NLM southern edge boundary volume conservation.
F VolCons(3) NLM eastern edge boundary volume conservation.
F VolCons(4) NLM northern edge boundary volume conservation.
14.000 T0 Background potential temperature (C) constant.
35.000 S0 Background salinity (PSU) constant.
1027.000 R0 Background density (kg/m3) used in linear Equation
of State.
1.7000E-04 Tcoef Thermal expansion coefficient (1/Celsius).
0.0000E+00 Scoef Saline contraction coefficient (1/PSU).
1.000 gamma2 Slipperiness variable: free-slip (1.0) or
no-slip (-1.0).
F LuvSrc Turning OFF momentum point Sources/Sinks.
F LwSrc Turning OFF volume influx point Sources/Sinks.
F LtracerSrc(01) Turning OFF point Sources/Sinks on tracer 01: temp
F LtracerSrc(02) Turning OFF point Sources/Sinks on tracer 02:
F LsshCLM Turning OFF processing of SSH climatology.
F Lm2CLM Turning OFF processing of 2D momentum climatology.
F Lm3CLM Turning OFF processing of 3D momentum climatology.
F LtracerCLM(01) Turning OFF processing of climatology tracer 01: temp
F LtracerCLM(02) Turning OFF processing of climatology tracer 02:
F LnudgeM2CLM Turning OFF nudging of 2D momentum climatology.
F LnudgeM3CLM Turning OFF nudging of 3D momentum climatology.
F LnudgeTCLM(01) Turning OFF nudging of climatology tracer 01: temp
F LnudgeTCLM(02) Turning OFF nudging of climatology tracer 02:

T Hout(idFsur) Write out free-surface.
T Hout(idUbar) Write out 2D U-momentum component.
T Hout(idVbar) Write out 2D V-momentum component.
T Hout(idUvel) Write out 3D U-momentum component.
T Hout(idVvel) Write out 3D V-momentum component.
T Hout(idWvel) Write out W-momentum component.
T Hout(idOvel) Write out omega vertical velocity.

1 inp_lib Using standard NetCDF library for input files.
1 out_lib Using standard NetCDF library for output files.

Output/Input Files:

Output Restart File:
Output History File: tidebore_his.nc
Input Grid File: tidebore_grd.nc
Tidal Forcing File: tidebore_tideforcing.nc
Input Forcing File 01: tidebore_frc.nc
Input Lateral Boundary File 01: tidebore_bry.nc
ROMS I/O variables Metadata File: /mnt/e/numerical_model/roms_1053/ROMS/External/varinfo.dat
Input/Output USER File: MyFile.dat

Generic User Parameters:

1.0000E+00 user(01) User parameter 01.

Tile partition information for Grid 01: 698x178x10 tiling: 2x2

tile Istr Iend Jstr Jend Npts

0 1 349 1 89 310610
1 350 698 1 89 310610
2 1 349 90 178 310610
3 350 698 90 178 310610

Tile minimum and maximum fractional coordinates for Grid 01:
(interior points only)

tile Xmin Xmax Ymin Ymax grid

0 0.50 349.50 0.50 89.50 RHO-points
1 349.50 698.50 0.50 89.50 RHO-points
2 0.50 349.50 89.50 178.50 RHO-points
3 349.50 698.50 89.50 178.50 RHO-points

0 1.00 349.50 0.50 89.50 U-points
1 349.50 698.00 0.50 89.50 U-points
2 1.00 349.50 89.50 178.50 U-points
3 349.50 698.00 89.50 178.50 U-points

0 0.50 349.50 1.00 89.50 V-points
1 349.50 698.50 1.00 89.50 V-points
2 0.50 349.50 89.50 178.00 V-points
3 349.50 698.50 89.50 178.00 V-points

Maximum halo size in XI and ETA directions:

HaloSizeI(1) = 722
HaloSizeJ(1) = 202
TileSide(1) = 355
TileSize(1) = 33725


Tracer Advection Scheme: NLM
========================

Variable Grid Horizontal Vertical
--------- ---- ------------ ------------

temp 1 Upstream3 Centered4

1 Upstream3 Centered4

Akima4 Fourth-order Akima advection
Centered2 Second-order centered differences advection
Centered4 Fourth-order centered differences advection
HSIMT Third High-order Spatial Inteporlation at Middle Time Advection with TVD limiter
MPDATA Multidimensional Positive Definite Advection Algorithm, recursive method
Splines Conservative Parabolic Splines Reconstruction Advection (only vertical; not recommended)
Split_U3 Split third-order Upstream Advection
Upstream3 Third-order Upstream-biased Advection (only horizontal)


Lateral Boundary Conditions: NLM
============================

Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------

zeta 1 Closed Closed Chapman Imp Closed

ubar 1 Closed Closed Flather Closed

vbar 1 Closed Closed Flather Closed

u 1 Closed Closed Radiation Closed

v 1 Closed Closed Radiation Closed

temp 1 Closed Closed Radiation Closed

1 Closed Closed Radiation Closed

tke 1 Closed Closed Radiation Closed

Activated C-preprocessing Options:

tide_bore tide_bore
ANA_BSFLUX Analytical kinematic bottom salinity flux
ANA_BTFLUX Analytical kinematic bottom temperature flux
ANA_INITIAL Analytical initial conditions
ASSUMED_SHAPE Using assumed-shape arrays
BOUNDARY_ALLREDUCE Using mpi_allreduce in mp_boundary routine
COLLECT_ALLREDUCE Using mpi_allreduce in mp_collect routine
DEBUGGING Internal debugging switch activated
DOUBLE_PRECISION Double precision arithmetic numerical kernel.
GLS_MIXING Generic Length-Scale turbulence closure
LIMIT_BSTRESS Limit bottom stress to maintain bottom velocity direction
KANTHA_CLAYSON Kantha and Clayson stability function formulation
MASKING Land/Sea masking
MPI MPI distributed-memory configuration
NONLINEAR Nonlinear Model
!NONLIN_EOS Linear Equation of State for seawater
N2S2_HORAVG Horizontal smoothing of buoyancy and shear
POWER_LAW Power-law shape time-averaging barotropic filter
PRSGRD31 Standard density Jacobian formulation (Song, 1998)
PROFILE Time profiling activated
K_GSCHEME Third-order upstream advection of TKE fields
RAMP_TIDES Ramping tidal forcing for one day
REDUCE_ALLREDUCE Using mpi_allreduce in mp_reduce routine
RI_SPLINES Parabolic Spline Reconstruction for Richardson Number
RHO_SURF Include difference between rho0 and surface density
RST_SINGLE Single precision fields in restart NetCDF file
SOLVE3D Solving 3D Primitive Equations
SPLINES_VDIFF Parabolic Spline Reconstruction for Vertical Diffusion
SPLINES_VVISC Parabolic Spline Reconstruction for Vertical Viscosity
SSH_TIDES Add tidal elevation to SSH climatology
TS_FIXED Diagnostic configuration, no evolution of tracer
UV_ADV Advection of momentum
UV_COR Coriolis term
UV_U3HADVECTION Third-order upstream horizontal advection of 3D momentum
UV_C4VADVECTION Fourth-order centered vertical advection of momentum
UV_QDRAG Quadratic bottom stress
UV_TIDES Add tidal currents to 2D momentum climatologies
VAR_RHO_2D Variable density barotropic mode
WET_DRY Wetting and drying activated

Process Information:

Node # 0 (pid= 17943) is active.
Node # 1 (pid= 17944) is active.
Node # 2 (pid= 17945) is active.
Node # 3 (pid= 17946) is active.

INITIAL: Configuring and initializing forward nonlinear model ...
*******

GET_GRID_NF90 - angle between XI-axis and EAST: angler
(Grid = 01, File: tidebore_grd.nc)
(Min = 3.26211878E-01 Max = 3.33637132E-01)
GET_GRID_NF90 - bathymetry at RHO-points: h
(Grid = 01, File: tidebore_grd.nc)
(Min = 1.00000000E+00 Max = 3.24228767E+01)
GET_GRID_NF90 - Coriolis parameter at RHO-points: f
(Grid = 01, File: tidebore_grd.nc)
(Min = 7.26867133E-05 Max = 7.42329861E-05)
GET_GRID_NF90 - reciprocal XI-grid spacing: pm
(Grid = 01, File: tidebore_grd.nc)
(Min = 4.92013714E-03 Max = 4.94395978E-03)
GET_GRID_NF90 - reciprocal ETA-grid spacing: pn
(Grid = 01, File: tidebore_grd.nc)
(Min = 5.24281944E-03 Max = 5.33272306E-03)
GET_GRID_NF90 - longitude of RHO-points: lon_rho
(Grid = 01, File: tidebore_grd.nc)
(Min = 1.19707389E+02 Max = 1.21209230E+02)
GET_GRID_NF90 - longitude of U-points: lon_u
(Grid = 01, File: tidebore_grd.nc)
(Min = 1.19708388E+02 Max = 1.21208232E+02)
GET_GRID_NF90 - longitude of V-points: lon_v
(Grid = 01, File: tidebore_grd.nc)
(Min = 1.19707686E+02 Max = 1.21208934E+02)
GET_GRID_NF90 - longitude of PSI-points: lon_psi
(Grid = 01, File: tidebore_grd.nc)
(Min = 1.19708684E+02 Max = 1.21207936E+02)
GET_GRID_NF90 - latitude of RHO-points lat_rho
(Grid = 01, File: tidebore_grd.nc)
(Min = 2.98942272E+01 Max = 3.05974376E+01)
GET_GRID_NF90 - latitude of U-points: lat_u
(Grid = 01, File: tidebore_grd.nc)
(Min = 2.98945205E+01 Max = 3.05971403E+01)
GET_GRID_NF90 - latitude of V-points: lat_v
(Grid = 01, File: tidebore_grd.nc)
(Min = 2.98950304E+01 Max = 3.05966189E+01)
GET_GRID_NF90 - latitude of PSI-points lat_psi
(Grid = 01, File: tidebore_grd.nc)
(Min = 2.98953238E+01 Max = 3.05963216E+01)
GET_GRID_NF90 - mask on RHO-points: mask_rho
(Grid = 01, File: tidebore_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID_NF90 - mask on U-points: mask_u
(Grid = 01, File: tidebore_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID_NF90 - mask on V-points: mask_v
(Grid = 01, File: tidebore_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)
GET_GRID_NF90 - mask on PSI-points: mask_psi
(Grid = 01, File: tidebore_grd.nc)
(Min = 0.00000000E+00 Max = 1.00000000E+00)

Vertical S-coordinate System, Grid 01:

level S-coord Cs-curve Z at hmin at hc half way at hmax

10 0.0000000 0.0000000 0.000 0.000 0.000 0.000
9 -0.1000000 -0.0889041 -0.100 -0.100 -1.497 -2.894
8 -0.2000000 -0.1848065 -0.200 -0.200 -3.104 -6.007
7 -0.3000000 -0.2864445 -0.300 -0.300 -4.800 -9.301
6 -0.4000000 -0.3921616 -0.400 -0.400 -6.561 -12.723
5 -0.5000000 -0.5000000 -0.500 -0.500 -8.356 -16.211
4 -0.6000000 -0.6078384 -0.600 -0.600 -10.150 -19.700
3 -0.7000000 -0.7135555 -0.700 -0.700 -11.911 -23.122
2 -0.8000000 -0.8151935 -0.800 -0.800 -13.608 -26.416
1 -0.9000000 -0.9110959 -0.900 -0.900 -15.215 -29.529
0 -1.0000000 -1.0000000 -1.000 -1.000 -16.711 -32.423

Time Splitting Weights for Grid 01: ndtfast = 5 nfast = 7
==================================

Primary Secondary Accumulated to Current Step

1 0.0119306869089341 0.2000000000000000 0.0119306869089341 0.2000000000000000
2 0.0475049776008206 0.1976138626182132 0.0594356645097548 0.3976138626182132
3 0.1047629417620785 0.1881128670980491 0.1641986062718333 0.5857267297162623
4 0.1760826281731956 0.1671602787456334 0.3402812344450289 0.7528870084618957
5 0.2416469636634479 0.1319437531109943 0.5819281981084767 0.8848307615728900
6 0.2602974116474963 0.0836143603783047 0.8422256097559730 0.9684451219511947
7 0.1577743902440272 0.0315548780488054 1.0000000000000002 1.0000000000000002

ndtfast, nfast = 5 7 nfast/ndtfast = 1.40000

Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):

1.000000000000 1.084198606272 0.542099303136 1.000000000000 1.000000000000

Power filter parameters, Fgamma, gamma = 0.28400 0.00000

Metrics information for Grid 01:
===============================

Minimum X-grid spacing, DXmin = 2.02267018E-01 km Water points = 2.02270761E-01 km
Maximum X-grid spacing, DXmax = 2.03246367E-01 km Water points = 2.03244460E-01 km
Minimum Y-grid spacing, DYmin = 1.87521457E-01 km Water points = 1.87577543E-01 km
Maximum Y-grid spacing, DYmax = 1.90737066E-01 km Water points = 1.90722487E-01 km
Minimum Z-grid spacing, DZmin = 1.00000000E-01 m Water points = 1.00000000E-01 m
Maximum Z-grid spacing, DZmax = 3.48859427E+00 m Water points = 3.48859427E+00 m

Minimum barotropic Courant Number = 2.70752784E-01
Maximum barotropic Courant Number = 1.54350222E+00
Maximum Coriolis Courant Number = 4.45333339E-03

ANA_INITIAL - vertically integrated u-momentum component: ubar
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - vertically integrated v-momentum component: vbar
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - free-surface: zeta
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - u-momentum component: u
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - v-momentum component: v
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
ANA_INITIAL - potential temperature: temp
(Grid = 01, Min = 1.40000000E+01 Max = 1.40000000E+01)
ANA_INITIAL - :
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD_NF90 - tidal period
(Grid = 01, Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD_NF90 - tidal elevation amplitude
(Grid = 01, Min = 0.00000000E+00 Max = 2.71615040E+00)
GET_2DFLD_NF90 - tidal elevation phase angle
(Grid = 01, Min = 0.00000000E+00 Max = 6.27040629E+00)
GET_2DFLD_NF90 - tidal current inclination angle
(Grid = 01, Min = 0.00000000E+00 Max = 6.22587082E+00)
GET_2DFLD_NF90 - tidal current phase angle
(Grid = 01, Min = 0.00000000E+00 Max = 6.22587082E+00)
GET_2DFLD_NF90 - maximum tidal current, ellipse major axis
(Grid = 01, Min = 0.00000000E+00 Max = 1.88859173E+00)
GET_2DFLD_NF90 - minimum tidal current, ellipse minor axis
(Grid = 01, Min = -3.80372082E-01 Max = 1.31444068E-01)
GET_2DFLD_NF90 - surface u-momentum stress, 0001-01-01 00:00:00.00
(Grid=01, Rec=1, Index=1, File: tidebore_frc.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 0.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00) regrid = F
GET_2DFLD_NF90 - surface v-momentum stress, 0001-01-01 00:00:00.00
(Grid=01, Rec=1, Index=1, File: tidebore_frc.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 0.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00) regrid = F
GET_2DFLD_NF90 - surface net heat flux, 0001-01-01 00:00:00.00
(Grid=01, Rec=1, Index=1, File: tidebore_frc.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 0.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00) regrid = F
GET_NGFLD_NF90 - free-surface eastern boundary condition, 0001-01-01 00:00:00.00
(Grid= 01, Rec=1, Index=2, File: tidebore_bry.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 0.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD_NF90 - 2D u-momentum eastern boundary condition, 0001-01-01 00:00:00.00
(Grid= 01, Rec=1, Index=2, File: tidebore_bry.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 0.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD_NF90 - 2D v-momentum eastern boundary condition, 0001-01-01 00:00:00.00
(Grid= 01, Rec=1, Index=2, File: tidebore_bry.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 0.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00)

Basin information for Grid 01:

Maximum grid stiffness ratios: rx0 = 2.000000E-01 (Beckmann and Haidvogel)
rx1 = 4.273600E+00 (Haney)

Initial domain volumes: TotVolume = 6.0672671833E+09 m3
MinCellVol = 3.8235622845E+03 m3
MaxCellVol = 1.3439476784E+05 m3
Max/Min = 3.5149098627E+01


NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000001 - 000000025000)

GET_2DFLD_NF90 - surface u-momentum stress, 0001-01-16 00:00:00.00
(Grid=01, Rec=2, Index=2, File: tidebore_frc.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 15.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00) regrid = F
GET_2DFLD_NF90 - surface v-momentum stress, 0001-01-16 00:00:00.00
(Grid=01, Rec=2, Index=2, File: tidebore_frc.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 15.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00) regrid = F
GET_2DFLD_NF90 - surface net heat flux, 0001-01-16 00:00:00.00
(Grid=01, Rec=2, Index=2, File: tidebore_frc.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 15.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00) regrid = F
GET_NGFLD_NF90 - free-surface eastern boundary condition, 0001-01-16 00:00:00.00
(Grid= 01, Rec=2, Index=1, File: tidebore_bry.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 15.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD_NF90 - 2D u-momentum eastern boundary condition, 0001-01-16 00:00:00.00
(Grid= 01, Rec=2, Index=1, File: tidebore_bry.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 15.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD_NF90 - 2D v-momentum eastern boundary condition, 0001-01-16 00:00:00.00
(Grid= 01, Rec=2, Index=1, File: tidebore_bry.nc)
(Tmin= 0.0000 Tmax= 180.0000) t = 15.0000
(Min = 0.00000000E+00 Max = 0.00000000E+00)

TIME-STEP YYYY-MM-DD hh:mm:ss.ss KINETIC_ENRG POTEN_ENRG TOTAL_ENRG NET_VOLUME
C => (i,j,k) Cu Cv Cw Max Speed

0 0001-01-01 00:00:00.00 0.000000E+00 2.666835E+01 2.666835E+01 9.713888E+09
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS_NF90 - creating history file, Grid 01: tidebore_his.nc
WRT_HIS_NF90 - wrote history fields (Index=1,1) in record = 1
DEF_RST - Illegal output type, io_type = 0
Found Error: 03 Line: 85 Source: ROMS/Utility/def_rst.F
Found Error: 03 Line: 576 Source: ROMS/Nonlinear/output.F
Found Error: 03 Line: 518 Source: ROMS/Nonlinear/main3d.F
Found Error: 03 Line: 301 Source: ROMS/Drivers/nl_roms.h, ROMS_run

Elapsed wall CPU time for each process (seconds):

Node # 0 CPU: 0.938
Total: 3.969
Average: 0.992
Minimum: 0.938
Maximum: 1.031

Nonlinear model elapsed CPU time profile, Grid: 01

Allocation and array initialization .............. 0.922 (23.2283 %)
Ocean state initialization ....................... 0.062 ( 1.5748 %)
Reading of input data ............................ 0.734 (18.5039 %)
Computation of global information integrals ...... 0.062 ( 1.5748 %)
2D/3D coupling, vertical metrics ................. 0.172 ( 4.3307 %)
Omega vertical velocity .......................... 0.062 ( 1.5748 %)
Total: 2.016 50.7874 %

Unique kernel(s) regions profiled ................ 2.016 50.7874 %
Residual, non-profiled code ...................... 1.953 49.2126 %


All percentages are with respect to total time = 3.969


MPI communications profile, Grid: 01

Message Passage: 2D halo exchanges ............... 0.016 ( 0.3937 %)
Message Passage: 3D halo exchanges ............... 0.016 ( 0.3937 %)
Message Passage: data broadcast .................. 0.938 (23.6220 %)
Message Passage: data reduction .................. 0.047 ( 1.1811 %)
Message Passage: data gathering .................. 0.375 ( 9.4488 %)
Message Passage: data scattering.................. 0.516 (12.9921 %)
Total: 1.906 48.0315 %

>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>

Dynamic and Automatic memory (MB) usage for Grid 01: 698x178x10 tiling: 2x2

tile Dynamic Automatic USAGE MPI-Buffers

0 142.73 22.18 164.91 11.87
1 143.14 22.18 165.31 11.87
2 144.27 22.18 166.45 11.87
3 144.68 22.18 166.86 11.87

TOTAL 574.83 88.70 663.53 47.48

<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
Found Error: 03 Line: 85 Source: ROMS/Utility/close_io.F, close_file_nf90
Found Error: 03 Line: 121 Source: ROMS/Utility/close_io.F, close_inp
Found Error: 03 Line: 85 Source: ROMS/Utility/close_io.F, close_file_nf90
Found Error: 03 Line: 85 Source: ROMS/Utility/close_io.F, close_file_nf90
Found Error: 03 Line: 85 Source: ROMS/Utility/close_io.F, close_file_nf90

ROMS/TOMS - Output NetCDF summary for Grid 01:
number of time records written in HISTORY file = 1

Analytical header files used:

ROMS/Functionals/ana_btflux.h
ROMS/Functionals/ana_initial.h
Found Error: 03 Line: 381 Source: ROMS/Utility/close_io.F, close_out

ROMS/TOMS - Output error ............ exit_flag: 3


Node # 1 CPU: 1.031
Node # 2 CPU: 1.016
Node # 3 CPU: 0.984
ERROR: Abnormal termination: NetCDF OUTPUT.
REASON: No error

The .in file I used is :
! Application title.

TITLE = tide_bore

! C-preprocessing Flag.

MyAppCPP = tide_bore

! Input variable information file name. This file needs to be processed
! first so all information arrays can be initialized properly.

VARNAME = /mnt/e/numerical_model/roms_1053/ROMS/External/varinfo.dat

! Number of nested grids.

Ngrids = 1

! Number of grid nesting layers. This parameter is used to allow refinement
! and composite grid combinations.

NestLayers = 1

! Number of grids in each nesting layer [1:NestLayers].

GridsInLayer = 1

! Grid dimension parameters. See notes below in the Glossary for how to set
! these parameters correctly.

Lm == 698 ! Number of I-direction INTERIOR RHO-points
Mm == 178 ! Number of J-direction INTERIOR RHO-points
N == 10 ! Number of vertical levels

Nbed = 0 ! Number of sediment bed layers

NAT = 2 ! Number of active tracers (usually, 2)
NPT = 0 ! Number of inactive passive tracers
NCS = 0 ! Number of cohesive (mud) sediment tracers
NNS = 0 ! Number of non-cohesive (sand) sediment tracers

! Domain decomposition parameters for serial, distributed-memory or
! shared-memory configurations used to determine tile horizontal range
! indices (Istr,Iend) and (Jstr,Jend), [1:Ngrids].

NtileI == 2 ! I-direction partition
NtileJ == 2 ! J-direction partition

! Set horizontal and vertical advection schemes for active and inert
! tracers. A different advection scheme is allowed for each tracer.
! For example, a positive-definite (monotonic) algorithm can be activated
! for salinity and inert tracers, while a different one is set for
! temperature. [1:NAT+NPT,Ngrids] values are expected.
!
! Keyword Advection Algorithm
!
! A4 4th-order Akima (horizontal/vertical)
! C2 2nd-order centered differences (horizontal/vertical)
! C4 4th-order centered differences (horizontal/vertical)
! HSIMT 3th-order HSIMT-TVD (horizontal/vertical)
! MPDATA recursive flux corrected MPDATA (horizontal/vertical)
! SPLINES parabolic splines (only vertical)
! SU3 split third-order upstream (horizontal/vertical)
! U3 3rd-order upstream-biased (only horizontal)
!
! The user has the option of specifying the full Keyword or the first
! two letters, regardless if using uppercase or lowercase. If nested
! grids, specify values for each grid (see glossary below).

Hadvection == U3 \ ! temperature
U3 ! salinity

Vadvection == C4 \ ! temperature
C4 ! salinity

! Adjoint-based algorithms can have different horizontal and schemes
! for active and inert tracers.

ad_Hadvection == U3 \ ! temperature
U3 ! salinity

ad_Vadvection == C4 \ ! temperature
C4 ! salinity

LBC(isFsur) == Clo Clo Cha Clo ! free-surface
LBC(isUbar) == Clo Clo Fla Clo ! 2D U-momentum
LBC(isVbar) == Clo Clo Fla Clo ! 2D V-momentum
LBC(isUvel) == Clo Clo Rad Clo ! 3D U-momentum
LBC(isVvel) == Clo Clo Rad Clo ! 3D V-momentum
LBC(isMtke) == Clo Clo Rad Clo ! mixing TKE

LBC(isTvar) == Clo Clo Rad Clo \ ! temperature
Clo Clo Rad Clo ! salinity

! Adjoint-based algorithms can have different lateral boundary
! conditions keywords.

ad_LBC(isFsur) == Per Clo Per Clo ! free-surface
ad_LBC(isUbar) == Per Clo Per Clo ! 2D U-momentum
ad_LBC(isVbar) == Per Clo Per Clo ! 2D U-momentum
ad_LBC(isUvel) == Per Clo Per Clo ! 3D U-momentum
ad_LBC(isVvel) == Per Clo Per Clo ! 3D V-momentum
ad_LBC(isMtke) == Per Clo Per Clo ! mixing TKE

ad_LBC(isTvar) == Per Clo Per Clo \ ! temperature
Per Clo Per Clo ! salinity

! Set lateral open boundary edge volume conservation switch for
! nonlinear model and adjoint-based algorithms. Usually activated
! with radiation boundary conditions to enforce global mass
! conservation, except if tidal forcing is enabled. [1:Ngrids].

VolCons(west) == F ! western boundary
VolCons(east) == F ! eastern boundary
VolCons(south) == F ! southern boundary
VolCons(north) == F ! northern boundary

ad_VolCons(west) == F ! western boundary
ad_VolCons(east) == F ! eastern boundary
ad_VolCons(south) == F ! southern boundary
ad_VolCons(north) == F ! northern boundary

! Time-Stepping parameters.

NTIMES == 25000
DT == 60.00
NDTFAST == 5.00

! Number of timesteps for computing observation impacts during the
! analysis-forecast cycle.

NTIMES_ANA == 720 ! analysis interval
NTIMES_FCT == 720 ! forecast interval

! Model iteration loops parameters.

ERstr = 1
ERend = 1
Nouter = 1
Ninner = 1
Nintervals = 1

! Number of eigenvalues (NEV) and eigenvectors (NCV) to compute for the
! Lanczos/Arnoldi problem in the Generalized Stability Theory (GST)
! analysis. NCV must be greater than NEV (see documentation below).

NEV = 2 ! Number of eigenvalues
NCV = 10 ! Number of eigenvectors

! Input/Output parameters.

NRREC == 0
LcycleRST == T
NRST == 30
NSTA == 1
NFLT == 1
NINFO == 1

! Output history, quicksave, average, and diagnostic files parameters.

LDEFOUT == T
NHIS == 10
NDEFHIS == 0
NQCK == 0
NDEFQCK == 0
NTSAVG == 1
NAVG == 30
NDEFAVG == 0
NTSDIA == 5
NDIA == 0
NDEFDIA == 0

! Output tangent linear and adjoint models parameters.

LcycleTLM == F
NTLM == 720
NDEFTLM == 0
LcycleADJ == F
NADJ == 720
NDEFADJ == 0
NSFF == 720
NOBC == 720

! GST output and check pointing restart parameters.

LmultiGST = F ! one eigenvector per file
LrstGST = F ! GST restart switch
MaxIterGST = 500 ! maximum number of iterations
NGST = 10 ! check pointing interval

! Relative accuracy of the Ritz values computed in the GST analysis.

Ritz_tol = 1.0d-15

! Harmonic/biharmonic horizontal diffusion of tracer for nonlinear model
! and adjoint-based algorithms: [1:NAT+NPT,Ngrids].

TNU2 == 0.0d0 0.0d0 ! m2/s
TNU4 == 2*0.0d0 ! m4/s

ad_TNU2 == 0.0d0 0.0d0 ! m2/s
ad_TNU4 == 0.0d0 0.0d0 ! m4/s

! Harmonic/biharmonic, horizontal viscosity coefficient for nonlinear model
! and adjoint-based algorithms: [Ngrids].

VISC2 == 0.001d0 ! m2/s
VISC4 == 0.0d0 ! m4/s

ad_VISC2 == 0.0d0 ! m2/s
ad_VISC4 == 0.0d0 ! m4/s

! Logical switches (TRUE/FALSE) to increase/decrease horizontal viscosity
! and/or diffusivity in specific areas of the application domain (like
! sponge areas) for the desired application grid.

LuvSponge == F ! horizontal momentum
LtracerSponge == F F ! temperature, salinity, inert

! Vertical mixing coefficients for tracers in nonlinear model and
! basic state scale factor in adjoint-based algorithms: [1:NAT+NPT,Ngrids]

AKT_BAK == 5.0d-6 1.0d-6 ! m2/s

ad_AKT_fac == 1.0d0 1.0d0 ! nondimensional

! Vertical mixing coefficient for momentum for nonlinear model and
! basic state scale factor in adjoint-based algorithms: [Ngrids].

AKV_BAK == 5.0d-5 ! m2/s

ad_AKV_fac == 1.0d0 ! nondimensional

! Upper threshold values to limit vertical mixing coefficients computed
! from vertical mixing parameterizations. Although this is an engineering
! fix, the vertical mixing values inferred from ocean observations are
! rarely higher than this upper limit value.

AKT_LIMIT == 1.0d-3 1.0d-3 ! m2/s

AKV_LIMIT == 1.0d-3 ! m2/s

! Turbulent closure parameters.

AKK_BAK == 5.0d-6 ! m2/s
AKP_BAK == 5.0d-6 ! m2/s
TKENU2 == 0.0d0 ! m2/s
TKENU4 == 0.0d0 ! m4/s

! Generic length-scale turbulence closure parameters.

GLS_P == 3.0d0 ! K-epsilon
GLS_M == 1.5d0
GLS_N == -1.0d0
GLS_Kmin == 7.6d-6
GLS_Pmin == 1.0d-12

GLS_CMU0 == 0.5477d0
GLS_C1 == 1.44d0
GLS_C2 == 1.92d0
GLS_C3M == -0.4d0
GLS_C3P == 1.0d0
GLS_SIGK == 1.0d0
GLS_SIGP == 1.30d0

! Constants used in surface turbulent kinetic energy flux computation.

CHARNOK_ALPHA == 1400.0d0 ! Charnok surface roughness
ZOS_HSIG_ALPHA == 0.5d0 ! roughness from wave amplitude
SZ_ALPHA == 0.25d0 ! roughness from wave dissipation
CRGBAN_CW == 100.0d0 ! Craig and Banner wave breaking

! Constants used in momentum stress computation.

RDRG == 3.0d-04 ! m/s
RDRG2 == 1.0d-03 ! nondimensional
Zob == 0.02d0 ! m
Zos == 0.02d0 ! m

! Height (m) of atmospheric measurements for Bulk fluxes parameterization.

BLK_ZQ == 10.0d0 ! air humidity
BLK_ZT == 10.0d0 ! air temperature
BLK_ZW == 10.0d0 ! winds

! Minimum depth for wetting and drying.

DCRIT == 0.10d0 ! m

! Various parameters.

WTYPE == 1
LEVSFRC == 15
LEVBFRC == 1


Vtransform == 1 ! transformation equation
Vstretching == 1 ! stretching function


THETA_S == 1.0d0 ! surface stretching parameter
THETA_B == 1.0d0 ! bottom stretching parameter
TCLINE == 1.0d0 ! critical depth (m)

! Mean Density and Brunt-Vaisala frequency.

RHO0 = 1025.0d0 ! kg/m3
BVF_BAK = 1.0d-5 ! 1/s2

! Time-stamp assigned for model initialization, reference time
! origin for tidal forcing, and model reference time for output
! NetCDF units attribute.

DSTART = 0.0d0 ! days
TIDE_START = 0.0d0 ! days
TIME_REF = 0.0d0 ! yyyymmdd.dd

! Nudging/relaxation time scales, inverse scales will be computed
! internally, [1:Ngrids].

TNUDG == 2*0.0d0 ! days
ZNUDG == 0.0d0 ! days
M2NUDG == 0.0d0 ! days
M3NUDG == 0.0d0 ! days

! Factor between passive (outflow) and active (inflow) open boundary
! conditions, [1:Ngrids]. If OBCFAC > 1, nudging on inflow is stronger
! than on outflow (recommended).

OBCFAC == 0.0d0 ! nondimensional

! Linear equation of State parameters:

R0 == 1027.0d0 ! kg/m3
T0 == 14.0d0 ! Celsius
S0 == 35.0d0 ! nondimensional
TCOEF == 1.7d-4 ! 1/Celsius
SCOEF == 0.0d0 ! nondimensional

! Slipperiness parameter: 1.0 (free slip) or -1.0 (no slip)

GAMMA2 == 1.0d0

! Logical switches (TRUE/FALSE) to activate horizontal momentum transport
! point Sources/Sinks (like river runoff transport) and mass point
! Sources/Sinks (like volume vertical influx), [1:Ngrids].

LuvSrc == F ! horizontal momentum transport
LwSrc == F ! volume vertical influx

! Logical switches (TRUE/FALSE) to activate tracers point Sources/Sinks
! (like river runoff) and to specify which tracer variables to consider:
! [1:NAT+NPT,Ngrids]. See glossary below for details.

LtracerSrc == F F 12*F ! temperature, salinity, inert

! Logical switches (TRUE/FALSE) to read and process climatology fields.
! See glossary below for details.

LsshCLM == F ! sea-surface height
Lm2CLM == F ! 2D momentum
Lm3CLM == F ! 3D momentum

LtracerCLM == F F ! temperature, salinity, inert

! Logical switches (TRUE/FALSE) to nudge the desired climatology field(s).
! If not analytical climatology fields, users need to turn ON the logical
! switches above to process the fields from the climatology NetCDF file
! that are needed for nudging. See glossary below for details.

LnudgeM2CLM == F ! 2D momentum
LnudgeM3CLM == F ! 3D momentum

LnudgeTCLM == F F ! temperature, salinity, inert

! Starting (DstrS) and ending (DendS) day for adjoint sensitivity forcing.
! DstrS must be less or equal to DendS. If both values are zero, their
! values are reset internally to the full range of the adjoint integration.

DstrS == 0.0d0 ! starting day
DendS == 0.0d0 ! ending day

! Starting and ending vertical levels of the 3D adjoint state variables
! whose sensitivity is required.

KstrS == 1 ! starting level
KendS == 1 ! ending level

! Logical switches (TRUE/FALSE) to specify the adjoint state variables
! whose sensitivity is required.

Lstate(isFsur) == F ! free-surface
Lstate(isUbar) == F ! 2D U-momentum
Lstate(isVbar) == F ! 2D V-momentum
Lstate(isUvel) == F ! 3D U-momentum
Lstate(isVvel) == F ! 3D V-momentum
Lstate(isWvel) == F ! 3D W-momentum

Lstate(isTvar) == F F ! NT tracers

! Logical switches (TRUE/FALSE) to specify the state variables for
! which Forcing Singular Vectors or Stochastic Optimals is required.

Fstate(isFsur) == F ! free-surface
Fstate(isUbar) == F ! 2D U-momentum
Fstate(isVbar) == F ! 2D V-momentum
Fstate(isUvel) == F ! 3D U-momentum
Fstate(isVvel) == F ! 3D V-momentum
Fstate(isTvar) == F F ! NT tracers

Fstate(isUstr) == F ! surface U-stress
Fstate(isVstr) == F ! surface V-stress
Fstate(isTsur) == F F ! NT surface tracers flux

! Stochastic Optimals time decorrelation scale (days) assumed for
! red noise processes.

SO_decay == 2.0d0 ! days

! Stochastic Optimals surface forcing standard deviation for
! dimensionalization.

SO_sdev(isFsur) == 1.0d0 ! free-surface
SO_sdev(isUbar) == 1.0d0 ! 2D U-momentum
SO_sdev(isVbar) == 1.0d0 ! 2D V-momentum
SO_sdev(isUvel) == 1.0d0 ! 3D U-momentum
SO_sdev(isVvel) == 1.0d0 ! 3D V-momentum
SO_sdev(isTvar) == 1.0d0 1.0d0 ! NT tracers

SO_sdev(isUstr) == 1.0d0 ! surface U-stress
SO_sdev(isVstr) == 1.0d0 ! surface V-stress
SO_sdev(isTsur) == 1.0d0 1.0d0 ! NT surface tracers flux

! Logical switches (TRUE/FALSE) to activate writing of fields into
! HISTORY output file.

Hout(idUvel) == T ! u 3D U-velocity
Hout(idVvel) == T ! v 3D V-velocity
Hout(idu3dE) == F ! u_eastward 3D U-eastward at RHO-points
Hout(idv3dN) == F ! v_northward 3D V-northward at RHO-points
Hout(idWvel) == T ! w 3D W-velocity
Hout(idOvel) == T ! omega omega vertical velocity
Hout(idUbar) == T ! ubar 2D U-velocity
Hout(idVbar) == T ! vbar 2D V-velocity
Hout(idu2dE) == F ! ubar_eastward 2D U-eastward at RHO-points
Hout(idv2dN) == F ! vbar_northward 2D V-northward at RHO-points
Hout(idFsur) == T ! zeta free-surface
Hout(idBath) == T ! bath time-dependent bathymetry

Hout(idTvar) == F F ! temp, salt temperature and salinity

Hout(idpthR) == F ! z_rho time-varying depths of RHO-points
Hout(idpthU) == F ! z_u time-varying depths of U-points
Hout(idpthV) == F ! z_v time-varying depths of V-points
Hout(idpthW) == F ! z_w time-varying depths of W-points

Hout(idUsms) == F ! sustr surface U-stress
Hout(idVsms) == F ! svstr surface V-stress
Hout(idUbms) == F ! bustr bottom U-stress
Hout(idVbms) == F ! bvstr bottom V-stress

Hout(idUbrs) == F ! bustrc bottom U-current stress
Hout(idVbrs) == F ! bvstrc bottom V-current stress
Hout(idUbws) == F ! bustrw bottom U-wave stress
Hout(idVbws) == F ! bvstrw bottom V-wave stress
Hout(idUbcs) == F ! bustrcwmax bottom max wave-current U-stress
Hout(idVbcs) == F ! bvstrcwmax bottom max wave-current V-stress

Hout(idUbot) == F ! Ubot bed wave orbital U-velocity
Hout(idVbot) == F ! Vbot bed wave orbital V-velocity
Hout(idUbur) == F ! Ur bottom U-velocity above bed
Hout(idVbvr) == F ! Vr bottom V-velocity above bed

Hout(idW2xx) == F ! Sxx_bar 2D radiation stress, Sxx component
Hout(idW2xy) == F ! Sxy_bar 2D radiation stress, Sxy component
Hout(idW2yy) == F ! Syy_bar 2D radiation stress, Syy component
Hout(idU2rs) == F ! Ubar_Rstress 2D radiation U-stress
Hout(idV2rs) == F ! Vbar_Rstress 2D radiation V-stress
Hout(idU2Sd) == F ! ubar_stokes 2D U-Stokes velocity
Hout(idV2Sd) == F ! vbar_stokes 2D V-Stokes velocity

Hout(idW3xx) == F ! Sxx 3D radiation stress, Sxx component
Hout(idW3xy) == F ! Sxy 3D radiation stress, Sxy component
Hout(idW3yy) == F ! Syy 3D radiation stress, Syy component
Hout(idW3zx) == F ! Szx 3D radiation stress, Szx component
Hout(idW3zy) == F ! Szy 3D radiation stress, Szy component
Hout(idU3rs) == F ! u_Rstress 3D U-radiation stress
Hout(idV3rs) == F ! v_Rstress 3D V-radiation stress
Hout(idU3Sd) == F ! u_stokes 3D U-Stokes velocity
Hout(idV3Sd) == F ! v_stokes 3D V-Stokes velocity

Hout(idWamp) == F ! Hwave wave height
Hout(idWlen) == F ! Lwave wave length
Hout(idWdir) == F ! Dwave wave direction
Hout(idWptp) == F ! Pwave_top wave surface period
Hout(idWpbt) == F ! Pwave_bot wave bottom period
Hout(idWorb) == F ! Ub_swan wave bottom orbital velocity
Hout(idWdis) == F ! Wave_dissip wave dissipation

Hout(idPair) == F ! Pair surface air pressure
Hout(idTair) == F ! Tair surface air temperature
Hout(idUair) == F ! Uair surface U-wind component
Hout(idVair) == F ! Vair surface V-wind component

Hout(idTsur) == F F ! shflux, ssflux surface net heat and salt flux
Hout(idLhea) == F ! latent latent heat flux
Hout(idShea) == F ! sensible sensible heat flux
Hout(idLrad) == F ! lwrad longwave radiation flux
Hout(idSrad) == F ! swrad shortwave radiation flux
Hout(idEmPf) == F ! EminusP E-P flux
Hout(idevap) == F ! evaporation evaporation rate
Hout(idrain) == F ! rain precipitation rate

Hout(idDano) == F ! rho density anomaly
Hout(idVvis) == F ! AKv vertical viscosity
Hout(idTdif) == F ! AKt vertical T-diffusion
Hout(idSdif) == F ! AKs vertical Salinity diffusion
Hout(idHsbl) == F ! Hsbl depth of surface boundary layer
Hout(idHbbl) == F ! Hbbl depth of bottom boundary layer
Hout(idMtke) == F ! tke turbulent kinetic energy
Hout(idMtls) == F ! gls turbulent length scale

! Logical switches (TRUE/FALSE) to activate writing of extra inert passive
! tracers other than biological and sediment tracers into the HISTORY
! output file. An inert passive tracer is one that it is only advected and
! diffused. Other processes are ignored. These tracers include, for example,
! dyes, pollutants, oil spills, etc. NPT values are expected. However, these
! switches can be activated using compact parameter specification.

Hout(inert) == F ! dye_01, ... inert passive tracers


! Generic User parameters, [1:NUSER].

NUSER = 1
USER = 1.d0


! Input and Output files processing library to use:
!
! [1] Standard NetCDF-3 or NetCDF-4 library
! [2] Serial or Parallel I/O with Parallel-IO (PIO) library (MPI only)

INP_LIB = 1
OUT_LIB = 1

! PIO library methods for reading/writing NetCDF files:
!
! [0] parallel read and write of PnetCDF (CDF-5, not recommended)
! [1] parallel read and write of NetCDF3 (64-bit offset)
! [2] serial read and write of NetCDF3 (64-bit offset)
! [3] parallel read and serial write of NetCDF4/HDF5
! [4] parallel read and write of NETCDF4/HDF5

PIO_METHOD = 2

! PIO library MPI processes set-up:

PIO_IOTASKS = 2 ! number of I/O tasks to define
PIO_STRIDE = 1 ! stride in the MPI-ran between I/O tasks
PIO_BASE = 0 ! offset for the first I/O task
PIO_AGGREG = 1 ! number of MPI-aggregators to use

! PIO library rearranger methods for moving data between computational and I/O
! processes:
!
! [1] Box rearrangement
! [2] Subset rearrangement

PIO_REARR = 1

! PIO library rearranger flag for MPI communications between computational
! and I/O processes:
!
! [0] Point-to-Point (low-level communications)
! [1] Collective (high-level grouped communications)

PIO_REARRCOM = 0

! PIO library rearranger flow control direction flag for MPI communications
! between computational and I/O processes:
!
! [0] Enable computational to I/O processes, and vice versa
! [2] Enable computational to I/O processes only
! [3] Enable I/O to computational processes only
! [4] Disable flow control

PIO_REARRDIR = 2

! PIO rearranger options for computational to I/O processes (C2I):

PIO_C2I_HS = T ! Enable C2I handshake (T/F)
PIO_C2I_Send = T ! Enable C2I Isends (T/F)
PIO_C2I_Preq = 64 ! Maximum pending C2I requests

! PIO rearranger options for I/O to computational processes (I2C):

PIO_I2C_HS = T ! Enable I2C handshake (T/F)
PIO_I2C_Send = T ! Enable I2C Isends (T/F)
PIO_I2C_Preq = 65 ! Maximum pending I2C requests



! Input NetCDF file names, [1:Ngrids].

GRDNAME == tidebore_grd.nc
ININAME == tidebore_ini.nc
ITLNAME ==
IRPNAME ==
IADNAME ==
FWDNAME ==
ADSNAME ==

FOInameA ==
FOInameB ==

FCTnameA ==
FCTnameB ==

NGCNAME =


NBCFILES == 1 ! number of boundary files

BRYNAME == tidebore_bry.nc
CLMNAME ==


NUDNAME ==


SSFNAME ==


TIDENAME == tidebore_tideforcing.nc


NFFILES == 1 ! number of unique forcing files

FRCNAME == tidebore_frc.nc ! forcing file 1, grid 1


DAINAME == tidebore_dai.nc
GSTNAME == tidebore_gst.nc
RSTNAME ==
HISNAME == tidebore_his.nc
QCKNAME == hzec_qck.nc
TLMNAME == hzec_tlm.nc
TLFNAME == hzec_tlf.nc
ADJNAME == hzec_adj.nc
AVGNAME == hzec_avg.nc
HARNAME == hzec_har.nc
DIANAME == hzec_dia.nc
STANAME == hzec_sta.nc
FLTNAME == hzec_flt.nc


APARNAM = s4dvar.in
SPOSNAM = stations.in
FPOSNAM = floats.in
BPARNAM = bio_Fennel.in
SPARNAM = sediment.in
USRNAME = MyFile.dat

[smchen]

Try giving RSTNAME a filename.

[zhqianj]

Thanks， I named rst filen as tidebore_rst.nc . and it still does not work.

the results says:


0 0001-01-01 00:00:00.00 0.000000E+00 2.666835E+01 2.666835E+01 9.713888E+09
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00
DEF_HIS_NF90 - creating history file, Grid 01: tidebore_his.nc
WRT_HIS_NF90 - wrote history fields (Index=1,1) in record = 1
DEF_RST_NF90 - creating restart file, Grid 01: tidebore_rst.nc
Found Error: ** Line: 164 Source: ROMS/Utility/def_var.F, def_var_nf90

DEF_VAR_NF90 - Grid 01, unable to define variable:
in NetCDF file: tidebore_rst.nc
Found Error: 03 Line: 1050 Source: ROMS/Utility/def_rst.F, def_rst_nf90
Found Error: 03 Line: 85 Source: ROMS/Utility/def_rst.F
Found Error: 03 Line: 576 Source: ROMS/Nonlinear/output.F
Found Error: 03 Line: 518 Source: ROMS/Nonlinear/main3d.F
Found Error: 03 Line: 301 Source: ROMS/Drivers/nl_roms.h, ROMS_run

Elapsed wall CPU time for each process (seconds):

Node # 0 CPU: 1.469
Total: 5.984
Average: 1.496
Minimum: 1.469
Maximum: 1.562

[smchen]

I have no idea why your rst file has problems. Or just switch off the rst by setting NRST=0 to see if the model runs well to make sure the problems are rst files.

-- Shi-Ming