Thank you so much, Kate.
It seems that we can set both the DSTART and Time_ref to 0 because I ran a model with them set to 0 succesfully.
My own model is just based on that model. It seems there is no report shown it is a velocity or density problem.
Yes, the NINFO is 1. I have three output files: ***_rst.nc, ***_his.nc and ***_avg.nc. The restart file is the ***_rst.nc file, right? what do you mean by the 'last record'? last record of the value of the parameter?
The last records of the log file is:
Code: Select all
2221 1 03:45:45 1.763320E-02 8.136258E+01 8.138021E+01 3.613824E+11
2222 1 03:46:30 1.770285E-02 8.134668E+01 8.136438E+01 3.613458E+11
2223 1 03:47:15 1.777491E-02 8.133072E+01 8.134850E+01 3.613089E+11
2224 1 03:48:00 1.784953E-02 8.131470E+01 8.133255E+01 3.612717E+11
2225 1 03:48:45 1.792612E-02 8.129861E+01 8.131654E+01 3.612343E+11
2226 1 03:49:30 1.800896E-02 8.128246E+01 8.130047E+01 3.611966E+11
2227 1 03:50:15 1.808512E-02 8.126624E+01 8.128433E+01 3.611587E+11
Blowing-up: Saving latest model state into RESTART file
WRT_RST - wrote re-start fields (Index=1,2) into time record = 0000001
I have tried a smaller timestep with 30 s (45 s before), but after 2 days, also it blew up.
I attached the .h file and .in file as follows, any more suggestion? looking forward!
Code: Select all
/* Basic physics options */
#define UV_COR
#define UV_ADV
#define UV_VIS2
#define TS_DIF2
#define SOLVE3D
#define SALINITY
#define MIX_GEO_UV /* mixing on geopotential (constant Z) surfaces */
#define MIX_GEO_TS /* mixing on geopotential (constant Z) surfaces */
#define NONLIN_EOS
/* Basic numerics options */
/* #define INLINE_2DIO */
#define CURVGRID
#define MASKING
#define SPLINES
#define TS_MPDATA
#undef TS_U3HADVECTION /* define if 3rd-order upstream horiz. advection */
#undef TS_SVADVECTION /* define if splines vertical advection */
#define UV_SADVECTION /* turn ON or OFF splines vertical advection */
/* average options */
#define AVERAGES
#define AVERAGES_AKV
#define AVERAGES_AKT
#define AVERAGES_AKS
#define AVERAGES_FLUXES
/* Basic Diagnostic options */
#undef DIAGNOSTICS_TS /*compute and write diagnostic terms for tracer equations balances*/
#ifdef DIAGNOSTICS_TS
# define DIAGNOSTICS_TS_split/*compute and write diagnostic terms in XI- and ETA- direction*/
# define DIAGNOSTICS_StoZcor /*Compute and write correction terms to horizontal flux diagnositcs*/
# define DIAGNOSTICS_salt /*Compute and write salinity diagnositc terms*/
# undef DIAGNOSTICS_UV /*Compute and write momentum diagnostic terms*/
#endif
/* Outputs */
#define NO_WRITE_GRID /* define if not writing grid arrays */
#undef WRITE_WATER /* use if only writing water points data*/
#define INLINE_2DIO
/* Surface and bottom boundary conditions */
#define ANA_BSFLUX /* analytical bottom salinity flux */
#define ANA_BTFLUX /* analytical bottom temperature flux */
#define ANA_SSFLUX /* analytical surface salinity flux */
#define ANA_STFLUX /* analytical surface temperature flux */
#define UV_LOGDRAG /* turn ON or OFF logarithmic bottom friction */
#undef QCORRECTION /* define if net heat flux correction */
#undef SST_RELAXATION
#define BULK_FLUXES
#ifdef BULK_FLUXES
# define LONGWAVE_OUT
# define ANA_CLOUD
# define ANA_RAIN
#else
# define ANA_SMFLUX
# define ANA_STFLUX
#endif
#define ATM_PRESS
#define EMINUSP
/* Vertical subgridscale turbulence closure */
#ifdef SOLVE3D /*!!!! Masked by Yun Li, Jan.11, 2008!!!!*/
# define GLS_MIXING /* Activate Generic Length-Scale mixing */
# undef LMD_MIXING
# ifdef MY25_MIXING
# undef N2S2_HORAVG /* Activate horizontal smoothing of buoyancy/shear */
# undef KANTHA_CLAYSON
# endif
# ifdef LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_BKPP
# define LMD_SKPP
# undef LMD_NONLOCAL
# endif
# ifdef GLS_MIXING
# define KANTHA_CLAYSON /* Kantha and Clayson stability function formulation */
# undef N2S2_HORAVG /* Activate horizontal smoothing of buoyancy/shear */
# undef CANUTO_A /* Canuto A-stability function formulation */
# undef CANUTO_B /* Canuto B-stability function formulation */
# endif
#endif
/* Open boundary conditions */
#define EAST_FSCHAPMAN
#define EAST_M2FLATHER
#define EAST_M3RADIATION
#define EAST_TRADIATION
#define EAST_TNUDGING
#define WESTERN_WALL
#define SOUTHERN_WALL
#define NORTHERN_WALL
#define OBC_DATA
#ifdef OBC_DATA
# define EAST_FSOBC
# define EAST_M2OBC
# define EAST_TOBC
#endif
#define TS_PSOURCE
#define UV_PSOURCE
#ifdef UV_PSOURCE
# define PSOURCE_FSCHAPMAN
#endif
/* For Restart */
#define PERFECT_RESTART
#ifdef PERFECT_RESTART
# undef ANA_INITIAL
#endif
/* For sediment */
#define SEDIMENT
#ifdef SEDIMENT
# define SUSPLOAD
# define RIVER_SEDIMENT
# undef SED_DENS
# define ANA_SPFLUX
# define ANA_BPFLUX
#endif
Code: Select all
! Grid dimension parameters. See notes below in the Glossary for how to set
! these parameters correctly.
Lm == 160 ! Number of I-direction INTERIOR RHO-points
Mm == 240 ! Number of J-direction INTERIOR RHO-points
N == 20 ! Number of vertical levels
Nbed = 1 ! 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 = 1 ! 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 == 4 ! I-direction partition
NtileJ == 6 ! J-direction partition
! Time-Stepping parameters.
NTIMES == 19200 ! 1985280 ! 1034(2006-1-1 to 2008-10-31)x1920x1 baroclinic steps/day
DT == 45.0d0
NDTFAST == 30
! 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 ! no Restart
LcycleRST == T ! the latest two restart saved every NRST times
NRST == 19200 ! 10 day x 1920 steps/day rst
NSTA == 1
NFLT == 1
NINFO == 1 ! display on the screen
! Output history, average, diagnostic files parameters.
LDEFOUT == T
NHIS == 320 ! 4 hours
NDEFHIS == 0!57600 ! 30 days
NTSAVG == 1
NAVG == 1920 ! 24 hr
NDEFAVG == 0!57600 ! 30 days
NTSDIA == 1
NDIA == 1920 ! 24 hours
NDEFDIA == 0!57600 ! 10 days
! Output tangent linear and adjoint models parameters.
LcycleTLM == F
NTLM == 72
NDEFTLM == 0
LcycleADJ == F
NADJ == 72
NDEFADJ == 0
NSFF == 72 ! new variables from ROMS 333
NOBC == 72 ! new variables from ROMS 333