Model crashed when read river runoff suspended sediment

[xkhmyx]

I was trying to run the longbay (SC) model.
When sediment is undefined (only hydro), there was no problem for the whole year 2008.
When sediment is defined (a total of 6 sed classes, two from sea bed and four from rivers), the model crashed in the middle of the year 2008.

Details are below.
I suspect this is a issue in sediment part, but am not sure.
Checked the river input nc file, but did not find any problem.

Kevin Xu
kxu@coastal.edu

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
361434 14129.99583 5.448856E-03 9.651963E+01 9.652508E+01 1.994789E+11 0
361435 14129.99653 5.452989E-03 9.651949E+01 9.652494E+01 1.994786E+11 0
361436 14129.99722 5.457445E-03 9.651935E+01 9.652481E+01 1.994784E+11 0
361437 14129.99792 5.462220E-03 9.651921E+01 9.652467E+01 1.994781E+11 0
361438 14129.99861 5.467310E-03 9.651907E+01 9.652454E+01 1.994778E+11 0
361439 14129.99931 5.472707E-03 9.651894E+01 9.652441E+01 1.994775E+11 0
GET_2DFLD - surface u-wind component, t = 14130.0042
(File: longbay_atmo.nc, Rec=****, Index=2)
(Tmin= 13514.0000 Tmax= 14244.9917)
(Min = 2.11503053E+00 Max = 2.11503053E+00)
GET_2DFLD - surface v-wind component, t = 14130.0042
(File: longbay_atmo.nc, Rec=****, Index=2)
(Tmin= 13514.0000 Tmax= 14244.9917)
(Min = 4.75043631E+00 Max = 4.75043631E+00)
361440 14130.00000 5.478408E-03 9.651880E+01 9.652428E+01 1.994773E+11 0
WRT_HIS - wrote history fields (Index=1,1) into time record = 0000002
GET_NGFLD - river runoff mass transport, t = 14131.0000
(File: longbay_river.nc, Rec=0618, Index=2)
(Tmin= 13514.0000 Tmax= 14598.0000)
(Min = -3.00158500E+02 Max = -1.35071500E+02)
GET_NGFLD - river runoff potential temperature, t = 14131.0000
(File: longbay_river.nc, Rec=0618, Index=2)
(Tmin= 13514.0000 Tmax= 14598.0000)
(Min = 1.00000000E+01 Max = 1.00000000E+01)
GET_NGFLD - river runoff salinity, t = 14131.0000
(File: longbay_river.nc, Rec=0618, Index=2)
(Tmin= 13514.0000 Tmax= 14598.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - river runoff suspended sediment concentrationt = 14131.0000
(File: longbay_river.nc, Rec=0618, Index=2)
(Tmin= 13514.0000 Tmax= 14598.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - river runoff suspended sediment concentrationt = 14131.0000
(File: longbay_river.nc, Rec=0618, Index=2)
(Tmin= 13514.0000 Tmax= 14598.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_NGFLD - river runoff suspended sediment concentrationt = 14131.0000
(File: longbay_river.nc, Rec=0618, Index=2)
(Tmin= 13514.0000 Tmax= 14598.0000)
(Min = 0.00000000E+00 Max = 6.20000000E-02)
GET_NGFLD - river runoff suspended sediment concentration, size class 04At line 221 of file diag.f90
Fortran runtime error: End of file
rank 0 in job 1 longbay.coastal.edu_43857 caused collective abort of all ranks
exit status of rank 0: killed by signal 9

[jcwarner]

GET_NGFLD - river runoff suspended sediment concentration, size class 04At line 221 of file diag.f90
Fortran runtime error: End of file
rank 0 in job 1 longbay.coastal.edu_43857 caused collective abort of all ranks
exit status of rank 0: killed by signal 9


1) What is on line 221 of file diag.f90?

2) "Fortran runtime error: End of file"
Are there enough time steps in the file?

[xkhmyx]

1) What is on line 221 of file diag.f90?
It is IF (Master) WRITE (stdout,20) iic(ng)-1, tdays(ng), avgke, &
& avgpe, avgkp, volume, trd


More info is:
IF (Master) WRITE (stdout,10) 'STEP', 'time[DAYS]', &
& 'KINETIC_ENRG', 'POTEN_ENRG', &
& 'TOTAL_ENRG', 'NET_VOLUME', &
& 'trd'
10 FORMAT (/,3x,a,2x,a,2x,a,4x,a,4x,a,3x,a,2x,a,/)
END IF
IF (Master) WRITE (stdout,20) iic(ng)-1, tdays(ng), avgke, &
& avgpe, avgkp, volume, trd
20 FORMAT (i7,f12.5,4(1pe14.6),i4)


2) "Fortran runtime error: End of file"
Are there enough time steps in the file?

Model time step is 60s, and the river runoff is daily.
Wind is every 20 min.

[arango]

It crashed because of temperature and salinity. Nothing to do with the sediment. If you are using a recent version of ROMS, you need to be sure that you activate input parameter LtracerSrc:

Code: Select all

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

  LtracerSrc == T T                        ! temperature, salinity, inert

Check release notes for the latest version or the following trac ticket.

[xkhmyx]

I am using ROMS 3.0
svn: $Id: Version 110 2007-09-20 04:31:40Z arango $

All Rivers' temperature is 10 degree C, and salinity is 0.
If it is a problem due to temperature and salinity, why did the model run the whole year (for hydrodynamics only) without any problem?