stations.in

From WikiROMS
Jump to navigationJump to search
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.


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

Note 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(idBath) == T  ! Write out time-dependent bathymetry
Sout(idFsur) == T  ! Write out free-surface
Sout(idUbar) == T  ! Write out 2D U-velocity
Sout(idVbar) == T  ! Write out 2D V-velocity
Sout(idUvel) == T  ! Write out 3D U-velocity
Sout(idVvel) == T  ! Write out 3D V-velocity
Sout(idWvel) == T  ! Write out 3D W-velocity
Sout(idOvel) == T  ! Write out 3D omega vertical velocity

Sout(idTvar) == T T  ! Write out tracer type variables, (NT,Ngrids) values are expected

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

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