I am running a 14 year simulation over the Aegean sea but for some reason,

which I am struggling to understand, the model blows up after 4 years. To

be honest I am surprised because the current configuration of the model

appeared to be ok. I have run the same region for 6 years with different

atmospheric forcing (6hr NCEP) and the model run smoothly without any problem.

At the moment I am using 3hr ECMWF era-interim data. The issue that i can't

understand is that the last record of restart file has only NaNs. From

log's file last records it doesn't seem anything unusual.

Code: Select all

```
...
175331 1095 19:39:00 2.216486E-03 2.547863E+03 2.547865E+03 5.064988E+13
(133,142,30) 3.073807E-02 0.000000E+00 2.351532E+01 6.851973E-01
175332 1095 19:48:00 2.218229E-03 2.547863E+03 2.547865E+03 5.064989E+13
(133,142,30) 3.073789E-02 0.000000E+00 2.351516E+01 6.855414E-01
175333 1095 19:57:00 2.219928E-03 2.547863E+03 2.547865E+03 5.064989E+13
(133,142,30) 3.073771E-02 0.000000E+00 2.351500E+01 6.873933E-01
175334 1095 20:06:00 2.221593E-03 2.547863E+03 2.547865E+03 5.064990E+13
(133,142,30) 3.073753E-02 0.000000E+00 2.351484E+01 6.891778E-01
175335 1095 20:15:00 2.223204E-03 2.547863E+03 2.547865E+03 5.064991E+13
(133,142,30) 3.073734E-02 0.000000E+00 2.351466E+01 6.909007E-01
175336 1095 20:24:00 2.224774E-03 2.547863E+03 2.547865E+03 5.064991E+13
(133,142,30) 3.073714E-02 0.000000E+00 2.351449E+01 6.925654E-01
175337 1095 20:33:00 2.226283E-03 2.547863E+03 2.547865E+03 5.064992E+13
(133,142,30) 3.073694E-02 0.000000E+00 2.351431E+01 6.941484E-01
175338 1095 20:42:00 2.227745E-03 2.547863E+03 2.547865E+03 5.064993E+13
(133,142,30) 3.073672E-02 0.000000E+00 2.351412E+01 6.956356E-01
175339 1095 20:51:00 2.229138E-03 2.547863E+03 2.547865E+03 5.064994E+13
(133,142,30) 3.073651E-02 0.000000E+00 2.351392E+01 6.970196E-01
175340 1095 21:00:00 2.230479E-03 2.547863E+03 2.547865E+03 5.064995E+13
(133,142,30) 3.073630E-02 0.000000E+00 2.351371E+01 6.983050E-01
GET_2DFLD - surface u-wind component, t = 1096 00:00:00
(Rec=0002920, Index=1, File: aeg_wind_era_2002.nc)
(Tmin= 731.1250 Tmax= 1096.0000)
(Min = -2.31089458E+00 Max = 9.14213168E+00)
GET_2DFLD - surface v-wind component, t = 1096 00:00:00
(Rec=0002920, Index=1, File: aeg_wind_era_2002.nc)
(Tmin= 731.1250 Tmax= 1096.0000)
(Min = -1.24266469E+34 Max = 5.11613192E+33)
GET_2DFLD - surface air pressure, t = 1096 00:00:00
(Rec=0002920, Index=1, File: aeg_Pair_era_2002.nc)
(Tmin= 731.1250 Tmax= 1096.0000)
(Min = 1.00522447E+03 Max = 1.01613286E+03)
GET_2DFLD - cloud fraction, t = 1096 00:00:00
(Rec=0002920, Index=1, File: aeg_cloud_era_2002.nc)
(Tmin= 731.1250 Tmax= 1096.0000)
(Min = 1.54924877E-01 Max = 9.99684666E-01)
GET_2DFLD - surface air temperature, t = 1096 00:00:00
(Rec=0002920, Index=1, File: aeg_Tair_era_2002.nc)
(Tmin= 731.1250 Tmax= 1096.0000)
(Min = 3.44416372E+00 Max = 1.65264554E+01)
GET_2DFLD - surface air relative humidity, t = 1096 00:00:00
(Rec=0002920, Index=1, File: aeg_Qair_era_2002.nc)
(Tmin= 731.1250 Tmax= 1096.0000)
(Min = 7.52880407E-01 Max = 9.95894587E-01)
175341 1095 21:09:00 2.231745E-03 2.547863E+03 2.547865E+03 5.064995E+13
(133,142,30) 3.073609E-02 0.000000E+00 2.351352E+01 6.995120E-01
175342 1095 21:18:00 NaN NaN NaN NaN
(051,062,01) NaN NaN NaN 0.000000E+00
Blowing-up: Saving latest model state into RESTART file
WRT_RST - wrote re-start fields (Index=1,1) into time record = 0000003
```

