!
!  ROMS/TOMS IO NetCDF variables.
!
!========================================================= Hernan G. Arango ===
!  Copyright (c) 2002-2020 The ROMS/TOMS Group                                !
!    Licensed under a MIT/X style license                                     !
!    See License_ROMS.txt                                                     !
!==============================================================================
!
! Note! This file contains some Tclm entries for historic reasons.
! These will no longer be read, but contain information on the
! climatology nudging used in CoSine. Many fields are using the same
! values since there is so little data.
! For instance, idTclm(iS1_N) points to NH4.
!
!==============================================================================
!
!  This file contains all ROMS/TOMS Input/Output variables information.  It   !
!  is used to manage NetCDF files.  The following information is needed for   !
!  any IO variable:                                                           !
!                                                                             !
!       Vinfo(1)      Field variable name.                                    !
!       Vinfo(2)      Long-name attribute.                                    !
!       Vinfo(3)      Units attribute.                                        !
!       Vinfo(4)      Field type attribute.                                   !
!       Vinfo(5)      Associated time variable name.                          !
!       Vinfo(6)      Index variable name used in information arrays.         !
!       Vinfo(7)      Staggered C-grid variable type:                         !
!                       'nulvar' => non-grided variable.                      !
!                       'p2dvar' => 2D PHI-variable.                          !
!                       'r2dvar' => 2D RHO-variable.                          !
!                       'u2dvar' => 2D U-variable.                            !
!                       'v2dvar' => 2D V-variable.                            !
!                       'p3dvar' => 3D PHI-variable.                          !
!                       'r3dvar' => 3D RHO-variable.                          !
!                       'u3dvar' => 3D U-variable.                            !
!                       'v3dvar' => 3D V-variable.                            !
!                       'w3dvar' => 3D W-variable.                            !
!                       'b3dvar' => 3D BED-sediment.                          !
!                       'l3dvar' => 3D spectral light variable.               !
!                       'l4dvar' => 4D spectral light variable.               !
!       scale         Scale to convert input data to model units.             !
!                                                                             !
!  Enclose Vinfo variables with single quotes so this file can be read with   !
!  a free format statement.  The  Vinfo(6) variable must be unique and case   !
!  sensitive.  It is checked in "mod_ncparm" to set the appropiate ID index   !
!  for the information arrays.                                                !
!                                                                             !
!  Notice that there is commented information for each variable  including    !
!  input/output type and internal model units in brackets.  There are some    !
!  variables that are special and require a particular units for Input and    !
!  and output.                                                                !
!                                                                             !
!==============================================================================
!
! The value in Vinfo(5) here is used to name the unlimited dimension in all
! output NetCDF files.  By default it is set to 'time'.  Alternatively, you
! can use 'ocean_time' as dimension name so both variable and dimension
! names are the same.  This facilitates processing in some NetCDF packages.
! Notice that we use the prefix 'ocean' in the variable name to differentiate
! between oceanic and atmospheric time clocks in fully coupled applications.
!
! Notice that salinity does not have physical units. A commented [PSU] can
! be found below to indicate that the practical salinity scale was used to
! determine conductivity. See,
!
! http://marine.rutgers.edu/po/users/forum/viewtopic.php?t=294
!
!------------------------------------------------------------------------------
!  Git Repository
!------------------------------------------------------------------------------
!

'$URL$'

!
!------------------------------------------------------------------------------
!  Model state variables.
!------------------------------------------------------------------------------
!

'ocean_time'                                       ! Input/Output
  'time since initialization'
  'second'                                         ! [s]
  'time, scalar, series'
  'ocean_time'
  'idtime'
  'nulvar'
  1.0d0

'bath'                                             ! Input/Output
  'bathymetry'
  'meter'                                          ! [m]
  'bathymetry, scalar, series'
  'ocean_time'
  'idbath'
  'r2dvar'
  1.0d0

'zeta'                                             ! Input/Output
  'free-surface'
  'meter'                                          ! [m]
  'free-surface, scalar, series'
  'ocean_time'
  'idFsur'
  'r2dvar'
  1.0d0

'ubar'                                             ! Input/Ouput
  'vertically integrated u-momentum component'
  'meter second-1'                                 ! [m/s]
  'ubar-velocity, scalar, series'
  'ocean_time'
  'idUbar'
  'u2dvar'
  1.0d0

'vbar'                                             ! Input/Output
  'vertically integrated v-momentum component'
  'meter second-1'                                 ! [m/s]
  'vbar-velocity, scalar, series'
  'ocean_time'
  'idVbar'
  'v2dvar'
  1.0d0

'ubar_eastward'                                    ! Ouput
  'eastward vertically integrated momentum component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'ubar_eastward, scalar, series'
  'ocean_time'
  'idu2dE'
  'r2dvar'
  1.0d0

'vbar_northward'                                   ! Output
  'northward vertically integrated momentum component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'vbar_northward, scalar, series'
  'ocean_time'
  'idv2dN'
  'r2dvar'
  1.0d0

'u'                                                ! Input/Output
  'u-momentum component'
  'meter second-1'                                 ! [m/s]
  'u-velocity, scalar, series'
  'ocean_time'
  'idUvel'
  'u3dvar'
  1.0d0

'v'                                                ! Input/Output
  'v-momentum component'
  'meter second-1'                                 ! [m/s]
  'v-velocity, scalar, series'
  'ocean_time'
  'idVvel'
  'v3dvar'
  1.0d0

'u_eastward'                                       ! Output
  'eastward momentum component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'u_eastward, scalar, series'
  'ocean_time'
  'idu3dE'
  'r3dvar'
  1.0d0

'v_northward'                                      ! Output
  'northward momentum component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'v_northward, scalar, series'
  'ocean_time'
  'idv3dN'
  'r3dvar'
  1.0d0

'u_sur'                                            ! Output
  'model surface u-momentum component'
  'meter second-1'                                 ! [m/s]
  'u_sur, scalar, series'
  'ocean_time'
  'idUsur'
  'u2dvar'
  1.0d0

'v_sur'                                            ! Output
  'model surface v-momentum component'
  'meter second-1'                                 ! [m/s]
  'v_sur, scalar, series'
  'ocean_time'
  'idVsur'
  'v2dvar'
  1.0d0

'u_sur_eastward'                                   ! Output
  'model surface eastward momentum component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'u_sur_eastward, scalar, series'
  'ocean_time'
  'idUsuE'
  'r2dvar'
  1.0d0

'v_sur_northward'                                  ! Output
  'model surface northward momentum component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'v_sur_northward, scalar, series'
  'ocean_time'
  'idVsuN'
  'r2dvar'
  1.0d0

'w'                                                ! Output
  'vertical momentum component'
  'meter second-1'                                 ! [m/s]
  'w-velocity, scalar, series'
  'ocean_time'
  'idWvel'
  'w3dvar'
  1.0d0

'omega'                                            ! Output
  'S-coordinate vertical momentum component'
  'meter3 second-1'                                ! [m3/s]
  'omega, scalar, series'
  'ocean_time'
  'idOvel'
  'w3dvar'
  1.0d0

'omega'                                            ! Output
  'S-coordinate vertical momentum component'
  'meter3 second-1'                                ! [m3/s]
  'omega, scalar, series'
  'ocean_time'
  'idOclm'
  'w3dvar'
  1.0d0

'rho'                                              ! Output
  'density anomaly'
  'kilogram meter-3'                               ! [kg/m3]
  'density, scalar, series'
  'ocean_time'
  'idDano'
  'r3dvar'
  1.0d0

'temp'                                             ! Input/Output
  'potential temperature'
  'Celsius'                                        ! [Celsius]
  'temperature, scalar, series'
  'ocean_time'
  'idTvar(itemp)'
  'r3dvar'
  1.0d0

'salt'                                             ! Input/Output
  'salinity'
  'nondimensional'                                 ! [PSU]
  'salinity, scalar, series'
  'ocean_time'
  'idTvar(isalt)'
  'r3dvar'
  1.0d0

'temp_sur'                                         ! Output
  'model surface potential temperature'
  'Celsius'                                        ! [Celsius]
  'temp_sur, scalar, series'
  'ocean_time'
  'idsurT(itemp)'
  'r2dvar'
  1.0d0

'salt_sur'                                         ! Output
  'model surface salinity'
  'nondimensional'                                 ! [PSU]
  'salinity, scalar, series'
  'ocean_time'
  'idsurT(isalt)'
  'r2dvar'
  1.0d0

'z_rho'                                            ! Output
  'depth of RHO-points'
  'meter'                                          ! [meter]
  'z_rho, scalar, series'
  'ocean_time'
  'idpthR'
  'r3dvar'
  1.0d0

'z_u'                                              ! Output
  'depth of U-points'
  'meter'                                          ! [meter]
  'z_u, scalar, series'
  'ocean_time'
  'idpthU'
  'u3dvar'
  1.0d0

'z_v'                                              ! Output
  'depth of V-points'
  'meter'                                          ! [meter]
  'z_v, scalar, series'
  'ocean_time'
  'idpthV'
  'v3dvar'
  1.0d0

'z_w'                                              ! Output
  'depth of W-points'
  'meter'                                          ! [meter]
  'z_w, scalar, series'
  'ocean_time'
  'idpthW'
  'w3dvar'
  1.0d0

'EBF'                                           ! Output
  'Ekman buoyancy flux at PSI-points'
  'Watt kg-1'                                   ! [m/s]
  'EBF, scalar, series'
  'ocean_time'
  'idEBF'
  'p2dvar'
  1.0d0

'strain'                                           ! Output
  'strain'
  'second-1'                                       ! [1/s]
  'strain, scalar, series'
  'ocean_time'
  'idStr'
  'p3dvar'
  1.0d0

'strain_sur'                                       ! Output
  'surface strain'
  'second-1'                                       ! [1/s]
  'strain_sur, scalar, series'
  'ocean_time'
  'idStr2'
  'p2dvar'
  1.0d0

'rvort_sur'                                        ! Output
  'surface relative vorticity'
  'second-1'                                       ! [1/s]
  'rvort_sur, scalar, series'
  'ocean_time'
  'idRVsu'
  'p2dvar'
  1.0d0

'pvort_sur'                                        ! Output
  'surface potential vorticity'
  'meter-1 second-1'                               ! [1/m/s]
  'pvort_sur, scalar, series'
  'ocean_time'
  'idPVsu'
  'p2dvar'
  1.0d0

'pvortk_sur'                                        ! Output
  'surface potential vorticity, k component'
  'meter-1 second-1'                               ! [1/m/s]
  'pvortk_sur, scalar, series'
  'ocean_time'
  'idPVKsu'
  'p2dvar'
  1.0d0

'pvortk'                                       ! Input/Output
  'potential vorticity, k component'
  'meter-1 second-1'                               ! [1/m/s]
  'pvorticity k component, scalar, series'
  'ocean_time'
  'id3dPVK'
  'p3dvar'
  1.0d0

'pvorti'                                       ! Input/Output
  'potential vorticity, i component'
  'meter-1 second-1'                               ! [1/m/s]
  'pvorticity, i component, scalar, series'
  'ocean_time'
  'id3dPVI'
  'p3dvar'
  1.0d0

'pvortj'                                       ! Input/Output
  'potential vorticity, j component'
  'meter-1 second-1'                               ! [1/m/s]
  'pvorticity, j component, scalar, series'
  'ocean_time'
  'id3dPVJ'
  'p3dvar'
  1.0d0

'pvorti_sur'                                        ! Output
  'surface potential vorticity, i component'
  'meter-1 second-1'                               ! [1/m/s]
  'pvorti_sur, scalar, series'
  'ocean_time'
  'idPVIsu'
  'p2dvar'
  1.0d0

'pvortj_sur'                                        ! Output
  'surface potential vorticity, j component'
  'meter-1 second-1'                               ! [1/m/s]
  'pvortj_sur, scalar, series'
  'ocean_time'
  'idPVJsu'
  'p2dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Wet/dry masks.
!------------------------------------------------------------------------------
!

'wetdry_mask_psi'                                  ! Input/Output
  'wet/dry mask on PSI-points'
  'nondimensional'
  'wetdry_mask_psi, scalar, series'
  'ocean_time'
  'idPwet'
  'p2dvar'
  1.0d0

'wetdry_mask_rho'                                  ! Input/Output
  'wet/dry mask on RHO-points'
  'nondimensional'
  'wetdry_mask_rho, scalar, series'
  'ocean_time'
  'idRwet'
  'r2dvar'
  1.0d0

'wetdry_mask_u'                                    ! Input/Output
  'wet/dry mask on U-points'
  'nondimensional'
  'wetdry_mask_u, scalar, series'
  'ocean_time'
  'idUwet'
  'u2dvar'
  1.0d0

'wetdry_mask_v'                                    ! Input/Output
  'wet/dry mask on V-points'
  'nondimensional'
  'wetdry_mask_v, scalar, series'
  'ocean_time'
  'idVwet'
  'v2dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Detided variables
!------------------------------------------------------------------------------
!

'CosW'                                             ! Input/Output
  'time-accumulated COS(omega(k)*t)'
  'radians'                                        ! [radians]
  'CosW, scalar'
  'tide_period'
  'idCosW'
  'nulvar'
  1.0d0

'SinW'                                             ! Input/Output
  'time-accumulated SIN(omega(k)*t)'
  'radians'                                        ! [radians]
  'SinW, scalar'
  'tide_period'
  'idSinW'
  'nulvar'
  1.0d0

'CosWCosW'                                         ! Input/Output
  'time-accumulated COS(omega(k)*t)*COS(omega(l)*t) matrix'
  'radians'                                        ! [radians]
  'CosWCosW, scalar'
  'tide_period'
  'idCos2'
  'nulvar'
  1.0d0

'SinWSinW'                                         ! Input/Output
  'time-accumulated SIN(omega(k)*t)*SIN(omega(l)*t) matrix'
  'radians'                                        ! [radians]
  'SinWSinW, scalar'
  'tide_period'
  'idSin2'
  'nulvar'
  1.0d0

'SinWCosW'                                         ! Input/Output
  'time-accumulated SIN(omega(k)*t)*COS(omega(l)*t) matrix'
  'radians'                                        ! [radians]
  'SinWCosW, scalar'
  'tide_period'
  'idSWCW'
  'nulvar'
  1.0d0

'zeta_tide'                                        ! Input/Output
  'time-accumulated free-surface tide harmonics'
  'meter'                                          ! [m]
  'zeta_tide, scalar'
  'ocean_time'
  'idFsuH'
  'r2dvar'
  1.0d0

'zeta_detided'                                     ! Output
  'detided free-surface'
  'meter'                                          ! [m]
  'zeta_detided, scalar, series'
  'ocean_time'
  'idFsuD'
  'r2dvar'
  1.0d0

'ubar_tide'                                        ! Input/Ouput
  'time-accumulated 2D u-momentum tide harmonics'
  'meter second-1'                                 ! [m/s]
  'ubar_tide, scalar'
  'ocean_time'
  'idu2dH'
  'u2dvar'
  1.0d0

'ubar_detided'                                     ! Ouput
  'detided 2D u-momentum'
  'meter second-1'                                 ! [m/s]
  'ubar_detided, scalar, series'
  'ocean_time'
  'idu2dD'
  'u2dvar'
  1.0d0

'vbar_tide'                                        ! Input/Output
  'time-accumulated 2D v-momentum tide harmonics'
  'meter second-1'                                 ! [m/s]
  'vbar_tide, scalar'
  'ocean_time'
  'idv2dH'
  'v2dvar'
  1.0d0

'vbar_detided'                                     ! Output
  'detided 2D v-momentum'
  'meter second-1'                                 ! [m/s]
  'vbar_detided, scalar, series'
  'ocean_time'
  'idv2dD'
  'v2dvar'
  1.0d0

'u_tide'                                           ! Input/Output
  'time-accumulated 3D u-momentum tide harmonics'
  'meter second-1'                                 ! [m/s]
  'u_tide, scalar'
  'ocean_time'
  'idu3dH'
  'u3dvar'
  1.0d0

'u_detided'                                        ! Output
  'detided 3D u-momentum'
  'meter second-1'                                 ! [m/s]
  'u_detided, scalar, series'
  'ocean_time'
  'idu3dD'
  'u3dvar'
  1.0d0

'v_tide'                                           ! Input/Output
  'time-accumulated 3D v-momentum tide harmonics'
  'meter second-1'                                 ! [m/s]
  'v_tide, scalar'
  'ocean_time'
  'idv3dH'
  'v3dvar'
  1.0d0

'v_detided'                                        ! Output
  'detided 3D v-momentum'
  'meter second-1'                                 ! [m/s]
  'v_detided, scalar, series'
  'ocean_time'
  'idv3dD'
  'v3dvar'
  1.0d0

'temp_tide'                                        ! Input/Output
  'time-accumulated potential temperature tide harmonics'
  'Celsius'                                        ! [Celsius]
  'temp_tide, scalar'
  'ocean_time'
  'idTrcH(itemp)'
  'r3dvar'
  1.0d0

'temp_detided'                                     ! Output
  'detided potential temperature'
  'Celsius'                                        ! [Celsius]
  'v_detided, scalar, series'
  'ocean_time'
  'idTrcD(itemp)'
  'r3dvar'
  1.0d0

'temp_tide'                                        ! Input/Output
  'time-accumulated potential temperature tide harmonics'
  'Celsius'                                        ! [Celsius]
  'temp_tide, scalar'
  'ocean_time'
  'idTrcH(itemp)'
  'r3dvar'
  1.0d0

'temp_detided'                                     ! Output
  'detided potential temperature'
  'Celsius'                                        ! [Celsius]
  'temp_detided, scalar, series'
  'ocean_time'
  'idTrcD(itemp)'
  'r3dvar'
  1.0d0

'salt_tide'                                        ! Input/Output
  'time-accumulated salinity tide harmonics'
  'nondimensional'                                 ! [PSU]
  'salt_tide, scalar'
  'ocean_time'
  'idTrcH(isalt)'
  'r3dvar'
  1.0d0

'salt_detided'                                     ! Output
  'detided potential temperature'
  'nondimensional'                                 ! [PSU]
  'salt_detided, scalar, series'
  'ocean_time'
  'idTrcD(isalt)'
  'r3dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Momentum stress variables.
!------------------------------------------------------------------------------
!

'sustr'                                            ! Input/Output
  'surface u-momentum stress'
  'newton meter-2'                                 ! Input:  [N/m2]
  'surface u-momentum stress, scalar, series'      ! [m2/s2]
  'sms_time'                                       ! Output: [N/m2]
  'idUsms'
  'u2dvar'
  1.0d0

'svstr'                                            ! Input/Output
  'surface v-momentum stress'
  'newton meter-2'                                 ! Input:  [N/m2]
  'surface v-momentum stress, scalar, series'      ! [m2/s2]
  'sms_time'                                       ! Output: [N/m2]
  'idVsms'
  'v2dvar'
  1.0d0

'bustr'                                            ! Output
  'bottom u-momentum stress'
  'newton meter-2'                                 ! [m2/s2]
  'bottom u-momentum stress, scalar, series'       ! Output: [N/m2]
  'bms_time'
  'idUbms'
  'u2dvar'
  1.0d0

'bvstr'                                            ! Output
  'bottom v-momentum stress'
  'newton meter-2'                                 ! [m2/s2]
  'bottom v-momentum stress, scalar, series'       ! Output: [N/m2]
  'bms_time'
  'idVbms'
  'v2dvar'
  1.0d0

'bustrc'                                           ! Output
  'current-induced, bottom u-momentum stress'
  'newton meter-2'                                 ! [m2/s2]
  'bottom u-momentum stress, scalar, series'       ! Output: [N/m2]
  'bms_time'
  'idUbrs'
  'r2dvar'
  1.0d0

'bvstrc'                                           ! Ouput
  'current-induced, bottom v-momentum stress'
  'newton meter-2'                                 ! [m2/s2]
  'bottom v-momentum stress, scalar, series'       ! Output: [N/m2]
  'bms_time'
  'idVbrs'
  'r2dvar'
  1.0d0

'bustrw'                                           ! Output
  'wind-induced, bottom u-momentum stress'
  'newton meter-2'                                 ! [m2/s2]
  'bottom u-wave stress, scalar, series'           ! Output: [N/m2]
  'ocean_time'
  'idUbws'
  'r2dvar'
  1.0d0

'bvstrw'                                           ! Output
  'wind-induced, bottom v-momentum stress'
  'newton meter-2'                                 ! [m2/s2]
  'bottom v-wave stress, scalar, series'           ! Output: [N/m2]
  'ocean_time'
  'idVbws'
  'r2dvar'
  1.0d0

'bustrcwmax'                                       ! Output
  'maximum wind and current, bottom u-wave stress'
  'newton meter-2'                                 ! [m2/s2]
  'bustrcwmax, scalar, series'                     ! Output: [N/m2]
  'ocean_time'
  'idUbcs'
  'u2dvar'
  1.0d0

'bvstrcwmax'                                       ! Output
  'maximum wind and current, bottom v-wave stress'
  'newton meter-2'                                 ! [m2/s2]
  'bustrcwmax, scalar, series'                     ! Output: [N/m2]
  'ocean_time'
  'idVbcs'
  'v2dvar'
  1.0d0

'Ubot'                                             ! Output
  'wind-induced, bed wave orbital u-velocity'
  'meter second-1'                                 ! [m/s]
  'bed wave u-velocity, scalar, series'
  'ocean_time'
  'idUbot'
  'r2dvar'
  1.0d0

'Vbot'                                             ! Output
  'wind-induced, bed wave orbital v-velocity'
  'meter second-1'                                 ! [m/s]
  'bed wave v-velocity, scalar, series'
  'ocean_time'
  'idVbot'
  'r2dvar'
  1.0d0

'Ur'                                               ! Output
  'bottom u-momentum above bed'
  'meter second-1'                                 ! [m/s]
  'bottom u-velocity, scalar, series'
  'ocean_time'
  'idUbur'
  'r2dvar'
  1.0d0

'Vr'                                               ! Output
  'bottom v-momentum above bed'
  'meter second-1'                                 ! [m/s]
  'bottom v-velocity, scalar, series'
  'ocean_time'
  'idVbvr'
  'r2dvar'
  1.0d0
!
!------------------------------------------------------------------------------
!  Tracer flux variables.
!------------------------------------------------------------------------------
!

'shflux'                                           ! Input/Output
  'surface net heat flux'
  'watt meter-2'                                   ! Input:  [Watt/m2]
  'surface heat flux, scalar, series'              ! [Celsius m/s]
  'shf_time'                                       ! Output: [Watt/m2]
  'idTsur(itemp)'
  'r2dvar'
  1.0d0

'dQdSST'                                           ! Input
  'surface net heat flux sensitivity to SST'
  'watt meter-2 Celsius-1'                         ! Input:  [Watt/m2]
  'dQdSST, scalar, series'                         ! [Celsius m/s]
  'sst_time'
  'iddQdT'
  'r2dvar'
  1.0d0

'ssflux'                                           ! Output
  'surface net salt flux, (E-P)*SALT'
  'meter second-1'                                 ! [PSU m/s]
  'surface net salt flux, scalar, series'
  'ssf_time'
  'idTsur(isalt)'
  'r2dvar'
  1.0d0

'swflux'                                           ! Input
  'surface net freswater flux, (E-P)'
  'centimeter day-1'                               ! Input:  [m/s]
  'surface net salt flux, scalar, series'          ! [PSU m/s]
  'swf_time'
  'idsfwf'
  'r2dvar'
  1.157407d-7                                      ! 0.01/86400

'EminusP'                                          ! Input/Output
  'bulk_flux surface net freswater flux, (E-P)'
  'meter second-1'                                 ! computed by NLM ROMS
  'EminusP, scalar, series'                        ! bulk_flux.F and
  'ocean_time'                                     ! needed in adjoint-based
  'idEmPf'                                         ! applications
  'r2dvar'
  1.0d0

'SSSf'                                             ! Input/Output
  'surface net freswater flux from salinity nudging'
  'meter second-1'                                 ! computed by NLM ROMS
  'SSSf, scalar, series'                           ! bulk_flux.F and
  'ocean_time'                                     ! needed in adjoint-based
  'idSSSf'                                         ! applications
  'r2dvar'
  1.0d0

'fastice'                                          ! Input/Output
  'landfast ice extent'
  '1'                                              !
  'FastIce, scalar, series'                        !
  'ocean_time'                                     !
  'idFastIce'                                      !
  'r2dvar'
  1.0d0

'bhflux'                                           ! Input
  'bottom net heat flux'
  'watt meter-2'                                   ! Input:  [Watt/m2]
  'bottom heat flux, scalar, series'               ! [Celsius m/s]
  'bhf_time'
  'idTbot(itemp)'
  'r2dvar'
  1.0d0

'bwflux'                                           ! Input
  'bottom net freshwater flux'
  'centimeter day-1'                               ! Input:  [m/s]
  'bottom water flux, scalar, series'              ! [PSU m/s]
  'bwf_time'
  'idTbot(isalt)'
  'r2dvar'
  1.157407d-7                                      ! 0.01/86400

!
!------------------------------------------------------------------------------
! Ice-related variables
!------------------------------------------------------------------------------
!

'prate'
  'precipitation rate'
  'kilogram meter-2 second-1'
  'precipitation rate, scalar, series'
  'prate_time'
  'idPrate'
  'r2dvar'
  1.0d0

'icec'
  'observed ice concentration'
  'nondimensional'
  'ice concentration, scalar series'
  'icec_time'
  'idIcec'
  'r2dvar'
  1.0d0

'skt'
  'skin temperature'
  'degrees Kelvin'
  'skin temperature, scalar, series'
  'skt_time'
  'idSkt'
  'r2dvar'
  1.0d0

'uice'
  'u-component of ice velocity'
  'meter second-1'
  'u-component of ice velocity, scalar, series'
  'ocean_time'
  'idUice'
  'u2dvar'
  1.0d0

'uice'
  'u-component of ice velocity'
  'meter second-1'
  'u-component of ice velocity, scalar, series'
  'ocean_time'
  'idUIclm'
  'u2dvar'
  1.0d0

'vice'
  'v-component of ice velocity'
  'meter second-1'
  'v-component of ice velocity, scalar, series'
  'ocean_time'
  'idVice'
  'v2dvar'
  1.0d0

'vice'
  'v-component of ice velocity'
  'meter second-1'
  'v-component of ice velocity, scalar, series'
  'ocean_time'
  'idVIclm'
  'v2dvar'
  1.0d0

'uice_eastward'                                    ! Output
  'ice eastward momentum component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'uice_eastward, scalar, series'
  'ocean_time'
  'idUiceE'
  'r2dvar'
  1.0d0

'vice_northward'                                   ! Output
  'ice northward momentum component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'vice_northward, scalar, series'
  'ocean_time'
  'idViceN'
  'r2dvar'
  1.0d0

'aice'
  'fraction of cell covered by ice'
  'nondimensional'
  'ice concentration, scalar, series'
  'ocean_time'
  'idAice'
  'r2dvar'
  1.0d0

'aice'
  'fraction of cell covered by ice'
  'nondimensional'
  'ice concentration, scalar, series'
  'ocean_time'
  'idAIclm'
  'r2dvar'
  1.0d0

'hice'
  'average ice thickness in cell'
  'meter'
  'ice thickness, scalar, series'
  'ocean_time'
  'idHice'
  'r2dvar'
  1.0d0

'hice'
  'average ice thickness in cell'
  'meter'
  'ice thickness, scalar, series'
  'ocean_time'
  'idHIclm'
  'r2dvar'
  1.0d0

'snow_thick'
  'thickness of snow cover'
  'meter'
  'snow thickness, scalar, series'
  'ocean_time'
  'idHsno'
  'r2dvar'
  1.0d0

'snow_thick'
  'thickness of snow cover'
  'meter'
  'snow thickness, scalar, series'
  'ocean_time'
  'idHSNclm'
  'r2dvar'
  1.0d0

'tisrf'
  'temperature of ice surface'
  'degrees Celcius'
  'surface temperature, scalar, series'
  'ocean_time'
  'idTice'
  'r2dvar'
  1.0d0

'ti'
  'interior ice temperature'
  'degrees Celcius'
  'interior temperature, scalar, series'
  'ocean_time'
  'idTimid'
  'r2dvar'
  1.0d0

't2'
  'top of ice temperature'
  'degrees Celcius'
  'top of ice temperature, scalar, series'
  'ocean_time'
  'idT2ice'
  'r2dvar'
  1.0d0

'tau_iw'
  'ice-water friction velocity'
  'meter second-1'
  'friction velocity, scalar, series'
  'ocean_time'
  'idTauiw'
  'r2dvar'
  1.0d0

'chu_iw'
  'ice-water momentum transfer coefficient'
  'meter second-1'
  'transfer coefficient, scalar, series'
  'ocean_time'
  'idChuiw'
  'r2dvar'
  1.0d0

's0mk'
  'salinity of molecular sub-layer under ice'
  'nondimensional'
  'salinity, scalar, series'
  'ocean_time'
  'idS0mk'
  'r2dvar'
  1.0d0

't0mk'
  'temperature of molecular sub-layer under ice'
  'degrees Celsius'
  'temperature, scalar, series'
  'ocean_time'
  'idT0mk'
  'r2dvar'
  1.0d0

'wfr'
  'rate of frazil ice growth'
  'm^3 s-1'
  'freeze rate, scalar, series'
  'ocean_time'
  'idWfr'
  'r2dvar'
  1.0d0

'wai'
  'rate of ice melt on top'
  'm^3 s-1'
  'freeze/melt rate, scalar, series'
  'ocean_time'
  'idWai'
  'r2dvar'
  1.0d0

'wao'
  'rate of ice growth in leads'
  'm^3 s-1'
  'freeze/melt rate, scalar, series'
  'ocean_time'
  'idWao'
  'r2dvar'
  1.0d0

