!  Station Extraction Parameters.
!
!svn $Id: stations_bl_test.in 62 2007-05-16 01:36:06Z arango $
!========================================================= Hernan G. Arango ===
!  Copyright (c) 2002-2007 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.  All comments will 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(idBath) == F                          ! time-dependent bathymetry
Sout(idFsur) == T                          ! free-surface
Sout(idUbar) == T                          ! 2D U-velocity
Sout(idVbar) == T                          ! 2D V-velocity
Sout(idUvel) == T                          ! 3D U-velocity
Sout(idVvel) == T                          ! 3D V-velocity
Sout(idWvel) == T                          ! 3D W-velocity
Sout(idOvel) == T                          ! 3D omega vertical velocity

Sout(idTvar) == T T T                       ! all (NT) tracers

Sout(idUsms) == T                          ! surface U-stress
Sout(idVsms) == T                          ! surface V-stress
Sout(idUbms) == T                          ! bottom U-stress
Sout(idVbms) == T                          ! bottom V-stress
Sout(idUbws) == F                          ! bottom U-wave stress
Sout(idVbws) == F                          ! bottom V-wave stress

Sout(idUbed) == F                          ! bed wave orbital U-velocity
Sout(idVbed) == F                          ! bed wave orbital V-velocity
Sout(idUbot) == F                          ! bottom U-momentum above bed
Sout(idVbot) == F                          ! bottom V-momentum above bed

Sout(idTsur) == T F                        ! surface net heat and salt flux
Sout(idLhea) == F                          ! latent heat flux
Sout(idShea) == F                          ! sensible heat flux
Sout(idLrad) == F                          ! longwave radiation flux
Sout(idSrad) == F                          ! shortwave radiation flux
Sout(idevap) == F                          ! evaporation rate
Sout(idrain) == F                          ! precipitation rate

Sout(idDano) == T                          ! density anomaly
Sout(idVvis) == T                          ! vertical viscosity
Sout(idTdif) == T                          ! vertical T-diffusion
Sout(idSdif) == T                          ! vertical Salinity diffusion
Sout(idHsbl) == F                          ! depth of surface boundary layer
Sout(idHbbl) == F                          ! depth of bottom boundary layer
Sout(idMtke) == F                          ! turbulent kinetic energy
Sout(idMtls) == F                          ! 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:
!
!   Sout(idBott(isd50)),  isd50 = 1        ! mean grain diameter
!   Sout(idBott(idens)),  idens = 2        ! mean grain density
!   Sout(idBott(iwsed)),  iwsed = 3        ! mean settling velocity
!   Sout(idBott(itauc)),  itauc = 4        ! critical erosion stress
!   Sout(idBott(irlen)),  irlen = 5        ! ripple length
!   Sout(idBott(irhgt)),  irhgt = 6        ! ripple height
!   Sout(idBott(ibwav)),  ibwav = 7        ! wave excursion amplitude
!   Sout(idBott(izdef)),  izdef = 8        ! default bottom roughness
!   Sout(idBott(izapp)),  izapp = 9        ! apparent bottom roughness
!   Sout(idBott(izNik)),  izNik = 10       ! Nikuradse bottom roughness
!   Sout(idBott(izbio)),  izbio = 11       ! biological bottom roughness
!   Sout(idBott(izbfm)),  izbfm = 12       ! bed form bottom roughness
!   Sout(idBott(izbld)),  izbld = 13       ! bed load bottom roughness
!   Sout(idBott(izwbl)),  izwbl = 14       ! wave bottom roughness
!   Sout(idBott(iactv)),  iactv = 15       ! active layer thickness
!   Sout(idBott(ishgt)),  ishgt = 16       ! 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 == 124 

