!  Station Extraction Parameters.
!
!svn $Id: stations.in 294 2009-01-09 21:37:26Z arango $
!========================================================= Hernan G. Arango ===
!  Copyright (c) 2002-2009 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(idBath) == T                          ! 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                        ! all (NT) tracers

Sout(idUsms) == F                          ! surface U-stress
Sout(idVsms) == F                          ! surface V-stress
Sout(idUbms) == F                          ! bottom U-stress
Sout(idVbms) == F                          ! 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) == F 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) == T                          ! 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 ==  25

! 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      12.5416     44.8691     ! real position: 12.50673 44.86875
         1    1      14.26613    42.46277
         1    1      12.41580    44.36130
         1    1      13.7797     43.2927     ! real position: 13.75660 43.29690
         1    1      13.87950    43.3331
         1    1      12.46140    44.79600
         1    1      12.49320    44.76940
         1    1      12.85512    44.46100    !  8  CP2
         1    1      13.12453    44.54018    !  9  CP3
         1    1      14.02125    44.75067    ! 10  KB1
         1    1      13.30660    43.83513    ! 11  SS2
         1    1      13.42610    43.93067    ! 12  SS5
         1    1      13.90525    44.25672    ! 13  SS8
         1    1      14.29040    44.48123    ! 14  SS0 (aka SS10)
         1    1      13.36667    43.88355    ! 15  SS4
         1    1      13.50438    43.99562    ! 16  SS6
         1    1      14.17478    44.41022    ! 17 SS9
         1    1      12.63697    45.27892    ! 18 VR2
         1    1      13.02810    45.18778    ! 19 VR4
         1    1      13.28367    45.12493    ! 20 VR5
         1    1      13.53603    45.05807    ! 21 VR6
         1    1      12.51440    45.31663    ! 22 VR1
         1    1      12.49370    44.76980    ! 23 Po 20-m
         1    1      13.87970    43.33290    ! 24 Chienti 20-m
         1    1      14.29800    42.48060    ! 25 Pescara 20-m

!
!  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(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 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.
!                
!