'wio'
  'rate of ice growth under ice'
  'm^3 s-1'
  'freeze rate, scalar, series'
  'ocean_time'
  'idWio'
  'r2dvar'
  1.0d0

'wro'
  'rate of ice melt runoff'
  'm^3 s-1'
  'melt runoff rate, scalar, series'
  'ocean_time'
  'idWro'
  'r2dvar'
  1.0d0

'wdiv'
  'rate of ice divergence'
  'm^3 s-1'
  'ice divergence rate, scalar, series'
  'ocean_time'
  'idWdiv'
  'r2dvar'
  1.0d0

'iomflx'
  'ice-ocean mass flux'
  'meter second-1'
  'ice-ocean mass flux, scalar, series'
  'ocean_time'
  'idIomflx'
  'r2dvar'
  1.0d0

'ageice'
  'age of the ice'
  'sec'
  'ice age, scalar, series'
  'ocean_time'
  'idAgeice'
  'r2dvar'
  1.0d0

'sig11'
  'internal ice stress 11 component'
  'Newton meter-1'
  'ice stress 11, scalar, series'
  'ocean_time'
  'idSig11'
  'r2dvar'
  1.0d0

'sig22'
  'internal ice stress 22 component'
  'Newton meter-1'
  'ice stress 22, scalar, series'
  'ocean_time'
  'idSig22'
  'r2dvar'
  1.0d0

'sig12'
  'internal ice stress 12 component'
  'Newton meter-1'
  'ice stress 12, scalar, series'
  'ocean_time'
  'idSig12'
  'r2dvar'
  1.0d0

'sig11'
  'internal ice stress 11 component'
  'Newton meter-1'
  'ice stress 11, scalar, series'
  'ocean_time'
  'idS11clm'
  'r2dvar'
  1.0d0

'sig22'
  'internal ice stress 22 component'
  'Newton meter-1'
  'ice stress 22, scalar, series'
  'ocean_time'
  'idS22clm'
  'r2dvar'
  1.0d0

'sig12'
  'internal ice stress 12 component'
  'Newton meter-1'
  'ice stress 12, scalar, series'
  'ocean_time'
  'idS12clm'
  'r2dvar'
  1.0d0

'albedo_ice'
  'ice albedo'
  'fraction'
  'ice albedo, scalar, series'
  'ocean_time'
  'idAlbI'
  'r2dvar'
  1.0d0

'LW_down'
  'downward longwave radiation'
  'W m-2'
  'LW_down, scalar, series'
  'ocean_time'
  'idLWdn'
  'r2dvar'
  1.0d0

'SW_down'
  'downward shortwave radiation'
  'W m-2'
  'SW_down, scalar, series'
  'ocean_time'
  'idSWdn'
  'r2dvar'
  1.0d0

'lat_ice'
  'latent heat flux over ice'
  'W m-2'
  'lat_ice, scalar, series'
  'ocean_time'
  'idLatI'
  'r2dvar'
  1.0d0

'sens_ice'
  'sensible heat flux over ice'
  'W m-2'
  'sens_ice, scalar, series'
  'ocean_time'
  'idSenI'
  'r2dvar'
  1.0d0

'LW_up_ice'
  'upward longwave radiation over ice'
  'W m-2'
  'LW_up_ice, scalar, series'
  'ocean_time'
  'idLWupI'
  'r2dvar'
  1.0d0

'SW_up_ice'
  'upward shortwave radiation over ice'
  'W m-2'
  'SW_up_ice, scalar, series'
  'ocean_time'
  'idSWupI'
  'r2dvar'
  1.0d0

'SW_in_ice'
  'shortwave radiation into ice'
  'W m-2'
  'SW_in_ice, scalar, series'
  'ocean_time'
  'idSWinI'
  'r2dvar'
  1.0d0

'SW_thru_ice'
  'shortwave radiation through ice'
  'W m-2'
  'SW_thru_ice, scalar, series'
  'ocean_time'
  'idSWthruI'
  'r2dvar'
  1.0d0

'saltflux_ice'
  'salt flux under ice'
  'meter second-1'                                 ! [PSU m/s]
  'saltflux_ice, scalar, series'
  'ocean_time'
  'idSaltI'
  'r2dvar'
  1.0d0

'saltflux_leads'
  'salt flux in leads'
  'meter second-1'                                 ! [PSU m/s]
  'saltflux_leads, scalar, series'
  'ocean_time'
  'idSaltL'
  'r2dvar'
  1.0d0

'qao'
  'air-ocean heat flux'
  'W m-2'
  'qao, scalar, series'
  'ocean_time'
  'idQao'
  'r2dvar'
  1.0d0

'qai'
  'ice-air heat flux'
  'W m-2'
  'qai, scalar, series'
  'ocean_time'
  'idQai'
  'r2dvar'
  1.0d0

'qio'
  'ice-ocean heat flux'
  'W m-2'
  'qio, scalar, series'
  'ocean_time'
  'idQio'
  'r2dvar'
  1.0d0

'qi2'
  'ice heat flux'
  'W m-2'
  'qi2, scalar, series'
  'ocean_time'
  'idQi2'
  'r2dvar'
  1.0d0

'snoice'
  'snow-ice conversion'
  'm s-1'
  'snoice, scalar, series'
  'ocean_time'
  'idSnoI'
  'r2dvar'
  1.0d0

'frswpen'
  'shortwave penetration into ice'
  'm s-1'
  'frswpen, scalar, series'
  'ocean_time'
  'idfrswp'
  'r2dvar'
  1.0d0

'rfaci'
  'reduction factor for near-IR fraction of SW under ice'
  'm s-1'
  'rfaci, scalar, series'
  'ocean_time'
  'idrfaci'
  'r2dvar'
  1.0d0

'Runoff'
  'Fresh water runoff from land'
  'kg/s/m^2'
  'runoff, scalar, series'
  'runoff_time'
  'idRunoff'
  'r2dvar'
  1.0d0

'wg2d'
  'surface wind gustiness squared from NCEP fluxes'
  'meter2 second-2'
  'gustiness squared, scalar, series'
  'ocean_time'
  'idWg2d'
  'r2dvar'
  1.0d0

'cdd'
  'momentum transfer coefficient from NCEP fluxes'
  'nondimensional'
  'transfer coefficient, scalar, series'
  'ocean_time'
  'idCdd'
  'r2dvar'
  1.0d0

'chd'
  'sensible heat transfer coefficient from NCEP'
  'nondimensional'
  'transfer coefficient, scalar, series'
  'ocean_time'
  'idChd'
  'r2dvar'
  1.0d0

'ced'
  'latent heat transfer coefficient from NCEP'
  'nondimensional'
  'transfer coefficient, scalar, series'
  'ocean_time'
  'idCed'
  'r2dvar'
  1.0d0

'wg2m'
  'surface wind gustiness squared from model'
  'meter2 second-2'
  'gustiness squared, scalar, series'
  'ocean_time'
  'idWg2m'
  'r2dvar'
  1.0d0

'cdm'
  'momentum transfer coefficient from model'
  'nondimensional'
  'transfer coefficient, scalar, series'
  'ocean_time'
  'idCdm'
  'r2dvar'
  1.0d0

'chm'
  'sensible heat transfer coefficient from model'
  'nondimensional'
  'transfer coefficient, scalar, series'
  'ocean_time'
  'idChm'
  'r2dvar'
  1.0d0

'cem'
  'latent heat transfer coefficient from model'
  'nondimensional'
  'transfer coefficient, scalar, series'
  'ocean_time'
  'idCem'
  'r2dvar'
  1.0d0

'rhoa'
  'near-surface air density from NCEP fluxes'
  'kg meter-3'
  'air density, scalar, series'
  'ocean_time'
  'idRhoa'
  'r2dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Mixing variables.
!------------------------------------------------------------------------------
!

'tke'                                              ! Output
  'turbulent kinetic energy'
  'meter2 second-2'                                ! [m2/s2]
  'tke, scalar, series'
  'ocean_time'
  'idMtke'
  'w3dvar'
  1.0d0

'gls'                                              ! Output
  'turbulent generic length scale'
  'meter3 second-2'                                ! [m3/s2]
  'gls, scalar, series'
  'ocean_time'
  'idMtls'
  'w3dvar'
  1.0d0

'Lscale'                                           ! Output
  'vertical mixing turbulent length scale'
  'meter'                                          ! [m]
  'Lscale, scalar, series'
  'ocean_time'
  'idVmLS'
  'w3dvar'
  1.0d0

'AKv'                                              ! Input/Output
  'vertical viscosity coefficient'
  'meter2 second-1'                                ! [m2/s]
  'AKv, scalar, series'
  'ocean_time'
  'idVvis'
  'w3dvar'
  1.0d0

'AKt'                                              ! Input/Output
  'temperature vertical diffusion coefficient'
  'meter2 second-1'                                ! [m2/s]
  'AKt, scalar, series'
  'ocean_time'
  'idTdif'
  'w3dvar'
  1.0d0

'AKs'                                              ! Input/Output
  'salinity vertical diffusion coefficient'
  'meter2 second-1'                                ! [m2/s]
  'AKs, scalar, series'
  'ocean_time'
  'idSdif'
  'w3dvar'
  1.0d0

'AKk'                                              ! Input/Output
  'Turbulent kinetic energy vertical diffusion coefficient'
  'meter2 second-1'                                ! [m2/s]
  'AKk, scalar, series'
  'ocean_time'
  'idVmKK'
  'w3dvar'
  1.0d0

'AKp'                                              ! Input/Output
  'Turbulent length scale vertical diffusion coefficient'
  'meter2 second-1'                                ! [m2/s]
  'AKp, scalar, series'
  'ocean_time'
  'idVmKP'
  'w3dvar'
  1.0d0

'ghat_temp'                                        ! Input/Output
  'potential temperature nonlocal transport'
  'nondimensional'                                 ! [nondimensional]
  'ghat_temp, scalar, series'
  'ocean_time'
  'idGhat(itemp)'
  'w3dvar'
  1.0d0

'ghat_salt'                                        ! Input/Output
  'salinity nonlocal transport'
  'nondimensional'                                 ! {nondimensional]
  'ghat_salt, scalar, series'
  'ocean_time'
  'idGhat(isalt)'
  'w3dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Boundary variables.
!------------------------------------------------------------------------------
!

'zeta_west'                                        ! Input
  'free-surface western boundary condition'
  'meter'                                          ! [m]
  'zeta_west, scalar, series'
  'ocean_time'
  'idZbry(iwest)'
  'nulvar'
  1.0d0

'zeta_east'                                        ! Input
  'free-surface eastern boundary condition'
  'meter'                                          ! [m]
  'zeta_east, scalar, series'
  'ocean_time'
  'idZbry(ieast)'
  'nulvar'
  1.0d0

'zeta_south'                                       ! Input
  'free-surface southern boundary condition'
  'meter'                                          ! [m]
  'zeta_south, scalar, series'
  'ocean_time'
  'idZbry(isouth)'
  'nulvar'
  1.0d0

'zeta_north'                                       ! Input
  'free-surface northern boundary condition'
  'meter'                                          ! [m]
  'zeta_north, scalar, series'
  'ocean_time'
  'idZbry(inorth)'
  'nulvar'
  1.0d0

'ubar_west'                                        ! Input
  '2D u-momentum western boundary condition'
  'meter second-1'                                 ! [m/s]
  'ubar_west, scalar, series'
  'ocean_time'
  'idU2bc(iwest)'
  'nulvar'
  1.0d0

'ubar_east'                                        ! Input
  '2D u-momentum eastern boundary condition'
  'meter second-1'                                 ! [m/s]
  'ubar_east, scalar, series'
  'ocean_time'
  'idU2bc(ieast)'
  'nulvar'
  1.0d0

'ubar_south'                                       ! Input
  '2D u-momentum southern boundary condition'
  'meter second-1'                                 ! [m/s]
  'ubar_south, scalar, series'
  'ocean_time'
  'idU2bc(isouth)'
  'nulvar'
  1.0d0

'ubar_north'                                       ! Input
  '2D u-momentum northern boundary condition'
  'meter second-1'                                 ! [m/s]
  'ubar_north, scalar, series'
  'ocean_time'
  'idU2bc(inorth)'
  'nulvar'
  1.0d0

'vbar_west'                                        ! Input
  '2D v-momentum western boundary condition'
  'meter second-1'                                 ! [m/s]
  'vbar_west, scalar, series'
  'ocean_time'
  'idV2bc(iwest)'
  'nulvar'
  1.0d0

'vbar_east'                                        ! Input
  '2D v-momentum eastern boundary condition'
  'meter second-1'                                 ! [m/s]
  'vbar_east, scalar, series'
  'ocean_time'
  'idV2bc(ieast)'
  'nulvar'
  1.0d0

'vbar_south'                                       ! Input
  '2D v-momentum southern boundary condition'
  'meter second-1'                                 ! [m/s]
  'vbar_south, scalar, series'
  'ocean_time'
  'idV2bc(isouth)'
  'nulvar'
  1.0d0

'vbar_north'                                       ! Input
  '2D v-momentum northern boundary condition'
  'meter second-1'                                 ! [m/s]
  'vbar_north, scalar, series'
  'ocean_time'
  'idV2bc(inorth)'
  'nulvar'
  1.0d0

'u_west'                                           ! Input
  '3D u-momentum western boundary condition'
  'meter second-1'                                 ! [m/s]
  'u_west, scalar, series'
  'ocean_time'
  'idU3bc(iwest)'
  'nulvar'
  1.0d0

'u_east'                                           ! Input
  '3D u-momentum eastern boundary condition'
  'meter second-1'                                 ! [m/s]
  'u_east, scalar, series'
  'ocean_time'
  'idU3bc(ieast)'
  'nulvar'
  1.0d0

'u_south'                                          ! Input
  '3D u-momentum southern boundary condition'
  'meter second-1'                                 ! [m/s]
  'u_south, scalar, series'
  'ocean_time'
  'idU3bc(isouth)'
  'nulvar'
  1.0d0

'u_north'                                          ! Input
  '3D u-momentum northern boundary condition'
  'meter second-1'                                 ! [m/s]
  'u_north, scalar, series'
  'ocean_time'
  'idU3bc(inorth)'
  'nulvar'
  1.0d0

'v_west'                                           ! Input
  '3D v-momentum western boundary condition'
  'meter second-1'                                 ! [m/s]
  'v_west, scalar, series'
  'ocean_time'
  'idV3bc(iwest)'
  'nulvar'
  1.0d0

'v_east'                                           ! Input
  '3D v-momentum eastern boundary condition'
  'meter second-1'                                 ! [m/s]
  'v_east, scalar, series'
  'ocean_time'
  'idV3bc(ieast)'
  'nulvar'
  1.0d0

'v_south'                                          ! Input
  '3D v-momentum southern boundary condition'
  'meter second-1'                                 ! [m/s]
  'v_south, scalar, series'
  'ocean_time'
  'idV3bc(isouth)'
  'nulvar'
  1.0d0

'v_north'                                          ! Input
  '3D v-momentum northern boundary condition'
  'meter second-1'                                 ! [m/s]
  'v_north, scalar, series'
  'ocean_time'
  'idV3bc(inorth)'
  'nulvar'
  1.0d0

'temp_west'                                        ! Input
  'temperature western boundary condition'
  'Celsius'                                        ! [Celsius]
  'temp_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,itemp)'
  'nulvar'
  1.0d0

'temp_east'                                        ! Input
  'temperature eastern boundary condition'
  'Celsius'                                        ! [Celsius]
  'temp_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,itemp)'
  'nulvar'
  1.0d0

'temp_south'                                       ! Input
  'temperature southern boundary condition'
  'Celsius'                                        ! [Celsius]
  'temp_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,itemp)'
  'nulvar'
  1.0d0

'temp_north'                                       ! Input
  'temperature northern boundary condition'
  'Celsius'                                        ! [Celsius]
  'temp_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,itemp)'
  'nulvar'
  1.0d0

'salt_west'                                        ! Input
  'salinity western boundary condition'
  'nondimensional'                                 ! [PSU]
  'salt_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,isalt)'
  'nulvar'
  1.0d0

'salt_east'                                        ! Input
  'salinity eastern boundary condition'
  'nondimensional'                                 ! [PSU]
  'salt_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,isalt)'
  'nulvar'
  1.0d0

'salt_south'                                       ! Input
  'salinity southern boundary condition'
  'nondimensional'                                 ! [PSU]
  'salt_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,isalt)'
  'nulvar'
  1.0d0

'salt_north'                                       ! Input
  'salinity northern boundary condition'
  'nondimensional'                                 ! [PSU]
  'salt_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,isalt)'
  'nulvar'
  1.0d0

'aice_west'
  'ice concentration western boundary condition'
  ' '
  'aice_west, scalar, series'
  'bry_time'
  'idAibc(iwest)'
  'nulvar'
  1.0d0

'aice_east'
  'ice concentration eastern boundary condition'
  ' '
  'aice_east, scalar, series'
  'bry_time'
  'idAibc(ieast)'
  'nulvar'
  1.0d0

'aice_north'
  'ice concentration northern boundary condition'
  ' '
  'aice_north, scalar, series'
  'bry_time'
  'idAibc(inorth)'
  'nulvar'
  1.0d0

'aice_south'
  'ice concentration southern boundary condition'
  ' '
  'aice_south, scalar, series'
  'bry_time'
  'idAibc(isouth)'
  'nulvar'
  1.0d0

'hice_west'
  'ice thickness western boundary condition'
  'meter'
  'hice_west, scalar, series'
  'bry_time'
  'idHibc(iwest)'
  'nulvar'
  1.0d0

'hice_east'
  'ice thickness eastern boundary condition'
  'meter'
  'hice_east, scalar, series'
  'bry_time'
  'idHibc(ieast)'
  'nulvar'
  1.0d0

'hice_north'
  'ice thickness northern boundary condition'
  'meter'
  'hice_north, scalar, series'
  'bry_time'
  'idHibc(inorth)'
  'nulvar'
  1.0d0

'hice_south'
  'ice thickness southern boundary condition'
  'meter'
  'hice_south, scalar, series'
  'bry_time'
  'idHibc(isouth)'
  'nulvar'
  1.0d0

'snow_thick_west'
  'snow thickness western boundary condition'
  'meter'
  'snow_thick_west, scalar, series'
  'bry_time'
  'idHsnbc(iwest)'
  'nulvar'
  1.0d0

'snow_thick_east'
  'snow thickness eastern boundary condition'
  'meter'
  'snow_thick_east, scalar, series'
  'bry_time'
  'idHsnbc(ieast)'
  'nulvar'
  1.0d0

'snow_thick_north'
  'snow thickness northern boundary condition'
  'meter'
  'snow_thick_north, scalar, series'
  'bry_time'
  'idHsnbc(inorth)'
  'nulvar'
  1.0d0

'snow_thick_south'
  'snow thickness southern boundary condition'
  'meter'
  'snow_thick_south, scalar, series'
  'bry_time'
  'idHsnbc(isouth)'
  'nulvar'
  1.0d0

'ti_west'
  'ice interior temperature western boundary condition'
  'Celsius'
  'ti_west, scalar, series'
  'bry_time'
  'idTibc(iwest)'
  'nulvar'
  1.0d0

'ti_east'
  'ice interior temperature eastern boundary condition'
  'Celsius'
  'ti_east, scalar, series'
  'bry_time'
  'idTibc(ieast)'
  'nulvar'
  1.0d0

'ti_north'
  'ice interior temperature northern boundary condition'
  'Celsius'
  'ti_north, scalar, series'
  'bry_time'
  'idTibc(inorth)'
  'nulvar'
  1.0d0

'ti_south'
  'ice interior temperature southern boundary condition'
  'Celsius'
  'ti_south, scalar, series'
  'bry_time'
  'idTibc(isouth)'
  'nulvar'
  1.0d0

'sig11_west'
  'ice stress component 11 western boundary condition'
  'N meter-1'
  'sig11_west, scalar, series'
  'bry_time'
  'idS11bc(iwest)'
  'nulvar'
  1.0d0

'sig11_east'
  'ice stress component 11 eastern boundary condition'
  'N meter-1'
  'sig11_east, scalar, series'
  'bry_time'
  'idS11bc(ieast)'
  'nulvar'
  1.0d0

'sig11_north'
  'ice stress component 11 northern boundary condition'
  'N meter-1'
  'sig11_north, scalar, series'
  'bry_time'
  'idS11bc(inorth)'
  'nulvar'
  1.0d0

'sig11_south'
  'ice stress component 22 southern boundary condition'
  'N meter-1'
  'sig11_south, scalar, series'
  'bry_time'
  'idS11bc(isouth)'
  'nulvar'
  1.0d0

'sig22_west'
  'ice stress component 22 western boundary condition'
  'N meter-1'
  'sig22_west, scalar, series'
  'bry_time'
  'idS22bc(iwest)'
  'nulvar'
  1.0d0

'sig22_east'
  'ice stress component 22 eastern boundary condition'
  'N meter-1'
  'sig22_east, scalar, series'
  'bry_time'
  'idS22bc(ieast)'
  'nulvar'
  1.0d0

'sig22_north'
  'ice stress component 22 northern boundary condition'
  'N meter-1'
  'sig22_north, scalar, series'
  'bry_time'
  'idS22bc(inorth)'
  'nulvar'
  1.0d0

'sig22_south'
  'ice stress component 12 southern boundary condition'
  'N meter-1'
  'sig22_south, scalar, series'
  'bry_time'
  'idS22bc(isouth)'
  'nulvar'
  1.0d0

'sig12_west'
  'ice stress component 12 western boundary condition'
  'N meter-1'
  'sig12_west, scalar, series'
  'bry_time'
  'idS12bc(iwest)'
  'nulvar'
  1.0d0

'sig12_east'
  'ice stress component 12 eastern boundary condition'
  'N meter-1'
  'sig12_east, scalar, series'
  'bry_time'
  'idS12bc(ieast)'
  'nulvar'
  1.0d0

'sig12_north'
  'ice stress component 12 northern boundary condition'
  'N meter-1'
  'sig12_north, scalar, series'
  'bry_time'
  'idS12bc(inorth)'
  'nulvar'
  1.0d0

'sig12_south'
  'ice stress component 12 southern boundary condition'
  'N meter-1'
  'sig12_south, scalar, series'
  'bry_time'
  'idS12bc(isouth)'
  'nulvar'
  1.0d0

'uice_west'
  'ice u-momentum western boundary condition'
  'meter second-1'
  'uice_west, scalar, series'
  'bry_time'
  'idUibc(iwest)'
  'nulvar'
  1.0d0

'uice_east'
  'ice u-momentum eastern boundary condition'
  'meter second-1'
  'uice_east, scalar, series'
  'bry_time'
  'idUibc(ieast)'
  'nulvar'
  1.0d0

'uice_north'
  'ice u-momentum northern boundary condition'
  'meter second-1'
  'uice_north, scalar, series'
  'bry_time'
  'idUibc(inorth)'
  'nulvar'
  1.0d0

'uice_south'
  'ice u-momentum southern boundary condition'
  'meter second-1'
  'uice_south, scalar, series'
  'bry_time'
  'idUibc(isouth)'
  'nulvar'
  1.0d0

'vice_west'
  'ice v-momentum western boundary condition'
  'meter second-1'
  'vice_west, scalar, series'
  'bry_time'
  'idVibc(iwest)'
  'nulvar'
  1.0d0

'vice_east'
  'ice v-momentum eastern boundary condition'
  'meter second-1'
  'vice_east, scalar, series'
  'bry_time'
  'idVibc(ieast)'
  'nulvar'
  1.0d0

'vice_north'
  'ice v-momentum northern boundary condition'
  'meter second-1'
  'vice_north, scalar, series'
  'bry_time'
  'idVibc(inorth)'
  'nulvar'
  1.0d0

'vice_south'
  'ice v-momentum southern boundary condition'
  'meter second-1'
  'vice_south, scalar, series'
  'bry_time'
  'idVibc(isouth)'
  'nulvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Atmospheric forcing variables.
!------------------------------------------------------------------------------
!

'Pair'                                             ! Input
  'surface air pressure'
  'Pascal'                                       ! [ROMS wants millibar]
  'Pair, scalar, series'
  'pair_time'
  'idPair'
  'r2dvar'
  0.01d0

'Tair'                                             ! Input
  'surface air temperature'
  'Celsius'                                        ! [Kelvin]
  'Tair, scalar, series'
  'tair_time'
  'idTair'
  'r2dvar'
  1.0d0

'Qair'                                             ! Input
  'surface air relative humidity'
  'kg/kg'                                          ! [kg/kg]
  'Qair, scalar, series'
  'qair_time'
  'idQair'
  'r2dvar'
  1.00d0                                           ! was percent

'cloud'                                            ! Input
  'cloud fraction'
  'nondimensional'                                 ! [nondimensional]
  'cloud, scalar, series'
  'cloud_time'
  'idCfra'
  'r2dvar'
  1.0d0

'swrad'                                            ! Input/Output
  'solar shortwave radiation flux'
  'watt meter-2'                                   ! Input:  [Watt/m2]
  'shortwave radiation, scalar, series'            ! [Celsius m/s]
  'srf_time'                                       ! Output: [Watt/m2]
  'idSrad'
  'r2dvar'
  1.0d0

'albedo'                                           ! Input/Output
  'albedo'
  'fraction'                                       ! Input: [fraction]
  'albedo, scalar, series'                         ! [Celsius m/s]
  'albedo_time'                                    ! Output: [fraction]
  'idAlbe'
  'r2dvar'
  1.0d0

'lwrad_down'                                       ! Input
  'downwelling longwave radiation flux'
  'watt meter-2'                                   ! Input:  [Watt/m2]
  'downwelling longwave radiation, scalar, series' ! [Celsius m/s]
  'lrf_time'
  'idLdwn'
  'r2dvar'
  1.0d0

'lwrad'                                            ! Input/Output
  'net longwave radiation flux'
  'watt meter-2'                                   ! Input:  [Watt/m2]
  'longwave radiation, scalar, series'             ! [Celsius m/s]
  'lrf_time'                                       ! Output: [Watt/m2]
  'idLrad'
  'r2dvar'
  1.0d0

'latent'                                           ! Input/Output
  'net latent heat flux'
  'watt meter-2'                                   ! Input:  [Watt/m2]
  'latent heat flux, scalar, series'               ! [Celsius m/s]
  'lhf_time'                                       ! Output: [Watt/m2]
  'idLhea'
  'r2dvar'
  1.0d0

'sensible'                                         ! Input/Output
  'net sensible heat flux'
  'watt meter-2'                                   ! Input:  [Watt/m2]
  'sensible heat flux, scalar, series'             ! [Celsius m/s]
  'shf_time'                                       ! Output: [Watt/m2]
  'idShea'
  'r2dvar'
  1.0d0

'rain'                                             ! Input
  'rain fall rate'
  'kilogram meter-2 second-1'                      ! [kg m2/s]
  'rain, scalar, series'
  'rain_time'
  'idrain'
  'r2dvar'
  1.0d0                                            ! 1/rhow (0.001 m3/kg) if m/s

'snow'                                             ! Input
  'snow fall rate'
  'kilogram meter-2 second-1'                      ! [kg m2/s]
  'snow, scalar, series'
  'snow_time'
  'idsnow'
  'r2dvar'
  1.0d0                                            ! 1/rhow (0.001 m3/kg) if m/s

'evaporation'                                      ! Input
  'evaporation rate'
  'kilogram meter-2 second-1'                      ! [kg m2/s]
  'evaporation, scalar, series'
  'evap_time'
  'idevap'
  'r2dvar'
  1.0d0                                            ! 1/rhow (0.001 m3/kg) if m/s

'Uwind'                                            ! Input
  'surface u-wind component'
  'meter second-1'                                 ! [m/s]
  'u-wind, scalar, series'
  'wind_time'
  'idUair'
  'r2dvar'
  1.0d0

'Vwind'                                            ! Input
  'surface v-wind component'
  'meter second-1'                                 ! [m/s]
  'v-wind, scalar, series'
  'wind_time'
  'idVair'
  'r2dvar'
  1.0d0

'Uwind_eastward'                                   ! Output
  'eastward surface wind component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'Uwind_eastward, scalar, series'
  'ocean_time'
  'idUairE'
  'r2dvar'
  1.0d0

'Vwind_northward'                                  ! Output
  'northward surface wind component at RHO-points'
  'meter second-1'                                 ! [m/s]
  'Vwind_northward, scalar, series'
  'ocean_time'
  'idVairN'
  'r2dvar'
  1.0d0

'Dwave'                                            ! Input
  'wind-induced wave direction'
  'degrees'                                        ! [radians]
  'Dwave, scalar, series'
  'wave_time'
  'idWdir'
  'r2dvar'
  0.017453292519943295                             ! pi/180

'Hwave'                                            ! Input
  'wind-induced significant wave height'
  'meter'                                          ! [m]
  'Hwave, scalar, series'
  'wave_time'
  'idWamp'
  'r2dvar'
  1.0d0

'Lwave'                                            ! Input
  'wind-induced mean wavelength'
  'meter'                                          ! [m]
  'Lwave, scalar, series'
  'wave_time'
  'idWlen'
  'r2dvar'
  1.0

'Pwave_top'                                        ! Input
  'wind-induced peak surface wave Period'
  'second'                                         ! [s]
  'Pwave_top, scalar, series'
  'wave_time'
  'idWptp'
  'r2dvar'
  1.0d0

'Pwave_bot'                                        ! Input
  'wind-induced bottom wave Period'
  'second'                                         ! [s]
  'Pwave_bot, scalar, series'
  'wave_time'
  'idWpbt'
  'r2dvar'
  1.0d0

'Ub_swan'                                          ! Input
  'wind-induced bottom orbital velocity'
  'meter second-1'                                 ! [m/s]
  'Ub_swan, scalar, series'
  'wave_time'
  'idWorb'
  'r2dvar'
  1.0d0