! 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         268           1      ! WL station 8570283 Ocean City, MD
         1    0         258         172      ! WL station 8571892 Cambridge, MD
         1    0         186         222      ! WL station 8573364 Tolchester Beach, MD
         1    0         147         217      ! WL station 8574680 Baltimore, MD
         1    0         166         196      ! WL station 8575512 Annapolis, MD
         1    0         158         158      ! WL station 8577330 Solomons Island, MD
         1    0         198          56      ! WL station 8632200 Kiptopeke, VA
         1    0         168          44      ! WL station 8638863 CBBT, VA
         1    0         153          50      ! WL station 8638610 Sewells Point, VA
         1    0         160         129      ! WL station 8635750 Lewisetta, VA
         1    0         163           1      ! WL station 8651370 Duck, NC
         1    0         228         279      ! WL station 8573927 Chesapeake City, MD
         1    0         164         101      ! WL station 8636580 Windmill Pt., VA
         1    0         113         129      ! WL station 8635150 Colonial Beach, VA
         1    0           1         132      ! WL station 8594900 Washington, DC
         1    0         231         147      ! WL station 8571559 Mc Creadys Creek, MD
         1    0         160         253      ! WL station 8574070 Havre De Grace, MD
         1    0         218         279      ! Current CMIST station cb1301 (Chesapeake City)
         1    0         184         171      ! Current CMIST station CHB0304/0305 (Taylors Island)
         1    0         169         126      ! Current CMIST station CHB0303 (Smith Point)
         1    0         155          51      ! Current CMIST station cb0401 (Nav stn Norfolk)
         1    0         150          46      ! Current CMIST station cb0501 (South Craney Island)
         1    0         180         106      ! Current CMIST station cb0801 (Rappahannock Shoal Channel)
         1    0         171         164      ! Current CMIST station cb1001 (Cove Point LNG Pier)
         1    0         180         240      ! Current CBOS North Bay Buoy station
         1    0         178          42      ! Bay Program T/S station CB7.4
         1    0         171          67      ! Bay Program T/S station CB6.4
         1    0         178          85      ! Bay Program T/S station CB6.3
         1    0         178          99      ! Bay Program T/S station CB6.1
         1    0         178         114      ! Bay Program T/S station CB5.4
         1    0         181         138      ! Bay Program T/S station CB5.2
         1    0         179         164      ! Bay Program T/S station CB4.4
         1    0         183         180      ! Bay Program T/S station CB4.2C
         1    0         186         189      ! Bay Program T/S station CB4.1C
         1    0         176         217      ! Bay Program T/S station CB3.2
         1    0         175         234      ! Bay Program T/S station CB2.2
         1    0         174         253      ! Bay Program T/S station CB2.1
         1    0         146          51      ! Bay Program T/S station LE5.4
         1    0          91          51      ! Bay Program T/S station RET5.2
         1    0         150          69      ! Bay Program T/S station LE4.3
         1    0         162          76      ! Bay Program T/S station WE4.1
         1    0         145         102      ! Bay Program T/S station LE3.4
         1    0         115         102      ! Bay Program T/S station LE3.1
         1    0         163         131      ! Bay Program T/S station LE2.3
         1    0         143         131      ! Bay Program T/S station LE2.2
         1    0          97         133      ! Bay Program T/S station RET2.4
         1    0         162         156      ! Bay Program T/S station LE1.4
         1    0         146         157      ! Bay Program T/S station LE1.3
         1    0         114         157      ! Bay Program T/S station LE1.1
         1    0         150         216      ! Bay Program T/S station WT5.1
         1    0         158         230      ! Bay Program T/S station WT2.1
         1    0         159         255      ! Bay Program T/S station CB1.1
         1    0         214         202      ! Bay Program T/S station ET4.2
         1    0         205         192      ! Bay Program T/S station EE1.1
         1    0         205         179      ! Bay Program T/S station EE2.1
         1    0         259         174      ! Bay Program T/S station ET5.2
         1    0         227         138      ! Bay Program T/S station EE3.1
         1    0         227         125      ! Bay Program T/S station EE3.2
         1    0         242         119      ! Bay Program T/S station EE3.4
         1    0         234         111      ! Bay Program T/S station EE3.5
         1    0         168         180      ! Bay Program T/S station CB4.2W
         1    0         167         176      ! Bay Program T/S station CB4.3W
         1    0         179         152      ! Bay Program T/S station CB5.1
         1    0         245         145      ! Bay Program T/S station ET6.2
         1    0         181         122      ! Bay Program T/S station CB5.3
         1    0         173         107      ! Bay Program T/S station CB5.5
         1    1   -75.998841    36.996644    ! Bay axis transect points
         1    1   -76.068393    37.110738    ! Bay axis transect points
         1    1   -76.122488    37.238255    ! Bay axis transect points
         1    1   -76.145672    37.419463    ! Bay axis transect points
         1    1   -76.134080    37.607383    ! Bay axis transect points
         1    1   -76.130216    37.802013    ! Bay axis transect points
         1    1   -76.161128    38.010067    ! Bay axis transect points
         1    1   -76.234544    38.151007    ! Bay axis transect points
         1    1   -76.307960    38.291946    ! Bay axis transect points
         1    1   -76.362056    38.432886    ! Bay axis transect points
         1    1   -76.404560    38.607383    ! Bay axis transect points
         1    1   -76.423879    38.775168    ! Bay axis transect points
         1    1   -76.392968    38.942953    ! Bay axis transect points
         1    1   -76.342736    39.104027    ! Bay axis transect points
         1    1   -76.277048    39.238255    ! Bay axis transect points
         1    1   -76.195904    39.338926    ! Bay axis transect points
         1    1   -76.095440    39.399329    ! Bay axis transect points
         1    1   -76.025889    39.439597    ! Bay axis transect points
         1    1   -75.935954    39.487529    ! Bay axis transect points
         1    1   -75.872992    39.526615    ! Bay axis transect points
         1    1   -75.996173    36.936878    ! Bay across transect-1 points
         1    1   -75.992694    36.960640    ! Bay across transect-1 points
         1    1   -75.986898    36.981762    ! Bay across transect-1 points
         1    1   -75.982260    37.005524    ! Bay across transect-1 points
         1    1   -75.977623    37.029286    ! Bay across transect-1 points
         1    1   -75.972985    37.050408    ! Bay across transect-1 points
         1    1   -75.969507    37.071530    ! Bay across transect-1 points
         1    1   -76.236652    37.379437    ! Bay across transect-2 points
         1    1   -76.199992    37.358041    ! Bay across transect-2 points
         1    1   -76.165950    37.336646    ! Bay across transect-2 points
         1    1   -76.131909    37.309902    ! Bay across transect-2 points
         1    1   -76.100486    37.283157    ! Bay across transect-2 points
         1    1   -76.066445    37.261762    ! Bay across transect-2 points
         1    1   -76.032403    37.240366    ! Bay across transect-2 points
         1    1   -76.225081    37.892927    ! Bay across transect-3 points
         1    1   -76.192781    37.914323    ! Bay across transect-3 points
         1    1   -76.167403    37.935718    ! Bay across transect-3 points
         1    1   -76.139717    37.962462    ! Bay across transect-3 points
         1    1   -76.112031    37.989207    ! Bay across transect-3 points
         1    1   -76.086653    38.015951    ! Bay across transect-3 points
         1    1   -76.386320    38.400840    ! Bay across transect-4 points
         1    1   -76.372864    38.410640    ! Bay across transect-4 points
         1    1   -76.362398    38.417990    ! Bay across transect-4 points
         1    1   -76.353427    38.427790    ! Bay across transect-4 points
         1    1   -76.344456    38.437591    ! Bay across transect-4 points
         1    1   -76.335485    38.447391    ! Bay across transect-4 points
         1    1   -76.326514    38.457191    ! Bay across transect-4 points
         1    1   -76.389391    39.006112    ! Bay across transect-5 points
         1    1   -76.379659    38.997023    ! Bay across transect-5 points
         1    1   -76.370517    38.985169    ! Bay across transect-5 points
         1    1   -76.362850    38.974698    ! Bay across transect-5 points
         1    1   -76.354592    38.964227    ! Bay across transect-5 points
         1    1   -76.137387    39.403041    ! Bay across transect-6 points
         1    1   -76.131508    39.396797    ! Bay across transect-6 points
         1    1   -76.126609    39.390553    ! Bay across transect-6 points
         1    1   -76.120730    39.384308    ! Bay across transect-6 points
         1    1   -76.115831    39.378956    ! Bay across transect-6 points
         1    1   -76.110931    39.372712    ! Bay across transect-6 points

!
!  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 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(idBath)  Write out time-dependent bathymetry.
!  Sout(idFsur)  Write out free-surface.
!  Sout(idUbar)  Write out 2D U-velocity component.
!  Sout(idVbar)  Write out 2D V-velocity component.
!  Sout(idUvel)  Write out 3D U-velocity component.
!  Sout(idVvel)  Write out 3D V-velocity component.
!  Sout(idWvel)  Write out 3D W-velocity component.
!  Sout(idOvel)  Write out 3D omega vertical velocity.
!
!  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(idUbws)  Write out wind-induced, bottom U-wave stress.
!  Sout(idVbws)  Write out wind-induced, bottom V-wave stress.
!
!  Sout(idUbed)  Write out wind-induced, bed wave orbital U-velocity.
!  Sout(idVbed)  Write out wind-induced, bed wave orbital V-velocity.
!  Sout(idUbot)  Write out bottom U-momentum above bed.
!  Sout(idVbot)  Write out bottom V-momentum above bed.
!
!  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(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 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.
!                
!