For bathymetry smoothing i used Mathieu Dutour Sikiric LP_bathymetry package

Other values concerning Barotropic Courant number etc are

Code: Select all

```
Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):
1.000000000000 1.033396903681 0.516698451840 1.000000000000 1.000000000000
Power filter parameters, Fgamma, gamma = 0.28400 0.22720
Metrics information for Grid 01:
===============================
Minimum X-grid spacing, DXmin = 2.38738876E+00 km
Maximum X-grid spacing, DXmax = 2.54891064E+00 km
Minimum Y-grid spacing, DYmin = 2.68266913E+00 km
Maximum Y-grid spacing, DYmax = 2.86416860E+00 km
Minimum Z-grid spacing, DZmin = 6.60088890E-02 m
Maximum Z-grid spacing, DZmax = 1.47386036E+02 m
Minimum barotropic Courant Number = 5.94298262E-02
Maximum barotropic Courant Number = 6.98174811E-01
Maximum Coriolis Courant Number = 5.16545468E-02
...
Basin information for Grid 01:
Maximum grid stiffness ratios: rx0 = 6.000344E-02 (Beckmann and Haidvogel)
rx1 = 3.279927E+00 (Haney)
Initial basin volumes: TotVolume = 4.7390984047E+13 m3
MinVolume = 2.3233030571E+06 m3
MaxVolume = 8.6057460458E+08 m3
Max/Min = 3.7040996523E+02
```

Code: Select all

```
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Closed Chapman Exp Closed Closed
ubar 1 Closed Shchepetkin Closed Closed
vbar 1 Closed Shchepetkin Closed Closed
u 1 Closed Clamped Closed Closed
v 1 Closed Clamped Closed Closed
temp 1 Closed Clamped Closed Closed
salt 1 Closed Clamped Closed Closed
tke 1 Closed Gradient Closed Closed
Activated C-preprocessing Options:
AEGEAN Aegean Sea 14 years experiment
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ASSUMED_SHAPE Using assumed-shape arrays.
AVERAGES Writing out time-averaged nonlinear model fields.
BULK_FLUXES Surface bulk fluxes parameterization.
CURVGRID Orthogonal curvilinear grid.
DIFF_3DCOEF Horizontal, time-dependent 3D diffusion coefficient.
DJ_GRADPS Parabolic Splines density Jacobian (Shchepetkin, 2002).
DOUBLE_PRECISION Double precision arithmetic.
EMINUSP Compute Salt Flux using E-P.
FLOATS Simulated Lagrangian drifters.
KANTHA_CLAYSON Kantha and Clayson stability function formulation.
LONGWAVE_OUT Compute outgoing longwave radiation internally.
MASKING Land/Sea masking.
MIX_ISO_TS Mixing of tracers along isopycnal surfaces.
MIX_GEO_UV Mixing of momentum along geopotential surfaces.
MY25_MIXING Mellor/Yamada Level-2.5 mixing closure.
NONLINEAR Nonlinear Model.
NONLIN_EOS Nonlinear Equation of State for seawater.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
_OPENMP OpenMP parallel shared-memory directives.
POWER_LAW Power-law shape time-averaging barotropic filter.
PROFILE Time profiling activated .
QCORRECTION Surface net heat flux correction.
K_C4ADVECTION Fourth-order centered differences advection of TKE fields.
RADIATION_2D Use tangential phase speed in radiation conditions.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SPLINES Conservative parabolic spline reconstruction.
STATIONS Writing out station data.
TS_U3HADVECTION Third-order upstream horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
TS_DIF2 Harmonic mixing of tracers.
TS_SMAGORINSKY Smagorinksy-like time-dependent diffusion coefficients.
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_VIS2 Harmonic mixing of momentum.
UV_SMAGORINSKY Smagorinksy-like time-dependent viscosity coefficients.
VAR_RHO_2D Variable density barotropic mode.
VISC_3DCOEF Horizontal, time-dependent 3D viscosity coefficient.
```

of Centers of gravity and integrals. Is it possible that these parameters create the problem. For the 6 year run these values where the same and the model run flawlessly. I am also attaching the output from the records

of restart file. Any help will be highly appreciated.

Giannis