'Wave_dissip'                                      ! Input
  'wave dissipation'
  'Watts meter-2'                                  ! [W/m2]
  'Wave_dissip, scalar, series'
  'wave_time'
  'idWdis'
  'r2dvar'
  0.00097561                                       !1/rho0

'Wave_break'                                       ! Input
  'percent wave breaking'
  'nondimensionsal'                                ! []
  'Wave_break, scalar, series'
  'wave_time'
  'idWbrk'
  'r2dvar'
  1.0d0

'Sxx_bar'                                          ! Output
  'vertically-integrated radiation stress, Sxx-component'
  'meter2 second2'                                 ! [s]
  'Sxx, scalar, series'
  'ocean_time'
  'idW2xx'
  'r2dvar'
  1.0d0

'Sxy_bar'                                          ! Output
  'vertically-integrated radiation stress, Sxy-component'
  'meter2 second2'                                 ! [s]
  'Sxy_bar, scalar, series'
  'ocean_time'
  'idW2xy'
  'r2dvar'
  1.0d0

'Syy_bar'                                          ! Output
  'vertically-integrated radiation stress, Syy-component'
  'meter2 second2'                                 ! [s]
  'Syy_bar, scalar, series'
  'ocean_time'
  'idW2yy'
  'r2dvar'
  1.0d0

'Sxx'                                              ! Output
  '3D horizontal radiation stress, Sxx-component'
  'meter2 second2'                                 ! [s]
  'Sxx, scalar, series'
  'ocean_time'
  'idW3xx'
  'r3dvar'
  1.0d0

'Sxy'                                              ! Output
  '3D horizontal radiation stress, Sxy-component'
  'meter2 second2'                                 ! [s]
  'Sxy, scalar, series'
  'ocean_time'
  'idW3xy'
  'r3dvar'
  1.0d0

'Syy'                                              ! Output
  '3D radiation stress, Syy-component'
  'meter2 second2'                                 ! [s]
  'Syy, scalar, series'
  'ocean_time'
  'idW3yy'
  'r3dvar'
  1.0d0

'Szx'                                              ! Output
  '3D vertical radiation stress, Szx-component'
  'meter2 second2'                                 ! [s]
  'Szx, scalar, series'
  'ocean_time'
  'idW3zx'
  'r3dvar'
  1.0d0

'Szy'                                              ! Output
  '3D vertical radiation stress, Szy-component'
  'meter2 second2'                                 ! [s]
  'Szy, scalar, series'
  'ocean_time'
  'idW3zy'
  'r3dvar'
  1.0d0

'ubar_Rstress'                                     ! Output
  'vertically-integrated u-radiation stress'
  'newton meter-2'                                 ! [N/m2]
  'ubar_Rstress, scalar, series'
  'ocean_time'
  'idU2rs'
  'u2dvar'
  1.0d0

'ubar_stokes'                                      ! Output
  'vertically-integrated u-Stokes drift velocity'
  'meter second-1'                                 ! [m/s]
  'ubar_stokes-velocity, scalar, series'
  'ocean_time'
  'idU2Sd'
  'u2dvar'
  1.0d0

'vbar_Rstress'                                      ! Output
  'vertically-integrated v-radiation stress'
  'newton meter-2'                                  ! [N/m2]
  'ubar_Rstress, scalar, series'
  'ocean_time'
  'idV2rs'
  'v2dvar'
  1.0d0

'vbar_stokes'                                      ! Output
  'vertically-integrated v-Stokes drift velocity'
  'meter second-1'                                 ! [m/s]
  'vbar_stokes, scalar, series'
  'ocean_time'
  'idV2Sd'
  'v2dvar'
  1.0d0

'u_Rstress'                                         ! Output
  'u-radiation stress'
  'newton meter-2'                                  ! [N/m2]
  'u_Rstress, scalar, series'
  'ocean_time'
  'idU3rs'
  'u3dvar'
  1.0d0

'u_stokes'                                         ! Output
  'u-Stokes drift velocity'
  'meter second-1'                                 ! [m/s]
  'u_stokes, scalar, series'
  'wave_time'
  'idU3Sd'
  'u3dvar'
  1.0d0

'v_Rstress'                                         ! Output
  'v-radiation stress'
  'newton meter-2'                                  ! [N/m2]
  'v_Rstress, scalar, series'
  'ocean_time'
  'idV3rs'
  'v3dvar'
  1.0d0

'v_stokes'                                         ! Output
  'v-Stokes drift velocity'
  'meter second-1'                                 ! [m/s]
  'v_stokes, scalar, series'
  'ocean_time'
  'idV3Sd'
  'v3dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Tidal forcing variables.
!------------------------------------------------------------------------------
!

'tide_period'                                      ! Input
  'tidal period'
  'hours'                                          ! [s]
  'tide_period, scalar, series'
  'tide_period'
  'idTper'
  'nulvar'
  3600.0d0                                         ! 60*60

'tide_Eamp'                                        ! Input
  'tidal elevation amplitude'
  'meter'                                          ! [m]
  'SSH_Tamp, scalar, series'
  'tide_period'
  'idTzam'
  'r2dvar'
  1.0d0

'tide_Ephase'                                      ! Input
  'tidal elevation phase angle'
  'degrees'                                        ! [radians]
  'SSH_Tphase, scalar, series'
  'tide_period'
  'idTzph'
  'r2dvar'
  0.017453292519943295                             ! pi/180

'tide_Cphase'                                      ! Input
  'tidal current phase angle'
  'degrees'                                        ! [radians]
  'tide_Cphase, scalar, series'
  'tide_period'
  'idTvph'
  'r2dvar'
  0.017453292519943295                             ! pi/180

'tide_Cangle'                                      ! Input
  'tidal current inclination angle'
  'degrees'                                        ! [radians]
  'tide_Cangle, scalar, series'
  'tide_period'
  'idTvan'
  'r2dvar'
  0.017453292519943295                             ! pi/180

'tide_Cmax'                                        ! Input
  'maximum tidal current, ellipse major axis'
  'meter second-1'                                 ! [m/s]
  'tide_Cmax, scalar, series'
  'tide_period'
  'idTvma'
  'r2dvar'
  1.0d0

'tide_Cmin'                                        ! Input
  'minimum tidal current, ellipse minor axis'
  'meter second-1'                                 ! [m/s]
  'tide_Cmin, scalar, series'
  'tide_period'
  'idTvmi'
  'r2dvar'
  1.0d0

'tide_Pamp'                                        ! Input
  'tidal potential elevation amplitude'
  'meter'                                          ! [m]
  'SSH_Pamp, scalar, series'
  'tide_period'
  'idTpam'
  'r2dvar'
  1.0d0

'tide_Pphase'                                      ! Input
  'tidal potential elevation phase angle'
  'degrees'                                        ! [radians]
  'SSH_Pphase, scalar, series'
  'tide_period'
  'idTpph'
  'r2dvar'
  0.017453292519943295                             ! pi/180

'tide_m2amp'                                       ! Input
  'tidal velocity amplitude'
  'meter second-1'                                 ! m/s
  'tide_m2amp, scalar'
  'tide_period'
  'idM2am'
  'r2dvar'
  1.0

!
!------------------------------------------------------------------------------
!  River forcing variables.
!------------------------------------------------------------------------------
!

'river_Xposition'                                  ! Input
  'river runoff  XI-positions at RHO-points'
  'nondimensional'                                 ! [nondimensional]
  'river_Xposition, scalar'
  'river'
  'idRxpo'
  'nulvar'
  1.0d0

'river_Eposition'                                  ! Input
  'river runoff ETA-positions at RHO-points'
  'nondimensional'                                 ! [nondimensional]
  'river_Eposition'
  'river'
  'idRepo'
  'nulvar'
  1.0d0

'river_direction'                                  ! Input
  'river runoff direction'
  'nondimensional'                                 ! [nondimensional]
  'river_direction, scalar'
  'river'
  'idRdir'
  'nulvar'
  1.0d0

'river_Vshape'                                     ! Input
  'river runoff mass transport vertical profile'
  'nondimensional'                                 ! [nondimensional]
  'river_Vshape, scalar'
  'river'
  'idRvsh'
  'nulvar'
  1.0d0

'river_transport'                                  ! Input
  'river runoff mass transport'
  'meter3 second-1'                                ! [m3/s]
  'river_transport, scalar, series'
  'river_time'
  'idRtra'
  'nulvar'
  1.0d0

'river_flag'                                       ! Input
  'river runoff flag'
  'nondimensional'                                 ! [nondimensional]
  'river_flag, scalar'
  'river'
  'idRflg'
  'nulvar'
  1.0d0

'river_temp'                                       ! Input
  'river runoff potential temperature'
  'Celsius'                                        ! [Celsius]
  'river_temp, scalar, series'
  'river_time'
  'idRtrc(itemp)'
  'nulvar'
  1.0d0

'river_salt'                                       ! Input
  'river runoff salinity'
  'nondimensional'                                 ! [PSU]
  'river_salt, scalar, series'
  'river_time'
  'idRtrc(isalt)'
  'nulvar'
  1.0d0

'river_dye_'                                       ! Input
  'river runoff dye concentration'
  'kilogram meter-3'                               ! [kg/m3]
  'river_dye, scalar, series'
  'river_time'
  'idRtrc(inert(i))'
  'nulvar'
  1.0d0

'river_mud_'                                       ! Input
  'river runoff suspended sediment concentration'
  'kilogram meter-3'                               ! [kg/m3]
  'river_sediment, scalar, series'
  'river_time'
  'idRtrc(idmud)'
  'nulvar'
  1.0d0

'river_sand_'                                      ! Input
  'river runoff suspended sediment concentration'
  'kilogram meter-3'                               ! [kg/m3]
  'river_sediment, scalar, series'
  'river_time'
  'idRtrc(isand)'
  'nulvar'
  1.0d0

'river_nutrient'                                   ! Input
  'river runoff nutrient'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_nutrient, scalar, series'
  'river_time'
  'idRtrc(iNutr)'
  'nulvar'
  1.0d0

'river_NO3'                                        ! Input
  'river runoff NO3'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_NO3, scalar, series'
  'river_time'
  'idRtrc(iNO3_)'
  'nulvar'
  1.0d0

'river_nanophytoplankton'                          ! Input
  'river runoff nanophytoplankton'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_nanophytoplankton, scalar, series'
  'river_time'
  'idRtrc(iSphy)'
  'nulvar'
  1.0d0

'river_diatom'                                     ! Input
  'river runoff diatom'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_diatom, scalar, series'
  'river_time'
  'idRtrc(iLphy)'
  'nulvar'
  1.0d0

'river_microzooplankton'                           ! Input
  'river runoff microzooplankton'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_microzooplankton, scalar, series'
  'river_time'
  'idRtrc(iSzoo)'
  'nulvar'
  1.0d0

'river_mesozooplankton'                            ! Input
  'river runoff mesozooplankton'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_mesozooplankton, scalar, series'
  'river_time'
  'idRtrc(iLzoo)'
  'nulvar'
  1.0d0

'river_SiOH4'                                      ! Input
  'river runoff SiOH4'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'river_SiOH4, scalar, series'
  'river_time'
  'idRtrc(iSiOH)'
  'nulvar'
  1.0d0

'river_opal'                                       ! Input
  'river runoff opal'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'river_opal, scalar, series'
  'river_time'
  'idRtrc(iopal)'
  'nulvar'
  1.0d0

'river_PO4'                                        ! Input
  'river runoff PO4'
  'millimole_P meter-3'                            ! [millimole/m3]
  'river_PO4, scalar, series'
  'river_time'
  'idRtrc(iPO4_)'
  'nulvar'
  1.0d0

'river_NH4'                                        ! Input
  'river runoff NH4'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_NH4, scalar, series'
  'river_time'
  'idRtrc(iNH4_)'
  'nulvar'
  1.0d0

'river_DON'                                        ! Input
  'river runoff dissolved organic nitrogen'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_DON, scalar, series'
  'river_time'
  'idRtrc(iDON_)'
  'nulvar'
  1.0d0

'river_PON'                                        ! Input
  'river runoff particulate organic nitrogen'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_DON, scalar, series'
  'river_time'
  'idRtrc(iDON_)'
  'nulvar'
  1.0d0

'river_detritus'                                   ! Input
  'river runoff detritus'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_detritus, scalar, series'
  'river_time'
  'idRtrc(iSDet)'
  'nulvar'
  1.0d0

'river_SDeN'                                       ! Input
  'river runoff SDeN'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_SDeN, scalar, series'
  'river_time'
  'idRtrc(iSDeN)'
  'nulvar'
  1.0d0

'river_SDeC'                                       ! Input
  'river runoff SDeC'
  'millimole meter-3'                              ! [millimole/m3]
  'river_SDeC, scalar, series'
  'river_time'
  'idRtrc(iSDeC)'
  'nulvar'
  1.0d0

'river_TIC'                                        ! Input
  'river runoff TIC'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'river_TIC, scalar, series'
  'river_time'
  'idRtrc(iTIC_)'
  'nulvar'
  1.0d0

'river_alkalinity'                                 ! Input
  'river runoff alkalinity'
  'milliequivalents meter-3'                       ! [milliequivalens/m3]
  'river_alkalinity, scalar, series'
  'river_time'
  'idRtrc(iTAlk)'
  'nulvar'
  1.0d0

'river_Oxyg'                                       ! Input
  'river runoff oxygen'
  'millimole_oxygen meter-3'                       ! [millimole/m3]
  'river_Oxyg, scalar, series'
  'river_time'
  'idRtrc(iOxyg)'
  'nulvar'
  1.0d0

!--------------------------UMAINE variables
'river_NO3_'                                        ! Input
  'river runoff NO3'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_NO3, scalar, series'
  'river_time'
  'idRtrc(iNO3_)'
  'nulvar'
  1.0d0

'river_Nh4_'                                        ! Input
  'river runoff NH4'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'river_NH4, scalar, series'
  'river_time'
  'idRtrc(iNH4_)'
  'nulvar'
  1.0d0

'river_SiOH'                                        ! Input
  'river runoff SiOH'
  'millimole_silica meter-3'                     ! [millimole/m3]
  'river_SiOH, scalar, series'
  'river_time'
  'idRtrc(iSiOH)'
  'nulvar'
  1.0d0

'river_PO4_'                                        ! Input
  'river runoff PO4'
  'millimole_phosphorus meter-3'                     ! [millimole/m3]
  'river_PO4, scalar, series'
  'river_time'
  'idRtrc(iPO4_)'
  'nulvar'
  1.0d0

'river_S1_N'                                        ! Input
  'river runoff S1_N'
  'millimole_N meter-3'                     ! [millimole/m3]
  'river_S1_N, scalar, series'
  'river_time'
  'idRtrc(iS1_N)'
  'nulvar'
  1.0d0

'river_S1_C'                                        ! Input
  'river runoff S1_C'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_S1_C, scalar, series'
  'river_time'
  'idRtrc(iS1_C)'
  'nulvar'
  1.0d0

'river_S1CH'                                        ! Input
  'river runoff S1CH'
  'milligram_CHL meter-3'                     ! [millimole/m3]
  'river_S1CH, scalar, series'
  'river_time'
  'idRtrc(iS1CH)'
  'nulvar'
  1.0d0

'river_S2_N'                                        ! Input
  'river runoff S2_N'
  'millimole_N meter-3'                     ! [millimole/m3]
  'river_S2_N, scalar, series'
  'river_time'
  'idRtrc(iS2_N)'
  'nulvar'
  1.0d0

'river_S2_C'                                        ! Input
  'river runoff S2_C'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_S2_C, scalar, series'
  'river_time'
  'idRtrc(iS2_C)'
  'nulvar'
  1.0d0

'river_S2CH'                                        ! Input
  'river runoff S2CH'
  'milligram_CHL meter-3'                     ! [millimole/m3]
  'river_S2CH, scalar, series'
  'river_time'
  'idRtrc(iS2CH)'
  'nulvar'
  1.0d0

'river_S3_N'                                        ! Input
  'river runoff S3_N'
  'millimole_N meter-3'                     ! [millimole/m3]
  'river_S3_N, scalar, series'
  'river_time'
  'idRtrc(iS3_N)'
  'nulvar'
  1.0d0

'river_S3_C'                                        ! Input
  'river runoff S3_C'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_S3_C, scalar, series'
  'river_time'
  'idRtrc(iS3_C)'
  'nulvar'
  1.0d0

'river_S3CH'                                        ! Input
  'river runoff S3CH'
  'milligram_CHL meter-3'                     ! [millimole/m3]
  'river_S3CH, scalar, series'
  'river_time'
  'idRtrc(iS3CH)'
  'nulvar'
  1.0d0

'river_Z1_N'                                        ! Input
  'river runoff Z1_N'
  'millimole_N meter-3'                     ! [millimole/m3]
  'river_Z1_N, scalar, series'
  'river_time'
  'idRtrc(iZ1_N)'
  'nulvar'
  1.0d0

'river_Z1_C'                                        ! Input
  'river runoff Z1_C'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_Z1_C, scalar, series'
  'river_time'
  'idRtrc(iZ1_C)'
  'nulvar'
  1.0d0

'river_Z2_C'                                        ! Input
  'river runoff Z2_C'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_Z2_C, scalar, series'
  'river_time'
  'idRtrc(iZ2_C)'
  'nulvar'
  1.0d0

'river_Z2_N'                                        ! Input
  'river runoff Z2_N'
  'millimole_N meter-3'                     ! [millimole/m3]
  'river_Z2_N, scalar, series'
  'river_time'
  'idRtrc(iZ2_N)'
  'nulvar'
  1.0d0

'river_BAC_'                                        ! Input
  'river runoff BAC'
  'millimole_N meter-3'                     ! [millimole/m3]
  'river_BAC, scalar, series'
  'river_time'
  'idRtrc(iBAC_)'
  'nulvar'
  1.0d0

'river_DD_N'                                        ! Input
  'river runoff DD_N'
  'millimole_N meter-3'                     ! [millimole/m3]
  'river_DD_N, scalar, series'
  'river_time'
  'idRtrc(iDD_N)'
  'nulvar'
  1.0d0

'river_DD_C'                                        ! Input
  'river runoff DD_C'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_DD_C, scalar, series'
  'river_time'
  'idRtrc(iDD_C)'
  'nulvar'
  1.0d0

'river_DDCA'                                        ! Input
  'river runoff DDCA'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_DDCA, scalar, series'
  'river_time'
  'idRtrc(iDDCA)'
  'nulvar'
  1.0d0

'river_DDSi'                                        ! Input
  'river runoff DDSi'
  'millimole_Si meter-3'                     ! [millimole/m3]
  'river_DDSi, scalar, series'
  'river_time'
  'idRtrc(iDDSi)'
  'nulvar'
  1.0d0

'river_LDON'                                        ! Input
  'river runoff LDON'
  'millimole_N meter-3'                     ! [millimole/m3]
  'river_LDON, scalar, series'
  'river_time'
  'idRtrc(iLDON)'
  'nulvar'
  1.0d0

'river_LDOC'                                        ! Input
  'river runoff LDOC'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_LDOC, scalar, series'
  'river_time'
  'idRtrc(iLDOC)'
  'nulvar'
  1.0d0

'river_SDON'                                        ! Input
  'river runoff SDON'
  'millimole_N meter-3'                     ! [millimole/m3]
  'river_SDON, scalar, series'
  'river_time'
  'idRtrc(iSDON)'
  'nulvar'
  1.0d0

'river_SDOC'                                        ! Input
  'river runoff SDOC'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_SDOC, scalar, series'
  'river_time'
  'idRtrc(iSDOC)'
  'nulvar'
  1.0d0

'river_CLDC'                                        ! Input
  'river runoff CLDC'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_CLDC, scalar, series'
  'river_time'
  'idRtrc(iCLDC)'
  'nulvar'
  1.0d0

'river_CSDC'                                        ! Input
  'river runoff CSDC'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_CSDC, scalar, series'
  'river_time'
  'idRtrc(iCSDC)'
  'nulvar'
  1.0d0

'river_Oxyg'                                        ! Input
  'river runoff Oxyg'
  'millimole meter-3'                     ! [millimole/m3]
  'river_Oxyg, scalar, series'
  'river_time'
  'idRtrc(iOxyg)'
  'nulvar'
  1.0d0

  'river_Talk'                                        ! Input
  'river runoff Talk'
  'millimole meter-3'                     ! [millimole/m3]
  'river_Talk, scalar, series'
  'river_time'
  'idRtrc(iTAlk)'
  'nulvar'
  1.0d0

  'river_TIC'                                        ! Input
  'river runoff TIC'
  'millimole_C meter-3'                     ! [millimole/m3]
  'river_TIC, scalar, series'
  'river_time'
  'idRtrc(iTIC_)'
  'nulvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Boundary layer variables.
!------------------------------------------------------------------------------
!

'Hsbl'                                             ! Input/Output
  'depth of oceanic surface boundary layer'
  'meter'                                          ! [m]
  'Hsbl, scalar, series'
  'ocean_time'
  'idHsbl'
  'r2dvar'
  1.0d0

'Hbbl'                                             ! Input/Output
  'depth of oceanic bottom boundary layer'
  'meter'                                          ! [m]
  'Hbbl, scalar, series'
  'ocean_time'
  'idHbbl'
  'r2dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Climatotogy variables.
!------------------------------------------------------------------------------
!

'zeta'                                             ! Input
  'sea surface height climatology'
  'meter'                                          ! [m]
  'zeta, scalar, series'
  'ocean_time'
  'idSSHc'
  'r2dvar'
  1.0d0

'ubar'                                             ! Input
  'vertically integrated u-momentum climatology'
  'meter second-1'                                 ! [m/s]
  'ubar-velocity, scalar, series'
  'ocean_time'
  'idUbcl'
  'u2dvar'
  1.0d0

'vbar'                                             ! Input
  'vertically integrated v-momentum climatology'
  'meter second-1'                                 ! [m/s]
  'vbar-velocity, scalar, series'
  'ocean_time'
  'idVbcl'
  'v2dvar'
  1.0d0

'u'                                                ! Input
  'u-momentum component climatology'
  'meter second-1'                                 ! [m/s]
  'u-velocity, scalar, series'
  'ocean_time'
  'idUclm'
  'u3dvar'
  1.0d0

'v'                                                ! Input
  'v-momentum component climatology'
  'meter second-1'                                 ! [m/s]
  'v-velocity, scalar, series'
  'ocean_time'
  'idVclm'
  'v3dvar'
  1.0d0

'omega'
  'S-coordinate vertical momentum climatology'
  'meter3 second-1'
  'omega, scalar, series'
  'ocean_time'
  'idOclm'
  'w3dvar'
  1.0d0

'SSS'                                              ! Input
  'sea surface salinity climatology'
  'nondimensional'                                 ! [PSU]
  'SSS, scalar, series'
  'sss_time'
  'idSSSc'
  'r2dvar'
  1.0d0

'SST'                                              ! Input
  'sea surface temperature climatology'
  'Celsius'                                        ! [Celsius]
  'SST, scalar, series'
  'sst_time'
  'idSSTc'
  'r2dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Adjoint sensitivy variables.  The units below are the same as the state.
!  However, the input functional can have squared, energy, nondimensional,
!  or other units.
!------------------------------------------------------------------------------
!

'zeta'                                             ! Input
  'sea surface height functional'
  'meter'                                          ! [m]
  'zeta, scalar, series'
  'zeta_time'
  'idZads'
  'r2dvar'
  1.0d0

'ubar'                                             ! Input
  'vertically integrated u-momentum functional'
  'meter second-1'                                 ! [m/s]
  'ubar-velocity, scalar, series'
  'v2d_time'
  'idUbas'
  'u2dvar'
  1.0d0

'vbar'                                             ! Input
  'vertically integrated v-momentum functional'
  'meter second-1'                                 ! [m/s]
  'vbar-velocity, scalar, series'
  'v2d_time'
  'idVbas'
  'v2dvar'
  1.0d0

'u'                                                ! Input
  'u-momentum component functional'
  'meter second-1'                                 ! [m/s]
  'u-velocity, scalar, series'
  'v3d_time'
  'idUads'
  'u3dvar'
  1.0d0

'v'                                                ! Input
  'v-momentum component functional'
  'meter second-1'                                 ! [m/s]
  'v-velocity, scalar, series'
  'v3d_time'
  'idVads'
  'v3dvar'
  1.0d0

'temp'                                             ! Input
  'potential temperature functional'
  'Celsius'                                        ! [Celsius]
  'temp_time, scalar, series'
  'temp_time'
  'idTads(itemp)'
  'r3dvar'
  1.0d0

'salt'                                             ! Input
  'salinity functional'
  'nondimensional'                                 ! [PSU]
  'salt_time, scalar, series'
  'salt_time'
  'idTads(isalt)'
  'r3dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Tangent linear model impulse forcing variables.
!
!  Notice that additional metadata is need here because with need different
!  indices to avoid conflict with the forward basic state variables. They
!  are used at the same time. The passive tracers (NAT+1:NT) are assigned
!  internally to the same values as the basic state.
!------------------------------------------------------------------------------
!

'zeta'                                             ! Input
  'free-surface impulse forcing'
  'meter'                                          ! [m]
  'zeta, scalar, series'
  'zeta_time'
  'idZtlf'
  'r2dvar'
  1.0d0

'ubar'                                             ! Input
  'vertically integrated u-momentum impulse forcing'
  'meter second-1'                                 ! [m/s]
  'ubar-velocity, scalar, series'
  'v2d_time'
  'idUbtf'
  'u2dvar'
  1.0d0

'vbar'                                             ! Input
  'vertically integrated v-momentum impulse forcing'
  'meter second-1'                                 ! [m/s]
  'vbar-velocity, scalar, series'
  'v2d_time'
  'idVbtf'
  'v2dvar'
  1.0d0

'u'                                                ! Input
  'u-momentum component impulse forcing'
  'meter second-1'                                 ! [m/s]
  'u-velocity, scalar, series'
  'v3d_time'
  'idUtlf'
  'u3dvar'
  1.0d0

'v'                                                ! Input
  'v-momentum component impulse forcing'
  'meter second-1'                                 ! [m/s]
  'v-velocity, scalar, series'
  'v3d_time'
  'idVtlf'
  'v3dvar'
  1.0d0

'temp'                                             ! Input
  'potential temperature impulse forcing'
  'Celsius'                                        ! [Celsius]
  'temp_time, scalar, series'
  'temp_time'
  'idTtlf(itemp)'
  'r3dvar'
  1.0d0

'salt'                                             ! Input
  'salinity functional'
  'nondimensional'                                 ! [PSU]
  'salt_time, scalar, series'
  'salt_time'
  'idTtlf(isalt)'
  'r3dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Observation variables: Optimal Interpolation and Nudging.
!------------------------------------------------------------------------------
!

'SSH'                                              ! Input
  'sea surface height observations'
  'meter'                                          ! [m]
  'SSHobs, scalar, series'
  'SSHobs_time'
  'idSSHo'
  'r2dvar'
  1.0d0

'SSH_err'                                          ! Input
  'sea surface height error variance'
  'nondimensional'                                 ! [nondimenisional]
  'SSH_err, scalar, series'
  'SSHobs_time'
  'idSSHe'
  'r2dvar'
  1.0d0

'u'                                                ! Input
  'u-momentum component observations'
  'meter second-1'                                 ! [m/s]
  'Uobs, scalar, series'
  'Vobs_time'
  'idUobs'
  'u3dvar'
  1.0d0

'v'                                                ! Input
  'v-momentum component observations'
  'meter second-1'                                 ! [m/s]
  'Vobs, scalar, series'
  'Vobs_time'
  'idVobs'
  'v3dvar'
  1.0d0

'Evel'                                             ! Input
  'horizontal velocity error variance'
  'nondimensional'                                 ! [nondimensional]
  'Evel, scalar, series'
  'Vobs_time'
  'idUVer'
  'r3dvar'
  1.0d0

'Usur'                                             ! Input
  'surface u-momentum component observations'
  'meter second-1'                                 ! [m/s]
  'Usur, scalar, series'
  'Vsur_time'
  'idUsur'
  'u2dvar'
  1.0d0

'Vsur'                                             ! Input
  'surface v-momentum component observations'
  'meter second-1'                                 ! [m/s]
  'Vsur, scalar, series'
  'Vsur_time'
  'idVsur'
  'v2dvar'
  1.0d0

'Esur'                                             ! Input
  'surface velocity error variance'
  'nondimensional'                                 ! [nondimensional]
  'Esur, scalar, series'
  'Vsur_time'
  'idUVse'
  'r2dvar'
  1.0d0

'SST'                                              ! Input
  'sea surface temperature observations'
  'Celsius'                                        ! [Celsius]
  'SSTobs, scalar, series'
  'SSTobs_time'
  'idSSTo'
  'r2dvar'
  1.0d0

'SST_err'                                          ! Input
  'sea surface temperature error variance'
  'nondimensional'                                 ! [nondimensional]
  'SST_err, scalar, series'
  'SSTobs_time'
  'idSSTe'
  'r2dvar'
  1.0d0

'temp'                                             ! Input
  'potential temperature observations'
  'Celsius'                                        ! [nondimensional]
  'Tobs, scalar, series'
  'Tobs_time'
  'idTobs(itemp)'
  'r3dvar'
  1.0d0

'temp_err'                                         ! Input
  'potential temperature error variance'
  'nondimensional'                                 ! [nondimensional]
  'temp_err, scalar, series'
  'Tobs_time'
  'idTerr(itemp)'
  'r3dvar'
  1.0d0

'salt'                                             ! Input
  'salinity observations'
  'nondimensional'                                 ! [PSU]
  'Sobs, scalar, series'
  'Sobs_time'
  'idTobs(isalt)'
  'r3dvar'
  1.0d0

'salt_err'                                         ! Input
  'salinity error variance'
  'nondimensional'                                 ! [nondimensional]
  'salt_err, scalar, series'
  'Sobs_time'
  'idTerr(isalt)'
  'r3dvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Inert passive tracer variables.
