3*************************************************** hernan g. arango ***
4** copyright(c) 2002-2025 the roms group **
5** licensed under a mit/x style license **
6** see license_roms.md **
7************************************************************************
9** defines sediment model input parameters in
output netcdf files. **
10** it is included in routine
"def_info.F". **
12************************************************************************
18 vinfo( 1)=
'minlayer_thick'
19 vinfo( 2)=
'depositional bed layer minimum thickness'
21 status=def_var(ng, model, piofile, piovar, pio_type, &
22 & 1, (/0/), aval, vinfo, ncname, &
23 & setparaccess = .false.)
24 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
26 vinfo( 1)=
'newlayer_thick'
27 vinfo( 2)=
'depositional bed thickness criteria to crate a '// &
30 status=def_var(ng, model, piofile, piovar, pio_type, &
31 & 1, (/0/), aval, vinfo, ncname, &
32 & setparaccess = .false.)
33 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
36 vinfo( 1)=
'bedload_coeff'
37 vinfo( 2)=
'bedload transport rate coefficient'
38 status=def_var(ng, model, piofile, piovar, pio_type, &
39 & 1, (/0/), aval, vinfo, ncname, &
40 & setparaccess = .false.)
41 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
46 vinfo( 2)=
'median sediment grain diameter used in '// &
47 &
'uniform initial conditions'
49 status=def_var(ng, model, piofile, piovar, pio_type, &
50 & 1, (/seddim/), aval, vinfo, ncname, &
51 & setparaccess = .false.)
52 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
55 vinfo( 2)=
'sediment grain density used in '// &
56 &
'uniform initial conditions'
57 vinfo( 3)=
'kilogram meter-3'
58 status=def_var(ng, model, piofile, piovar, pio_type, &
59 & 1, (/seddim/), aval, vinfo, ncname, &
60 & setparaccess = .false.)
61 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
64 vinfo( 2)=
'sediment concentration used in '// &
65 &
'uniform initial conditions'
66 vinfo( 3)=
'kilogram meter-3'
67 status=def_var(ng, model, piofile, piovar, pio_type, &
68 & 1, (/seddim/), aval, vinfo, ncname, &
69 & setparaccess = .false.)
70 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
74 vinfo( 2)=
'sediment particle settling velocity'
75 vinfo( 3)=
'meter second-1'
76 status=def_var(ng, model, piofile, piovar, pio_type, &
77 & 1, (/seddim/), aval, vinfo, ncname, &
78 & setparaccess = .false.)
79 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
82 vinfo( 2)=
'sediment surface layer erosing rate'
83 vinfo( 3)=
'kilogram meter-2 second-1'
84 status=def_var(ng, model, piofile, piovar, pio_type, &
85 & 1, (/seddim/), aval, vinfo, ncname, &
86 & setparaccess = .false.)
87 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
90 vinfo( 2)=
'sediment critical shear for erosion'
91 vinfo( 3)=
'meter-2 second-2'
92 status=def_var(ng, model, piofile, piovar, pio_type, &
93 & 1, (/seddim/), aval, vinfo, ncname, &
94 & setparaccess = .false.)
95 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
98 vinfo( 2)=
'sediment critical shear for deposition'
99 vinfo( 3)=
'meter-2 second-2'
100 status=def_var(ng, model, piofile, piovar, pio_type, &
101 & 1, (/seddim/), aval, vinfo, ncname, &
102 & setparaccess = .false.)
103 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN
106 vinfo( 2)=
'sediment porosity'
107 status=def_var(ng, model, piofile, piovar, pio_type, &
108 & 1, (/seddim/), aval, vinfo, ncname, &
109 & setparaccess = .false.)
110 IF (founderror(exit_flag, noerror, __line__, myfile))
RETURN