! Station Extraction Parameters.
!
!svn $Id: stations_test_head.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) == 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) == F ! 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) == T ! bath time-dependent bathymetry
Sout(idTvar) == T T ! temp, salt, ... all (NT) tracers
Sout(idUsms) == F ! sustr surface U-stress
Sout(idVsms) == F ! svstr surface V-stress
Sout(idUbms) == F ! bustr bottom U-stress
Sout(idVbms) == F ! 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) == F F ! shflux, ssflux surface net heat and salt flux
Sout(idLhea) == F ! latent latent heat flux
Sout(idShea) == F ! sensible sensible heat flux
Sout(idLrad) == F ! lwrad longwave radiation flux
Sout(idSrad) == F ! swrad shortwave radiation flux
Sout(idEmPf) == F ! EminusP E-P flux
Sout(idevap) == F ! evaporation evaporation rate
Sout(idrain) == F ! rain precipitation rate
Sout(idDano) == F ! rho density anomaly
Sout(idVvis) == F ! AKv vertical viscosity
Sout(idTdif) == F ! AKt vertical T-diffusion
Sout(idSdif) == F ! AKs vertical Salinity diffusion
Sout(idHsbl) == F ! Hsbl depth of surface boundary layer
Sout(idHbbl) == F ! 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.
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 == 71
! 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 0 120 262 "PORTS 8726694 -TPA Cruise Terminal 2"
1 0 138 251 "PORTS 8726674 -East Bay"
1 0 120 248 "PORTS 8726673 -Seabulk"
1 0 123 244 "PORTS 8726669 -Berth 223"
1 0 125 58 "NWLON 8726217 Cortez "
1 0 115 59 "NWLON 8726243 Anna Maria Out "
1 0 120 85 "NWLON 8726273 Desoto Point "
1 0 89 72 "NWLON 8726347 Egmont Key "
1 0 88 86 "NWLON 8726364 Mulet Key "
1 0 127 128 "NWLON 8726384 Port Manatee "
1 0 64 101 "NWLON 8726428 Tierra Verde "
1 0 51 149 "NWLON 8726520 St. Petersburg "
1 0 165 182 "NWLON 8726537 Apollo Beach "
1 0 86 189 "NWLON 8726607 Old Tampa Bay "
1 0 99 203 "NWLON 8726641 Gandy Bridge "
1 0 122 238 "NWLON 8726657 Davis Island "
1 0 138 246 "NWLON 8726667 Mckay Bay Ent "
1 0 2 249 "NWLON 8726689 Bay Aristocrat "
1 0 1 1 "NWLON 8726724 Clearwater Bea "
1 0 18 276 "NWLON 8726738 Safety Harbor "
1 0 174 1 "NWLON 8725858 Venice Pier "
1 0 80 43 "C_01 Outer Egmont "
1 0 86 71 "C_02 Inner Egmont "
1 0 97 108 "PORTS t01010 -Sunshine Skyway Brid-C_03"
1 0 99 130 "PORTS t03010 -Port Manatee Current Predictions-C_04"
1 0 85 189 "PORTS t02010 -Old Port Tampa Currents C_05 "
1 0 73 1 "C_06 Gulf of America "
1 0 57 152 "C_10 t1002s--St Petersb "
1 0 82 139 "C_11 t1102s--Little Bay "
1 0 80 112 "C_12 t1201s--Tarpon Key "
1 0 95 109 "C_13 t1302s--Sun Bridge "
1 0 107 139 "C_14 Beacon Key "
1 0 120 155 "C_15 Mangrove Pt., Cut E "
1 0 95 60 "C_20 Southwest Channel "
1 0 102 62 "C_21 t2101s--PassageKey "
1 0 118 66 "C_22 t2201s--Anna Maria "
1 0 94 94 "C_23 Mullet Key "
1 0 62 104 "C_24 t2401s--PinellasBu "
1 0 131 89 "C_26 t2601s--ManateeRiv "
1 0 133 230 "C_27 t2701s--PinellasPt "
1 0 74 170 "C_30 Cut J Channel "
1 0 140 200 "C_31 Cut C Channel "
1 0 94 202 "C_32 t3201s--GandyBrEst "
1 0 92 196 "C_33 t3301s--GandyBrWst "
1 0 80 221 "C_34 t3401s--FranklndBr "
1 0 51 254 "C_35 t3501s--CourtneyBr "
1 0 145 187 "C_36 Gadsden Point "
1 0 76 177 "C_40 t4001s--Xmas isl "
1 0 84 181 "C_41 Pinic I, Cut K "
1 0 79 161 "C_42 t4201s--Venebian I "
1 0 63 187 "C_43 4301s--Old TB Str "
1 0 118 168 "C_44 Intersec Y-cut "
1 0 130 198 "C_46 t4601s--AlphiaRv "
1 0 88 128 "C_50 t5001s--Point Pine "
1 0 93 130 "C_51 t5101s--Outer Mana "
1 0 106 130 "C_52 Manatee Channel "
1 0 115 129 "C_53 t5301s--Inner Mana "
1 0 104 86 "C_54 t5401s--Waterway L "
1 0 99 130 "C_55 Manatee Channel "
1 0 85 72 "C_56 Egmont Chan 2 "
1 0 80 42 "C_60 Outer Egmont 2 "
1 0 82 70 "S1 outer Egmont Channel "
1 0 92 126 "S2 Port Manatee Channel PORTS 8726413 -C-CUT "
1 0 138 199 "S3 Lower Hillsboro. Bay "
1 0 92 167 "HIST h00001 -Cut G Channel Buoy 6G"
1 0 123 158 "HIST h00002 -Cut E Channel Buoy 4E"
1 0 62 91 "HIST h000003 -Pass A Grille Channel"
1 1 -82.445243 27.846195 "Alifia Entrance Near Hillsbourough C-Cut"
1 1 -82.556397 27.693264 "Near C-Cut"
1 1 -82.983074 27.590473 "Sea Buoy"
1 1 -82.516667 27.773333 "PORTS 8726524 -Gadsden Cut, Tampa Bay, FL"
!
! 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(isd50) Write out mean grain diameter.
! Sout(idens) Write out mean grain density.
! Sout(iwsed) Write out mean settling velocity.
! Sout(itauc) Write out critical erosion stress.
! Sout(irlen) Write out ripple length.
! Sout(irhgt) Write out ripple height.
! Sout(ibwav) Write out wave excursion amplitude.
! Sout(izdef) Write out default bottom roughness.
! Sout(izapp) Write out apparent bottom roughness.
! Sout(izNik) Write out Nikuradse bottom roughness.
! Sout(izbio) Write out biological bottom roughness.
! Sout(izbfm) Write out bed form bottom roughness.
! Sout(izbld) Write out bed load bottom roughness.
! Sout(izwbl) Write out wave bottom roughness.
! Sout(iactv) Write out active layer thickness.
! Sout(ishgt) Write out saltation height.
!
!------------------------------------------------------------------------------
! 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.
!
!