!------------------------------------------------------------------------------
!
!  Only one entry should be specified here.  This entry contain all the
!  information necessary to process all NPT fields.  The variable name is
!  arbitrary and you may choose any unique name for Vinfo(1) that is not
!  already defined. The identifier for Vinfo(6) must be "idTvar(inert(i))".
!

'dye_'                                             ! Input/Output
  'dye concentration'
  'kilogram meter-3'                               ! [kg/m3]
  'dye_, scalar, series'
  'ocean_time'
  'idTvar(inert(i))'
  'r3dvar'
  1.0d0

'dye_west_'                                        ! Input
  'dye concentration western boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'dye_west_, scalar, series'
  'dye_time'
  'idTbry(iwest,inert(i))'
  'nulvar'
  1.0d0

'dye_east_'                                        ! Input
  'dye concentration eastern boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'dye_east_, scalar, series'
  'dye_time'
  'idTbry(ieast,inert(i))'
  'nulvar'
  1.0d0

'dye_south_'                                       ! Input
  'dye concentration southern boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'dye_south_, scalar, series'
  'dye_time'
  'idTbry(isouth,inert(i))'
  'nulvar'
  1.0d0

'dye_north_'                                       ! Input
  'dye concentration northern boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'dye_north_, scalar, series'
  'dye_time'
  'idTbry(inorth,inert(i))'
  'nulvar'
  1.0d0

!
!------------------------------------------------------------------------------
!  Sediment variables:  The actual variable name for some of the sediment
!                       tracers are set-up internally in "mod_ncparam.F".
!  It depends on the mumber of cohesive and noncohesive sediment tracers
!  for a particular application. For example if you have 1 cohesive and
!  2 noncohesive sediment classes, the variable names are:
!
!    "mud_"          becomes   "mud_01"
!    "mud_west_"     becomes   "mud_west_01"
!    "mud_east_"     becomes   "mud_east_01"
!    "mud_south_"    becomes   "mud_south_01"
!    "mud_north_"    becomes   "mud_north_01"
!
!    "sand_"         becomes   "sand_01",       "sand_02"
!    "sand_west"     becomes   "sand_west_01",  "sand_west_02"
!    "sand_east"     becomes   "sand_east_01",  "sand_east_02"
!    "sand_south"    becomes   "sand_south_01", "sand_south_02"
!    "sand_north"    becomes   "sand_north_01", "sand_north_02"
!------------------------------------------------------------------------------
!

'mud_'                                             ! Input/Output
  'suspended cohesive sediment'
  'kilogram meter-3'                               ! [kg/m3]
  'mud_, scalar, series'
  'ocean_time'
  'idTvar(idmud(i))'
  'r3dvar'
  1.0d0

'mud_west_'                                        ! Input
  'suspended cohesive sediment western boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'mud_west_, scalar, series'
  'mud_time'
  'idTbry(iwest,idmud(i))'
  'r3dvar'
  1.0d0

'mud_east_'                                        ! Input
  'suspended cohesive sediment eastern boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'mud_east_, scalar, series'
  'mud_time'
  'idTbry(ieast,idmud(i))'
  'nulvar'
  1.0d0

'mud_south_'                                       ! Input
  'suspended cohesive sediment southern boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'mud_south_, scalar, series'
  'mud_time'
  'idTbry(isouth,idmud(i))'
  'nulvar'
  1.0d0

'mud_north_'                                       ! Input
  'suspended cohesive sediment northern boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'mud_north_, scalar, series'
  'mud_time'
  'idTbry(inorth,idmud(i))'
  'nulvar'
  1.0d0

'sand_'                                            ! Input/Output
  'suspended noncohesive sediment'
  'kilogram meter-3'                               ! [kg/m3]
  'sand_, scalar, series'
  'ocean_time'
  'idTvar(isand(i))'
  'r3dvar'
  1.0d0

'sand_west_'                                       ! Input
  'suspended noncohesive sediment western boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'sand_west_, scalar, series'
  'sand_time'
  'idTbry(iwest,isand(i))'
  'nulvar'
  1.0d0

'sand_east_'                                       ! Input
  'suspended noncohesive sediment eastern boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'sand_east_, scalar, series'
  'sand_time'
  'idTbry(ieast,isand(i))'
  'nulvar'
  1.0d0

'sand_south_'                                      ! Input
  'suspended noncohesive sediment southern boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'sand_south_, scalar, series'
  'sand_time'
  'idTbry(isouth,isand(i))'
  'nulvar'
  1.0d0

'sand_north_'                                      ! Input
  'suspended noncohesive sediment northern boundary condition'
  'kilogram meter-3'                               ! [kg/m3]
  'sand_north_, scalar, series'
  'sand_time'
  'idTbry(inorth,isand(i))'
  'nulvar'
  1.0d0

'mudfrac_'                                         ! Input/Output
  'cohesive sediment fraction'
  'nondimensional'                                 ! [nondimensional]
  'mudfrac_, scalar, series'
  'ocean_time'
  'idfrac'
  'b3dvar'
  1.0d0

'sandfrac_'                                        ! Input/Output
  'noncohesive sediment fraction'
  'nondimensional'                                 ! [nondimensional]
  'sandfrac_, scalar, series'
  'ocean_time'
  'idfrac'
  'b3dvar'
  1.0d0

'mudmass_'                                         ! Input/Output
  'cohesive sediment mass'
  'kilogram meter-2'                               ! [kg/m2]
  'mudmass_, scalar, series'
  'ocean_time'
  'idBmas'
  'b3dvar'
  1.0d0

'sandmass_'                                        ! Input/Output
  'noncohesive sediment mass'
  'kilogram meter-2'                               ! [kg/m2]
  'sandmass_, scalar, series'
  'ocean_time'
  'idBmas'
  'b3dvar'
  1.0d0

'bedload_Usand_'                                   ! Input/Output
  'bed load flux of sand in U-direction'
  'kilogram meter-1 s-1'                           ! [kg/m/s]
  'bedload_Usand_, scalar, series'
  'ocean_time'
  'idUbld'
  'u2dvar'
  1.0d0

'bedload_Umud_'                                    ! Input/Output
  'bed load flux of mud in U-direction'
  'kilogram meter-1 s-1'                           ! [kg/m/s]
  'bedload_Umud_, scalar, series'
  'ocean_time'
  'idUbld'
  'u2dvar'
  1.0d0

'bedload_Vsand_'                                   ! Input/Output
  'bed load flux of sand in V-direction'
  'kilogram meter-1 s-1'                           ! [kg/m/s]
  'bedload_Vsand_, scalar, series'
  'ocean_time'
  'idVbld'
  'v2dvar'
  1.0d0

'bedload_Vmud_'                                    ! Input/Output
  'bed load flux of mud in V-direction'
  'kilogram meter-1 s-1'                           ! [kg/m/s]
  'bedload_Vmud_, scalar, series'
  'ocean_time'
  'idVbld'
  'v2dvar'
  1.0d0

'ero_flux'                                         ! Output
  'erosion flux'
  'kilogram meter-2'                               ! [kg/m2]
  'ero_flux, scalar, series'
  'ocean_time'
  'idBott(idefx)'
  'r2dvar'
  1.0d0

'dep_net'                                          ! Output
  'erosion or deposition'
  'meter'                                          ! [m]
  'dep_net, scalar, series'
  'ocean_time'
  'idBott(idnet)'
  'r2dvar'
  1.0d0

'tau_crit_off'                                     ! Output
  'Offset for critical shear stress equation'
  'kg/Pa/s'                                        ! [kg/Pa/s]
  'tau_crit_off, scalar, series'
  'ocean_time'
  'idBott(idoff)'
  'r2dvar'
  1.0d0

'tau_crit_slope'                                   ! Output
  'Slope for critical shear stress equation'
  'Pa/s'                                           ! [Pa/s]
  'tau_crit_slope, scalar, series'
  'ocean_time'
  'idBott(idslp)'
  'r2dvar'
  1.0d0

'tau_crit_time'                                    ! Output
  'Nudging time for critical shear stress equation'
  's-1'                                            ! [s-1]
  'tau_crit_time, scalar, series'
  'ocean_time'
  'idBott(idtim)'
  'r2dvar'
  1.0d0

'bed_thickness'                                    ! Input/Output
  'sediment bed layer thickness'
  'meter'                                          ! [m]
  'bed_thickness, scalar, series'
  'ocean_time'
  'idSbed(ithck)'
  'b3dvar'
  1.0d0

'bed_age'                                          ! Input/Output
  'sediment layer age'
  'seconds'                                        ! [s]
  'bed_age, scalar, series'
  'ocean_time'
  'idSbed(iaged)'
  'b3dvar'
  1.0d0

'bed_porosity'                                     ! Input/Output
  'sediment layer porosity'
  'nondimensional'                                 ! [nondimensional]
  'bed_porosity, scalar, series'
  'ocean_time'
  'idSbed(iporo)'
  'b3dvar'
  1.0d0

'bed_biodiff'                                      ! Input/Output
  'biodiffusivity at bottom of each layer'
  'meter2 second-1'                                ! [m2/s]
  'bed_biodiff, scalar, series'
  'ocean_time'
  'idSbed(idiff)'
  'b3dvar'
  1.0d0

'bed_tau_crit'                                     ! Input/Output
  'tau critical in each layer'
  'Newtons meter-2'                                ! [N/m2]
  'bed_tau_crit, scalar, series'
  'ocean_time'
  'idSbed(ibtcr)'
  'b3dvar'
  1.0d0

'grain_diameter'                                   ! Input/Output
  'sediment median grain diameter size'
  'meter'                                          ! [m]
  'grain_diameter, scalar, series'
  'ocean_time'
  'idBott(isd50)'
  'r2dvar'
  1.0d0

'grain_density'                                    ! Input/Output
  'sediment median grain density'
  'kilogram meter3'                                ! [kg/m3]
  'grain_density, scalar, series'
  'ocean_time'
  'idBott(idens)'
  'r2dvar'
  1.0d0

'settling_vel'                                     ! Input/Output
  'sediment median grain settling velocity'
  'meter second-1'                                 ! [m/s]
  'settling_vel, scalar, series'
  'ocean_time'
  'idBott(iwsed)'
  'r2dvar'
  1.0d0

'erosion_stress'                                   ! Input/Output
  'sediment median critical erosion stress'
  'newton meter-2'                                 ! [N/m2]
  'erosion_stress, scalar, series'
  'ocean_time'
  'idBott(itauc)'
  'r2dvar'
  1.0d0

'ripple_length'                                    ! Input/Output
  'bottom ripple length'
  'meter'                                          ! [m]
  'ripple_length, scalar, series'
  'ocean_time'
  'idBott(irlen)'
  'r2dvar'
  1.0d0

'ripple_height'                                    ! Input/Output
  'bottom ripple height'
  'meter'                                          ! [m]
  'ripple_height, scalar, series'
  'ocean_time'
  'idBott(irhgt)'
  'r2dvar'
  1.0d0

'bed_wave_amp'                                     ! Input/Output
  'bed wave excursion amplitude'
  'meter'                                          ! [m]
  'bed_wave_amp, scalar, series'
  'ocean_time'
  'idBott(ibwav)'
  'r2dvar'
  1.0d0

'rdrag'                                            ! Input/Output
  'linear bottom drag coefficient'
  'meter second-1'                                 ! [m/s]
  'rdrag, scalar, series'
  'ocean_time'
  'idragL'
  'r2dvar'
  1.0d0

'rdrag2'                                           ! Input/Output
  'quadratic bottom drag coefficient'
  'nondimensional'                                 ! [nondimensional]
  'rdrag2, scalar, series'
  'ocean_time'
  'idragQ'
  'r2dvar'
  1.0d0

'wavedrag'                                         ! Input/Output
  'linear bottom drag coefficient'
  'meter second-1'                                 ! [m/s]
  'wavedrag, scalar, series'
  'ocean_time'
  'idragW'
  'r2dvar'
  1.0d0

'ZoBot'                                            ! Input/Output
  'time invariant, bottom roughness length'
  'meter'                                          ! [m]
  'ZoBot, scalar, series'
  'ocean_time'
  'idZoBL'
  'r2dvar'
  1.0d0

'Zo_def'                                           ! Input/Output
  'default bottom roughness length'
  'meter'                                          ! [m]
  'Zo_def, scalar, series'
  'ocean_time'
  'idBott(izdef)'
  'r2dvar'
  1.0d0

'Zo_app'                                           ! Input/Output
  'apparent bottom roughness length'
  'meter'                                          ! [m]
  'Zo_app, scalar, series'
  'ocean_time'
  'idBott(izapp)'
  'r2dvar'
  1.0d0

'Zo_Nik'                                           ! Input/Output
  'Nikaradse bottom roughness length'
  'meter'                                          ! [m]
  'Zo_Nik, scalar, series'
  'ocean_time'
  'idBott(izNik)'
  'r2dvar'
  1.0d0

'Zo_bio'                                           ! Input/Output
  'biological bottom roughness length'
  'meter'                                          ! [m]
  'Zo_bio, scalar, series'
  'ocean_time'
  'idBott(izbio)'
  'r2dvar'
  1.0d0

'Zo_bedform'                                       ! Input/Output
  'bed form bottom roughness length'
  'meter'                                          ! [m]
  'Zo_bedform, scalar, series'
  'ocean_time'
  'idBott(izbfm)'
  'r2dvar'
  1.0d0

'Zo_bedload'                                       ! Input/Output
  'bed load bottom roughness length'
  'meter'                                          ! [m]
  'Zo_bedload, scalar, series'
  'ocean_time'
  'idBott(izbld)'
  'r2dvar'
  1.0d0

'Zo_wbl'                                           ! Input/Output
  'wave bottom boundary layer roughness length'
  'meter'                                          ! [m]
  'Zo_wbl, scalar, series'
  'ocean_time'
  'idBott(izwbl)'
  'r2dvar'
  1.0d0

'active_layer_thickness'                           ! Input/Output
  'active layer thickness'
  'meter'                                          ! [m]
  'active_layer_thickness, scalar, series'
  'ocean_time'
  'idBott(iactv)'
  'r2dvar'
  1.0d0

'cohesive_behavior'                                ! Output
  'cohesive behavior'
  '1'                                              ! []
  'cohesive_behavior, scalar, series'
  'ocean_time'
  'idBott(idprp)'
  'r2dvar'
  1.0d0

'saltation'                                        ! Input/Output
  'bottom saltation height'
  'meter'                                          ! [m]
  'saltation, scalar, series'
  'ocean_time'
  'idBott(ishgt)'
  'r2dvar'
  1.0d0

'dmix_offset'                                      ! Input/Output
  'dmix erodibility profile offset'
  'meter'                                          ! [m]
  'dmix_offset, scalar, series'
  'ocean_time'
  'idBott(idoff)'
  'r2dvar'
  1.0d0

'dmix_slope'                                      ! Input/Output
  'dmix erodibility profile slope'
  '_'                                             ! [_]
  'dmix_slope, scalar, series'
  'ocean_time'
  'idBott(idslp)'
  'r2dvar'
  1.0d0

'dmix_time'                                        ! Input/Output
  'dmix erodibility profile time scale'
  'seconds'                                        ! [s]
  'dmix_time, scalar, series'
  'ocean_time'
  'idBott(idtim)'
  'r2dvar'
  1.0d0

'tau_crit_slope_c'                                 ! Input/Output
  'tau_crit_slope_c'
  'N/m2'                                           ! [N/m2]
  'tau_crit_slope_c, scalar, series'
  'ocean_time'
  'idBott(itcsc)'
  'r2dvar'
  1.0d0

'tau_crit_offset_c'                                ! Input/Output
  'tau_crit_offset_c'
  'N/m2'                                           ! [N/m2]
  'tau_crit_offset_c, scalar, series'
  'ocean_time'
  'idBott(itcoc)'
  'r2dvar'
  1.0d0

'tau_crit_rtime'                                   ! Input/Output
  'tau_crit_rtime'
  's'                                              ! [s]
  'tau_crit_rtime, scalar, series'
  'ocean_time'
  'idBott(itcrt)'
  'r2dvar'
  1.0d0

'Db_max'                                           ! Input/Output
  'Db_max'
  '?'                                              ! [?]
  'Db_max, scalar, series'
  'ocean_time'
  'idBott(idbmx)'
  'r2dvar'
  1.0d0

'Db_mm'                                            ! Input/Output
  'Db_mm'
  '?'                                              ! [?]
  'Db_mm, scalar, series'
  'ocean_time'
  'idBott(idbmm)'
  'r2dvar'
  1.0d0

'Db_zs'                                            ! Input/Output
  'Db_zs'
  '?'                                              ! [?]
  'Db_zs, scalar, series'
  'ocean_time'
  'idBott(idbzs)'
  'r2dvar'
  1.0d0

'Db_zm'                                            ! Input/Output
  'Db_zm'
  '?'                                              ! [?]
  'Db_zm, scalar, series'
  'ocean_time'
  'idBott(idbzm)'
  'r2dvar'
  1.0d0

'Db_zphi'                                          ! Input/Output
  'Db_zphi'
  '?'                                              ! [?]
  'Db_zphi, scalar, series'
  'ocean_time'
  'idBott(idbzp)'
  'r2dvar'
  1.0d0

'wet_dry_masking'                                  ! Input/Output
  'wet_dry mask'
  '1'                                              ! []
  'wet masking, scalar, series'
  'ocean_time'
  'idWmsk'
  'r2dvar'
  1.0d0

!
!---------------------------------------------------------------------
!  Define names of variables associated with Fennel et al (2006),
!  NEMURO, NPZD or UMaine ecosystem models.
!---------------------------------------------------------------------
!

'NO3'                                              ! Input/Output
  'nitrate concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3, scalar, series'
  'ocean_time'
  'idTvar(iNO3_)'
  'r3dvar'
  1.0d0

'NO3'                                              ! Input/Output
  'nitrate concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3, scalar, series'
  'ocean_time'
  'idTclm(iNO3_)'
  'r3dvar'
  1.0d0

'NO3'                                              ! Input
  'nitrate concentration functional'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3, scalar, series'
  'ocean_time'
  'idTads(iNO3_)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'ammonium concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4, scalar, series'
  'ocean_time'
  'idTvar(iNH4_)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'ammonium concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4, scalar, series'
  'ocean_time'
  'idTclm(iNH4_)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input
  'ammonium concentration functional'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4, scalar, series'
  'ocean_time'
  'idTads(iNH4_)'
  'r3dvar'
  1.0d0

'PON'                                              ! Input/Output
  'particulate organic nitrogen concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'PON, scalar, series'
  'ocean_time'
  'idTvar(iPON_)'
  'r3dvar'
  1.0d0

'PON'                                              ! Input/Output
  'particulate organic nitrogen concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'PON, scalar, series'
  'ocean_time'
  'idTclm(iPON_)'
  'r3dvar'
  1.0d0

'DON'                                              ! Input/Output
  'dissolved organic nitrogen concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'DON, scalar, series'
  'ocean_time'
  'idTvar(iDON_)'
  'r3dvar'
  1.0d0

'DON'                                              ! Input/Output
  'dissolved organic nitrogen concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'DON, scalar, series'
  'ocean_time'
  'idTclm(iDON_)'
  'r3dvar'
  1.0d0

'SiOH4'                                            ! Input/Output
  'silicate concentration'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH4, scalar, series'
  'ocean_time'
  'idTvar(iSiOH)'
  'r3dvar'
  1.0d0

'SiOH4'                                            ! Input/Output
  'silicate concentration'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH4, scalar, series'
  'ocean_time'
  'idTclm(iSiOH)'
  'r3dvar'
  1.0d0

'opal'                                             ! Input/Output
  'particulate organic silica concentration'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'silica, scalar, series'
  'ocean_time'
  'idTvar(iopal)'
  'r3dvar'
  1.0d0

'opal'                                             ! Input/Output
  'particulate organic silica concentration'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'silica, scalar, series'
  'ocean_time'
  'idTclm(iopal)'
  'r3dvar'
  1.0d0

'FeSp'                                             ! Input/Output
  'small phytop. iron concentration'
  'millimole_iron  meter-3'                        ! [millimole/m3]
  'iron, scalar, series'
  'ocean_time'
  'idTvar(iFeSp)'
  'r3dvar'
  1.0d0

'FeSp'                                             ! Input/Output
  'small phytop. iron concentration'
  'millimole_iron  meter-3'                        ! [millimole/m3]
  'iron, scalar, series'
  'ocean_time'
  'idTclm(iFeSp)'
  'r3dvar'
  1.0d0

'FeLp'                                            ! Input/Output
  'large phytop. iron concentration'
  'millimole_iron  meter-3'                       ! [millimole/m3]
  'iron, scalar, series'
  'ocean_time'
  'idTvar(iFeLp)'
  'r3dvar'
  1.0d0

'FeLp'                                            ! Input/Output
  'large phytop. iron concentration'
  'millimole_iron  meter-3'                       ! [millimole/m3]
  'iron, scalar, series'
  'ocean_time'
  'idTclm(iFeLp)'
  'r3dvar'
  1.0d0

'FeD'                                             ! Input/Output
  'disolved iron concentration'
  'millimole_iron  meter-3'                       ! [millimole/m3]
  'iron, scalar, series'
  'ocean_time'
  'idTvar(iFeD_)'
  'r3dvar'
  1.0d0

'FeD'                                             ! Input/Output
  'disolved iron concentration'
  'millimole_iron  meter-3'                       ! [millimole/m3]
  'iron, scalar, series'
  'ocean_time'
  'idTclm(iFeD_)'
  'r3dvar'
  1.0d0

'phytoplankton'                                    ! Input/Output
  'phytoplankton concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'phytoplankton, scalar, series'
  'ocean_time'
  'idTvar(iPhyt)'
  'r3dvar'
  1.0d0

'phytoplankton'                                    ! Input
  'phytoplankton concentration functional'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'phytoplankton, scalar, series'
  'ocean_time'
  'idTads(iPhyt)'
  'r3dvar'
  1.0d0

'nanophytoplankton'                                ! Input/Output
  'nanophytoplankton biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'nanophytoplankton, scalar, series'
  'ocean_time'
  'idTvar(iSphy)'
  'r3dvar'
  1.0d0

'nanophytoplankton'                                ! Input/Output
  'nanophytoplankton biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'nanophytoplankton, scalar, series'
  'ocean_time'
  'idTclm(iSphy)'
  'r3dvar'
  1.0d0

'diatom'                                           ! Input/Output
  'diatom biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'diatom, scalar, series'
  'ocean_time'
  'idTvar(iLphy)'
  'r3dvar'
  1.0d0

'diatom'                                           ! Input/Output
  'diatom biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'diatom, scalar, series'
  'ocean_time'
  'idTclm(iLphy)'
  'r3dvar'
  1.0d0

'zooplankton'                                      ! Input/Output
  'zooplankton concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'zooplankton, scalar, series'
  'ocean_time'
  'idTvar(iZoop)'
  'r3dvar'
  1.0d0

'zooplankton'                                      ! Input
  'zooplankton concentration functional'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'zooplankton, scalar, series'
  'ocean_time'
  'idTads(iZoop)'
  'r3dvar'
  1.0d0

'microzooplankton'                                 ! Input/Output
  'microzooplankton biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'microzooplankton, scalar, series'
  'ocean_time'
  'idTvar(iSzoo)'
  'r3dvar'
  1.0d0

'microzooplankton'                                 ! Input/Output
  'microzooplankton biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'microzooplankton, scalar, series'
  'ocean_time'
  'idTclm(iSzoo)'
  'r3dvar'
  1.0d0

'mesozooplankton'                                  ! Input/Output
  'mesozooplankton biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'mesozooplankton, scalar, series'
  'ocean_time'
  'idTvar(iLzoo)'
  'r3dvar'
  1.0d0

'mesozooplankton'                                  ! Input/Output
  'mesozooplankton biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'mesozooplankton, scalar, series'
  'ocean_time'
  'idTclm(iLzoo)'
  'r3dvar'
  1.0d0

'Pzooplankton'                                     ! Input/Output
  'predator-zooplankton biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'Pzooplankton, scalar, series'
  'ocean_time'
  'idTvar(iPzoo)'
  'r3dvar'
  1.0d0

'Pzooplankton'                                     ! Input/Output
  'predator-zooplankton biomass'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'Pzooplankton, scalar, series'
  'ocean_time'
  'idTclm(iPzoo)'
  'r3dvar'
  1.0d0

'chlorophyll'                                      ! Input/Output
  'chlorophyll concentration'
  'milligrams_chlorophyll meter-3'                 ! [milligrams/m3]
  'chlorophyll, scalar, series'
  'ocean_time'
  'idTvar(iChlo)'
  'r3dvar'
  1.0d0

'chlorophyll'                                      ! Input
  'chlorophyll concentration functional'
  'milligrams_chlorophyll meter-3'                 ! [milligrams/m3]
  'chlorophyll, scalar, series'
  'ocean_time'
  'idTads(iChlo)'
  'r3dvar'
  1.0d0

'detritus'                                         ! Input/Output
  'detritus concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'detritus, scalar, series'
  'ocean_time'
  'idTvar(iSDet)'
  'r3dvar'
  1.0d0

'detritus'                                         ! Input
  'detritus concentration functional'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'detritus, scalar, series'
  'ocean_time'
  'idTads(iSDet)'
  'r3dvar'
  1.0d0

'LdetritusN'                                       ! Input/Output
  'large fraction nitrogen detritus concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'LdetritusN, scalar, series'
  'ocean_time'
  'idTvar(iLDeN)'
  'r3dvar'
  1.0d0

'LdetritusN'                                       ! Input
  'large fraction nitrogen detritus concentration functional'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'LdetritusN, scalar, series'
  'ocean_time'
  'idTads(iLDeN)'
  'r3dvar'
  1.0d0

'SdetritusN'                                       ! Input/Output
  'small fraction nitrogen detritus concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'SdetritusN, scalar, series'
  'ocean_time'
  'idTvar(iSDeN)'
  'r3dvar'
  1.0d0

'SdetritusN'                                       ! Input
  'small fraction nitrogen detritus concentration functional'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'SdetritusN, scalar, series'
  'ocean_time'
  'idTads(iSDeN)'
  'r3dvar'
  1.0d0

'FeSp'                                             ! Input/Output
  'small phytop. iron concentration'
  'millimole_iron  meter-3'                        ! [millimole/m3]
  'iron, scalar, series'
  'ocean_time'
  'idTvar(iFeSp)'
  'r3dvar'
  1.0d0

'FeLp'                                            ! Input/Output
  'large phytop. iron concentration'
  'millimole_iron  meter-3'                       ! [millimole/m3]
  'iron, scalar, series'
  'ocean_time'
  'idTvar(iFeLp)'
  'r3dvar'
  1.0d0

'FeD'                                             ! Input/Output
  'dissolved iron concentration'
  'millimole_iron  meter-3'                       ! [millimole/m3]
  'iron, scalar, series'
  'ocean_time'
  'idTvar(iFeD_)'
  'r3dvar'
  1.0d0

'PONsed'                                          ! Input/Output
  'Sediment PON'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'PONsed, scalar, series'
  'ocean_time'
  'idPONsed'
  'r2dvar'
  1.0d0

'OPALsed'                                         ! Input/Output
  'Sediment OPAL'
  'millimole_silica meter-3'                      ! [millimole/m3]
  'OPALsed, scalar, series'
  'ocean_time'
  'idOPALsed'
  'r2dvar'
  1.0d0

'DENITsed'                                        ! Input/Output
  'Sediment DENIT'
  'millimole_nitrogen meter-3'                    ! [millimole/m3]
  'DENITsed, scalar, series'
  'ocean_time'
  'idDENITsed'
  'r2dvar'
  1.0d0

'PON_burial'                                       ! Input/Output
  'Buried PON'
  'millimole_nitrogen meter-3'                      ! [millimole/m3]
  'PON_burial, scalar, series'
  'ocean_time'
  'idPONbur'
  'r2dvar'
  1.0d0

'OPAL_burial'                                     ! Input/Output
  'Buried OPAL'
  'millimole_silica meter-3'                      ! [millimole/m3]
  'OPAL_burial, scalar, series'
  'ocean_time'
  'idOPALbur'
  'r2dvar'
  1.0d0

'primary_prod'                                    ! Input/Output
  'Net primary production'
  'millimole_nitrate meter-3'                      ! [millimole/m3]
  'primary_prod, scalar, series'
  'ocean_time'
  'idNPP'
  'r2dvar'
  1.0d0

'TIC'                                              ! Input/Output
  'total inorganic carbon'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC, scalar, series'
  'ocean_time'
  'idTvar(iTIC_)'
  'r3dvar'
  1.0d0

'TIC'                                              ! Input
  'total inorganic carbon functional'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC, scalar, series'
  'ocean_time'
  'idTads(iTIC_)'
  'r3dvar'
  1.0d0

'alkalinity'                                       ! Input/Output
  'total alkalinity'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity, scalar, series'
  'ocean_time'
  'idTvar(iTAlk)'
  'r3dvar'
  1.0d0

'alkalinity'                                       ! Input
  'total alkalinity functional'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity, scalar, series'
  'ocean_time'
  'idTads(iTAlk)'
  'r3dvar'
  1.0d0

'LdetritusC'                                       ! Input/Output
  'large fraction carbon detritus concentration'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'LdetritusC, scalar, series'
  'ocean_time'
  'idTvar(iLDeC)'
  'r3dvar'
  1.0d0

'LdetritusC'                                       ! Input
  'large fraction carbon detritus concentration functional'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'LdetritusC, scalar, series'
  'ocean_time'
  'idTads(iLDeC)'
  'r3dvar'
  1.0d0

'SdetritusC'                                       ! Input/Output
  'small fraction carbon detritus concentration'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'SdetritusC, scalar, series'
  'ocean_time'
  'idTvar(iSDeC)'
  'r3dvar'
  1.0d0

