stations.in

Standard Input Script - stations.in

The stations.in file sets the parameters and locations for stations output. The name of this file is set by the SPOSNAM keyword in the ocean.in file. This standard input ASCII file is organized in several sections as shown below, with links to more detailed explanation where required.

Notice: Detailed information about ROMS input script file syntax can be found here.

Notice: A default stations.in input script is provided in the User/External subdirectory. Also there are several stations input scripts in the ROMS/External subdirectory which are used in the distributed test cases. They are usually named stations_app.in where app is the lowercase of the test case CPP option.

Introduction

This input file provides the positions of model station data to be extracted and written into STATIONS output NetCDF file. This file is used to monitor model fields at frequent time intervals. It is like putting a moored hydrographic station inside the model. It may be used to sample model at observation locations for inter-comparison analysis. It is also a good tool for debugging; a station can be put at a trouble spot to help find algorithm or set-up problems. It is like putting a print statement in the code to display the values of all variables at a particular location.

Stations Output Model Control Switches

• This switch is used to control the writing of stations data within nested and/or multiple connected grids, [1:Ngrids] values are expected. By default this switch is set to TRUE in mod_scalars.F for all grids when the CPP option STATIONS is activated. The user can control which grid to process by turning on/off this switch.

Lstations == T

• This set of switches determine what fields are written in the stations output file, [Sout(:,ng), ng=1, Ngrids]. Except for tracer type variables, Ngrids values are expected for all switches.

Sout(idUvel) == T  ! u 3D U-velocity
Sout(idVvel) == T  ! v 3D V-velocity
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(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(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(idEmPf) == F  ! E-P 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

Sediment Switches

• This set of switches 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. There are currently 16 properties for the exposed sediment layer.
  isd50 = 1  ! mean grain diameter
  idens = 2  ! mean grain density
  iwsed = 3  ! mean settling velocity
  itauc = 4  ! critical erosion stress
  irlen = 5  ! ripple length
  irhgt = 6  ! ripple height
  ibwav = 7  ! wave excursion amplitude
  izdef = 8  ! default bottom roughness
  izapp = 9  ! apparent bottom roughness
  izNik = 10  ! Nikuradse bottom roughness
  izbio = 11  ! biological bottom roughness
  izbfm = 12  ! bed form bottom roughness
  izbld = 13  ! bed load bottom roughness
  izwbl = 14  ! wave bottom roughness
  iactv = 15  ! active layer thickness
  ishgt = 16  ! saltation height
  Since there are 16 properties Sout(idBott) must have 16 switches:
  ! 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

Station Locations

• 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.

NSTATION == 8

• 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. These locations are used to interpolate model data before writing to output stations file (STANAME). 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.

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