ROMS
Loading...
Searching...
No Matches
sediment_def.h
Go to the documentation of this file.
1/*
2** git $id$
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************************************************************************
8** **
9** defines sediment model input parameters in output netcdf files. **
10** it is included in routine "def_info.F". **
11** **
12************************************************************************
13*/
14
15!
16! Define sediment model parameters.
17!
18 vinfo( 1)='minlayer_thick'
19 vinfo( 2)='depositional bed layer minimum thickness'
20 vinfo( 3)='meter'
21 status=def_var(ng, model, ncid, varid, nf_type, &
22 & 1, (/0/), aval, vinfo, ncname, &
23 & setparaccess = .false.)
24 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
25
26 vinfo( 1)='newlayer_thick'
27 vinfo( 2)='depositional bed thickness criteria to crate a '// &
28 & 'new layer'
29 vinfo( 3)='meter'
30 status=def_var(ng, model, ncid, varid, nf_type, &
31 & 1, (/0/), aval, vinfo, ncname, &
32 & setparaccess = .false.)
33 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
34
35#ifdef BEDLOAD
36 vinfo( 1)='bedload_coeff'
37 vinfo( 2)='bedload transport rate coefficient'
38 status=def_var(ng, model, ncid, varid, nf_type, &
39 & 1, (/0/), aval, vinfo, ncname, &
40 & setparaccess = .false.)
41 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
42#endif
43
44#ifdef ANA_SEDIMENT
45 vinfo( 1)='Sd50'
46 vinfo( 2)='median sediment grain diameter used in '// &
47 & 'uniform initial conditions'
48 vinfo( 3)='meter'
49 status=def_var(ng, model, ncid, varid, nf_type, &
50 & 1, (/seddim/), aval, vinfo, ncname, &
51 & setparaccess = .false.)
52 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
53
54 vinfo( 1)='Srho'
55 vinfo( 2)='sediment grain density used in '// &
56 & 'uniform initial conditions'
57 vinfo( 3)='kilogram meter-3'
58 status=def_var(ng, model, ncid, varid, nf_type, &
59 & 1, (/seddim/), aval, vinfo, ncname, &
60 & setparaccess = .false.)
61 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
62
63 vinfo( 1)='Csed'
64 vinfo( 2)='sediment concentration used in '// &
65 & 'uniform initial conditions'
66 vinfo( 3)='kilogram meter-3'
67 status=def_var(ng, model, ncid, varid, nf_type, &
68 & 1, (/seddim/), aval, vinfo, ncname, &
69 & setparaccess = .false.)
70 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
71#endif
72
73 vinfo( 1)='Wsed'
74 vinfo( 2)='sediment particle settling velocity'
75 vinfo( 3)='meter second-1'
76 status=def_var(ng, model, ncid, varid, nf_type, &
77 & 1, (/seddim/), aval, vinfo, ncname, &
78 & setparaccess = .false.)
79 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
80
81 vinfo( 1)='Erate'
82 vinfo( 2)='sediment surface layer erosing rate'
83 vinfo( 3)='kilogram meter-2 second-1'
84 status=def_var(ng, model, ncid, varid, nf_type, &
85 & 1, (/seddim/), aval, vinfo, ncname, &
86 & setparaccess = .false.)
87 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
88
89 vinfo( 1)='tau_ce'
90 vinfo( 2)='sediment critical shear for erosion'
91 vinfo( 3)='meter-2 second-2'
92 status=def_var(ng, model, ncid, varid, nf_type, &
93 & 1, (/seddim/), aval, vinfo, ncname, &
94 & setparaccess = .false.)
95 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
96
97 vinfo( 1)='tau_cd'
98 vinfo( 2)='sediment critical shear for deposition'
99 vinfo( 3)='meter-2 second-2'
100 status=def_var(ng, model, ncid, varid, nf_type, &
101 & 1, (/seddim/), aval, vinfo, ncname, &
102 & setparaccess = .false.)
103 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
104
105 vinfo( 1)='poros'
106 vinfo( 2)='sediment porosity'
107 status=def_var(ng, model, ncid, varid, nf_type, &
108 & 1, (/seddim/), aval, vinfo, ncname, &
109 & setparaccess = .false.)
110 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
subroutine output(ng)
Definition output.F:4