'SdetritusC'                                       ! Input
  'small fraction carbon detritus concentration functional'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'SdetritusC, scalar, series'
  'ocean_time'
  'idTads(iSDeC)'
  'r3dvar'
  1.0d0

'oxygen'                                           ! Input/Output
  'dissolved oxygen concentration'
  'millimole_oxygen meter-3'                       ! [millimole/m3]
  'Oxygen, scalar, series'
  'ocean_time'
  'idTvar(iOxyg)'
  'r3dvar'
  1.0d0

'oxygen'                                           ! Input
  'dissolved oxygen concentration functional'
  'millimole_oxygen meter-3'                       ! [millimole/m3]
  'Oxygen, scalar, series'
  'ocean_time'
  'idTads(iOxyg)'
  'r3dvar'
  1.0d0

!
!  Diagnostic variables.
!

'CO2_airsea'                                       ! Input/Output
  'CO2 air-sea flux'
  'millimole_carbon meter-2 day-1'                 ! [millimole/m2/day]
  'CO2_airsea, scalar, series'
  'ocean_time'
  'iDbio2(iCOfx)'
  'r2dvar'
  1.0d0

'denitrification'                                  ! Input/Output
  'denitrification'
  'millimole_nitrogen meter-2 day-1'               ! [millimole/m2/day]
  'denitrification, scalar, series'
  'ocean_time'
  'iDbio2(iDNIT)'
  'r2dvar'
  1.0d0

'pCO2'                                             ! Input/Output
  'partial pressure of CO2'
  'atm'                                            ! [atm]
  'pCO2, scalar, series'
  'ocean_time'
  'iDbio2(ipCO2)'
  'r2dvar'
  1.0d0

'O2_airsea'                                        ! Input/Output
  'O2 air-sea flux'
  'millimole_oxygen meter-2 day-1'                 ! [millimole/m2/day]
  'O2_airsea, scalar, series'
  'ocean_time'
  'iDbio2(iO2fx)'
  'r2dvar'
  1.0d0

'P_Production'                                     ! Input/Output
  'primary production'
  'millimole_nitrogen meter-3 day-1'               ! [millimole/m3/day]
  'P_Production, scalar, series'
  'ocean_time'
  'iDbio3(iPPro)'
  'r3dvar'
  1.0d0

'NO3_uptake'                                       ! Input/Output
  'NO3 uptake'
  'millimole_nitrogen meter-3 day-1'               ! [millimole/m3/day]
  'NO3_uptake, scalar, series'
  'ocean_time'
  'iDbio3(iNO3u)'
  'r3dvar'
  1.0d0

!--------------UMAINE variables
'NO3'                                              ! Input
  'nitrate concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3, scalar, series'
  'ocean_time'
  'idTvar(iNO3_)'
  'r3dvar'
  1.0d0

'NO3'                                              ! Input
  'nitrate concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3, scalar, series'
  'ocean_time'
  'idTclm(iNO3_)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'ammonium concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4, scalar, series'
  'ocean_time'
  'idTvar(iNH4_)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'ammonium concentration'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4, scalar, series'
  'ocean_time'
  'idTclm(iNH4_)'
  'r3dvar'
  1.0d0

'SiOH'                                             ! Input/Output
  'silicate concentration'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH, scalar, series'
  'ocean_time'
  'idTvar(iSiOH)'
  'r3dvar'
  1.0d0

'SiOH'                                             ! Input/Output
  'silicate concentration'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH, scalar, series'
  'ocean_time'
  'idTclm(iSiOH)'
  'r3dvar'
  1.0d0

'PO4'                                              ! Input/Output
  'phosphate concentration'
  'millimole_PO4 meter-3'                          ! [millimole/m3]
  'PO4, scalar, series'
  'ocean_time'
  'idTvar(iPO4_)'
  'r3dvar'
  1.0d0

'PO4'                                              ! Input/Output
  'phosphate concentration'
  'millimole_PO4 meter-3'                          ! [millimole/m3]
  'PO4, scalar, series'
  'ocean_time'
  'idTclm(iPO4_)'
  'r3dvar'
  1.0d0

'S1_N'                                              ! Input/Output
  'S1_N'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_N, scalar, series'
  'ocean_time'
  'idTvar(iS1_N)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'S1_N'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_N, scalar, series'
  'ocean_time'
  'idTclm(iS1_N)'
  'r3dvar'
  1.0d0

'S1_C'                                              ! Input/Output
  'S1_C'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_C, scalar, series'
  'ocean_time'
  'idTvar(iS1_C)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'S1_C'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_C, scalar, series'
  'ocean_time'
  'idTclm(iS1_C)'
  'r3dvar'
  5.6d0

'S1CH'                                              ! Input/Output
  'S1CH'
  'milligram meter-3'                     ! [millimole/m3]
  'S1CH, scalar, series'
  'ocean_time'
  'idTvar(iS1CH)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'S1CH'
  'milligram meter-3'                     ! [millimole/m3]
  'S1CH, scalar, series'
  'ocean_time'
  'idTclm(iS1CH)'
  'r3dvar'
  3.4d0

'S2_N'                                              ! Input/Output
  'S2_N'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_N, scalar, series'
  'ocean_time'
  'idTvar(iS2_N)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'S2_N'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_N, scalar, series'
  'ocean_time'
  'idTclm(iS2_N)'
  'r3dvar'
  1.0d0

'S2_C'                                              ! Input/Output
  'S2_C'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_C, scalar, series'
  'ocean_time'
  'idTvar(iS2_C)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'S2_C'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_C, scalar, series'
  'ocean_time'
  'idTclm(iS2_C)'
  'r3dvar'
  5.6d0

'S2CH'                                              ! Input/Output
  'S2CH'
  'milligram meter-3'                     ! [millimole/m3]
  'S2CH, scalar, series'
  'ocean_time'
  'idTvar(iS2CH)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'S2CH'
  'milligram meter-3'                     ! [millimole/m3]
  'S2CH, scalar, series'
  'ocean_time'
  'idTclm(iS2CH)'
  'r3dvar'
  3.4d0

'S3_C'                                              ! Input/Output
  'S3_C'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_C, scalar, series'
  'ocean_time'
  'idTvar(iS3_C)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'S3_C'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_C, scalar, series'
  'ocean_time'
  'idTclm(iS3_C)'
  'r3dvar'
  5.6d0

'S3_N'                                              ! Input/Output
  'S3_N'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_N, scalar, series'
  'ocean_time'
  'idTvar(iS3_N)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'S3_N'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_N, scalar, series'
  'ocean_time'
  'idTclm(iS3_N)'
  'r3dvar'
  1.0d0

'S3CH'                                              ! Input/Output
  'S3CH'
  'milligram meter-3'                     ! [millimole/m3]
  'S3CH, scalar, series'
  'ocean_time'
  'idTvar(iS3CH)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'S3CH'
  'milligram meter-3'                     ! [millimole/m3]
  'S3CH, scalar, series'
  'ocean_time'
  'idTclm(iS3CH)'
  'r3dvar'
  3.4d0

'Z1_N'                                              ! Input/Output
  'Z1_N'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_N, scalar, series'
  'ocean_time'
  'idTvar(iZ1_N)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'Z1_N'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_N, scalar, series'
  'ocean_time'
  'idTclm(iZ1_N)'
  'r3dvar'
  1.0d0

'Z1_C'                                              ! Input/Output
  'Z1_C'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_C, scalar, series'
  'ocean_time'
  'idTvar(iZ1_C)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'Z1_C'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_C, scalar, series'
  'ocean_time'
  'idTclm(iZ1_C)'
  'r3dvar'
  5.6d0

'Z2_N'                                              ! Input/Output
  'Z2_N'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_N, scalar, series'
  'ocean_time'
  'idTvar(iZ2_N)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'Z2_N'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_N, scalar, series'
  'ocean_time'
  'idTclm(iZ2_N)'
  'r3dvar'
  1.0d0

'Z2_C'                                              ! Input/Output
  'Z2_C'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_C, scalar, series'
  'ocean_time'
  'idTvar(iZ2_C)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'Z2_C'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_C, scalar, series'
  'ocean_time'
  'idTclm(iZ2_C)'
  'r3dvar'
  5.6d0

'BAC_'                                              ! Input/Output
  'BAC'
  'millimole meter-3'                     ! [millimole/m3]
  'BAC, scalar, series'
  'ocean_time'
  'idTvar(iBAC_)'
  'r3dvar'
  1.0d0

'BAC'                                              ! Input/Output
  'BAC'
  'millimole meter-3'                     ! [millimole/m3]
  'BAC, scalar, series'
  'ocean_time'
  'idTclm(iBAC_)'
  'r3dvar'
  1.0d0

'DD_N'                                              ! Input/Output
  'DD_N'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_N, scalar, series'
  'ocean_time'
  'idTvar(iDD_N)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'DD_N'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_N, scalar, series'
  'ocean_time'
  'idTclm(iDD_N)'
  'r3dvar'
  1.0d0

'DD_C'                                              ! Input/Output
  'DD_C'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_C, scalar, series'
  'ocean_time'
  'idTvar(iDD_C)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'DD_C'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_C, scalar, series'
  'ocean_time'
  'idTclm(iDD_C)'
  'r3dvar'
  5.6d0

'DDSi'                                              ! Input/Output
  'DDSi'
  'millimole meter-3'                     ! [millimole/m3]
  'DDSi, scalar, series'
  'ocean_time'
  'idTvar(iDDSi)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'DDSi'
  'millimole meter-3'                     ! [millimole/m3]
  'DDSi, scalar, series'
  'ocean_time'
  'idTclm(iDDSi)'
  'r3dvar'
  1.0d0

'LDON'                                              ! Input/Output
  'LDON'
  'millimole meter-3'                     ! [millimole/m3]
  'LDON, scalar, series'
  'ocean_time'
  'idTvar(iLDON)'
  'r3dvar'
  1.0d0

'LDOC'                                              ! Input/Output
  'LDON'
  'millimole meter-3'                     ! [millimole/m3]
  'LDON, scalar, series'
  'ocean_time'
  'idTclm(iLDON)'
  'r3dvar'
  0.1d0

'LDOC'                                              ! Input/Output
  'LDOC'
  'millimole meter-3'                     ! [millimole/m3]
  'LDOC, scalar, series'
  'ocean_time'
  'idTvar(iLDOC)'
  'r3dvar'
  1.0d0

'LDOC'                                              ! Input/Output
  'LDOC'
  'millimole meter-3'                     ! [millimole/m3]
  'LDOC, scalar, series'
  'ocean_time'
  'idTclm(iLDOC)'
  'r3dvar'
  1.0d0

'SDON'                                              ! Input/Output
  'SDON'
  'millimole meter-3'                     ! [millimole/m3]
  'SDON, scalar, series'
  'ocean_time'
  'idTvar(iSDON)'
  'r3dvar'
  1.0d0

'SDOC'                                              ! Input/Output
  'SDON'
  'millimole meter-3'                     ! [millimole/m3]
  'SDON, scalar, series'
  'ocean_time'
  'idTclm(iSDON)'
  'r3dvar'
  0.065d0

'SDOC'                                              ! Input/Output
  'SDOC'
  'millimole meter-3'                     ! [millimole/m3]
  'SDOC, scalar, series'
  'ocean_time'
  'idTvar(iSDOC)'
  'r3dvar'
  1.0d0

'SDOC'                                              ! Input/Output
  'SDOC'
  'millimole meter-3'                     ! [millimole/m3]
  'SDOC, scalar, series'
  'ocean_time'
  'idTclm(iSDOC)'
  'r3dvar'
  1.0d0

'DDCA'                                              ! Input/Output
  'DDCA'
  'millimole meter-3'                     ! [millimole/m3]
  'DDCA, scalar, series'
  'ocean_time'
  'idTvar(iDDCA)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'DDCA'
  'millimole meter-3'                     ! [millimole/m3]
  'DDCA, scalar, series'
  'ocean_time'
  'idTclm(iDDCA)'
  'r3dvar'
  5.6d0

'CLDC'                                              ! Input/Output
  'CLDC'
  'millimole meter-3'                     ! [millimole/m3]
  'CLDC, scalar, series'
  'ocean_time'
  'idTvar(iCLDC)'
  'r3dvar'
  1.0d0

'LDOC'                                              ! Input/Output
  'CLDC'
  'millimole meter-3'                     ! [millimole/m3]
  'CLDC, scalar, series'
  'ocean_time'
  'idTclm(iCLDC)'
  'r3dvar'
  1.4d0

'CSDC'                                              ! Input/Output
  'CSDC'
  'millimole meter-3'                     ! [millimole/m3]
  'CSDC, scalar, series'
  'ocean_time'
  'idTvar(iCSDC)'
  'r3dvar'
  1.0d0

'SDOC'                                              ! Input/Output
  'CSDC'
  'millimole meter-3'                     ! [millimole/m3]
  'CSDC, scalar, series'
  'ocean_time'
  'idTclm(iCSDC)'
  'r3dvar'
  0.0088d0

'oxyg'                                           ! Input/Output
  'dissolved oxygen concentration'
  'millimole_O2 meter-3'                       ! [millimole/m3]
  'Oxygen, scalar, series'
  'ocean_time'
  'idTvar(iOxyg)'
  'r3dvar'
  1.0d0

'Oxyg'                                           ! Input/Output
  'dissolved oxygen concentration'
  'millimole_O2 meter-3'                       ! [millimole/m3]
  'Oxygen, scalar, series'
  'ocean_time'
  'idTclm(iOxyg)'
  'r3dvar'
  1.0d0

'TIC'                                              ! Input
  'total inorganic carbon'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC, scalar, series'
  'ocean_time'
  'idTvar(iTIC_)'
  'r3dvar'
  1.0d0

'TIC'                                              ! Input
  'total inorganic carbon'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC, scalar, series'
  'ocean_time'
  'idTclm(iTIC_)'
  'r3dvar'
  1.0d0

'Talk'                                       ! Input/Output
  'total alkalinity'
  'millimole meter-3'                       ! [millimole/m3]
  'alkalinity, scalar, series'
  'ocean_time'
  'idTvar(iTAlk)'
  'r3dvar'
  1.0d0

'Talk'                                       ! Input/Output
  'total alkalinity'
  'millimole meter-3'                       ! [millimole/m3]
  'alkalinity, scalar, series'
  'ocean_time'
  'idTclm(iTAlk)'
  'r3dvar'
  1.0d0

!--------UMAINE variables

!
!  Boundary variables.
!

'NO3_east'                                         ! Input
  'nitrate eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iNO3_)'
  'nulvar'
  1.0d0

'NO3_west'                                         ! Input
  'nitrate western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iNO3_)'
  'nulvar'
  1.0d0

'NO3_south'                                        ! Input
  'nitrate southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iNO3_)'
  'nulvar'
  1.0d0

'NO3_north'                                        ! Input
  'nitrate northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iNO3_)'
  'nulvar'
  1.0d0

'PO4_east'                                         ! Input
  'phosphate eastern boundary condition'
  'millimole_P meter-3'                     ! [millimole/m3]
  'PO4_east, scalar, series'
  'PO4_time'
  'idTbry(ieast,iPO4_)'
  'nulvar'
  1.0d0

'PO4_west'                                         ! Input
  'phosphate western boundary condition'
  'millimole_P meter-3'                     ! [millimole/m3]
  'PO4_west, scalar, series'
  'PO4_time'
  'idTbry(iwest,iPO4_)'
  'nulvar'
  1.0d0

'PO4_south'                                        ! Input
  'phosphate southern boundary condition'
  'millimole_P meter-3'                     ! [millimole/m3]
  'PO4_south, scalar, series'
  'PO4_time'
  'idTbry(isouth,iPO4_)'
  'nulvar'
  1.0d0

'PO4_north'                                        ! Input
  'phosphate northern boundary condition'
  'millimole_P meter-3'                     ! [millimole/m3]
  'PO4_north, scalar, series'
  'PO4_time'
  'idTbry(inorth,iPO4_)'
  'nulvar'
  1.0d0

'TIC_east'                                         ! Input
  'TIC eastern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC_east, scalar, series'
  'TIC_time'
  'idTbry(ieast,iTIC_)'
  'nulvar'
  1.0d0

'TIC_west'                                         ! Input
  'TIC western boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC_west, scalar, series'
  'TIC_time'
  'idTbry(iwest,iTIC_)'
  'nulvar'
  1.0d0

'TIC_north'                                        ! Input
  'TIC northern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC_north, scalar, series'
  'TIC_time'
  'idTbry(inorth,iTIC_)'
  'nulvar'
  1.0d0

'TIC_south'                                        ! Input
  'TIC southern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC_south, scalar, series'
  'TIC_time'
  'idTbry(isouth,iTIC_)'
  'nulvar'
  1.0d0

'oxygen_east'                                      ! Input
  'oxygen eastern boundary condition'
  'millimole_oxygen meter-3'                       ! [millimole/m3]
  'oxygen_east, scalar, series'
  'oxygen_time'
  'idTbry(ieast,iOxyg)'
  'nulvar'
  1.0d0

'oxygen_west'                                      ! Input
  'oxygen western boundary condition'
  'millimole_oxygen meter-3'                       ! [millimole/m3]
  'oxygen_west, scalar, series'
  'oxygen_time'
  'idTbry(iwest,iOxyg)'
  'nulvar'
  1.0d0

'oxygen_north'                                     ! Input
  'oxygen northern boundary condition'
  'millimole_oxygen meter-3'                       ! [millimole/m3]
  'oxygen_north, scalar, series'
  'oxygen_time'
  'idTbry(inorth,iOxyg)'
  'nulvar'
  1.0d0

'oxygen_south'                                     ! Input
  'oxygen southern boundary condition'
  'millimole_oxygen meter-3'                       ! [millimole/m3]
  'oxygen_south, scalar, series'
  'oxygen_time'
  'idTbry(isouth,iOxyg)'
  'nulvar'
  1.0d0

'alkalinity_east'                                  ! Input
  'alkalinity eastern boundary condition'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity_east, scalar, series'
  'alkalinity_time'
  'idTbry(ieast,iTAlk)'
  'nulvar'
  1.0d0

'alkalinity_west'                                  ! Input
  'alkalinity western boundary condition'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity_west, scalar, series'
  'alkalinity_time'
  'idTbry(iwest,iTAlk)'
  'nulvar'
  1.0d0

'alkalinity_north'                                 ! Input
  'alkalinity northern boundary condition'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity_north, scalar, series'
  'alkalinity_time'
  'idTbry(inorth,iTAlk)'
  'nulvar'
  1.0d0

'alkalinity_south'                                 ! Input
  'alkalinity southern boundary condition'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity_south, scalar, series'
  'alkalinity_time'
  'idTbry(isouth,iTAlk)'
  'nulvar'
  1.0d0

'NH4_east'                                         ! Input
  'ammonium eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iNH4_)'
  'nulvar'
  1.0d0

'NH4_west'                                         ! Input
  'ammonium western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iNH4_)'
  'nulvar'
  1.0d0

'NH4_south'                                        ! Input
  'ammonium southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iNH4_)'
  'nulvar'
  1.0d0

'NH4_north'                                        ! Input
  'ammonium northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iNH4_)'
  'nulvar'
  1.0d0

'DON_east'                                         ! Input
  'dissolved organic nitrogen eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'DON_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iDON_)'
  'nulvar'
  1.0d0

'DON_west'                                         ! Input
  'dissolved organic nitrogen western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'DON_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iDON_)'
  'nulvar'
  1.0d0

'DON_south'                                        ! Input
  'dissolved organic nitrogen southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'DON_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iDON_)'
  'nulvar'
  1.0d0

'DON_north'                                        ! Input
  'dissolved organic nitrogen northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'DON_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iDON_)'
  'nulvar'
  1.0d0

'PON_east'                                         ! Input
  'particulate organic nitrogen eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'PON_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iPON_)'
  'nulvar'
  1.0d0

'PON_west'                                         ! Input
  'particulate organic nitrogen western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'PON_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iPON_)'
  'nulvar'
  1.0d0

'PON_south'                                        ! Input
  'particulate organic nitrogen southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'PON_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iPON_)'
  'nulvar'
  1.0d0

'PON_north'                                        ! Input
  'particulate organic nitrogen northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'PON_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iPON_)'
  'nulvar'
  1.0d0

'SiOH4_east'                                       ! Input
  'silicate eastern boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH4_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iSiOH)'
  'nulvar'
  1.0d0

'SiOH4_west'                                       ! Input
  'silicate western boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH4_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iSiOH)'
  'nulvar'
  1.0d0

'SiOH4_south'                                      ! Input
  'silicate southern boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH4_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iSiOH)'
  'nulvar'
  1.0d0

'SiOH4_north'                                      ! Input
  'silicate northern boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH4_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iSiOH)'
  'nulvar'
  1.0d0

'opal_east'                                        ! Input
  'particulate organic silica eastern boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'opal_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iopal)'
  'nulvar'
  1.0d0

'opal_west'                                        ! Input
  'particulate organic silica western boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'opal_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iopal)'
  'nulvar'
  1.0d0

'opal_south'                                       ! Input
  'particulate organic silica southern boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'opal_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iopal)'
  'nulvar'
  1.0d0

'opal_north'                                       ! Input
  'particulate organic silica northern boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'opal_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iopal)'
  'nulvar'
  1.0d0

'FeSp_east'                                        ! Input
  'small phytop. iron eastern boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeSp_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iFeSp)'
  'nulvar'
  1.0d0

'FeSp_west'                                        ! Input
  'small phytop. iron western boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeSp_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iFeSp)'
  'nulvar'
  1.0d0

'FeSp_north'                                       ! Input
  'small phytop. iron northern boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeSp_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iFeSp)'
  'nulvar'
  1.0d0

'FeSp_south'                                       ! Input
  'small phytop. iron southern boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeSp_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iFeSp)'
  'nulvar'
  1.0d0

'FeLp_east'                                        ! Input
  'large phytop. iron eastern boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeLp_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iFeLp)'
  'nulvar'
  1.0d0

'FeLp_west'                                        ! Input
  'large phytop. iron western boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeLp_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iFeLp)'
  'nulvar'
  1.0d0

'FeLp_north'                                       ! Input
  'large phytop. iron northern boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeLp_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iFeLp)'
  'nulvar'
  1.0d0

'FeLp_south'                                       ! Input
  'large phytop. iron southern boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeLp_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iFeLp)'
  'nulvar'
  1.0d0

'FeD_east'                                         ! Input
  'dissolved iron eastern boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeD_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iFeD_)'
  'nulvar'
  1.0d0

'FeD_west'                                         ! Input
  'dissolved iron western boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeD_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iFeD_)'
  'nulvar'
  1.0d0

'FeD_north'                                        ! Input
  'dissolved iron northern boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeD_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iFeD_)'
  'nulvar'
  1.0d0

'FeD_south'                                        ! Input
  'dissolved iron southern boundary condition'
  'millimole_iron meter-3'                         ! [millimole/m3]
  'FeD_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iFeD_)'
  'nulvar'
  1.0d0

'phyt_east'                                        ! Input
  'phytoplankton eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'phyt_east, scalar, series'
  'phyt_time'
  'idTbry(ieast,iPhyt)'
  'nulvar'
  1.0d0

'phyt_west'                                        ! Input
  'phytoplankton western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'phyt_west, scalar, series'
  'phyt_time'
  'idTbry(iwest,iPhyt)'
  'nulvar'
  1.0d0

'phyt_south'                                       ! Input
  'phytoplankton southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'phyt_south, scalar, series'
  'phyt_time'
  'idTbry(isouth,iPhyt)'
  'nulvar'
  1.0d0

'phyt_north'                                       ! Input
  'phytoplankton northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'phyt_north, scalar, series'
  'phyt_time'
  'idTbry(inorth,iPhyt)'
  'nulvar'
  1.0d0

'nanophy_east'                                     ! Input
  'nanophytoplankton eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'nanophy_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iSphy)'
  'nulvar'
  1.0d0

'nanophy_west'                                     ! Input
  'nanophytoplankton western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'nanophy_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iSphy)'
  'nulvar'
  1.0d0

'nanophy_south'                                    ! Input
  'nanophytoplankton southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'nanophy_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iSphy)'
  'nulvar'
  1.0d0

'nanophy_north'                                    ! Input
  'nanophytoplankton northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'nanophy_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iSphy)'
  'nulvar'
  1.0d0

'diatom_east'                                      ! Input
  'diatom eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'diatom_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iLphy)'
  'nulvar'
  1.0d0

'diatom_west'                                      ! Input
  'diatom western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'diatom_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iLphy)'
  'nulvar'
  1.0d0

'diatom_south'                                     ! Input
  'diatom southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'diatom_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iLphy)'
  'nulvar'
  1.0d0

'diatom_north'                                     ! Input
  'diatom northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'diatom_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iLphy)'
  'nulvar'
  1.0d0

'zoop_east'                                        ! Input
  'zooplankton eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'zoop_east, scalar, series'
  'zoop_time'
  'idTbry(ieast,iZoop)'
  'nulvar'
  1.0d0

'zoop_west'                                        ! Input
  'zooplankton western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'zoop_west, scalar, series'
  'zoop_time'
  'idTbry(iwest,iZoop)'
  'nulvar'
  1.0d0

'zoop_south'                                       ! Input
  'zooplankton southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'zoop_south, scalar, series'
  'zoop_time'
  'idTbry(isouth,iZoop)'
  'nulvar'
  1.0d0

'zoop_north'                                       ! Input
  'zooplankton northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'zoop_north, scalar, series'
  'zoop_time'
  'idTbry(inorth,iZoop)'
  'nulvar'
  1.0d0

'microzoo_east'                                    ! Input
  'microzooplankton eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'microzoo_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iSzoo)'
  'nulvar'
  1.0d0

'microzoo_west'                                    ! Input
  'microzooplankton western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'microzoo_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iSzoo)'
  'nulvar'
  1.0d0

'microzoo_south'                                   ! Input
  'microzooplankton southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'microzoo_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iSzoo)'
  'nulvar'
  1.0d0

'microzoo_north'                                   ! Input
  'microzooplankton northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'microzoo_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iSzoo)'
  'nulvar'
  1.0d0

'mesozoo_east'                                     ! Input
  'mesozooplankton eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'mesozoo_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iLzoo)'
  'nulvar'
  1.0d0

'mesozoo_west'                                     ! Input
  'mesozooplankton western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'mesozoo_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iLzoo)'
  'nulvar'
  1.0d0

'mesozoo_south'                                    ! Input
  'mesozooplankton southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'mesozoo_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iLzoo)'
  'nulvar'
  1.0d0

'mesozoo_north'                                    ! Input
  'mesozooplankton northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'mesozoo_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iLzoo)'
  'nulvar'
  1.0d0

'Pzoo_east'                                        ! Input
  'predator-zooplankton eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'Pzoo_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iPzoo)'
  'nulvar'
  1.0d0

'Pzoo_west'                                        ! Input
  'predator-zooplankton western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'Pzoo_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iPzoo)'
  'nulvar'
  1.0d0

'Pzoo_south'                                       ! Input
  'predator-zooplankton southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'Pzoo_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iPzoo)'
  'nulvar'
  1.0d0

'Pzoo_north'                                       ! Input
  'predator-zooplankton northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'Pzoo_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iPzoo)'
  'nulvar'
  1.0d0

'chlo_east'                                        ! Input
  'chlorophyll eastern boundary condition'
  'milligram meter-3'                              ! [milligram/m3]
  'chlo_east, scalar, series'
  'chlo_time'
  'idTbry(ieast,iChlo)'
  'nulvar'
  1.0d0

'chlo_west'                                        ! Input
  'chlorophyll western boundary condition'
  'milligram meter-3'                              ! [milligram/m3]
  'chlo_west, scalar, series'
  'chlo_time'
  'idTbry(iwest,iChlo)'
  'nulvar'
  1.0d0

'chlo_south'                                       ! Input
  'chlorophyll southern boundary condition'
  'milligram meter-3'                              ! [milligram/m3]
  'chlo_south, scalar, series'
  'chlo_time'
  'idTbry(isouth,iChlo)'
  'nulvar'
  1.0d0

'chlo_north'                                       ! Input
  'chlorophyll northern boundary condition'
  'milligram meter-3'                              ! [milligram/m3]
  'chlo_north, scalar, series'
  'chlo_time'
  'idTbry(inorth,iChlo)'
  'nulvar'
  1.0d0

'detritus_east'                                    ! Input
  'detritus eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'detritus_east, scalar, series'
  'detritus_time'
  'idTbry(ieast,iSDet)'
  'nulvar'
  1.0d0

'detritus_west'                                    ! Input
  'detritus western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'detritus_west, scalar, series'
  'detritus_time'
  'idTbry(iwest,iSDet)'
  'nulvar'
  1.0d0

'detritus_south'                                   ! Input
  'detritus southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'detritus_south, scalar, series'
  'detritus_time'
  'idTbry(isouth,iSDet)'
  'nulvar'
  1.0d0

'detritus_north'                                   ! Input
  'detritus northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'detritus_north, scalar, series'
  'detritus_time'
  'idTbry(inorth,iSDet)'
  'nulvar'
  1.0d0

'LDeN_east'                                        ! Input
  'large fraction nitrogen detritus eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'LDeN_east, scalar, series'
  'LDeN_time'
  'idTbry(ieast,iLDeN)'
  'nulvar'
  1.0d0

'LDeN_west'                                        ! Input
  'large fraction nitrogen detritus western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'LDeN_west, scalar, series'
  'LDeN_time'
  'idTbry(iwest,iLDeN)'
  'nulvar'
  1.0d0

'LDeN_south'                                       ! Input
  'large fraction nitrogen detritus southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'LDeN_south, scalar, series'
  'LDeN_time'
  'idTbry(isouth,iLDeN)'
  'nulvar'
  1.0d0

