first step nan by using 'roms2roms'

[XMUZSM]

Hi, I have a parent restart file (theta_s=8 and theta_b=0.1). And I want to modified stretching, so I change parameters(theta_s=7, theta_b=2) by using code "roms2roms". After that, I spin up model with the modified version restart files.

First step like:

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 2012-07-01 21:00:00.00 5.550321E-03 2.136531E+04 2.136532E+04 2.709732E+16 0.00
(356,327,29) 4.685782E-02 1.507398E-02 0.000000E+00 2.364141E+00
DEF_HIS - creating history file, Grid 01: res/his.seddy2.s2b7_1w.nc
WRT_HIS - wrote history fields (Index=1,1) in record = 0000001
DEF_QUICK - creating quicksave file, Grid 01: res/qck.seddy2.s2b7_1w.nc
WRT_QUICK - wrote quicksave fields (Index=1,1) in record = 0000001
DEF_AVG - creating average file, Grid 01: res/avg.seddy2.s2b7_1w.nc
DEF_DIAGS - creating diagnostics file, Grid 01: res/dia.seddy2.s2b7_1w.nc
1 2012-07-01 21:03:20.00 NaN NaN NaN NaN 19.37
(000,000,00) 0.000000E+00 0.000000E+00 0.000000E+00 NaN
NaN error
NaN error
Found Error: 01 Line: 298 Source: ROMS/Nonlinear/main3d.F
Found Error: 01 Line: 298 Source: ROMS/Drivers/nl_ocean.h

Blowing-up: Saving latest model state into RESTART file

NaN error
NaN error
NaN error


the model can run if I change little (theta_s=8,theta_b=0.2), that's weird

[wilkin]

I suspect that in the vertical interpolation either the new surface or bottom level is outside the range of the old levels and extrapolation creates NANs. You might need to create a rule to catch this, or explicitly set the extrapolation method to 'nearest'

[XMUZSM]

I suspect that in the vertical interpolation either the new surface or bottom level is outside the range of the old levels and extrapolation creates NANs. You might need to create a rule to catch this, or explicitly set the extrapolation method to 'nearest'

Hi, it is weird that there is no NaN in the interpolated data. I switch from 'linear' to 'nearest '. but it didn't work.

[wilkin]

Not interpolation method but extrapolation method (assuming you are using Matlab interp … you don't say whose code you are using)