-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathstations.in
More file actions
executable file
·456 lines (435 loc) · 28.2 KB
/
stations.in
File metadata and controls
executable file
·456 lines (435 loc) · 28.2 KB
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
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
! Station Extraction Parameters.
!
!svn $Id: stations.in 889 2018-02-10 03:32:52Z arango $
!========================================================= Hernan G. Arango ===
! Copyright (c) 2002-2018 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) == F ! u_eastward 3D U-eastward at RHO-points
Sout(idv3dN) == F ! 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) == T ! ubar 2D U-velocity
Sout(idVbar) == T ! 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) == 22*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) == F ! 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) == F ! AKv vertical viscosity
Sout(idTdif) == T ! AKt vertical T-diffusion
Sout(idSdif) == F ! AKs vertical Salinity diffusion
Sout(idHsbl) == T ! Hsbl depth of surface boundary layer
Sout(idHbbl) == T ! Hbbl depth of bottom boundary layer
Sout(idMtke) == F ! tke turbulent kinetic energy
Sout(idMtls) == F ! gls turbulent length scale
! Logical switches (TRUE/FALSE) to activate writing of exposed sediment
! layer properties into STATIONS output file. Currently, MBOTP properties
! are expected for the bottom boundary layer and/or sediment models:
!
! idBott( 1=isd50) grain_diameter mean grain diameter
! idBott( 2=idens) grain_density mean grain density
! idBott( 3=iwsed) settling_vel mean settling velocity
! idBott( 4=itauc) erosion_stress critical erosion stress
! idBott( 5=irlen) ripple_length ripple length
! idBott( 6=irhgt) ripple_height ripple height
! idBott( 7=ibwav) bed_wave_amp wave excursion amplitude
! idBott( 8=izdef) Zo_def default bottom roughness
! idBott( 9=izapp) Zo_app apparent bottom roughness
! idBott(10=izNik) Zo_Nik Nikuradse bottom roughness
! idBott(11=izbio) Zo_bio biological bottom roughness
! idBott(12=izbfm) Zo_bedform bed form bottom roughness
! idBott(13=izbld) Zo_bedload bed load bottom roughness
! idBott(14=izwbl) Zo_wbl wave bottom roughness
! idBott(15=iactv) active_layer_thickness active layer thickness
! idBott(16=ishgt) saltation saltation height
!
! 1 1 1 1 1 1 1
! 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6
Sout(idBott) == F F F F F F F F F F F F F F F F
! 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 == 140 ! Added 12 new stations for A.Bever.
! NSTATION == 128
! 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 -76.0383 38.2200 !water level 1 1 Bishops Head
1 1 -76.0683 38.5733 ! 2 2 Cambridge
1 1 -76.5283 39.2200 ! 3 3 Francis Scott Key Bridge
1 1 -76.4800 38.9833 ! 4 4 Annapolis
1 1 -76.3867 38.4317 ! 5 5 Cove Point
1 1 -76.4517 38.3167 ! 6 6 Solomons Island
1 1 -76.5333 38.1333 ! 7 7 Piney Point
1 1 -75.9883 37.1667 ! 8 8 Kiptopeke
1 1 -76.0150 37.5383 ! 9 9 Rappahannoc Light
1 1 -76.4650 37.9950 ! 10 10 Lewisetta
1 1 -76.3333 37.2500 ! 11 11 York River East Rear Range Light
1 1 -76.4783 37.2267 ! 12 12 Yorktown USCG Training Center
1 1 -76.3383 36.8883 ! 13 13 South Craney Island
1 1 -76.3300 36.9467 ! 14 14 Sewells Point
1 1 -76.1133 36.9667 ! 15 15 Chesapeake Bay Bridge Tunnel
1 1 -75.6967 36.9050 ! 16 16 Chesapeake Light
1 1 -76.0067 36.9300 ! 17 17 Cape Henry
1 1 -75.7300 36.1918 ! 91.d0 2.d0 18 18 Duck
1 1 -75.5989 37.6204 ! 89.d0 33.d0 19 19 Wach
1 1 -76.2500 39.2133 ! 20 20 Tolchester
1 1 -76.5783 39.2667 ! 21 21 Baltimore
1 1 -76.38000 38.47333 ! CBOS stations 22 1 MB
1 1 -76.11333 39.37500 ! 23 2 NB
1 1 -76.15927 38.63498 ! 24 3 CHOP
1 1 -76.52650 38.38867 ! 25 4 PAX
1 1 -76.08134 39.54511 !bay program T,S 26 1 CB1.1
1 1 -76.02467 39.44011 ! 27 2 CB2.1
1 1 -76.17467 39.34678 ! 28 3 CB2.2
1 1 -76.23801 39.24844 ! 29 4 CB3.1
1 1 -76.30634 39.16345 ! 30 5 CB3.2
1 1 -76.35968 38.99511 ! 31 6 CB3.3C
1 1 -76.34634 39.00178 ! 32 7 CB3.3E
1 1 -76.38802 39.00344 ! 33 8 CB3.3W
1 1 -76.39440 38.92706 ! 34 9 CB4.0C
1 1 -76.38690 38.92706 ! 35 10 CB4.0E
1 1 -76.43273 38.92733 ! 36 11 CB4.0W
1 1 -76.39967 38.82511 ! 37 12 CB4.1C
1 1 -76.37106 38.81650 ! 38 13 CB4.1E
1 1 -76.46273 38.81345 ! 39 14 CB4.1W
1 1 -76.41773 38.64484 ! 40 15 CB4.2C
1 1 -76.39995 38.64484 ! 41 16 CB4.2E
1 1 -76.50134 38.64345 ! 42 17 CB4.2W
1 1 -76.43467 38.55651 ! 43 18 CB4.3C
1 1 -76.38967 38.55651 ! 44 19 CB4.3E
1 1 -76.49301 38.55651 ! 45 20 CB4.3W
1 1 -76.34300 38.41318 ! 46 21 CB4.4
1 1 -76.29300 38.31846 ! 47 22 CB5.1
1 1 -76.37550 38.32512 ! 48 23 CB5.1W
1 1 -76.22800 38.13679 ! 49 24 CB5.2
1 1 -76.16800 37.91179 ! 50 25 CB5.3
1 1 -76.17466 37.80013 ! 51 26 CB5.4
1 1 -76.29467 37.81347 ! 52 27 CB5.4W
1 1 -76.18967 37.69180 ! 53 28 CB5.5
1 1 -76.16216 37.58847 ! 54 29 CB6.1
1 1 -76.15633 37.48680 ! 55 30 CB6.2
1 1 -76.15966 37.41153 ! 56 31 CB6.3
1 1 -76.20799 37.23653 ! 57 32 CB6.4
1 1 -75.98966 37.68346 ! 58 33 CB7.1
1 1 -75.97466 37.77513 ! 59 34 CB7.1N
1 1 -76.05799 37.58124 ! 60 35 CB7.1S
1 1 -76.07966 37.41153 ! 61 36 CB7.2
1 1 -76.02466 37.41153 ! 62 37 CB7.2E
1 1 -76.12521 37.11681 ! 63 38 CB7.3
1 1 -76.05383 37.22875 ! 64 39 CB7.3E
1 1 -76.02048 36.99570 ! 65 40 CB7.4
1 1 -75.99937 37.06237 ! 66 41 CB7.4N
1 1 -76.16772 36.99542 ! 67 42 CB8.1
1 1 -76.03494 36.94737 ! 68 43 CB8.1E
1 1 -76.24967 38.88345 ! 69 44 EE1.1
1 1 -76.27467 38.65012 ! 70 45 EE2.1
1 1 -76.30801 38.53345 ! 71 46 EE2.2
1 1 -76.01633 38.28345 ! 72 47 EE3.0
1 1 -75.97465 38.20012 ! 73 48 EE3.1
1 1 -75.92467 37.98010 ! 74 49 EE3.2
1 1 -75.76632 37.94179 ! 75 50 EE3.3
1 1 -75.79132 37.90846 ! 76 51 EE3.4
1 1 -75.84438 37.79652 ! 77 52 EE3.5
1 1 -75.95800 39.57511 ! 78 53 ET1.1
1 1 -75.89967 39.50845 ! 79 54 ET2.3
1 1 -76.21634 38.99178 ! 80 55 ET4.2
1 1 -76.05800 38.58012 ! 81 56 ET5.2
1 1 -75.88299 38.33345 ! 82 57 ET6.2
1 1 -76.60162 38.42513 ! 83 58 LE1.1
1 1 -76.51106 38.37873 ! 84 59 LE1.2
1 1 -76.48801 38.34068 ! 85 60 LE1.3
1 1 -76.42134 38.31207 ! 86 61 LE1.4
1 1 -76.58302 38.16679 ! 87 62 LE2.2
1 1 -76.34772 38.02152 ! 88 63 LE2.3
1 1 -76.61592 37.75925 ! 89 64 LE3.1
1 1 -76.55047 37.66914 ! 90 65 LE3.2
1 1 -76.47442 37.68842 ! 91 66 LE3.3
1 1 -76.44486 37.63189 ! 92 67 LE3.4
1 1 -76.28467 37.59680 ! 93 68 LE3.6
1 1 -76.28300 37.60680 ! 94 69 LE3.6N
1 1 -76.29300 37.57263 ! 95 70 LE3.6S
1 1 -76.30661 37.53069 ! 96 71 LE3.7
1 1 -76.69125 37.41883 ! 97 72 LE4.1
1 1 -76.57811 37.29044 ! 98 73 LE4.2
1 1 -76.43094 37.23392 ! 99 74 LE4.3
1 1 -76.64833 37.20297 ! 100 75 LE5.1
1 1 -76.59306 37.05600 ! 101 76 LE5.2
1 1 -76.47544 36.99044 ! 102 77 LE5.3
1 1 -76.39275 36.95486 ! 103 78 LE5.4
1 1 -76.30300 36.99681 ! 104 79 LE5.5
1 1 -76.29050 36.98180 ! 105 80 LE5.5A
1 1 -76.28966 36.97458 ! 106 81 LE5.5B
1 1 -76.31328 36.99903 ! 107 82 LE5.5-W
1 1 -76.33836 36.90456 ! 108 83 LE5.6
1 1 -76.66413 38.49068 ! 109 84 RET1.1
1 1 -76.99053 38.36263 ! 110 85 RET2.4
1 1 -76.82220 37.91730 ! 111 86 RET3.1
1 1 -76.71195 37.81158 ! 112 87 RET3.2
1 1 -76.79715 37.57125 ! 113 88 RET4.2
1 1 -76.78889 37.50869 ! 114 89 RET4.3
1 1 -76.78219 37.20294 ! 115 90 RET5.2
1 1 -76.68024 38.58179 ! 116 91 TF1.7
1 1 -77.04859 38.70650 ! 117 92 TF2.1
1 1 -76.90928 38.01847 ! 118 93 TF3.3
1 1 -76.92329 37.22182 ! 119 94 TF5.6A
1 1 -76.34634 37.31181 ! 120 95 WE4.1
1 1 -76.38634 37.24181 ! 121 96 WE4.2
1 1 -76.39050 37.25181 ! 122 97 WE4.2N
1 1 -76.38634 37.23681 ! 123 98 WE4.2S
1 1 -76.37300 37.17681 ! 124 99 WE4.3
1 1 -76.29300 37.11015 ! 125 100 WE4.4
1 1 -76.34162 39.38344 ! 126 101 WT2.1
1 1 -76.52469 39.20844 ! 127 102 WT5.1
1 1 -76.47469 39.07511 ! 128 103 WT6.1
1 1 -76.073378 39.54025 ! 129 1 CBIBS_SUS
1 1 -76.39108 39.151892 ! 130 2 CBIBS_PAT
1 1 -76.446768 38.963602 ! 131 3 CBIBS_ANA
1 1 -77.035672 38.78765 ! 132 4 CBIBS_UPT
1 1 -76.414605 38.556375 ! 133 5 CBIBS_GOR
1 1 -76.335572 38.033345 ! 134 6 CBIBS_POT
1 1 -76.26037 37.568152 ! 135 7 CBIBS_STR
1 1 -76.777298 37.204095 ! 136 8 CBIBS_JAM
1 1 -76.297501 36.870227 ! 137 9 CBIBS_NOR_Moved_In_Grid
1 1 -76.045418 36.97414 ! 138 10 CBIBS_FSL
1 1 -77.114333 37.305068 ! 139 11 CBIBS_VCU
1 1 -76.39524 38.56424 ! 140 12 CBIBS_GOR2_WaterDepthMatch
!
! 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.
!
!