'LDeN_north'                                       ! Input
  'large fraction nitrogen detritus northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'LDeN_north, scalar, series'
  'LDeN_time'
  'idTbry(inorth,iLDeN)'
  'nulvar'
  1.0d0

'SDeN_east'                                        ! Input
  'small fraction nitrogen detritus eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'SDeN_east, scalar, series'
  'SDeN_time'
  'idTbry(ieast,iSDeN)'
  'nulvar'
  1.0d0

'SDeN_west'                                        ! Input
  'small fraction nitrogen detritus western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'SDeN_west, scalar, series'
  'SDeN_time'
  'idTbry(iwest,iSDeN)'
  'nulvar'
  1.0d0

'SDeN_north'                                       ! Input
  'small fraction nitrogen detritus northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'SDeN_north, scalar, series'
  'SDeN_time'
  'idTbry(inorth,iSDeN)'
  'nulvar'
  1.0d0

'SDeN_south'                                       ! Input
  'small fraction nitrogen detritus southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'SDeN_south, scalar, series'
  'SDeN_time'
  'idTbry(isouth,iSDeN)'
  'nulvar'
  1.0d0

'LdeC_east'                                        ! Input
  'large fraction carbon detritus eastern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'LdeC_east, scalar, series'
  'LdeC_time'
  'idTbry(ieast,iLDeC)'
  'nulvar'
  1.0d0

'LdeC_west'                                        ! Input
  'large fraction carbon detritus western boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'LdeC_west, scalar, series'
  'LdeC_time'
  'idTbry(iwest,iLDeC)'
  'nulvar'
  1.0d0

'LdeC_south'                                       ! Input
  'large fraction carbon detritus southern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'LdeC_south, scalar, series'
  'LdeC_time'
  'idTbry(isouth,iLDeC)'
  'nulvar'
  1.0d0

'LdeC_north'                                       ! Input
  'large fraction carbon detritus northern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'LdeC_north, scalar, series'
  'LdeC_time'
  'idTbry(inorth,iLDeC)'
  'nulvar'
  1.0d0

'SdeC_east'                                        ! Input
  'small fraction carbon detritus eastern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'SdeC_east, scalar, series'
  'SdeC_time'
  'idTbry(ieast,iSDeC)'
  'nulvar'
  1.0d0

'SdeC_west'                                        ! Input
  'small fraction carbon detritus western boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'SdeC_west, scalar, series'
  'SdeC_time'
  'idTbry(iwest,iSDeC)'
  'nulvar'
  1.0d0

'SdeC_south'                                       ! Input
  'small fraction carbon detritus southern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'SdeC_south, scalar, series'
  'SdeC_time'
  'idTbry(isouth,iSDeC)'
  'nulvar'
  1.0d0

'SdeC_north'                                       ! Input
  'small fraction carbon detritus northern boundary condition'
  'millimole_carbon meter-3'
  'SdeC_north, scalar, series'                     ! [millimole/m3]
  'SdeC_time'
  'idTbry(inorth,iSDeC)'
  'nulvar'
  1.0d0

!----------UMAINE variables
'S1_N_east'                                         ! Input
  'S1_N eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_N_east, scalar, series'
  'S1_N_time'
  'idTbry(ieast,iS1_N)'
  'nulvar'
  1.0d0

'S1_N_west'                                         ! Input
  'S1_N western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_N_west, scalar, series'
  'S1_N_time'
  'idTbry(iwest,iS1_N)'
  'nulvar'
  1.0d0

'S1_N_south'                                        ! Input
  'S1_N southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_N_south, scalar, series'
  'S1_N_time'
  'idTbry(isouth,iS1_N)'
  'nulvar'
  1.0d0

'S1_N_north'                                        ! Input
  'S1_N northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_N_north, scalar, series'
  'S1_N_time'
  'idTbry(inorth,iS1_N)'
  'nulvar'
  1.0d0

'S1_C_east'                                         ! Input
  'S1_C eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_C_east, scalar, series'
  'S1_C_time'
  'idTbry(ieast,iS1_C)'
  'nulvar'
  1.0d0

'S1_C_west'                                         ! Input
  'S1_C western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_C_west, scalar, series'
  'S1_C_time'
  'idTbry(iwest,iS1_C)'
  'nulvar'
  1.0d0

'S1_C_south'                                        ! Input
  'S1_C southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_C_south, scalar, series'
  'S1_C_time'
  'idTbry(isouth,iS1_C)'
  'nulvar'
  1.0d0

'S1_C_north'                                        ! Input
  'S1_C northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S1_C_north, scalar, series'
  'S1_C_time'
  'idTbry(inorth,iS1_C)'
  'nulvar'
  1.0d0

'S1CH_east'                                         ! Input
  'S1CH eastern boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S1CH_east, scalar, series'
  'S1CH_time'
  'idTbry(ieast,iS1CH)'
  'nulvar'
  1.0d0

'S1CH_west'                                         ! Input
  'S1CH western boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S1CH_west, scalar, series'
  'S1CH_time'
  'idTbry(iwest,iS1CH)'
  'nulvar'
  1.0d0

'S1CH_south'                                        ! Input
  'S1CH southern boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S1CH_south, scalar, series'
  'S1CH_time'
  'idTbry(isouth,iS1CH)'
  'nulvar'
  1.0d0

'S1CH_north'                                        ! Input
  'S1CH northern boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S1CH_north, scalar, series'
  'S1CH_time'
  'idTbry(inorth,iS1CH)'
  'nulvar'
  1.0d0

'S2_N_east'                                         ! Input
  'S2_N eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_N_east, scalar, series'
  'S2_N_time'
  'idTbry(ieast,iS2_N)'
  'nulvar'
  1.0d0

'S2_N_west'                                         ! Input
  'S2_N western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_N_west, scalar, series'
  'S2_N_time'
  'idTbry(iwest,iS2_N)'
  'nulvar'
  1.0d0

'S2_N_south'                                        ! Input
  'S2_N southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_N_south, scalar, series'
  'S2_N_time'
  'idTbry(isouth,iS2_N)'
  'nulvar'
  1.0d0

'S2_N_north'                                        ! Input
  'S2_N northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_N_north, scalar, series'
  'S2_N_time'
  'idTbry(inorth,iS2_N)'
  'nulvar'
  1.0d0

'S2_C_east'                                         ! Input
  'S2_C eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_C_east, scalar, series'
  'S2_C_time'
  'idTbry(ieast,iS2_C)'
  'nulvar'
  1.0d0

'S2_C_west'                                         ! Input
  'S2_C western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_C_west, scalar, series'
  'S2_C_time'
  'idTbry(iwest,iS2_C)'
  'nulvar'
  1.0d0

'S2_C_south'                                        ! Input
  'S2_C southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_C_south, scalar, series'
  'S2_C_time'
  'idTbry(isouth,iS2_C)'
  'nulvar'
  1.0d0

'S2_C_north'                                        ! Input
  'S2_C northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S2_C_north, scalar, series'
  'S2_C_time'
  'idTbry(inorth,iS2_C)'
  'nulvar'
  1.0d0

'S2CH_east'                                         ! Input
  'S2CH eastern boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S2CH_east, scalar, series'
  'S2CH_time'
  'idTbry(ieast,iS2CH)'
  'nulvar'
  1.0d0

'S2CH_west'                                         ! Input
  'S2CH western boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S2CH_west, scalar, series'
  'S2CH_time'
  'idTbry(iwest,iS2CH)'
  'nulvar'
  1.0d0

'S2CH_south'                                        ! Input
  'S2CH southern boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S2CH_south, scalar, series'
  'S2CH_time'
  'idTbry(isouth,iS2CH)'
  'nulvar'
  1.0d0

'S2CH_north'                                        ! Input
  'S2CH northern boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S2CH_north, scalar, series'
  'S2CH_time'
  'idTbry(inorth,iS2CH)'
  'nulvar'
  1.0d0

'S3_N_east'                                         ! Input
  'S3_N eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_N_east, scalar, series'
  'S3_N_time'
  'idTbry(ieast,iS3_N)'
  'nulvar'
  1.0d0

'S3_N_west'                                         ! Input
  'S3_N western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_N_west, scalar, series'
  'S3_N_time'
  'idTbry(iwest,iS3_N)'
  'nulvar'
  1.0d0

'S3_N_south'                                        ! Input
  'S3_N southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_N_south, scalar, series'
  'S3_N_time'
  'idTbry(isouth,iS3_N)'
  'nulvar'
  1.0d0

'S3_N_north'                                        ! Input
  'S3_N northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_N_north, scalar, series'
  'S3_N_time'
  'idTbry(inorth,iS3_N)'
  'nulvar'
  1.0d0

'S3_C_east'                                         ! Input
  'S3_C eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_C_east, scalar, series'
  'S3_C_time'
  'idTbry(ieast,iS3_C)'
  'nulvar'
  1.0d0

'S3_C_west'                                         ! Input
  'S3_C western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_C_west, scalar, series'
  'S3_C_time'
  'idTbry(iwest,iS3_C)'
  'nulvar'
  1.0d0

'S3_C_south'                                        ! Input
  'S3_C southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_C_south, scalar, series'
  'S3_C_time'
  'idTbry(isouth,iS3_C)'
  'nulvar'
  1.0d0

'S3_C_north'                                        ! Input
  'S3_C northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'S3_C_north, scalar, series'
  'S3_C_time'
  'idTbry(inorth,iS3_C)'
  'nulvar'
  1.0d0

'S3CH_east'                                         ! Input
  'S3CH eastern boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S3CH_east, scalar, series'
  'S3CH_time'
  'idTbry(ieast,iS3CH)'
  'nulvar'
  1.0d0

'S3CH_west'                                         ! Input
  'S3CH western boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S3CH_west, scalar, series'
  'S3CH_time'
  'idTbry(iwest,iS3CH)'
  'nulvar'
  1.0d0

'S3CH_south'                                        ! Input
  'S3CH southern boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S3CH_south, scalar, series'
  'S3CH_time'
  'idTbry(isouth,iS3CH)'
  'nulvar'
  1.0d0

'S3CH_north'                                        ! Input
  'S3CH northern boundary condition'
  'milligram meter-3'                     ! [milligram/m3]
  'S3CH_north, scalar, series'
  'S3CH_time'
  'idTbry(inorth,iS3CH)'
  'nulvar'
  1.0d0

'Z1_N_east'                                         ! Input
  'Z1_N eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_N_east, scalar, series'
  'Z1_N_time'
  'idTbry(ieast,iZ1_N)'
  'nulvar'
  1.0d0

'Z1_N_west'                                         ! Input
  'Z1_N western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_N_west, scalar, series'
  'Z1_N_time'
  'idTbry(iwest,iZ1_N)'
  'nulvar'
  1.0d0

'Z1_N_south'                                        ! Input
  'Z1_N southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_N_south, scalar, series'
  'Z1_N_time'
  'idTbry(isouth,iZ1_N)'
  'nulvar'
  1.0d0

'Z1_N_north'                                        ! Input
  'Z1_N northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_N_north, scalar, series'
  'Z1_N_time'
  'idTbry(inorth,iZ1_N)'
  'nulvar'
  1.0d0

'Z1_C_east'                                         ! Input
  'Z1_C eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_C_east, scalar, series'
  'Z1_C_time'
  'idTbry(ieast,iZ1_C)'
  'nulvar'
  1.0d0

'Z1_C_west'                                         ! Input
  'Z1_C western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_C_west, scalar, series'
  'Z1_C_time'
  'idTbry(iwest,iZ1_C)'
  'nulvar'
  1.0d0

'Z1_C_south'                                        ! Input
  'Z1_C southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_C_south, scalar, series'
  'Z1_C_time'
  'idTbry(isouth,iZ1_C)'
  'nulvar'
  1.0d0

'Z1_C_north'                                        ! Input
  'Z1_C northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z1_C_north, scalar, series'
  'Z1_C_time'
  'idTbry(inorth,iZ1_C)'
  'nulvar'
  1.0d0

'Z2_N_east'                                         ! Input
  'Z2_N eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_N_east, scalar, series'
  'Z2_N_time'
  'idTbry(ieast,iZ2_N)'
  'nulvar'
  1.0d0

'Z2_N_west'                                         ! Input
  'Z2_N western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_N_west, scalar, series'
  'Z2_N_time'
  'idTbry(iwest,iZ2_N)'
  'nulvar'
  1.0d0

'Z2_N_south'                                        ! Input
  'Z2_N southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_N_south, scalar, series'
  'Z2_N_time'
  'idTbry(isouth,iZ2_N)'
  'nulvar'
  1.0d0

'Z2_N_north'                                        ! Input
  'Z2_N northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_N_north, scalar, series'
  'Z2_N_time'
  'idTbry(inorth,iZ2_N)'
  'nulvar'
  1.0d0

'Z2_C_east'                                         ! Input
  'Z2_C eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_C_east, scalar, series'
  'Z2_C_time'
  'idTbry(ieast,iZ2_C)'
  'nulvar'
  1.0d0

'Z2_C_west'                                         ! Input
  'Z2_C western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_C_west, scalar, series'
  'Z2_C_time'
  'idTbry(iwest,iZ2_C)'
  'nulvar'
  1.0d0

'Z2_C_south'                                        ! Input
  'Z2_C southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_C_south, scalar, series'
  'Z2_C_time'
  'idTbry(isouth,iZ2_C)'
  'nulvar'
  1.0d0

'Z2_C_north'                                        ! Input
  'Z2_C northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'Z2_C_north, scalar, series'
  'Z2_C_time'
  'idTbry(inorth,iZ2_C)'
  'nulvar'
  1.0d0

'DD_N_east'                                         ! Input
  'DD_N eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_N_east, scalar, series'
  'DD_N_time'
  'idTbry(ieast,iDD_N)'
  'nulvar'
  1.0d0

'DD_N_west'                                         ! Input
  'DD_N western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_N_west, scalar, series'
  'DD_N_time'
  'idTbry(iwest,iDD_N)'
  'nulvar'
  1.0d0

'DD_N_south'                                        ! Input
  'DD_N southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_N_south, scalar, series'
  'DD_N_time'
  'idTbry(isouth,iDD_N)'
  'nulvar'
  1.0d0

'DD_N_north'                                        ! Input
  'DD_N northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_N_north, scalar, series'
  'DD_N_time'
  'idTbry(inorth,iDD_N)'
  'nulvar'
  1.0d0

'DD_C_east'                                         ! Input
  'DD_C eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_C_east, scalar, series'
  'DD_C_time'
  'idTbry(ieast,iDD_C)'
  'nulvar'
  1.0d0

'DD_C_west'                                         ! Input
  'DD_C western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_C_west, scalar, series'
  'DD_C_time'
  'idTbry(iwest,iDD_C)'
  'nulvar'
  1.0d0

'DD_C_south'                                        ! Input
  'DD_C southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_C_south, scalar, series'
  'DD_C_time'
  'idTbry(isouth,iDD_C)'
  'nulvar'
  1.0d0

'DD_C_north'                                        ! Input
  'DD_C northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DD_C_north, scalar, series'
  'DD_C_time'
  'idTbry(inorth,iDD_C)'
  'nulvar'
  1.0d0

'BAC_east'                                         ! Input
  'BAC eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'BAC_east, scalar, series'
  'BAC_time'
  'idTbry(ieast,iBAC_)'
  'nulvar'
  1.0d0

'BAC_west'                                         ! Input
  'BAC western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'BAC_west, scalar, series'
  'BAC_time'
  'idTbry(iwest,iBAC_)'
  'nulvar'
  1.0d0

'BAC_south'                                        ! Input
  'BAC southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'BAC_south, scalar, series'
  'BAC_time'
  'idTbry(isouth,iBAC_)'
  'nulvar'
  1.0d0

'BAC_north'                                        ! Input
  'BAC northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'BAC_north, scalar, series'
  'BAC_time'
  'idTbry(inorth,iBAC_)'
  'nulvar'
  1.0d0

'DDCA_east'                                         ! Input
  'DDCA eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DDCA_east, scalar, series'
  'DDCA_time'
  'idTbry(ieast,iDDCA)'
  'nulvar'
  1.0d0

'DDCA_west'                                         ! Input
  'DDCA western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DDCA_west, scalar, series'
  'DDCA_time'
  'idTbry(iwest,iDDCA)'
  'nulvar'
  1.0d0

'DDCA_south'                                        ! Input
  'DDCA southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DDCA_south, scalar, series'
  'DDCA_time'
  'idTbry(isouth,iDDCA)'
  'nulvar'
  1.0d0

'DDCA_north'                                        ! Input
  'DDCA northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DDCA_north, scalar, series'
  'DDCA_time'
  'idTbry(inorth,iDDCA)'
  'nulvar'
  1.0d0

'DDSi_east'                                         ! Input
  'DDSi eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DDSi_east, scalar, series'
  'DDSi_time'
  'idTbry(ieast,iDDSi)'
  'nulvar'
  1.0d0

'DDSi_west'                                         ! Input
  'DDSi western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DDSi_west, scalar, series'
  'DDSi_time'
  'idTbry(iwest,iDDSi)'
  'nulvar'
  1.0d0

'DDSi_south'                                        ! Input
  'DDSi southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DDSi_south, scalar, series'
  'DDSi_time'
  'idTbry(isouth,iDDSi)'
  'nulvar'
  1.0d0

'DDSi_north'                                        ! Input
  'DDSi northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'DDSi_north, scalar, series'
  'DDSi_time'
  'idTbry(inorth,iDDSi)'
  'nulvar'
  1.0d0

'LDON_east'                                         ! Input
  'LDON eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'LDON_east, scalar, series'
  'LDON_time'
  'idTbry(ieast,iLDON)'
  'nulvar'
  1.0d0

'LDON_west'                                         ! Input
  'LDON western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'LDON_west, scalar, series'
  'LDON_time'
  'idTbry(iwest,iLDON)'
  'nulvar'
  1.0d0

'LDON_south'                                        ! Input
  'LDON southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'LDON_south, scalar, series'
  'LDON_time'
  'idTbry(isouth,iLDON)'
  'nulvar'
  1.0d0

'LDON_north'                                        ! Input
  'LDON northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'LDON_north, scalar, series'
  'LDON_time'
  'idTbry(inorth,iLDON)'
  'nulvar'
  1.0d0

'LDOC_east'                                         ! Input
  'LDOC eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'LDOC_east, scalar, series'
  'LDOC_time'
  'idTbry(ieast,iLDOC)'
  'nulvar'
  1.0d0

'LDOC_west'                                         ! Input
  'LDOC western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'LDOC_west, scalar, series'
  'LDOC_time'
  'idTbry(iwest,iLDOC)'
  'nulvar'
  1.0d0

'LDOC_south'                                        ! Input
  'LDOC southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'LDOC_south, scalar, series'
  'LDOC_time'
  'idTbry(isouth,iLDOC)'
  'nulvar'
  1.0d0

'LDOC_north'                                        ! Input
  'LDOC northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'LDOC_north, scalar, series'
  'LDOC_time'
  'idTbry(inorth,iLDOC)'
  'nulvar'
  1.0d0

'SDON_east'                                         ! Input
  'SDON eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'SDON_east, scalar, series'
  'SDON_time'
  'idTbry(ieast,iSDON)'
  'nulvar'
  1.0d0

'SDON_west'                                         ! Input
  'SDON western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'SDON_west, scalar, series'
  'SDON_time'
  'idTbry(iwest,iSDON)'
  'nulvar'
  1.0d0

'SDON_south'                                        ! Input
  'SDON southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'SDON_south, scalar, series'
  'SDON_time'
  'idTbry(isouth,iSDON)'
  'nulvar'
  1.0d0

'SDON_north'                                        ! Input
  'SDON northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'SDON_north, scalar, series'
  'SDON_time'
  'idTbry(inorth,iSDON)'
  'nulvar'
  1.0d0

'SDOC_east'                                         ! Input
  'SDOC eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'SDOC_east, scalar, series'
  'SDOC_time'
  'idTbry(ieast,iSDOC)'
  'nulvar'
  1.0d0

'SDOC_west'                                         ! Input
  'SDOC western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'SDOC_west, scalar, series'
  'SDOC_time'
  'idTbry(iwest,iSDOC)'
  'nulvar'
  1.0d0

'SDOC_south'                                        ! Input
  'SDOC southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'SDOC_south, scalar, series'
  'SDOC_time'
  'idTbry(isouth,iSDOC)'
  'nulvar'
  1.0d0

'SDOC_north'                                        ! Input
  'SDOC northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'SDOC_north, scalar, series'
  'SDOC_time'
  'idTbry(inorth,iSDOC)'
  'nulvar'
  1.0d0

'CLDC_east'                                         ! Input
  'CLDC eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'CLDC_east, scalar, series'
  'CLDC_time'
  'idTbry(ieast,iCLDC)'
  'nulvar'
  1.0d0

'CLDC_west'                                         ! Input
  'CLDC western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'CLDC_west, scalar, series'
  'CLDC_time'
  'idTbry(iwest,iCLDC)'
  'nulvar'
  1.0d0

'CLDC_south'                                        ! Input
  'CLDC southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'CLDC_south, scalar, series'
  'CLDC_time'
  'idTbry(isouth,iCLDC)'
  'nulvar'
  1.0d0

'CLDC_north'                                        ! Input
  'CLDC northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'CLDC_north, scalar, series'
  'CLDC_time'
  'idTbry(inorth,iCLDC)'
  'nulvar'
  1.0d0

'CSDC_east'                                         ! Input
  'CSDC eastern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'CSDC_east, scalar, series'
  'CSDC_time'
  'idTbry(ieast,iCSDC)'
  'nulvar'
  1.0d0

'CSDC_west'                                         ! Input
  'CSDC western boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'CSDC_west, scalar, series'
  'CSDC_time'
  'idTbry(iwest,iCSDC)'
  'nulvar'
  1.0d0

'CSDC_south'                                        ! Input
  'CSDC southern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'CSDC_south, scalar, series'
  'CSDC_time'
  'idTbry(isouth,iCSDC)'
  'nulvar'
  1.0d0

'CSDC_north'                                        ! Input
  'CSDC northern boundary condition'
  'millimole meter-3'                     ! [millimole/m3]
  'CSDC_north, scalar, series'
  'CSDC_time'
  'idTbry(inorth,iCSDC)'
  'nulvar'
  1.0d0

'NO3_east'                                         ! Input
  'nitrate eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3_east, scalar, series'
  'NO3_time'
  'idTbry(ieast,iNO3_)'
  'nulvar'
  1.0d0

'NO3_west'                                         ! Input
  'nitrate western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3_west, scalar, series'
  'NO3_time'
  'idTbry(iwest,iNO3_)'
  'nulvar'
  1.0d0

'NO3_south'                                        ! Input
  'nitrate southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3_south, scalar, series'
  'NO3_time'
  'idTbry(isouth,iNO3_)'
  'nulvar'
  1.0d0

'NO3_north'                                        ! Input
  'nitrate northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NO3_north, scalar, series'
  'NO3_time'
  'idTbry(inorth,iNO3_)'
  'nulvar'
  1.0d0

'PO4_east'                                         ! Input
  'phosphate eastern boundary condition'
  'millimole_P meter-3'                     ! [millimole/m3]
  'PO4_east, scalar, series'
  'PO4_time'
  'idTbry(ieast,iPO4_)'
  'nulvar'
  1.0d0

'PO4_west'                                         ! Input
  'phosphate western boundary condition'
  'millimole_P meter-3'                     ! [millimole/m3]
  'PO4_west, scalar, series'
  'PO4_time'
  'idTbry(iwest,iPO4_)'
  'nulvar'
  1.0d0

'PO4_south'                                        ! Input
  'phosphate southern boundary condition'
  'millimole_P meter-3'                     ! [millimole/m3]
  'PO4_south, scalar, series'
  'PO4_time'
  'idTbry(isouth,iPO4_)'
  'nulvar'
  1.0d0

'PO4_north'                                        ! Input
  'phosphate northern boundary condition'
  'millimole_P meter-3'                     ! [millimole/m3]
  'PO4_north, scalar, series'
  'PO4_time'
  'idTbry(inorth,iPO4_)'
  'nulvar'
  1.0d0

'TIC_east'                                         ! Input
  'TIC eastern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC_east, scalar, series'
  'TIC_time'
  'idTbry(ieast,iTIC_)'
  'nulvar'
  1.0d0

'TIC_west'                                         ! Input
  'TIC western boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC_west, scalar, series'
  'TIC_time'
  'idTbry(iwest,iTIC_)'
  'nulvar'
  1.0d0

'TIC_north'                                        ! Input
  'TIC northern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC_north, scalar, series'
  'TIC_time'
  'idTbry(inorth,iTIC_)'
  'nulvar'
  1.0d0

'TIC_south'                                        ! Input
  'TIC southern boundary condition'
  'millimole_carbon meter-3'                       ! [millimole/m3]
  'TIC_south, scalar, series'
  'TIC_time'
  'idTbry(isouth,iTIC_)'
  'nulvar'
  1.0d0

'Oxyg_east'                                      ! Input
  'oxygen eastern boundary condition'
  'millimole_O2 meter-3'                       ! [millimole/m3]
  'oxygen_east, scalar, series'
  'oxygen_time'
  'idTbry(ieast,iOxyg)'
  'nulvar'
  1.0d0

'Oxyg_west'                                      ! Input
  'oxygen western boundary condition'
  'millimole_O2 meter-3'                       ! [millimole/m3]
  'oxygen_west, scalar, series'
  'oxygen_time'
  'idTbry(iwest,iOxyg)'
  'nulvar'
  1.0d0

'Oxyg_north'                                     ! Input
  'oxygen northern boundary condition'
  'millimole_O2 meter-3'                       ! [millimole/m3]
  'oxygen_north, scalar, series'
  'oxygen_time'
  'idTbry(inorth,iOxyg)'
  'nulvar'
  1.0d0

'Oxyg_south'                                     ! Input
  'oxygen southern boundary condition'
  'millimole_O2 meter-3'                       ! [millimole/m3]
  'oxygen_south, scalar, series'
  'oxygen_time'
  'idTbry(isouth,iOxyg)'
  'nulvar'
  1.0d0

'Talk_east'                                  ! Input
  'alkalinity eastern boundary condition'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity_east, scalar, series'
  'alkalinity_time'
  'idTbry(ieast,iTAlk)'
  'nulvar'
  1.0d0

'Talk_west'                                  ! Input
  'alkalinity western boundary condition'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity_west, scalar, series'
  'alkalinity_time'
  'idTbry(iwest,iTAlk)'
  'nulvar'
  1.0d0

'Talk_north'                                 ! Input
  'alkalinity northern boundary condition'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity_north, scalar, series'
  'alkalinity_time'
  'idTbry(inorth,iTAlk)'
  'nulvar'
  1.0d0

'Talk_south'                                 ! Input
  'alkalinity southern boundary condition'
  'milliequivalents meter-3'                       ! [milliequivalents/m3]
  'alkalinity_south, scalar, series'
  'alkalinity_time'
  'idTbry(isouth,iTAlk)'
  'nulvar'
  1.0d0

'NH4_east'                                         ! Input
  'ammonium eastern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4_east, scalar, series'
  'NH4_time'
  'idTbry(ieast,iNH4_)'
  'nulvar'
  1.0d0

'NH4_west'                                         ! Input
  'ammonium western boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4_west, scalar, series'
  'NH4_time'
  'idTbry(iwest,iNH4_)'
  'nulvar'
  1.0d0

'NH4_south'                                        ! Input
  'ammonium southern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4_south, scalar, series'
  'NH4_time'
  'idTbry(isouth,iNH4_)'
  'nulvar'
  1.0d0

'NH4_north'                                        ! Input
  'ammonium northern boundary condition'
  'millimole_nitrogen meter-3'                     ! [millimole/m3]
  'NH4_north, scalar, series'
  'NH4_time'
  'idTbry(inorth,iNH4_)'
  'nulvar'
  1.0d0

  'SiOH_east'                                       ! Input
  'silicate eastern boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH_east, scalar, series'
  'SiOH_time'
  'idTbry(ieast,iSiOH)'
  'nulvar'
  1.0d0

'SiOH_west'                                       ! Input
  'silicate western boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH_west, scalar, series'
  'SiOH_time'
  'idTbry(iwest,iSiOH)'
  'nulvar'
  1.0d0

'SiOH_south'                                      ! Input
  'silicate southern boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH_south, scalar, series'
  'SiOH_time'
  'idTbry(isouth,iSiOH)'
  'nulvar'
  1.0d0

'SiOH_north'                                      ! Input
  'silicate northern boundary condition'
  'millimole_silica meter-3'                       ! [millimole/m3]
  'SiOH_north, scalar, series'
  'SiOH_time'
  'idTbry(inorth,iSiOH)'
  'nulvar'
  1.0d0
!---------- end UMAINE variables
!
!-----------------------------------------------------------------------
!  Define names of variables associated with EcoSim bio-optical module.
!-----------------------------------------------------------------------
!

'NO3'                                              ! Input/Output
  'nitrate concentration'
  'millimole_N03 meter-3'                          ! [millimole/m3]
  'NO3, scalar, series'
  'ocean_time'
  'idTvar(iNO3_)'
  'r3dvar'
  1.0d0

'NH4'                                              ! Input/Output
  'ammonium concentration'
  'millimole_NH4 meter-3'                          ! [millimole/m3]
  'NH4, scalar, series'
  'ocean_time'
  'idTvar(iNH4_)'
  'r3dvar'
  1.0d0

'PO4'                                              ! Input/Output
  'phosphate concentration'
  'millimole_PO4 meter-3'                          ! [millimole/m3]
  'PO4, scalar, series'
  'ocean_time'
  'idTvar(iPO4_)'
  'r3dvar'
  1.0d0

'FeO'                                              ! Input/Output
  'iron concentration'
  'millimole_FeO meter-3'                          ! [millimole/m3]
  'FeO, scalar, series'
  'ocean_time'
  'idTvar(iFeO_)'
  'r3dvar'
  1.0d0

