Model crashed when read river runoff suspended sediment

Sediment modeling collaborators: issues, applications, information exchange

Moderators: arango, robertson, rsignell

Post Reply
Message
Author
xkhmyx
Posts: 8
Joined: Fri May 11, 2007 7:59 pm
Location: Louisiana State University

Model crashed when read river runoff suspended sediment

#1 Unread post by 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
Posts: 1172
Joined: Wed Dec 31, 2003 6:16 pm
Location: USGS, USA

Re: Model crashed when read river runoff suspended sediment

#2 Unread post by 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
Posts: 8
Joined: Fri May 11, 2007 7:59 pm
Location: Louisiana State University

Re: Model crashed when read river runoff suspended sediment

#3 Unread post by 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.

User avatar
arango
Site Admin
Posts: 1347
Joined: Wed Feb 26, 2003 4:41 pm
Location: DMCS, Rutgers University
Contact:

Re: Model crashed when read river runoff suspended sediment

#4 Unread post by 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 :arrow: release notes for the latest version or the following :arrow: trac ticket.

xkhmyx
Posts: 8
Joined: Fri May 11, 2007 7:59 pm
Location: Louisiana State University

Re: Model crashed when read river runoff suspended sediment

#5 Unread post by 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?

Post Reply