-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathstations.in
353 lines (332 loc) · 21.8 KB
/
stations.in
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
! Station Extraction Parameters.
!
!svn $Id: stations.in 804 2016-06-06 23:02:36Z arango $
!========================================================= Hernan G. Arango ===
! Copyright (c) 2002-2016 The ROMS/TOMS Group !
! Licensed under a MIT/X style license !
! See License_ROMS.txt !
!==============================================================================
! !
! Input parameters can be entered in ANY order, provided that the parameter !
! KEYWORD (usually, upper case) is typed correctly followed by "=" or "==" !
! symbols. Any comment lines are allowed and must begin with an exclamation !
! mark (!) in column one. Comments may appear to the right of a parameter !
! specification to improve documentation. Comments will be ignored during !
! reading. Blank lines are also allowed and ignored. Continuation lines in !
! a parameter specification are allowed and must be preceded by a backslash !
! (\). In some instances, more than one value is required for a parameter. !
! If fewer values are provided, the last value is assigned for the entire !
! parameter array. The multiplication symbol (*), without blank spaces in !
! between, is allowed for a parameter specification. For example, in a two !
! grids nested application: !
! !
! AKT_BAK == 2*1.0d-6 2*5.0d-6 ! m2/s !
! !
! indicates that the first two entries of array AKT_BAK, in fortran column- !
! major order, will have the same value of "1.0d-6" for grid 1, whereas the !
! next two entries will have the same value of "5.0d-6" for grid 2. !
! !
! In multiple levels of nesting and/or multiple connected domains step-ups, !
! "Ngrids" entries are expected for some of these parameters. In such case, !
! the order of the entries for a parameter is extremely important. It must !
! follow the same order (1:Ngrids) as in the state variable declaration. The !
! USER may follow the above guidelines for specifying his/her values. These !
! parameters are marked by "==" plural symbol after the KEYWORD. !
! !
!==============================================================================
!
! Switch to control the writing of stations data within nested and/or multiple
! connected grids, [1:Ngrids].
Lstations == T
! Logical switches (TRUE/FALSE) to activate writing of fields in STATION
! output file, [Sout(:,ng), ng=1, Ngrids].
Sout(idUvel) == T ! u 3D U-velocity
Sout(idVvel) == T ! v 3D V-velocity
Sout(idu3dE) == T ! u_eastward 3D U-eastward at RHO-points
Sout(idv3dN) == T ! v_northward 3D V-nortward at RHO-points
Sout(idWvel) == T ! w 3D W-velocity
Sout(idOvel) == T ! omega 3D omega vertical velocity
Sout(idUbar) == F ! ubar 2D U-velocity
Sout(idVbar) == F ! vbar 2D V-velocity
Sout(idu2dE) == F ! ubar_eastward 2D U-eastward at RHO-points
Sout(idv2dN) == F ! vbar_northward 2D V-northward at RHO-points
Sout(idFsur) == T ! zeta free-surface
Sout(idBath) == F ! bath time-dependent bathymetry
Sout(idTvar) == T T ! temp, salt, ... all (NT) tracers
Sout(idUsms) == T ! sustr surface U-stress
Sout(idVsms) == T ! svstr surface V-stress
Sout(idUbms) == T ! bustr bottom U-stress
Sout(idVbms) == T ! bvstr bottom V-stress
Sout(idUbrs) == F ! bustrc bottom U-current stress
Sout(idVbrs) == F ! bvstrc bottom V-current stress
Sout(idUbws) == F ! bustrw bottom U-wave stress
Sout(idVbws) == F ! bvstrw bottom V-wave stress
Sout(idUbcs) == F ! bustrcwmax bottom max wave-current U-stress
Sout(idVbcs) == F ! bvstrcwmax bottom max wave-current V-stress
Sout(idUbot) == F ! Ubot bed wave orbital U-velocity
Sout(idVbot) == F ! Vbot bed wave orbital V-velocity
Sout(idUbur) == F ! Ur bottom U-velocity above bed
Sout(idVbvr) == F ! Vr bottom V-velocity above bed
Sout(idW2xx) == F ! Sxx_bar 2D radiation stress, Sxx component
Sout(idW2xy) == F ! Sxy_bar 2D radiation stress, Sxy component
Sout(idW2yy) == F ! Syy_bar 2D radiation stress, Syy component
Sout(idU2rs) == F ! Ubar_Rstress 2D radiation U-stress
Sout(idV2rs) == F ! Vbar_Rstress 2D radiation V-stress
Sout(idU2Sd) == F ! ubar_stokes 2D U-Stokes velocity
Sout(idV2Sd) == F ! vbar_stokes 2D V-Stokes velocity
Sout(idW3xx) == F ! Sxx 3D radiation stress, Sxx component
Sout(idW3xy) == F ! Sxy 3D radiation stress, Sxy component
Sout(idW3yy) == F ! Syy 3D radiation stress, Syy component
Sout(idW3zx) == F ! Szx 3D radiation stress, Szx component
Sout(idW3zy) == F ! Szy 3D radiation stress, Szy component
Sout(idU3rs) == F ! u_Rstress 3D U-radiation stress
Sout(idV3rs) == F ! v_Rstress 3D V-radiation stress
Sout(idU3Sd) == F ! u_stokes 3D U-Stokes velocity
Sout(idV3Sd) == F ! v_stokes 3D V-Stokes velocity
Sout(idWamp) == F ! Hwave wave height
Sout(idWlen) == F ! Lwave wave length
Sout(idWdir) == F ! Dwave wave direction
Sout(idWptp) == F ! Pwave_top wave surface period
Sout(idWpbt) == F ! Pwave_bot wave bottom period
Sout(idWorb) == F ! Ub_swan wave bottom orbital velocity
Sout(idWdis) == F ! Wave_dissip wave dissipation
Sout(idPair) == T ! Pair surface air pressure
Sout(idUair) == T ! Uair surface U-wind component
Sout(idVair) == T ! Vair surface V-wind component
Sout(idTsur) == T T ! shflux, ssflux surface net heat and salt flux
Sout(idLhea) == T ! latent latent heat flux
Sout(idShea) == T ! sensible sensible heat flux
Sout(idLrad) == T ! lwrad longwave radiation flux
Sout(idSrad) == T ! swrad shortwave radiation flux
Sout(idEmPf) == T ! EminusP E-P flux
Sout(idevap) == T ! evaporation evaporation rate
Sout(idrain) == T ! rain precipitation rate
Sout(idDano) == T ! rho density anomaly
Sout(idVvis) == T ! AKv vertical viscosity
Sout(idTdif) == T ! AKt vertical T-diffusion
Sout(idSdif) == T ! AKs vertical Salinity diffusion
Sout(idHsbl) == F ! Hsbl depth of surface boundary layer
Sout(idHbbl) == F ! Hbbl depth of bottom boundary layer
Sout(idMtke) == T ! tke turbulent kinetic energy
Sout(idMtls) == T ! gls turbulent length scale
! Logical switches (TRUE/FALSE) to activate writing of exposed sediment
! layer properties into STATIONS output file.
Sout(isd50) == F ! grain_diameter mean grain diameter
Sout(idens) == F ! grain_density mean grain density
Sout(iwsed) == F ! settling_vel mean settling velocity
Sout(itauc) == F ! erosion_stress critical erosion stress
Sout(irlen) == F ! ripple_length ripple length
Sout(irhgt) == F ! ripple_height ripple height
Sout(ibwav) == F ! bed_wave_amp wave excursion amplitude
Sout(izdef) == F ! Zo_def default bottom roughness
Sout(izapp) == F ! Zo_app apparent bottom roughness
Sout(izNik) == F ! Zo_Nik Nikuradse bottom roughness
Sout(izbio) == F ! Zo_bio biological bottom roughness
Sout(izbfm) == F ! Zo_bedform bed form bottom roughness
Sout(izbld) == F ! Zo_bedload bed load bottom roughness
Sout(izwbl) == F ! Zo_wbl wave bottom roughness
Sout(iactv) == F ! active_layer_t... active layer thickness
Sout(ishgt) == F ! saltation saltation height
! Number of stations to process in each nested grid. These values are
! essential because the station arrays are dynamically allocated using
! these values, [1:Ngrids].
NSTATION == 45
! Station locations for all grids in any desired order. The horizontal
! location for a particular station may be specified in terms of fractional
! (I,J) grid pairs (FLAG=0) or (longitude,latitude) grid pairs (FLAG=1).
! Here, FLAG is a special switch and may be used for multiple purposes.
! The GRID column indicates nested grid number to process. This value must
! be one in non-nested applications. The COMMENT section is ignored during
! reading and may be used to help documentation.
POS = GRID FLAG X-POS Y-POS COMMENT
1 1 9.618 59.104 ! BC-1 Frierfjorden
1 1 10.405 59.532 ! D-1 Ytre Drammensfjorden
1 1 10.421 59.628 ! D-2 Indre Drammensfjord
1 1 10.314 59.706 ! D-3 Indre Drammensfjord - Solumstranda
1 1 10.419 59.572 ! D-4 Solbergbukta
1 1 9.723 59.039 ! FG-1 Langesundsfjorden
1 1 9.797 59.023 ! GI-1 Håøyfjorden
1 1 11.002 59.109 ! I-1 Ramsø - Østerelva
1 1 11.159 59.085 ! I-2 Singlefjorden
1 1 10.962 59.147 ! I-3 Kjøkøy Øst
1 1 11.369 59.102 ! ID-1 Skysskaffern
1 1 11.385 59.075 ! ID-2 Kjellvik - Iddefjorden
1 1 10.473 59.434 ! HO-1 HortenIndre
1 1 10.769 59.31 ! KF-1 Krokstadfjorden
1 1 10.764 59.308 ! KF-1 Krokstadfjorden
1 1 10.706 59.307 ! KUR-1 Kurefjorden
1 1 10.052 59.019 ! LA-1 Larviksfjorden
1 1 10.374 59.485 ! LØ-1 Langøya
1 1 10.662 59.453 ! M-3 Kjellandsvik
1 1 10.675 59.468 ! M-4 Kambo
1 1 10.673 59.505 ! M-6 N-Jeløya
1 1 10.276 59.099 ! ME-1 Mefjorden-Indre
1 1 10.285 59.083 ! ME-2 Mefjorden-Ytre
1 1 10.66 59.443 ! MO-1 Indre Mossesundet
1 1 10.678 59.484 ! MO-2 Kippenes
1 1 10.754 59.041 ! OF-1 Torbjørnskjær
1 1 10.692 59.187 ! OF-2 Missingene
1 1 10.591 59.359 ! OF-4 Bastø
1 1 10.458 59.487 ! OF-5 Breiangen
1 1 10.64 59.59 ! OF-7 Filtvedt
1 1 11.314 59.112 ! R-5 Ringdalsfjorden
1 1 10.717 59.235 ! RA-1 Rauerfjorden
1 1 11.218 59.087 ! S-6 Kjeøy
1 1 11.162 59.114 ! S-9 Haslau - Singlefjorden
1 1 10.266 59.556 ! SAN-A Sandebukta
1 1 10.244 59.562 ! SAN-1 Sandebukta
1 1 10.257 59.551 ! SAN-3 Sandebukta
1 1 10.301 59.517 ! SAN-4 Kommersøya - Holmestrand
1 1 10.304 59.509 ! SAN-8 Sandebukta
1 1 10.246 59.077 ! SF-1 Sandefjord
1 1 10.233 59.1 ! SF-2 Sandefjord-Trangsholmen2 (Sandf-1)
1 1 10.226 59.117 ! SF-3 Sandefjord-Havn1 (Sandf-2)
1 1 10.834 59.137 ! Ø-1 Leira - Vesterelva
1 1 10.355 59.203 ! TØ-1 Vestfjorden
1 1 10.491 59.214 ! TØ-2 Vallø
! GLOSSARY:
! =========
!
! This input file provides the positions of stations data to be extracted
! and written into STATIONS output file. This file is use to monitor model
! fields at frequent time intervals. It is also used to sample the model at
! observation locations in order to build the penalty function used in
! variational data assimilation.
!
!------------------------------------------------------------------------------
! Control Switch.
!------------------------------------------------------------------------------
!
! Lstations Switch to control the writing of station data within nested
! and/or multiple connected grids. Ngrid values are expected.
! By default this switch is set to TRUE in "mod_scalars" for
! all grids when the CPP option STATIONS is activated. The
! USER can control which grid to process by turning on/off
! this switch.
!
!------------------------------------------------------------------------------
! Logical switches (TRUE/FALSE) to activate writing of fields in STATIONS
! output file. Except for tracer type variables, Ngrids values are expected
! for all switches.
!------------------------------------------------------------------------------
!
! Sout(idUvel) Write out 3D U-velocity component.
! Sout(idVvel) Write out 3D V-velocity component.
! Sout(idu3dE) Write out 3D Eastward velocity component at RHO-points.
! Sout(idv3dN) Write out 3D Northward velocity component at RHO-points.
! Sout(idWvel) Write out 3D W-velocity component.
! Sout(idOvel) Write out 3D omega vertical velocity.
! Sout(idUbar) Write out 2D U-velocity component.
! Sout(idVbar) Write out 2D V-velocity component.
! Sout(idu2dE) Write out 2D Eastward velocity component at RHO-points.
! Sout(idv2dN) Write out 2D Northward velocity component at RHO-points.
! Sout(idFsur) Write out free-surface.
! Sout(idBath) Write out time-dependent bathymetry.
!
! Sout(idTvar) Write out tracer type variables, (NT,Ngrids) values are
! Expected.
!
! Sout(idUsms) Write out surface U-momentum stress.
! Sout(idVsms) Write out surface V-momentum stress.
! Sout(idUbms) Write out bottom U-momentum stress.
! Sout(idVbms) Write out bottom V-momentum stress.
!
! Sout(idUbrs) Write out current-induced, U-momentum stress.
! Sout(idVbrs) Write out current-induced, V-momentum stress.
! Sout(idUbws) Write out wind-induced, bottom U-wave stress.
! Sout(idVbws) Write out wind-induced, bottom V-wave stress.
! Sout(idUbcs) Write out bottom maximum wave and current U-stress.
! Sout(idVbcs) Write out bottom maximum wave and current V-stress.
!
! Sout(idUbot) Write out wind-induced, bed wave orbital U-velocity.
! Sout(idVbot) Write out wind-induced, bed wave orbital V-velocity.
! Sout(idUbur) Write out bottom U-velocity above bed.
! Sout(idVbvr) Write out bottom V-velocity above bed.
!
! Sout(idW2xx) Write out 2D radiation stress, Sxx component.
! Sout(idW2xy) Write out 2D radiation stress, Sxy component.
! Sout(idW2yy) Write out 2D radiation stress, Syy component.
! Sout(idU2rs) Write out 2D U-radiation stress.
! Sout(idV2rs) Write out 2D V-radiation stress.
! Sout(idU2Sd) Write out 2D U-Stokes velocity.
! Sout(idV2Sd) Write out 2D V-Stokes velocity.
!
! Sout(idW3xx) Write out 3D radiation stress, Sxx component.
! Sout(idW3xy) Write out 3D radiation stress, Sxy component.
! Sout(idW3yy) Write out 3D radiation stress, Syy component.
! Sout(idW3zx) Write out 3D radiation stress, Szx component.
! Sout(idW3zy) Write out 3D radiation stress, Szy component.
! Sout(idU3rs) Write out 3D U-radiation stress.
! Sout(idV3rs) Write out 3D V-radiation stress.
! Sout(idU3Sd) Write out 3D U-Stokes velocity.
! Sout(idV3Sd) Write out 3D V-Stokes velocity.
!
! Sout(idWamp) Write out wave height.
! Sout(idWlen) Write out wave length.
! Sout(idWdir) Write out wave direction.
! Sout(idWptp) Write out wave surface period.
! Sout(idWpbt) Write out wave bottom period.
! Sout(idWorb) Write out wave bottom orbital velocity.
! Sout(idWdis) Write out wave dissipation.
!
! Sout(idPair) Write out surface air pressure.
! Sout(idUair) Write out surface U-wind component.
! Sout(idVair) Write out surface V-wind component.
!
! Sout(idTsur) Write out surface net heat and salt flux.
! Sout(idLhea) Write out latent heat flux.
! Sout(idShea) Write out sensible heat flux.
! Sout(idLrad) Write out long-wave radiation flux.
! Sout(idSrad) Write out short-wave radiation flux.
! Sout(idEmPf) Write out E-P flux.
! Sout(idevap) Write out evaporation rate.
! Sout(idrain) Write out precipitation rate.
!
! Sout(idDano) Write out density anomaly.
! Sout(idVvis) Write out vertical viscosity coefficient.
! Sout(idTdif) Write out vertical diffusion coefficient of temperature.
! Sout(idSdif) Write out vertical diffusion coefficient of salinity.
! Sout(idHsbl) Write out depth of oceanic surface boundary layer.
! Sout(idHbbl) Write out depth of oceanic bottom boundary layer.
! Sout(idMtke) Write out turbulent kinetic energy.
! Sout(idMtls) Write out turbulent kinetic energy times length scale.
!
! Sout(idBott) Write out exposed sediment layer properties, 1:MBOTP.
!
!------------------------------------------------------------------------------
! Station positions KEYWORDS.
!------------------------------------------------------------------------------
!
! NSTATION Number of stations to process in each nested grid. Ngrids
! values are expected. These values are essential because
! the station arrays are dynamically allocated using these
! values.
!
! POS Station locations for all grids in any desired order. The
! horizontal location for a particular station can be
! specified in terms of fractional (I,J) grid pairs (FLAG=0)
! or (longitude,latitude) grid pairs (FLAG=1). Here, FLAG is
! a special switch and may be used for multiple purposes.
! The GRID column indicates nested grid number to process.
! This value must be one in non-nested applications. The
! COMMENT section is ignored during reading and may be used
! to help documentation. This locations are used to
! interpolate model data before writting to output STATIONS
! file. The stations locations are read using the following
! statement:
!
! is(1:Ngrids)=0
! DO WHILE (.true.)
! read (unit,*,ERR=10) igrid, flag, Xpos, Ypos
! ng=MAX(1,ABS(ng))
! is(ng)=ic(ng)+1
! SCALARS(ng)%Sflag(is(ng))=flag
! SCALARS(ng)%SposX(is(ng))=Xpos
! SCALARS(ng)%SposY(is(ng))=Ypos
! END DO
!
! Since there is no decoding during the reading of the stations
! locations, the POS keyword must be the last one to process.
! That is, it should be at the end of the input specification.
!
!