'SiO'                                              ! Input/Output
  'silica concentration'
  'millimole_SiO meter-3'                          ! [millimole/m3]
  'SiO, scalar, series'
  'ocean_time'
  'idTvar(iSiO_)'
  'r3dvar'
  1.0d0

'DIC'                                              ! Input/Output
  'DIC concentration'
  'millimole_DIC meter-3'                          ! [millimole/m3]
  'DIC concentration, scalar, series'
  'ocean_time'
  'idTvar(iDIC_)'
  'r3dvar'
  1.0d0

'Bac_C'                                            ! Input/Output
  'bacteria, carbon group'
  'millimole_C meter-3'                            ! [millimole/m3]
  'Bac_C, scalar, series'
  'ocean_time'
  'idTvar(iBacC)'
  'r3dvar'
  1.0d0

'Bac_F'                                            ! Input/Output
  'bacteria, iron Group'
  'millimole_Fe meter-3'                           ! [millimole/m3]
  'Bac_F, scalar, series'
  'ocean_time'
  'idTvar(iBacF)'
  'r3dvar'
  1.0d0

'Bac_N'                                            ! Input/Output
  'bacteria, nitrogen group'
  'millimole_N meter-3'                            ! [millimole/m3]
  'Bac_N, scalar, series'
  'ocean_time'
  'idTvar(iBacN)'
  'r3dvar'
  1.0d0

'Bac_P'                                            ! Input/Output
  'bacteria, phosphorus group'
  'millimole_P meter-3'                            ! [millimole/m3]
  'Bac_P, scalar, series'
  'ocean_time'
  'idTvar(iBacP)'
  'r3dvar'
  1.0d0

'CDM_C'                                            ! Input/Output
  'color degradational matter'
  'millimole_CDMC meter-3'                         ! [millimole/m3]
  'CDM_C, scalar, series'
  'ocean_time'
  'idTvar(iCDMC)'
  'r3dvar'
  1.0d0

'DOM_C'                                            ! Input/Output
  'dissolved organic matter, carbon group'
  'millimole_C meter-3'                            ! [millimole/m3]
  'DOM_C, scalar, series'
  'ocean_time'
  'idTvar(iDOMC)'
  'r3dvar'
  1.0d0

'DOM_N'                                            ! Input/Output
  'dissolved organic matter, nitrogen group'
  'millimole_N meter-3'                            ! [millimole/m3]
  'DOM_N, scalar, series'
  'ocean_time'
  'idTvar(iDOMN)'
  'r3dvar'
  1.0d0

'DOM_P'                                            ! Input/Output
  'dissolved organic matter,) phosphorus group'
  'millimole_P meter-3'                            ! [millimole/m3]
  'DOM_P, scalar, series'
  'ocean_time'
  'idTvar(iDOMP)'
  'r3dvar'
  1.0d0

'Fec_C'                                            ! Input/Output
  'fecal matter, carbon group'
  'millimole_C meter-3'                            ! [millimole/m3]
  'Fec_C, scalar, series'
  'ocean_time'
  'idTvar(iFecC)'
  'r3dvar'
  1.0d0

'Fec_F'                                            ! Input/Output
  'fecal matter, iron group'
  'millimole_Fe meter-3'                           ! [millimole/m3]
  'Fec_F, scalar, series'
  'ocean_time'
  'idTvar(iFecF)'
  'r3dvar'
  1.0d0

'Fec_N'                                            ! Input/Output
  'fecal matter, nitrogen group'
  'millimole_N meter-3'                            ! [millimole/m3]
  'Fec_N, scalar, series'
  'ocean_time'
  'idTvar(iFecN)'
  'r3dvar'
  1.0d0

'Fec_P'                                            ! Input/Output
  'fecal matter, phosphorus group'
  'millimole_P meter-3'                            ! [millimole/m3]
  'Fec_P, scalar, series'
  'ocean_time'
  'idTvar(iFecP)'
  'r3dvar'
  1.0d0

'Fec_S'                                            ! Input/Output
  'fecal matter, silica group'
  'millimole_SiO meter-3'                          ! [millimole/m3]
  'Fec_S, scalar, series'
  'ocean_time'
  'idTvar(iFecS)'
  'r3dvar'
  1.0d0

'Phy_C'                                            ! Input/Output
  'phytoplankton, carbon group'
  'millimole_C meter-3'                            ! [millimole/m3]
  'Phy_C, scalar, series'
  'ocean_time'
  'idTvar(iPhyC)'
  'r3dvar'
  1.0d0

'Phy_F'                                            ! Input/Output
  'phytoplankton, iron group'
  'millimole_Fe meter-3'                           ! [millimole/m3]
  'Phy_F, scalar, series'
  'ocean_time'
  'idTvar(iPhyF)'
  'r3dvar'
  1.0d0

'Phy_N'                                            ! Input/Output
  'phytoplankton, nitrogen group'
  'millimole_N meter-3'                            ! [millimole/m3]
  'Phy_N, scalar, series'
  'ocean_time'
  'idTvar(iPhyN)'
  'r3dvar'
  1.0d0

'Phy_P'                                            ! Input/Output
  'phytoplankton, phosphorus group'
  'millimole_P meter-3'                            ! [millimole/m3]
  'Phy_P, scalar, series'
  'ocean_time'
  'idTvar(iPhyP)'
  'r3dvar'
  1.0d0

'Phy_S'                                            ! Input/Output
  'phytoplankton, silica group'
  'millimole_SiO meter-3'                          ! [millimole/m3]
  'Phy_S, scalar, series'
  'ocean_time'
  'idTvar(iPhyS)'
  'r3dvar'
  1.0d0

'Pigs_XY'                                          ! Input/Output
  'Phytoplankton X, pigment Y'
  'milligrams meter-3'                             ! [milligrams/m3]
  'Pigs_XY, scalar, series'
  'ocean_time'
  'idTvar(iPigs)'
  'r3dvar'
  1.0d0

!
!-----------------------------------------------------------------------
!  Define names of variables for time-averaged quadratic terms.
!-----------------------------------------------------------------------
!

'zeta2'                                            ! Output
  'time-averaged squared free-surface'
  'meter2'                                         ! [m2]
  'zeta2, scalar, series'
  'ocean_time'
  'idZZav'
  'r2dvar'
  1.0d0

'ubar2'                                            ! Output
  'time-averaged squared vertically integrated u-momentum'
  'meter2 second-2'                                ! [m2/s2]
  'ubar2, scalar, series'
  'ocean_time'
  'idU2av'
  'u2dvar'
  1.0d0

'vbar2'                                            ! Output
  'time-averaged squared vertically integrated v-momentum'
  'meter2 second-2'                                ! [m2/s2]
  'vbar2, scalar, series'
  'ocean_time'
  'idV2av'
  'v2dvar'
  1.0d0

'Huon'                                             ! Output
  'time-averaged u-volume flux'
  'meter3 second-1'                                ! [m3/s]
  'Huon, scalar, series'
  'ocean_time'
  'idHUav'
  'u3dvar'
  1.0d0

'Hvom'                                             ! Output
  'time-averaged v-volume flux'
  'meter3 second-1'                                ! [m3/s]
  'Hvom, scalar, series'
  'ocean_time'
  'idHVav'
  'v3dvar'
  1.0d0

'uu'                                               ! Output
  'time-averaged u-momentum times u-momentum'
  'meter2 second-2'                                ! [m2/s2]
  'uu, scalar, series'
  'ocean_time'
  'idUUav'
  'u3dvar'
  1.0d0

'uv'                                               ! Output
  'time-averaged u-momentum times v-momentum'
  'meter2 second-2'                                ! [m2/s2]
  'uv, scalar, series'
  'ocean_time'
  'idUVav'
  'r3dvar'
  1.0d0

'vv'                                               ! Output
  'time-averaged v-momentum times v-momentum'
  'meter2 second-2'                                ! [m2/s2]
  'uv, scalar, series'
  'ocean_time'
  'idVVav'
  'v3dvar'
  1.0d0

'speed'                                            ! Output
  'time-averaged speed'
  'meter2 second-2'                                ! [m2/s2]
  'speed, scalar, series'
  'ocean_time'
  'idSUVav'
  'r3dvar'
  1.0d0

'tracer2'                                          ! Output
  'time-averaged squared tracer'
  'tracer2'                                        ! [tracer2]
  'tracer2, scalar, series'
  'ocean_time'
  'idTTav(itrc)'
  'r3dvar'
  1.0d0

'Huontracer'                                       ! Output
  'time-averaged tracer u-volume flux'
  'meter3 second-1 tracer'                         ! [m3/s tracer]
  'Huontracer, scalar, series'
  'ocean_time'
  'iHUTav(itrc)'
  'u3dvar'
 1.0d0

'utracer'                                          ! Output
  'time-averaged u-momentum times tracer'
  'meter second-1 tracer'                          ! [m/s tracer]
  'utracer, scalar, series'
  'ocean_time'
  'idUTav(itrc)'
  'u3dvar'
  1.0d0

'Hvomtracer'                                       ! Output
  'time-averaged tracer v-volume flux'
  'meter3 second-1 tracer'                         ! [m3/s tracer]
  'Hvomtracer, scalar, series'
  'ocean_time'
  'iHVTav(itrc)'
  'v3dvar'
 1.0d0

'vtracer'                                          ! Output
  'time-averaged v-momentum times tracer'
  'meter second-1 tracer'                          ! [m/s tracer]
  'vtracer, scalar, series'
  'ocean_time'
  'idVTav(itrc)'
  'v3dvar'
  1.0d0

!
!-----------------------------------------------------------------------
!  Define names of variables for diagnostic fields.
!-----------------------------------------------------------------------
!
!  2D momentum diagnostic terms.
!

'pvorticity_bar'                                   ! Input/Output
  '2D potential vorticity'
  'meter-1 second-1'                               ! [1/m/s]
  'pvorticity_bar, scalar, series'
  'ocean_time'
  'id2dPV'
  'p2dvar'
  1.0d0

'rvorticity_bar'                                   ! Input/Output
  '2D relative vorticity'
  'second-1'                                       ! [1/s]
  'rvorticity_bar, scalar, series'
  'ocean_time'
  'id2dRV'
  'p2dvar'
  1.0d0

'ubar_accel'                                       ! Input/Output
  '2D u-momentum, acceleration term'
  'meter second-2'                                 ! [m/s2]
  'ubar_accel, scalar, series'
  'ocean_time'
  'idDu2d(M2rate)'
  'u2dvar'
  1.0d0

'vbar_accel'                                       ! Input/Output
  '2D v-momentum, acceleration term'
  'meter second-2'                                 ! [m/s2]
  'vbar_accel, scalar, series'
  'ocean_time'
  'idDv2d(M2rate)'
  'v2dvar'
  1.0d0

'ubar_prsgrd'                                      ! Input/Output
  '2D u-momentum, pressure gradient term'
  'meter second-2'                                 ! [m/s2]
  'ubar_prsgrd, scalar, series'
  'ocean_time'
  'idDu2d(M2pgrd)'
  'u2dvar'
  1.0d0

'vbar_prsgrd'                                      ! Input/Output
  '2D v-momentum, pressure gradient term'
  'meter second-2'                                 ! [m/s2]
  'vbar_prsgrd, scalar, series'
  'ocean_time'
  'idDv2d(M2pgrd)'
  'v2dvar'
  1.0d0

'ubar_cor'                                         ! Input/Output
  '2D u-momentum, Coriolis term'
  'meter second-2'                                 ! [m/s2]
  'ubar_cor, scalar, series'
  'ocean_time'
  'idDu2d(M2fcor)'
  'u2dvar'
  1.0d0

'vbar_cor'                                         ! Input/Output
  '2D v-momentum, Coriolis term'
  'meter second-2'                                 ! [m/s2]
  'vbar_cor, scalar, series'
  'ocean_time'
  'idDv2d(M2fcor)'
  'v2dvar'
  1.0d0

'ubar_hadv'                                        ! Input/Output
  '2D u-momentum, horizontal advection term'
  'meter second-2'                                 ! [m/s2]
  'ubar_hadv, scalar, series'
  'ocean_time'
  'idDu2d(M2hadv)'
  'u2dvar'
  1.0d0

'ubar_xadv'                                        ! Input/Output
  '2D u-momentum, horizontal XI-advection term'
  'meter second-2'                                 ! [m/s2]
  'ubar_xadv, scalar, series'
  'ocean_time'
  'idDu2d(M2xadv)'
  'u2dvar'
  1.0d0

'ubar_yadv'                                        ! Input/Output
  '2D u-momentum, horizontal ETA-advection term'
  'meter second-2'                                 ! [m/s2]
  'ubar_yadv, scalar, series'
  'ocean_time'
  'idDu2d(M2yadv)'
  'u2dvar'
  1.0d0

'vbar_hadv'                                        ! Input/Output
  '2D v-momentum, horizontal advection term'
  'meter second-2'                                 ! [m/s2]
  'vbar_hadv, scalar, series'
  'ocean_time'
  'idDv2d(M2hadv)'
  'v2dvar'
  1.0d0

'vbar_xadv'                                        ! Input/Output
  '2D v-momentum, horizontal XI-advection term'
  'meter second-2'                                 ! [m/s2]
  'vbar_xadv, scalar, series'
  'ocean_time'
  'idDv2d(M2xadv)'
  'v2dvar'
  1.0d0

'vbar_yadv'                                        ! Input/Output
  '2D v-momentum, horizontal ETA-advection term'
  'meter second-2'                                 ! [m/s2]
  'vbar_yadv, scalar, series'
  'ocean_time'
  'idDv2d(M2yadv)'
  'v2dvar'
  1.0d0

'ubar_hrad'                                        ! Input/Output
  '2D u-momentum, horizontal radiation stress term'
  'meter second-2'                                 ! [m/s2]
  'ubar_hrad, scalar, series'
  'ocean_time'
  'idDu2d(M2hrad)'
  'u2dvar'
  1.0d0

'vbar_hrad'                                        ! Input/Output
  '2D v-momentum, horizontal radiation stress term'
  'meter second-2'                                 ! [m/s2]
  'vbar_hrad, scalar, series'
  'ocean_time'
  'idDv2d(M2hrad)'
  'v2dvar'
  1.0d0

'ubar_hvisc'                                       ! Input/Output
  '2D u-momentum, horizontal viscosity term'
  'meter second-2'                                 ! [m/s2]
  'ubar_hvisc, scalar, series'
  'ocean_time'
  'idDu2d(M2hvis)'
  'u2dvar'
  1.0d0

'ubar_xvisc'                                       ! Input/Output
  '2D u-momentum, horizontal XI-viscosity term'
  'meter second-2'                                 ! [m/s2]
  'ubar_xvisc, scalar, series'
  'ocean_time'
  'idDu2d(M2xvis)'
  'u2dvar'
  1.0d0

'ubar_yvisc'                                       ! Input/Output
  '2D u-momentum, horizontal ETA-viscosity term'
  'meter second-2'                                 ! [m/s2]
  'ubar_yvisc, scalar, series'
  'ocean_time'
  'idDu2d(M2yvis)'
  'u2dvar'
  1.0d0

'vbar_hvisc'                                       ! Input/Output
  '2D v-momentum, horizontal viscosity term'
  'meter second-2'                                 ! [m/s2]
  'vbar_hvisc, scalar, series'
  'ocean_time'
  'idDv2d(M2hvis)'
  'v2dvar'
  1.0d0

'vbar_xvisc'                                       ! Input/Output
  '2D v-momentum, horizontal XI-viscosity term'
  'meter second-2'                                 ! [m/s2]
  'vbar_xvisc, scalar, series'
  'ocean_time'
  'idDv2d(M2xvis)'
  'v2dvar'
  1.0d0

'vbar_yvisc'                                       ! Input/Output
  '2D v-momentum, horizontal ETA-viscosity term'
  'meter second-2'                                 ! [m/s2]
  'vbar_yvisc, scalar, series'
  'ocean_time'
  'idDv2d(M2yvis)'
  'v2dvar'
  1.0d0

'ubar_sstr'                                        ! Input/Output
  '2D u-momentum, surface stress term'
  'meter second-2'                                 ! [m/s2]
  'ubar_sstr, scalar, series'
  'ocean_time'
  'idDu2d(M2sstr)'
  'u2dvar'
  1.0d0

'ubar_bstr'                                        ! Input/Output
  '2D u-momentum, bottom stress term'
  'meter second-2'                                 ! [m/s2]
  'ubar_bstr, scalar, series'
  'ocean_time'
  'idDu2d(M2bstr)'
  'u2dvar'
  1.0d0

'vbar_sstr'                                        ! Input/Output
  '2D v-momentum, surface stress term'
  'meter second-2'                                 ! [m/s2]
  'vbar_sstr, scalar, series'
  'ocean_time'
  'idDv2d(M2sstr)'
  'v2dvar'
  1.0d0

'vbar_bstr'                                        ! Input/Output
  '2D v-momentum, bottom stress term'
  'meter second-2'                                 ! [m/s2]
  'vbar_bstr, scalar, series'
  'ocean_time'
  'idDv2d(M2bstr)'
  'v2dvar'
  1.0d0

!
!  3D momentum diagnostic terms.
!

'pvorticity'                                       ! Input/Output
  'potential vorticity'
  'meter-1 second-1'                               ! [1/m/s]
  'pvorticity, scalar, series'
  'ocean_time'
  'id3dPV'
  'p3dvar'
  1.0d0

'rvorticity'                                       ! Input/Output
  'relative vorticity'
  'second-1'                                       ! [1/s]
  'rvorticity, scalar, series'
  'ocean_time'
  'id3dRV'
  'p3dvar'
  1.0d0

'u_accel'                                          ! Input/Output
  '3D u-momentum, acceleration term'
  'meter second-2'                                 ! [m/s2]
  'u_accel, scalar, series'
  'ocean_time'
  'idDu3d(M3rate)'
  'u3dvar'
  1.0d0

'v_accel'                                          ! Input/Output
  '3D v-momentum, acceleration term'
  'meter second-2'                                 ! [m/s2]
  'v_accel, scalar, series'
  'ocean_time'
  'idDv3d(M3rate)'
  'v3dvar'
  1.0d0

'u_prsgrd'                                         ! Input/Output
  '3D u-momentum, pressure gradient term'
  'meter second-2'                                 ! [m/s2]
  'u_prsgrd, scalar, series'
  'ocean_time'
  'idDu3d(M3pgrd)'
  'u3dvar'
  1.0d0

'v_prsgrd'                                         ! Input/Output
  '3D v-momentum, pressure gradient term'
  'meter second-2'                                 ! [m/s2]
  'v_prsgrd, scalar, series'
  'ocean_time'
  'idDv3d(M3pgrd)'
  'v3dvar'
  1.0d0

'u_cor'                                            ! Input/Output
  '3D u-momentum, Coriolis term'
  'meter second-2'                                 ! [m/s2]
  'u_cor, scalar, series'
  'ocean_time'
  'idDu3d(M3fcor)'
  'u3dvar'
  1.0d0

'v_cor'                                            ! Input/Output
  '3D v-momentum, Coriolis term'
  'meter second-2'                                 ! [m/s2]
  'v_cor, scalar, series'
  'ocean_time'
  'idDv3d(M3fcor)'
  'v3dvar'
  1.0d0

'u_hadv'                                           ! Input/Output
  '3D u-momentum, horizontal advection term'
  'meter second-2'                                 ! [m/s2]
  'u_hadv, scalar, series'
  'ocean_time'
  'idDu3d(M3hadv)'
  'u3dvar'
  1.0d0

'u_xadv'                                           ! Input/Output
  '3D u-momentum, horizontal XI-advection term'
  'meter second-2'                                 ! [m/s2]
  'u_xadv, scalar, series'
  'ocean_time'
  'idDu3d(M3xadv)'
  'u3dvar'
  1.0d0

'u_yadv'                                           ! Input/Output
  '3D u-momentum, horizontal ETA-advection term'
  'meter second-2'                                 ! [m/s2]
  'u_yadv, scalar, series'
  'ocean_time'
  'idDu3d(M3yadv)'
  'u3dvar'
  1.0d0

'v_hadv'                                           ! Input/Output
  '3D v-momentum, horizontal advection term'
  'meter second-2'                                 ! [m/s2]
  'v_hadv, scalar, series'
  'ocean_time'
  'idDv3d(M3hadv)'
  'v3dvar'
  1.0d0

'v_xadv'                                           ! Input/Output
  '3D v-momentum, horizontal XI-advection term'
  'meter second-2'                                 ! [m/s2]
  'v_xadv, scalar, series'
  'ocean_time'
  'idDv3d(M3xadv)'
  'v3dvar'
  1.0d0

'v_yadv'                                           ! Input/Output
  '3D v-momentum, horizontal ETA-advection term'
  'meter second-2'                                 ! [m/s2]
  'v_yadv, scalar, series'
  'ocean_time'
  'idDv3d(M3yadv)'
  'v3dvar'
  1.0d0

'u_vadv'                                           ! Input/Output
  '3D u-momentum, vertical advection term'
  'meter second-2'                                 ! [m/s2]
  'u_vadv, scalar, series'
  'ocean_time'
  'idDu3d(M3vadv)'
  'u3dvar'
  1.0d0

'v_vadv'                                           ! Input/Output
  '3D v-momentum, vertical advection term'
  'meter second-2'                                 ! [m/2]
  'v_vadv, scalar, series'
  'ocean_time'
  'idDv3d(M3vadv)'
  'v3dvar'
  1.0d0

'u_hrad'                                           ! Input/Output
  '3D u-momentum, horizontal radiation stress  term'
  'meter second-2'                                 ! [m/s2]
  'u_hrad, scalar, series'
  'ocean_time'
  'idDu3d(M3hrad)'
  'u3dvar'
  1.0d0

'v_hrad'                                           ! Input/Output
  '3D v-momentum, horizontal radiation stress term'
  'meter second-2'                                 ! [m/s2]
  'v_hrad, scalar, series'
  'ocean_time'
  'idDv3d(M3hrad)'
  'v3dvar'
  1.0d0

'u_vrad'                                           ! Input/Output
  '3D u-momentum, vertical radiation stress term'
  'meter second-2'                                 ! [m/s2]
  'u_vrad, scalar, series'
  'ocean_time'
  'idDu3d(M3vrad)'
  'u3dvar'
  1.0d0

'v_vrad'                                           ! Input/Output
  '3D v-momentum, vertical radiation stress term'
  'meter second-2'                                 ! [m/2]
  'v_vrad, scalar, series'
  'ocean_time'
  'idDv3d(M3vrad)'
  'v3dvar'
  1.0d0

'u_hvisc'                                          ! Input/Output
  '3D u-momentum, horizontal viscosity term'
  'meter second-2'                                 ! [m/s2]
  'u_hvisc, scalar, series'
  'ocean_time'
  'idDu3d(M3hvis)'
  'u3dvar'
  1.0d0

'u_xvisc'                                          ! Input/Output
  '3D u-momentum, horizontal XI-viscosity term'
  'meter second-2'                                 ! [m/s2]
  'u_xvisc, scalar, series'
  'ocean_time'
  'idDu3d(M3xvis)'
  'u3dvar'
  1.0d0

'u_yvisc'                                          ! Input/Output
  '3D u-momentum, horizontal ETA-viscosity term'
  'meter second-2'                                 ! [m/s2]
  'u_yvisc, scalar, series'
  'ocean_time'
  'idDu3d(M3yvis)'
  'u3dvar'
  1.0d0

'v_hvisc'                                          ! Input/Output
  '3D v-momentum, horizontal viscosity term'
  'meter second-2'                                 ! [m/s2]
  'v_hvisc, scalar, series'
  'ocean_time'
  'idDv3d(M3hvis)'
  'v3dvar'
  1.0d0

'v_xvisc'                                          ! Input/Output
  '3D v-momentum, horizontal XI-viscosity term'
  'meter second-2'                                 ! [m/s2]
  'v_xvisc, scalar, series'
  'ocean_time'
  'idDv3d(M3xvis)'
  'v3dvar'
  1.0d0

'v_yvisc'                                          ! Input/Output
  '3D v-momentum, horizontal ETA-viscosity term'
  'meter second-2'                                 ! [m/s2]
  'v_yvisc, scalar, series'
  'ocean_time'
  'idDv3d(M3yvis)'
  'v3dvar'
  1.0d0

'u_vvisc'                                          ! Input/Output
  '3D u-momentum, vertical viscosity term'
  'meter second-2'                                 ! [m/s2]
  'u_vvisc, scalar, series'
  'ocean_time'
  'idDu3d(M3vvis)'
  'u3dvar'
  1.0d0

'v_vvisc'                                          ! Input/Output
  '3D v-momentum, vertical viscosity term'
  'meter second-2'                                 ! [m/s2]
  'v_vvisc, scalar, series'
  'ocean_time'
  'idDv3d(M3vvis)'
  'v3dvar'
  1.0d0

!
!  Tracer diagnostics terms. Notice that the variable names are generic
!  to allow string concatenation with the actual tracer name. These
!  entries MUST be entered in this file after all the active and
!  passive tracer variable information.
!

'_rate'                                            ! Input/Output
  'time rate of change'
  'second-1'                                       ! [1/s T]
  '_rate, scalar, series'
  'ocean_time'
  'idDtrc(iTrate)'
  'r3dvar'
  1.0d0

'_hadv'                                            ! Input/Output
  'horizontal advection term'
  'second-1'                                       ! [1/s T]
  '_hadv, scalar, series'
  'ocean_time'
  'idDtrc(iThadv)'
  'r3dvar'
  1.0d0

'_xadv'                                            ! Input/Output
  'horizontal XI-advection term'
  'second-1'                                       ! [1/s T]
  '_xadv, scalar, series'
  'ocean_time'
  'idDtrc(iTxadv)'
  'r3dvar'
  1.0d0

'_yadv'                                            ! Input/Output
  'horizontal ETA-advection term'
  'second-1'                                       ! [1/s T]
  '_yadv, scalar, series'
  'ocean_time'
  'idDtrc(iTyadv)'
  'r3dvar'
  1.0d0

'_vadv'                                            ! Input/Output
  'vertical advection term'
  'second-1'                                       ! [1/s T]
  '_vadv, scalar, series'
  'ocean_time'
  'idDtrc(iTvadv)'
  'r3dvar'
  1.0d0

'_hdiff'                                           ! Input/Output
  'horizontal diffusion term'
  'second-1'                                       ! [1/s T]
  '_hdiff, scalar, series'
  'ocean_time'
  'idDtrc(iThdif)'
  'r3dvar'
  1.0d0

'_xdiff'                                           ! Input/Output
  'horizontal XI-diffusion term'
  'second-1'                                       ! [1/s T]
  '_xdiff, scalar, series'
  'ocean_time'
  'idDtrc(iTxdif)'
  'r3dvar'
  1.0d0

'_ydiff'                                           ! Input/Output
  'horizontal ETA-diffusion term'
  'second-1'                                       ! [1/s T]
  '_ydiff, scalar, series'
  'ocean_time'
  'idDtrc(iTydif)'
  'r3dvar'
  1.0d0

'_sdiff'                                           ! Input/Output
  'horizontal S-diffusion rotated tensor term'
  'second-1'                                       ! [1/s T]
  '_sdiff, scalar, series'
  'ocean_time'
  'idDtrc(iTsdif)'
  'r3dvar'
  1.0d0

'_vdiff'                                           ! Input/Output
  'vertical diffusion term'
  'second-1'
  '_vdiff, scalar, series'
  'ocean_time'
  'idDtrc(iTvdif)'
  'r3dvar'
  1.0d0

!
!-----------------------------------------------------------------------
!  Forward trajectory variables for tangent linear and adjoint models.
!-----------------------------------------------------------------------
!

'rzeta'                                            ! Input/Output
  'RHS of free-surface equation'
  'meter3 second-1'                                ! [m3/s]
  'rzeta, scalar, series'
  'ocean_time'
  'idRzet'
  'r2dvar'
  1.0d0

'rubar'                                            ! Input/Output
  'RHS of vertically integrated u-momentum'
  'meter4 second-2'                                ! [m4/s2]
  'rubar, scalar, series'
  'ocean_time'
  'idRu2d'
  'u2dvar'
  1.0d0

'rvbar'                                            ! Input/Output
  'RHS of vertically integrated v-momentum'
  'meter4 second-2'                                ! [m4/s2]
  'rvbar, scalar, series'
  'ocean_time'
  'idRv2d'
  'v2dvar'
  1.0d0

'rufrc'                                            ! Input/Output
  'RHS of vertically integrated u-coupling terms'
  'meter4 second-2'                                ! [m4/s2]
  'rufrc, scalar, series'
  'ocean_time'
  'idRuct'
  'u2dvar'
  1.0d0

'rvfrc'                                            ! Input/Output
  'RHS of vertically integrated v-coupling terms'
  'meter4 second-2'                                ! [m4/s2]
  'rvfrc, scalar, series'
  'ocean_time'
  'idRvct'
  'v2dvar'
  1.0d0

'DU_avg1'                                          ! Input/Output
  'time averaged u-flux for 2D equations'
  'meter3 second-1'                                ! [m3/s]
  'DU_avg1, scalar, series'
  'ocean_time'
  'idUfx1'
  'u2dvar'
  1.0d0

'DU_avg2'                                          ! Input/Output
  'time averaged u-flux for 3D equations coupling'
  'meter3 second-1'                                ! [m3/s]
  'DU_avg2, scalar, series'
  'ocean_time'
  'idUfx2'
  'u2dvar'
  1.0d0

'DV_avg1'                                          ! Input/Output
  'time averaged v-flux for 2D equations'
  'meter3 second-1'                                ! [m3/s]
  'DV_avg1, scalar, series'
  'ocean_time'
  'idVfx1'
  'v2dvar'
  1.0d0

'DV_avg2'                                          ! Input/Output
  'time averaged v-flux for 3D equations coupling'
  'meter3 second-1'                                ! [m3/s]
  'DV_avg2, scalar, series'
  'ocean_time'
  'idVfx2'
  'v2dvar'
  1.0d0

'ru'                                               ! Input/Output
  'RHS of total u-momentum'
  'meter4 second-2'                                ! [m4/s2]
  'ru, scalar, series'
  'ocean_time'
  'idRu3d'
  'u3dvar'
  1.0d0

'rv'                                               ! Input/Output
  'RHS of total v-momentum'
  'meter4 second-2'                                ! [m4/s2]
  'rv, scalar, series'
  'ocean_time'
  'idRv3d'
  'v3dvar'
  1.0d0

!
!-----------------------------------------------------------------------
!  Varational data assimilation variables.
!-----------------------------------------------------------------------
!

'Nobs'                                             ! Input
  'number of observations per survey time'
  'nondimensional'                                 ! [nondimensional]
  'nulvar'
  'nulvar'
  'idNobs'
  'nulvar'
  1.0d0

'survey_time'                                      ! Input
  'observations survey time'
  'day'                                            ! [day]
  'nulvar'
  'nulvar'
  'idOday'
  'nulvar'
  1.0d0

'obs_type'                                         ! Input
  'associated model state variable type'
  'nondimensional'                                 ! [nondimensional]
  'nulvar'
  'datum'
  'idOtyp'
  'nulvar'
  1.0d0

'obs_time'                                         ! Input
  'time of observation'
  'day'                                            ! [day]
  'nulvar'
  'datum'
  'idObsT'
  'nulvar'
  1.0d0

'obs_lon'                                         ! Input
  'observation longitude'
  'degrees_east'                                  ! [degrees_east]
  'nulvar'
  'datum'
  'idOlon'
  'nulvar'
  1.0d0

'obs_lat'                                         ! Input
  'observation latitude'
  'degrees_north'                                 ! [degrees_north]
  'nulvar'
  'datum'
  'idOlat'
  'nulvar'
  1.0d0

'obs_depth'                                        ! Input/Output
  'depth of observation'
  'meter'                                          ! [m]
  'nulvar'
  'datum'
  'idObsD'
  'nulvar'
  1.0d0

'obs_Xgrid'                                        ! Input
  'x-grid observation location'
  'nondimensional'                                 ! [nondimensional]
  'nulvar'
  'datum'
  'idObsX'
  'nulvar'
  1.0d0

'obs_Ygrid'                                        ! Input
  'y-grid observation location'
  'nondimensional'                                 ! [nondimensional]
  'nulvar'
  'datum'
  'idObsY'
  'nulvar'
  1.0d0

'obs_Zgrid'                                        ! Input/Output
  'z-grid observation location'
  'nondimensional'                                 ! [nondimensional]
  'nulvar'
  'datum'
  'idObsZ'
  'nulvar'
  1.0d0

'obs_value'                                        ! Input
  'observation value'
  'state variable units'                           ! [state units]
  'nulvar'
  'datum'
  'idOval'
  'nulvar'
  1.0d0

'obs_variance'                                     ! Input
  'observation global variance'
  'squared state variable units'                   ! [squared state units]
  'nulvar'
  'datum'
  'idOvar'
  'nulvar'
  1.0d0

'obs_error'                                        ! Input
  'observation error, assigned weight'
  'inverse squared state variable units'           ! [1/(squared state units)]
  'nulvar'
  'datum'
  'idOerr'
  'nulvar'
  1.0d0

'obs_scale'                                        ! Output
  'observation screening/normalization scale'
  'nondimensional'
  'nulvar'
  'datum'
  'idObsS'
  'nulvar'
  1.0d0

'Kh'                                               ! Input
  'convolution horizontal diffusion coefficient'
  'meter2 second-1'                                ! [m2/s]
  'Kh, scalar, series'
  'ocean_time'
  'idKhor'
  'r2dvar'
  1.0d0

'Kv'                                               ! Input
  'convolution vertical diffusion coefficient'
  'meter2 second-1'                                ! [m2/s]
  'Kv, scalar, series'
  'ocean_time'
  'idKver'
  'w3dvar'
  1.0d0

'NLmodel_initial'                                  ! Output
  'initial nonlinear model at observation locations'
  'nulvar'
  'nulvar'
  'datum'
  'idNLmi'
  'nulvar'
  1.0d0

'NLmodel_value'                                    ! Output
  'nonlinear model at observation locations'
  'nulvar'
  'nulvar'
  'datum'
  'idNLmo'
  'nulvar'
  1.0d0

'TLmodel_value'                                    ! Output
  'tangent linear model at observation locations'
  'nulvar'
  'nulvar'
  'datum'
  'idTLmo'
  'nulvar'
  1.0d0

'misfit_initial'                                   ! Output
  'initial model-observation misfit'
  'nulvar'
  'nulvar'
  'datum'
  'idMOMi'
  'nulvar'
  1.0d0

'misfit_final'                                     ! Output
  'final model-observation misfit'
  'nulvar'
  'nulvar'
  'datum'
  'idMOMf'
  'nulvar'
  1.0d0

!-------------------------------------------------------------
! GOANPZ variables
! BEST_NPZ variables
!-------------------------------------------------------------

'NO3'
  'nitrate concentration'
  'millimole nitrogen meter-3'
  'nitrate, scalar, series'
  'bio_time'
  'idTvar(iNO3)'
  'r3dvar'
  1.0d0

'no3_west'
  'nitrate western boundary condition'
  'mmol N per m3'
  'no3_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iNO3)'
  'nulvar'
  1.0d0

'no3_east'
  'nitrate eastern boundary condition'
  'mmol N per m3'
  'no3_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iNO3)'
  'nulvar'
  1.0d0

'no3_south'
  'nitrate southern boundary condition'
  'mmol N per m3'
  'no3_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iNO3)'
  'nulvar'
  1.0d0

'no3_north'
  'nitrate northern boundary condition'
  'mmol N per m3'
  'no3_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iNO3)'
  'nulvar'
  1.0d0
!---------------------------------------------------------------------------
'NH4'
  'ammonia concentration'
  'millimole nitrogen meter-3'
  'ammonia, scalar, series'
  'bio_time'
  'idTvar(iNH4)'
  'r3dvar'
  1.0d0

'nh4_west'
  'ammonia western boundary condition'
  'mmol N per m3'
  'nh4_west, scalar, series '
  'bio_time'
  'idTbry(iwest,iNH4)'
  'nulvar'
  1.0d0

'nh4_east'
  'ammonia eastern boundary condition'
  'mmol N per m3'
  'nh4_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iNH4)'
  'nulvar'
  1.0d0

'nh4_south'
  'ammonia southern boundary condition'
  'mmol N per m3'
  'nh4_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iNH4)'
  'nulvar'
  1.0d0

'nh4_north'
  'ammonia northern boundary condition'
  'mmol N per m3'
  'nh4_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iNH4)'
  'nulvar'
  1.0d0

!------------------------------------------------------------------
'PhS'
  'small phytoplankton concentration'
  'milligram carbon meter-3  '
  'small phytoplankton, scalar, series   '
  'bio_time'
  'idTvar(iPhS)'
  'r3dvar'
  1.0d0

'phs_west'
  'small phytoplankton western boundary condition'
  'mg C per m3'
  'phs_west, scalar, series '
  'bio_time'
  'idTbry(iwest,iPhS)'
  'nulvar'
  1.0d0

'phs_east'
  'small phytoplankton eastern boundary condition'
  'mg C per m3'
  'phs_east, scalar, series '
  'bio_time'
  'idTbry(ieast,iPhS)'
  'nulvar'
  1.0d0

'phs_south'
  'small phytoplankton southern boundary condition'
  'mg C per m3'
  'phs_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iPhS)'
  'nulvar'
  1.0d0

'phs_north'
  'small phytoplankton northern boundary condition'
  'mg C per m3'
  'phs_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iPhS)'
  'nulvar'
  1.0d0

!-------------------------------------------------------------
'PhL'
  'large phytoplankton concentration'
  'milligram carbon meter-3'
  'large phytoplankton, scalar, series'
  'bio_time'
  'idTvar(iPhL)'
  'r3dvar'
  1.0d0

'phl_west'
  'large phytoplankton western boundary condition'
  'mg C per m3'
  'phl_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iPhL)'
  'nulvar'
  1.0d0

'phl_east'
  'large phytoplankton eastern boundary condition'
  'mg C per m3'
  'phl_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iPhL)'
  'nulvar'
  1.0d0

'phl_south'
  'large phytoplankton southern boundary condition'
  'mg C per m3'
  'phl_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iPhL)'
  'nulvar'
  1.0d0

'phl_north'
  'large phytoplankton northern boundary condition'
  'mg C per m3'
  'phl_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iPhL)'
  'nulvar'
  1.0d0
!--------------------------------------------------------
'MZS'
  'small microzooplankton concentration'
  'milligram carbon meter-3'
  'small microzooplankton, scalar, series'
  'bio_time'
  'idTvar(iMZS)'
  'r3dvar'
  1.0d0

'mzs_west'
  'large microzooplankton western boundary condition'
  'mg C per m3'
  'mzl_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iMZS)'
  'nulvar'
  1.0d0

'mzs_east'
  'large microzooplankton eastern boundary condition'
  'mg C per m3'
  'mzl_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iMZS)'
  'nulvar'
  1.0d0

'mzs_south'
  'large microzooplankton southern boundary condition'
  'mg C per m3'
  'mzl_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iMZS)'
  'nulvar'
  1.0d0

'mzs_north'
  'large microzooplankton northern boundary condition'
  'mg C per m3'
  'mzl_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iMZS)'
  'nulvar'
  1.0d0
!---------------------------------------------------------
'MZL'
  'large microzooplankton concentration'
  'milligram carbon meter-3'
  'large microzooplankton, scalar, series'
  'bio_time'
  'idTvar(iMZL)'
  'r3dvar'
  1.0d0

'mzl_west'
  'large microzooplankton western boundary condition'
  'mg C per m3'
  'mzl_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iMZL)'
  'nulvar'
  1.0d0

'mzl_east'
  'large microzooplankton eastern boundary condition'
  'mg C per m3'
  'mzl_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iMZL)'
  'nulvar'
  1.0d0

'mzl_south'
  'large microzooplankton southern boundary condition'
  'mg C per m3'
  'mzl_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iMZL)'
  'nulvar'
  1.0d0

'mzl_north'
  'large microzooplankton northern boundary condition'
  'mg C per m3'
  'mzl_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iMZL)'
  'nulvar'
  1.0d0

!----------------------------------------------------------
'Cop'
  'small coastal copepod concentration'
  'milligram carbon meter-3'
  'copepod, scalar, series'
  'bio_time'
  'idTvar(iCop)'
  'r3dvar'
  1.0d0

'cop_west'
  'copepods western boundary condition'
  'mg C per m3'
  'cop_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iCop)'
  'nulvar'
  1.0d0

'cop_east'
  'copepods eastern boundary condition'
  'mg C per m3'
  'cop_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iCop)'
  'nulvar'
  1.0d0

'cop_south'
  'copepods southern boundary condition'
  'mg C per m3'
  'cop_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iCop)'
  'nulvar'
  1.0d0

'cop_north'
  'copepods northern boundary condition'
  'mg C per m3'
  'cop_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iCop)'
  'nulvar'
  1.0d0

!----------------------------------------------------------
'NCa'
  'neocalanus spp. concentration'
  'milligram carbon meter-3'
  'neocalanus, scalar, series'
  'bio_time'
  'idTvar(iNCa)'
  'r3dvar'
  1.0d0

'nca_west'
  'neocalanus western boundary condition'
  'mg C per m3'
  'nca_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iNCa)'
  'nulvar'
  1.0d0

'nca_east'
  'neocalanus eastern boundary condition'
  'mg C per m3'
  'nca_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iNCa)'
  'nulvar'
  1.0d0

'nca_south'
  'neocalanus southern boundary condition'
  'mg C per m3'
  'nca_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iNCa)'
  'nulvar'
  1.0d0

'nca_north'
  'neocalanus northern boundary condition'
  'mg C per m3'
  'nca_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iNCa)'
  'nulvar'
  1.0d0

!----------------------------------------------------------
'NCaS'
  'neocalanus spp. concentration'
  'milligram carbon meter-3'
  'neocalanus, scalar, series'
  'bio_time'
  'idTvar(iNCaS)'
  'r3dvar'
  1.0d0

'nca_westS'
  'neocalanus western boundary condition'
  'mg C per m3'
  'nca_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iNCaS)'
  'nulvar'
  1.0d0

'nca_eastS'
  'neocalanus eastern boundary condition'
  'mg C per m3'
  'nca_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iNCaS)'
  'nulvar'
  1.0d0

'nca_southS'
  'neocalanus southern boundary condition'
  'mg C per m3'
  'nca_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iNCaS)'
  'nulvar'
  1.0d0

'nca_northS'
  'neocalanus northern boundary condition'
  'mg C per m3'
  'nca_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iNCaS)'
  'nulvar'
  1.0d0

!----------------------------------------------------------
'NCaO'
  'Offshore neocalanus spp. concentration'
  'milligram carbon meter-3'
  'neocalanus, scalar, series'
  'bio_time'
  'idTvar(iNCaO)'
  'r3dvar'
  1.0d0

'nca_westO'
  'Offshore neocalanus western boundary condition'
  'mg C per m3'
  'nca_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iNCaO)'
  'nulvar'
  1.0d0

'nca_eastO'
  'Offshore neocalanus eastern boundary condition'
  'mg C per m3'
  'nca_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iNCaO)'
  'nulvar'
  1.0d0

'nca_southO'
  'Offshore neocalanus southern boundary condition'
  'mg C per m3'
  'nca_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iNCaO)'
  'nulvar'
  1.0d0

'nca_northO'
  'Offshore neocalanus northern boundary condition'
  'mg C per m3'
  'nca_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iNCaO)'
  'nulvar'
  1.0d0

!-----------------------------------------------------------
'Eup'
  'euphausiid concentration'
  'milligram carbon meter-3'
  'euphausiid, scalar, series'
  'bio_time'
  'idTvar(iEup)'
  'r3dvar'
  1.0d0

'eup_west'
  'euphausiids western boundary condition'
  'mg C per m3'
  'eup_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iEup)'
  'nulvar'
  1.0d0

'eup_east'
  'euphausiids eastern boundary condition'
  'mg C per m3'
  'eup_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iEup)'
  'nulvar'
  1.0d0

'eup_south'
  'euphausiids southern boundary condition'
  'mg C per m3'
  'eup_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iEup)'
  'nulvar'
  1.0d0

'eup_north'
  'euphausiids northern boundary condition'
  'mg C per m3'
  'eup_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iEup)'
  'nulvar'
  1.0d0

!-----------------------------------------------------------
'EupS'
  'euphausiid concentration'
  'milligram carbon meter-3'
  'euphausiid, scalar, series'
  'bio_time'
  'idTvar(iEupS)'
  'r3dvar'
  1.0d0

'eup_westS'
  'euphausiids western boundary condition'
  'mg C per m3'
  'eup_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iEupS)'
  'nulvar'
  1.0d0

'eup_eastS'
  'euphausiids eastern boundary condition'
  'mg C per m3'
  'eup_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iEupS)'
  'nulvar'
  1.0d0

'eup_southS'
  'euphausiids southern boundary condition'
  'mg C per m3'
  'eup_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iEupS)'
  'nulvar'
  1.0d0

'eup_northS'
  'euphausiids northern boundary condition'
  'mg C per m3'
  'eup_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iEupS)'
  'nulvar'
  1.0d0

!-----------------------------------------------------------
'EupO'
  'Offshore euphausiid concentration'
  'milligram carbon meter-3'
  'euphausiid, scalar, series'
  'bio_time'
  'idTvar(iEupO)'
  'r3dvar'
  1.0d0

'eup_westO'
  'Offshore euphausiids western boundary condition'
  'mg C per m3'
  'eup_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iEupO)'
  'nulvar'
  1.0d0

'eup_eastO'
  'Offshore euphausiids eastern boundary condition'
  'mg C per m3'
  'eup_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iEupO)'
  'nulvar'
  1.0d0

'eup_southO'
  'Offshore euphausiids southern boundary condition'
  'mg C per m3'
  'eup_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iEupO)'
  'nulvar'
  1.0d0

'eup_northO'
  'Offshore euphausiids northern boundary condition'
  'mg C per m3'
  'eup_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iEupO)'
  'nulvar'
  1.0d0

!-------------------------------------------------------
'Det'
  'detritus concentration'
  'milligram carbon meter-3'
  'detritus, scalar, series'
  'bio_time'
  'idTvar(iDet)'
  'r3dvar'
  1.0d0

'det_west'
  'detritus western boundary condition'
  'mg C per m3'
  'det_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iDet)'
  'nulvar'
  1.0d0

'det_east'
  'detritus eastern boundary condition'
  'mg C per m3'
  'det_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iDet)'
  'nulvar'
  1.0d0

'det_south'
  'detritus southern boundary condition'
  'mg C per m3'
  'det_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iDet)'
  'nulvar'
  1.0d0

'det_north'
  'detritus northern boundary condition'
  'mg C per m3'
  'det_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iDet)'
  'nulvar'
  1.0d0

!-------------------------------------------------------
'DetF'
  'Fast sinking detritus concentration'
  'milligram carbon meter-3'
  'detritus, scalar, series'
  'bio_time'
  'idTvar(iDetF)'
  'r3dvar'
  1.0d0

'detF_west'
  'Fast sinking detritus western boundary condition'
  'mg C per m3'
  'detF_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iDetF)'
  'nulvar'
  1.0d0

'detF_east'
  'Fast sinking detritus eastern boundary condition'
  'mg C per m3'
  'detF_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iDetF)'
  'nulvar'
  1.0d0

'detF_south'
  'Fast sinking detritus southern boundary condition'
  'mg C per m3'
  'detF_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iDetF)'
  'nulvar'
  1.0d0

'detF_north'
  'Fast sinking detritus northern boundary condition'
  'mg C per m3'
  'detF_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iDetF)'
  'nulvar'
  1.0d0

!--------------
'Iron'
  'iron concentration'
  'micromol Fe m-3'
  'iron, scalar, series'
  'bio_time'
  'idTvar(iFe)'
  'r3dvar'
  1.0d0

'iron_west'
  'iron western boundary condition'
  'micromol Fe m-3'
  'iron_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iFe)'
  'nulvar'
  1.0d0

'iron_east'
  'iron eastern boundary condition'
  'micromol Fe m-3'
  'iron_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iFe)'
  'nulvar'
  1.0d0

'iron_south'
  'iron southern boundary condition'
  'micromol Fe m-3'
  'iron_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iFe)'
  'nulvar'
  1.0d0

'iron_north'
  'iron northern boundary condition'
  'micromol Fe m-3'
  'iron_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iFe)'
  'nulvar'
  1.0d0

'Iron'
  'iron climatology'
  'micromol Fe m-3'
  'iron, scalar, series'
  'ocean_time'
  'idTclm(iFe)'
  'r3dvar'
  1.0d0

!ICE_BIO
'ice_thick'
  'ice thickness climatology'
  'meters'
  'ice_thick, scalar, series'
  'ocean_time'
  'idTclm(i1CI)'
  'r2dvar'
  1.0d0

!3DAKT
'AKt'
  'AKt from 3D run'
  'meter2 second-1 '
  'Akt, scalar, series'
  'ocean_time'
  'idTclm(iAKt3)'
  'r2dvar'
  1.0d0


!------------------------------------------------------------
 'Jel'
  'Jellyfish concentration'
  'milligram carbon meter-3'
  'jellyfish, scalar, series'
  'bio_time'
  'idTvar(iJel)'
  'r3dvar'
  1.0d0


'jel_west'
  'jellyfish western boundary condition'
  'mg C per m3'
  'jel_west, scalar, series'
  'bio_time'
  'idTbry(iwest,iJel)'
  'nulvar'
  1.0d0

'jel_east'
  'jellyfish eastern boundary condition'
  'mg C per m3'
  'jel_east, scalar, series'
  'bio_time'
  'idTbry(ieast,iJel)'
  'nulvar'
  1.0d0

'jel_south'
  'jellyfish southern boundary condition'
  'mg C per m3'
  'jel_south, scalar, series'
  'bio_time'
  'idTbry(isouth,iJel)'
  'nulvar'
  1.0d0

'jel_north'
  'jellyfish northern boundary condition'
  'mg C per m3'
  'jel_north, scalar, series'
  'bio_time'
  'idTbry(inorth,iJel)'
  'nulvar'
  1.0d0

!------------------------------------------------------------
!BENTHIC
!--------------------
'Ben'
  'benthos concentration'
  'milligram carbon meter-2'
  'benthos, scalar, series'
  'bio_time'
  'idBeTvar(iBen)'
  'r3dvar'
  1.0d0

'DetBen'
  'benthic detritus concentration'
  'milligram carbon meter-2'
  'benthic detritus, scalar, series'
  'bio_time'
  'idBeTvar(iBenDet)'
  'r3dvar'
  1.0d0

!-------
!ICE_BIO
!-------
'IcePhL'
  'Ice algae concentration'
  'milligram carbon meter-3'
  'Ice algae, scalar, series'
  'bio_time'
  'idIcePhL'
  'r2dvar'
  1.0d0

'Ice Algae_west'
  'Ice algae concentration western boundary condition'
  'meter'
  'Ice_algae_thick_west, scalar, series'
  'bry_time'
  'idIcePhLbc(iwest)'
  'nulvar'
  1.0d0

'Ice Algae_east'
  'Ice algae concentration eastern boundary condition'
  'meter'
  'Ice_algae_thick_east, scalar, series'
  'bry_time'
  'idIcePhLbc(ieast)'
  'nulvar'
  1.0d0

'Ice Algae_north'
  'Ice algae concentration northern boundary condition'
  'meter'
  'Ice_algae_thick_north, scalar, series'
  'bry_time'
  'idIcePhLbc(inorth)'
  'nulvar'
  1.0d0

'Ice Algae_south'
  'Ice algae concentration southern boundary condition'
  'meter'
  'Ice_algae_thick_south, scalar, series'
  'bry_time'
  'idIcePhLbc(isouth)'
  'nulvar'
  1.0d0

'IceNO3'
  'Ice nitrate concentration'
  'mgC/m3'
  'Ice nitrate, scalar, series'
  'ocean_time'
  'idIceNO3'
  'r2dvar'
  1.0d0

'Ice_Nitrate_west'
  'Ice nitrate concentration western boundary condition'
  'meter'
  'ice_nitrate_conc_west, scalar, series'
  'bry_time'
  'idIceNO3bc(iwest)'
  'nulvar'
  1.0d0

'Ice_Nitrate_east'
  'Ice nitrate concentration eastern boundary condition'
  'meter'
  'ice_nitrate_conc_east, scalar, series'
  'bry_time'
  'idIceNO3bc(ieast)'
  'nulvar'
  1.0d0

'Ice_Nitrate_north'
  'Ice nitrate concentration northern boundary condition'
  'meter'
  'ice_nitrate_conc_north, scalar, series'
  'bry_time'
  'idIceNO3bc(inorth)'
  'nulvar'
  1.0d0

'Ice_Nitrate_south'
  'Ice nitrate concentration southern boundary condition'
  'meter'
  'ice_nitrate_conc_south, scalar, series'
  'bry_time'
  'idIceNO3bc(isouth)'
  'nulvar'
  1.0d0

'IceNH4'
  'Ice ammonium concentration'
  'mgC/m3'
  'Ice ammonium, scalar, series'
  'ocean_time'
  'idIceNH4'
  'r2dvar'
  1.0d0

'Ice_Ammonium_west'
  'Ice Ammonium concentration western boundary condition'
  'meter'
  'ice_Ammonium_conc_west, scalar, series'
  'bry_time'
  'idIceNH4bc(iwest)'
  'nulvar'
  1.0d0

'Ice_Ammonium_east'
  'Ice Ammonium concentration eastern boundary condition'
  'meter'
  'ice_Ammonium_conc_east, scalar, series'
  'bry_time'
  'idIceNH4bc(ieast)'
  'nulvar'
  1.0d0

'Ice_Ammonium_north'
  'Ice Ammonium concentration northern boundary condition'
  'meter'
  'ice_Ammonium_conc_north, scalar, series'
  'bry_time'
  'idIceNH4bc(inorth)'
  'nulvar'
  1.0d0

'Ice_Ammonium_south'
  'Ice Ammonium concentration southern boundary condition'
  'meter'
  'ice_Ammonium_conc_south, scalar, series'
  'bry_time'
  'idIceNH4bc(isouth)'
  'nulvar'
  1.0d0

'IceZ'
  'Ice Thickness'
  'meters'
  'Ice thickness, series'
  'bio_time'
  'idIceBvar(iIceZ)'
  'r2dvar'
  1.0d0

'IceLog'
  'Logical Ice Counter'
  'unitless'
  'Logical Ice Counter, scalar, series'
  'bio_time'
  'idIceLog'
  'r2dvar'
  1.0d0

!------------------------------------------------------------
!STATIONARY
!-----------------------
'Stat1_3d'
  '1st 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat1)'
  'r3dvar'
  1.0d0

'Stat2_3d'
  '2nd 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat2)'
  'r3dvar'
  1.0d0

'Stat3_3d'
  '3rd 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat3)'
  'r3dvar'
  1.0d0

'Stat4_3d'
  '4th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat4)'
  'r3dvar'
  1.0d0

'Stat5_3d'
  '5th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat5)'
  'r3dvar'
  1.0d0

'Stat6_3d'
  '6th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat6)'
  'r3dvar'
  1.0d0

'Stat7_3d'
  '7th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat7)'
  'r3dvar'
  1.0d0

'Stat8_3d'
  '8th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat8)'
  'r3dvar'
  1.0d0


'Stat9_3d'
  '9th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat9)'
  'r3dvar'
  1.0d0


'Stat10_3d'
  '10th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat10)'
  'r3dvar'
  1.0d0


'Stat11_3d'
  '11th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat11)'
  'r3dvar'
  1.0d0

'Stat12_3d'
  '12th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat12)'
  'r3dvar'
  1.0d0

'Stat13_3d'
  '13th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat13)'
  'r3dvar'
  1.0d0

'Stat14_3d'
  '14th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat14)'
  'r3dvar'
  1.0d0

'Stat15_3d'
  '15th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat15)'
  'r3dvar'
  1.0d0

'Stat16_3d'
  '16th 3D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTSvar(i3Stat16)'
  'r3dvar'
  1.0d0
!-----------------------------
!PROD3
!---------------------------
'prod_PhS'
  'primary production small phytoplankton'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT3var(iPhSprd)'
  'r3dvar'
  1.0d0

'prod_PhL'
  'primary production large phytoplankton'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT3var(iPhLprd)'
  'r3dvar'
  1.0d0

'prod_MZS'
  'secondary production small microzooplankton'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT3var(iMZSprd)'
  'r3dvar'
  1.0d0

'prod_MZL'
  'secondary production large microzooplankton'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT3var(iMZLprd)'
  'r3dvar'
  1.0d0

'prod_Cop'
  'secondary production copepods'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT3var(iCopPrd)'
  'r3dvar'
  1.0d0

'prod_NCa'
  'secondary production Neocalanus'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT3var(iNCaPrd)'
  'r3dvar'
  1.0d0

'prod_Eup'
  'secondary production Euphausiids'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT3var(iEupPrd)'
  'r3dvar'
  1.0d0

'prod_Jel'
  'secondary production Jellyfish'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT3var(iJelPrd)'
  'r3dvar'
  1.0d0

!----------------------------
! PROD2
!----------------------------
'prod_Ben'
  'secondary production Benthos'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT2var(iBenPrd)'
  'r2dvar'
  1.0d0

'prod_IcePhL'
  'primary production ice algae'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT2var(iIAPrd)'
  'r2dvar'
  1.0d0

'prod_X'
  'temporary production array'
  'milligram carbon meter-3 d-1'
  'PrimProd scalar series'
  'bio_time'
  'idPT2var(iXPrd)'
  'r2dvar'
  1.0d0

!----------------------------
! BIOFLUX
!----------------------------
'bioflux'
  '2d bioflux variable'
  'flux mmol'
  'BioFlx scalar series'
  'bio_time'
  'idTBFvar(iBF)'
  'r2dvar'
  1.0d0

!------------------------------------------------------------
!STATIONARY2
!-----------------------
'Stat1_2d'
  '1st 2D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTS2var(i2Stat1)'
  'r2dvar'
  1.0d0

'Stat2_2d'
  '2nd 2D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTS2var(i2Stat2)'
  'r2dvar'
  1.0d0

'Stat3_2d'
  '3rd 2D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTS2var(i2Stat3)'
  'r2dvar'
  1.0d0

'Stat4_2d'
  '4th 2D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTS2var(i2Stat4)'
  'r2dvar'
  1.0d0

'Stat5_2d'
  '5th 2D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTS2var(i2Stat5)'
  'r2dvar'
  1.0d0

'Stat6_2d'
  '6th 2D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTS2var(i2Stat6)'
  'r2dvar'
  1.0d0

'Stat7_2d'
  '7th 2D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTS2var(i2Stat7)'
  'r2dvar'
  1.0d0

'Stat8_2d'
  '8th 2D stationary variable'
  'milligram carbon meter-3 d-1'
  'StatVar scalar series'
  'bio_time'
  'idTS2var(i2Stat8)'
  'r2dvar'
  1.0d0

'BgError_value'                                    ! Output
  'Background error at observation locations'
  'nulvar'
  'nulvar'
  'datum'
  'idBgEr'
  'nulvar'
  1.0d0

'BgThresh_value'                                   ! Output
  'Threshold for background quality control check of observations'
  'nulvar'
  'nulvar'
  'datum'
  'idBgTh'
  'nulvar'
  1.0d0