varinfo.dat.cobalt

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

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

'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

'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

! atmospheric forcing for Cobalt BGC model

'atmCO2'                                           ! Input
  'atmCO2'
  'ppmv'                                           ! Input: [ppmv]
  'atmCO2, scalar, series'                         ! [ppmv]
  'atmCO2_time'                                    ! Output: []
  'idAtmCO2'
  'r2dvar'
  1.0d0

'ironsed'                                          ! Input
  'ironsed'
  'mol.m-2.s-1'                                    ! Input: [mol.m-2.s-1]
  'ironsed, scalar, series'                        ! [mol.m-2.s-1]
  'ironsed_time'                                   ! Output: []
  'idIronSed'
  'r2dvar'
  1.0d0

'fecoast'                                          ! Input
  'fecoast'
  'mol Fe m kg-1 s-1'                              ! Input: [mol Fe m kg-1 s-1]
  'fecoast, scalar, series'                        ! [mol Fe m kg-1 s-1]
  'fecoast_time'                                   ! Output: []
  'idFeCoast'
  'r2dvar'
  1.0d0

'solublefe'                                        ! Input
  'solublefe'
  'mol.m-2.s-1'                                    ! Input: [ppmv]
  'soluble_fe, scalar, series'                     ! [ppmv]
  'solublefe_time'                                 ! Output: []
  'idSolubleFe'
  'r2dvar'
  1.0d0

'mineralfe'                                        ! Input
  'mineralfe'
  'g.m-2.s-1'                                      ! Input: [ppmv]
  'mineral_fe, scalar, series'                     ! [ppmv]
  'mineralfe_time'                                 ! Output: []
  'idMineralFe'
  'r2dvar'
  1.0d0

! end atmospheric forcing for Cobalt BGC model

'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

!------------------------------------------------------------------------------
! Cobalt river input
!------------------------------------------------------------------------------

'river_ldon'
  'ldon river source'
  'mol.m-2.s-1'
  'river_ldon, scalar, series'
  'runoff_time'
  'idriver_ldon'
  'r2dvar'
  1.0d0

'river_ldop'
  'ldop river source'
  'mol.m-2.s-1'
  'river_ldop, scalar, series'
  'runoff_time'
  'idriver_ldop'
  'r2dvar'
  1.0d0

'river_ndet'
  'ndet river source'
  'mol.m-2.s-1'
  'river_ndet, scalar, series'
  'runoff_time'
  'idriver_ndet'
  'r2dvar'
  1.0d0

'river_no3'
  'no3 river source'
  'mol.m-2.s-1'
  'river_no3, scalar, series'
  'runoff_time'
  'idriver_no3'
  'r2dvar'
  1.0d0

'river_po4'
  'po4 river source'
  'mol.m-2.s-1'
  'river_po4, scalar, series'
  'runoff_time'
  'idriver_po4'
  'r2dvar'
  1.0d0

'river_srdon'
  'srdon river source'
  'mol.m-2.s-1'
  'river_srdon, scalar, series'
  'runoff_time'
  'idriver_srdon'
  'r2dvar'
  1.0d0

'river_srdop'
  'srdop river source'
  'mol.m-2.s-1'
  'river_srdop, scalar, series'
  'runoff_time'
  'idriver_srdop'
  'r2dvar'
  1.0d0

'river_sldon'
  'sldon river source'
  'mol.m-2.s-1'
  'river_sldon, scalar, series'
  'runoff_time'
  'idriver_sldon'
  'r2dvar'
  1.0d0

'river_sldop'
  'sldop river source'
  'mol.m-2.s-1'
  'river_sldop, scalar, series'
  'runoff_time'
  'idriver_sldop'
  'r2dvar'
  1.0d0

'river_fed'
  'fed river source'
  'mol.m-2.s-1'
  'river_fed, scalar, series'
  'runoff_time'
  'idriver_fed'
  'r2dvar'
  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 Cobalt (idTvar, idRtrc)
!---------------------------------------------------------------------
!


'alk'
  'Alkalinity'
  'mol/kg'
  'alk, scalar, series'
  'ocean_time'
  'idTvar(ialk)'
  'r3dvar'
  1.0d0

'river_alk'
  'river runoff Alkalinity'
  'mol/kg'
  'alk, scalar, series'
  'ocean_time'
  'idRtrc(ialk)'
  'r3dvar'
  1.0d0


'cadet_arag'
  'Detrital CaCO3'
  'mol/kg'
  'cadet_arag, scalar, series'
  'ocean_time'
  'idTvar(icadet_arag)'
  'r3dvar'
  1.0d0

'river_cadet_arag'
  'river runoff Detrital CaCO3'
  'mol/kg'
  'cadet_arag, scalar, series'
  'ocean_time'
  'idRtrc(icadet_arag)'
  'r3dvar'
  1.0d0


'cadet_calc'
  'Detrital CaCO3'
  'mol/kg'
  'cadet_calc, scalar, series'
  'ocean_time'
  'idTvar(icadet_calc)'
  'r3dvar'
  1.0d0

'river_cadet_calc'
  'river runoff Detrital CaCO3'
  'mol/kg'
  'cadet_calc, scalar, series'
  'ocean_time'
  'idRtrc(icadet_calc)'
  'r3dvar'
  1.0d0


'dic'
  'Dissolved Inorganic Carbon'
  'mol/kg'
  'dic, scalar, series'
  'ocean_time'
  'idTvar(idic)'
  'r3dvar'
  1.0d0

'river_dic'
  'river runoff Dissolved Inorganic Carbon'
  'mol/kg'
  'dic, scalar, series'
  'ocean_time'
  'idRtrc(idic)'
  'r3dvar'
  1.0d0


'fed'
  'Dissolved Iron'
  'mol/kg'
  'fed, scalar, series'
  'ocean_time'
  'idTvar(ifed)'
  'r3dvar'
  1.0d0

'river_fed'
  'river runoff Dissolved Iron'
  'mol/kg'
  'fed, scalar, series'
  'ocean_time'
  'idRtrc(ifed)'
  'r3dvar'
  1.0d0


'fedet'
  'Detrital Iron'
  'mol/kg'
  'fedet, scalar, series'
  'ocean_time'
  'idTvar(ifedet)'
  'r3dvar'
  1.0d0

'river_fedet'
  'river runoff Detrital Iron'
  'mol/kg'
  'fedet, scalar, series'
  'ocean_time'
  'idRtrc(ifedet)'
  'r3dvar'
  1.0d0


'fedi'
  'Diazotroph Iron'
  'mol/kg'
  'fedi, scalar, series'
  'ocean_time'
  'idTvar(ifedi)'
  'r3dvar'
  1.0d0

'river_fedi'
  'river runoff Diazotroph Iron'
  'mol/kg'
  'fedi, scalar, series'
  'ocean_time'
  'idRtrc(ifedi)'
  'r3dvar'
  1.0d0


'felg'
  'Large Phytoplankton Iron'
  'mol/kg'
  'felg, scalar, series'
  'ocean_time'
  'idTvar(ifelg)'
  'r3dvar'
  1.0d0

'river_felg'
  'river runoff Large Phytoplankton Iron'
  'mol/kg'
  'felg, scalar, series'
  'ocean_time'
  'idRtrc(ifelg)'
  'r3dvar'
  1.0d0


'fesm'
  'Small Phytoplankton Iron'
  'mol/kg'
  'fesm, scalar, series'
  'ocean_time'
  'idTvar(ifesm)'
  'r3dvar'
  1.0d0

'river_fesm'
  'river runoff Small Phytoplankton Iron'
  'mol/kg'
  'fesm, scalar, series'
  'ocean_time'
  'idRtrc(ifesm)'
  'r3dvar'
  1.0d0


'ldon'
  'labile DON'
  'mol/kg'
  'ldon, scalar, series'
  'ocean_time'
  'idTvar(ildon)'
  'r3dvar'
  1.0d0

'river_ldon'
  'river runoff labile DON'
  'mol/kg'
  'ldon, scalar, series'
  'ocean_time'
  'idRtrc(ildon)'
  'r3dvar'
  1.0d0


'ldop'
  'labile DOP'
  'mol/kg'
  'ldop, scalar, series'
  'ocean_time'
  'idTvar(ildop)'
  'r3dvar'
  1.0d0

'river_ldop'
  'river runoff labile DOP'
  'mol/kg'
  'ldop, scalar, series'
  'ocean_time'
  'idRtrc(ildop)'
  'r3dvar'
  1.0d0


'lith'
  'Lithogenic Aluminosilicate'
  'g/kg'
  'lith, scalar, series'
  'ocean_time'
  'idTvar(ilith)'
  'r3dvar'
  1.0d0

'river_lith'
  'river runoff Lithogenic Aluminosilicate'
  'g/kg'
  'lith, scalar, series'
  'ocean_time'
  'idRtrc(ilith)'
  'r3dvar'
  1.0d0


'lithdet'
  'lithdet'
  'g/kg'
  'lithdet, scalar, series'
  'ocean_time'
  'idTvar(ilithdet)'
  'r3dvar'
  1.0d0

'river_lithdet'
  'river runoff lithdet'
  'g/kg'
  'lithdet, scalar, series'
  'ocean_time'
  'idRtrc(ilithdet)'
  'r3dvar'
  1.0d0


'nbact'
  'bacterial'
  'mol/kg'
  'nbact, scalar, series'
  'ocean_time'
  'idTvar(inbact)'
  'r3dvar'
  1.0d0

'river_nbact'
  'river runoff bacterial'
  'mol/kg'
  'nbact, scalar, series'
  'ocean_time'
  'idRtrc(inbact)'
  'r3dvar'
  1.0d0


'ndet'
  'ndet'
  'mol/kg'
  'ndet, scalar, series'
  'ocean_time'
  'idTvar(indet)'
  'r3dvar'
  1.0d0

'river_ndet'
  'river runoff ndet'
  'mol/kg'
  'ndet, scalar, series'
  'ocean_time'
  'idRtrc(indet)'
  'r3dvar'
  1.0d0


'ndi'
  'Diazotroph Nitrogen'
  'mol/kg'
  'ndi, scalar, series'
  'ocean_time'
  'idTvar(indi)'
  'r3dvar'
  1.0d0

'river_ndi'
  'river runoff Diazotroph Nitrogen'
  'mol/kg'
  'ndi, scalar, series'
  'ocean_time'
  'idRtrc(indi)'
  'r3dvar'
  1.0d0


'nlg'
  'Large Phytoplankton Nitrogen'
  'mol/kg'
  'nlg, scalar, series'
  'ocean_time'
  'idTvar(inlg)'
  'r3dvar'
  1.0d0

'river_nlg'
  'river runoff Large Phytoplankton Nitrogen'
  'mol/kg'
  'nlg, scalar, series'
  'ocean_time'
  'idRtrc(inlg)'
  'r3dvar'
  1.0d0


'nsm'
  'Small Phytoplankton Nitrogen'
  'mol/kg'
  'nsm, scalar, series'
  'ocean_time'
  'idTvar(insm)'
  'r3dvar'
  1.0d0

'river_nsm'
  'river runoff Small Phytoplankton Nitrogen'
  'mol/kg'
  'nsm, scalar, series'
  'ocean_time'
  'idRtrc(insm)'
  'r3dvar'
  1.0d0


'nh4'
  'Ammonia'
  'mol/kg'
  'nh4, scalar, series'
  'ocean_time'
  'idTvar(inh4)'
  'r3dvar'
  1.0d0

'river_nh4'
  'river runoff Ammonia'
  'mol/kg'
  'nh4, scalar, series'
  'ocean_time'
  'idRtrc(inh4)'
  'r3dvar'
  1.0d0


'no3'
  'Nitrate'
  'mol/kg'
  'no3, scalar, series'
  'ocean_time'
  'idTvar(ino3)'
  'r3dvar'
  1.0d0

'river_no3'
  'river runoff Nitrate'
  'mol/kg'
  'no3, scalar, series'
  'ocean_time'
  'idRtrc(ino3)'
  'r3dvar'
  1.0d0


'o2'
  'Oxygen'
  'mol/kg'
  'o2, scalar, series'
  'ocean_time'
  'idTvar(io2)'
  'r3dvar'
  1.0d0

'river_o2'
  'river runoff Oxygen'
  'mol/kg'
  'o2, scalar, series'
  'ocean_time'
  'idRtrc(io2)'
  'r3dvar'
  1.0d0


'pdet'
  'Detrital Phosphorus'
  'mol/kg'
  'pdet, scalar, series'
  'ocean_time'
  'idTvar(ipdet)'
  'r3dvar'
  1.0d0

'river_pdet'
  'river runoff Detrital Phosphorus'
  'mol/kg'
  'pdet, scalar, series'
  'ocean_time'
  'idRtrc(ipdet)'
  'r3dvar'
  1.0d0


'po4'
  'Phosphate'
  'mol/kg'
  'po4, scalar, series'
  'ocean_time'
  'idTvar(ipo4)'
  'r3dvar'
  1.0d0

'river_po4'
  'river runoff Phosphate'
  'mol/kg'
  'po4, scalar, series'
  'ocean_time'
  'idRtrc(ipo4)'
  'r3dvar'
  1.0d0


'srdon'
  'Semi-Refractory DON'
  'mol/kg'
  'srdon, scalar, series'
  'ocean_time'
  'idTvar(isrdon)'
  'r3dvar'
  1.0d0

'river_srdon'
  'river runoff Semi-Refractory DON'
  'mol/kg'
  'srdon, scalar, series'
  'ocean_time'
  'idRtrc(isrdon)'
  'r3dvar'
  1.0d0


'srdop'
  'Semi-Refractory DOP'
  'mol/kg'
  'srdop, scalar, series'
  'ocean_time'
  'idTvar(isrdop)'
  'r3dvar'
  1.0d0

'river_srdop'
  'river runoff Semi-Refractory DOP'
  'mol/kg'
  'srdop, scalar, series'
  'ocean_time'
  'idRtrc(isrdop)'
  'r3dvar'
  1.0d0


'sldon'
  'Semilabile DON'
  'mol/kg'
  'sldon, scalar, series'
  'ocean_time'
  'idTvar(isldon)'
  'r3dvar'
  1.0d0

'river_sldon'
  'river runoff Semilabile DON'
  'mol/kg'
  'sldon, scalar, series'
  'ocean_time'
  'idRtrc(isldon)'
  'r3dvar'
  1.0d0


'sldop'
  'Semilabile DOP'
  'mol/kg'
  'sldop, scalar, series'
  'ocean_time'
  'idTvar(isldop)'
  'r3dvar'
  1.0d0

'river_sldop'
  'river runoff Semilabile DOP'
  'mol/kg'
  'sldop, scalar, series'
  'ocean_time'
  'idRtrc(isldop)'
  'r3dvar'
  1.0d0


'sidet'
  'Detrital Silicon'
  'mol/kg'
  'sidet, scalar, series'
  'ocean_time'
  'idTvar(isidet)'
  'r3dvar'
  1.0d0

'river_sidet'
  'river runoff Detrital Silicon'
  'mol/kg'
  'sidet, scalar, series'
  'ocean_time'
  'idRtrc(isidet)'
  'r3dvar'
  1.0d0


'silg'
  'Large Phytoplankton Silicon'
  'mol/kg'
  'silg, scalar, series'
  'ocean_time'
  'idTvar(isilg)'
  'r3dvar'
  1.0d0

'river_silg'
  'river runoff Large Phytoplankton Silicon'
  'mol/kg'
  'silg, scalar, series'
  'ocean_time'
  'idRtrc(isilg)'
  'r3dvar'
  1.0d0


'sio4'
  'Silicate'
  'mol/kg'
  'sio4, scalar, series'
  'ocean_time'
  'idTvar(isio4)'
  'r3dvar'
  1.0d0

'river_sio4'
  'river runoff Silicate'
  'mol/kg'
  'sio4, scalar, series'
  'ocean_time'
  'idRtrc(isio4)'
  'r3dvar'
  1.0d0


'nsmz'
  'Small Zooplankton Nitrogen'
  'mol/kg'
  'nsmz, scalar, series'
  'ocean_time'
  'idTvar(insmz)'
  'r3dvar'
  1.0d0

'river_nsmz'
  'river runoff Small Zooplankton Nitrogen'
  'mol/kg'
  'nsmz, scalar, series'
  'ocean_time'
  'idRtrc(insmz)'
  'r3dvar'
  1.0d0


'nmdz'
  'Medium-sized zooplankton Nitrogen'
  'mol/kg'
  'nmdz, scalar, series'
  'ocean_time'
  'idTvar(inmdz)'
  'r3dvar'
  1.0d0

'river_nmdz'
  'river runoff Medium-sized zooplankton Nitrogen'
  'mol/kg'
  'nmdz, scalar, series'
  'ocean_time'
  'idRtrc(inmdz)'
  'r3dvar'
  1.0d0


'nlgz'
  'large Zooplankton Nitrogen'
  'mol/kg'
  'nlgz, scalar, series'
  'ocean_time'
  'idTvar(inlgz)'
  'r3dvar'
  1.0d0

'river_nlgz'
  'river runoff large Zooplankton Nitrogen'
  'mol/kg'
  'nlgz, scalar, series'
  'ocean_time'
  'idRtrc(inlgz)'
  'r3dvar'
  1.0d0

! coastal diatoms (NVO)

'nmd'
  'medium Phytoplankton Nitrogen'
  'mol/kg'
  'nmd, scalar, series'
  'ocean_time'
  'idTvar(inmd)'
  'r3dvar'
  1.0d0

'river_nmd'
  'river medium Phytoplankton Nitrogen'
  'mol/kg'
  'river_nmd, scalar, series'
  'ocean_time'
  'idRtrc(inmd)'
  'r3dvar'
  1.0d0

'simd'
  'medium Phytoplankton Silicon'
  'mol/kg'
  'simd, scalar, series'
  'ocean_time'
  'idTvar(isimd)'
  'r3dvar'
  1.0d0

'river_simd'
  'river medium Phytoplankton Silicon'
  'mol/kg'
  'river_simd, scalar, series'
  'ocean_time'
  'idRtrc(isimd)'
  'r3dvar'
  1.0d0

'femd'
  'medium Phytoplankton Iron'
  'mol/kg'
  'femd, scalar, series'
  'ocean_time'
  'idTvar(ifemd)'
  'r3dvar'
  1.0d0

'river_femd'
  'river medium Phytoplankton Iron'
  'mol/kg'
  'river_femd, scalar, series'
  'ocean_time'
  'idRtrc(ifemd)'
  'r3dvar'
  1.0d0

! Other BGC for Cobalt

'chl'
  'Chlorophyll'
  'ug kg-1'
  'chl, scalar, series'
  'ocean_time'
  'iDobgc(iochl)'
  'r3dvar'
  1.0d0

'co3_ion'
  'Carbonate ion'
  'mol/kg'
  'co3_ion, scalar, series'
  'ocean_time'
  'iDobgc(ioco3_ion)'
  'r3dvar'
  1.0d0

'htotal'
  'H+ ion concentration'
  'mol/kg'
  'htotal, scalar, series'
  'ocean_time'
  'iDobgc(iohtotal)'
  'r3dvar'
  1.0d0

'irr_mem'
  'Irradiance memory'
  'Watts/m^2'
  'irr_mem, scalar, series'
  'ocean_time'
  'iDobgc(ioirr_mem)'
  'r3dvar'
  1.0d0

'mu_mem_sm'
  'Aggreg memory small phyto'
  ''
  'mu_mem_sm, scalar, series'
  'ocean_time'
  'iDobgc(iomu_mem_sm)'
  'r3dvar'
  1.0d0

'mu_mem_di'
  'Aggreg memory diatoms'
  ''
  'mu_mem_di, scalar, series'
  'ocean_time'
  'iDobgc(iomu_mem_di)'
  'r3dvar'
  1.0d0

'mu_mem_lg'
  'Aggreg memory large phyto'
  ''
  'mu_mem_lg, scalar, series'
  'ocean_time'
  'iDobgc(iomu_mem_lg)'
  'r3dvar'
  1.0d0

'mu_mem_md'
  'Aggreg memory medium phyto'
  ''
  'mu_mem_md, scalar, series'
  'ocean_time'
  'iDobgc(iomu_mem_md)'
  'r3dvar'
  1.0d0

! Diagnostic variables for Cobalt

'cased'
  'Sediment CaCO3'
  'mol.m-3'
  'cased, scalar, series'
  'ocean_time'
  'iDbio2(icased)'
  'r2dvar'
  1.0d0

'cadet_arag_btf'
  'aragonite flux to Sediments'
  'mol m-2 s-1'
  'cadet_arag_btf, scalar, series'
  'ocean_time'
  'iDbio2(icadet_arag_btf)'
  'r2dvar'
  1.0d0

'cadet_calc_btf'
  'calcite flux to Sediments'
  'mol m-2 s-1'
  'cadet_calc_btf, scalar, series'
  'ocean_time'
  'iDbio2(icadet_calc_btf)'
  'r2dvar'
  1.0d0

'ndet_btf'
  'N flux to Sediments'
  'mol m-2 s-1'
  'ndet_btf, scalar, series'
  'ocean_time'
  'iDbio2(indet_btf)'
  'r2dvar'
  1.0d0

'pdet_btf'
  'P flux to Sediments'
  'mol m-2 s-1'
  'pdet_btf, scalar, series'
  'ocean_time'
  'iDbio2(ipdet_btf)'
  'r2dvar'
  1.0d0

'sidet_btf'
  'SiO2 flux to Sediments'
  'mol m-2 s-1'
  'sidet_btf, scalar, series'
  'ocean_time'
  'iDbio2(isidet_btf)'
  'r2dvar'
  1.0d0

'mxl_depth'
  'depth of the mixed layer hsbl'
  'm'
  'mxl_depth, scalar, series'
  'ocean_time'
  'iDbio2(imxl_depth)'
  'r2dvar'
  1.0d0

'mxl_level'
  'level of the mixed layer hsbl'
  'm'
  'mxl_level, scalar, series'
  'ocean_time'
  'iDbio2(imxl_level)'
  'r2dvar'
  1.0d0

'alk_btf'
  'Alkalinity bottom flux'
  'mol m-2 s-1'
  'alk_btf, scalar, series'
  'ocean_time'
  'iDbio2(ialk_btf)'
  'r2dvar'
  1.0d0

'dic_btf'
  'DIC bottom flux'
  'mol m-2 s-1'
  'dic_btf, scalar, series'
  'ocean_time'
  'iDbio2(idic_btf)'
  'r2dvar'
  1.0d0

'fed_btf'
  'fed bottom flux'
  'mol m-2 s-1'
  'fed_btf, scalar, series'
  'ocean_time'
  'iDbio2(ifed_btf)'
  'r2dvar'
  1.0d0

'nh4_btf'
  'Ammonia bottom flux'
  'mol m-2 s-1'
  'nh4_btf, scalar, series'
  'ocean_time'
  'iDbio2(inh4_btf)'
  'r2dvar'
  1.0d0

'no3_btf'
  'Nitrates bottom flux'
  'mol m-2 s-1'
  'no3_btf, scalar, series'
  'ocean_time'
  'iDbio2(ino3_btf)'
  'r2dvar'
  1.0d0

'o2_btf'
  'Oxygen bottom flux'
  'mol m-2 s-1'
  'o2_btf, scalar, series'
  'ocean_time'
  'iDbio2(io2_btf)'
  'r2dvar'
  1.0d0

'po4_btf'
  'Phosphate bottom flux'
  'mol m-2 s-1'
  'po4_btf, scalar, series'
  'ocean_time'
  'iDbio2(ipo4_btf)'
  'r2dvar'
  1.0d0

'sio4_btf'
  'Sillicate bottom flux'
  'mol m-2 s-1'
  'sio4_btf, scalar, series'
  'ocean_time'
  'iDbio2(isio4_btf)'
  'r2dvar'
  1.0d0

'jprod_ndet_100'
  'production of ndet in upper 100 meters'
  'check units'
  'jprod_ndet_100, scalar, series'
  'ocean_time'
  'iDbio2(ijprod_ndet_100)'
  'r2dvar'
  1.0d0

'jremin_ndet_100'
  'remineralization of ndet in upper 100 meters'
  'check units'
  'jremin_ndet_100, scalar, series'
  'ocean_time'
  'iDbio2(ijremin_ndet_100)'
  'r2dvar'
  1.0d0

'alpha'
  'Solubility of CO2 for air'
  'mol/kg/atm'
  'alpha, scalar, series'
  'ocean_time'
  'iDbio2(ialpha)'
  'r2dvar'
  1.0d0

'co2star'
  'CO2*water, or H2CO3 concentration'
  'mol/kg'
  'co2star, scalar, series'
  'ocean_time'
  'iDbio2(ico2star)'
  'r2dvar'
  1.0d0

'pco2surf'
  'oceanic pCO2'
  'ppmv'
  'pco2surf, scalar, series'
  'ocean_time'
  'iDbio2(ipco2surf)'
  'r2dvar'
  1.0d0

'pCO2'
  'oceanic pCO2'
  'ppmv'
  'pco2, scalar, series'
  'ocean_time'
  'iDbio3(ipCO2)'
  'r3dvar'
  1.0d0

'co2_flx'
  'air-sea fluxes in terms of contrib to DIC'
  'mol/kg'
  'co2_flx, scalar, series'
  'ocean_time'
  'iDbio2(ico2_flx)'
  'r2dvar'
  1.0d0

'o2_flx'
  'air-sea fluxes in terms of contrib to O2'
  'mol/kg'
  'o2_flx, scalar, series'
  'ocean_time'
  'iDbio2(io2_flx)'
  'r2dvar'
  1.0d0

'ironsed_flx'
  'iron fluxes from the sediments'
  'mol/kg'
  'ironsed_flx, scalar, series'
  'ocean_time'
  'iDbio2(iironsed_flx)'
  'r2dvar'
  1.0d0

'npp_100'
  'Net Primary production integral 100 meters'
  'mgC/m2/day'
  'npp_100, scalar, series'
  'ocean_time'
  'iDbio2(inpp_100)'
  'r2dvar'
  1.0d0

'jno3denit_wc_vint'
  'sedimentary denitrification'
  'mgC/m2/day'
  'jno3denit_wc_vint, scalar, series'
  'ocean_time'
  'iDbio2(ijno3denit_wc_vint)'
  'r2dvar'
  1.0d0

'mesozoo_200'
  'Mesozoo Biomass integral 200 meters'
  'mgC/m2'
  'mesozoo_200, scalar, series'
  'ocean_time'
  'iDbio2(imesozoo_200)'
  'r2dvar'
  1.0d0

'prod_n_100_sm'
  'nitrogen production for small phyto upper 100m'
  'mol.m-2.s-1'
  'prod_n_100_sm, scalar, series'
  'ocean_time'
  'iDbio2(iprod_n_100_sm)'
  'r2dvar'
  1.0d0

'aggloss_n_100_sm'
  'loss to aggregation for small phyto upper 100m'
  'mol.m-2.s-1'
  'aggloss_n_100_sm, scalar, series'
  'ocean_time'
  'iDbio2(iaggloss_n_100_sm)'
  'r2dvar'
  1.0d0

'zloss_n_100_sm'
  'loss to zooplankton for small phyto upper 100m'
  'mol.m-2.s-1'
  'zloss_n_100_sm, scalar, series'
  'ocean_time'
  'iDbio2(izloss_n_100_sm)'
  'r2dvar'
  1.0d0

'prod_n_100_lg'
  'nitrogen production for large phyto upper 100m'
  'mol.m-2.s-1'
  'prod_n_100_lg, scalar, series'
  'ocean_time'
  'iDbio2(iprod_n_100_lg)'
  'r2dvar'
  1.0d0

'aggloss_n_100_lg'
  'loss to aggregation for large phyto upper 100m'
  'mol.m-2.s-1'
  'aggloss_n_100_lg, scalar, series'
  'ocean_time'
  'iDbio2(iaggloss_n_100_lg)'
  'r2dvar'
  1.0d0

'zloss_n_100_lg'
  'loss to zooplankton for large phyto upper 100m'
  'mol.m-2.s-1'
  'zloss_n_100_lg, scalar, series'
  'ocean_time'
  'iDbio2(izloss_n_100_lg)'
  'r2dvar'
  1.0d0

'prod_n_100_di'
  'nitrogen production for diazo phyto upper 100m'
  'mol.m-2.s-1'
  'prod_n_100_di, scalar, series'
  'ocean_time'
  'iDbio2(iprod_n_100_di)'
  'r2dvar'
  1.0d0

'aggloss_n_100_di'
  'loss to aggregation for diazo phyto upper 100m'
  'mol.m-2.s-1'
  'aggloss_n_100_di, scalar, series'
  'ocean_time'
  'iDbio2(iaggloss_n_100_di)'
  'r2dvar'
  1.0d0

'zloss_n_100_di'
  'loss to zooplankton for diazo phyto upper 100m'
  'mol.m-2.s-1'
  'zloss_n_100_di, scalar, series'
  'ocean_time'
  'iDbio2(izloss_n_100_di)'
  'r2dvar'
  1.0d0

'prod_n_100_smz'
  'nitrogen production for small zoo upper 100m'
  'mol.m-2.s-1'
  'prod_n_100_smz, scalar, series'
  'ocean_time'
  'iDbio2(iprod_n_100_smz)'
  'r2dvar'
  1.0d0

'ingest_n_100_smz'
  'ingestion for small zoo upper 100m'
  'mol.m-2.s-1'
  'ingest_n_100_smz, scalar, series'
  'ocean_time'
  'iDbio2(iingest_n_100_smz)'
  'r2dvar'
  1.0d0

'zloss_n_100_smz'
  'loss to zooplankton for small zoo upper 100m'
  'mol.m-2.s-1'
  'zloss_n_100_smz, scalar, series'
  'ocean_time'
  'iDbio2(izloss_n_100_smz)'
  'r2dvar'
  1.0d0

'hploss_n_100_smz'
  'loss to highe predators for small zoo upper 100m'
  'mol.m-2.s-1'
  'hploss_n_100_smz, scalar, series'
  'ocean_time'
  'iDbio2(ihploss_n_100_smz)'
  'r2dvar'
  1.0d0

'prod_ndet_100_smz'
  'detrital production for small zoo upper 100m'
  'mol.m-2.s-1'
  'prod_ndet_100_smz, scalar, series'
  'ocean_time'
  'iDbio2(iprod_ndet_100_smz)'
  'r2dvar'
  1.0d0

'prod_n_100_mdz'
  'nitrogen production for medium zoo upper 100m'
  'mol.m-2.s-1'
  'prod_n_100_mdz, scalar, series'
  'ocean_time'
  'iDbio2(iprod_n_100_mdz)'
  'r2dvar'
  1.0d0

'ingest_n_100_mdz'
  'ingestion for medium zoo upper 100m'
  'mol.m-2.s-1'
  'ingest_n_100_mdz, scalar, series'
  'ocean_time'
  'iDbio2(iingest_n_100_mdz)'
  'r2dvar'
  1.0d0

'zloss_n_100_mdz'
  'loss to zooplankton for medium zoo upper 100m'
  'mol.m-2.s-1'
  'zloss_n_100_mdz, scalar, series'
  'ocean_time'
  'iDbio2(izloss_n_100_mdz)'
  'r2dvar'
  1.0d0

'hploss_n_100_mdz'
  'loss to highe predators for medium zoo upper 100m'
  'mol.m-2.s-1'
  'hploss_n_100_mdz, scalar, series'
  'ocean_time'
  'iDbio2(ihploss_n_100_mdz)'
  'r2dvar'
  1.0d0

'prod_ndet_100_mdz'
  'detrital production for medium zoo upper 100m'
  'mol.m-2.s-1'
  'prod_ndet_100_mdz, scalar, series'
  'ocean_time'
  'iDbio2(iprod_ndet_100_mdz)'
  'r2dvar'
  1.0d0

'prod_n_100_lgz'
  'nitrogen production for large zoo upper 100m'
  'mol.m-2.s-1'
  'prod_n_100_lgz, scalar, series'
  'ocean_time'
  'iDbio2(iprod_n_100_lgz)'
  'r2dvar'
  1.0d0

'ingest_n_100_lgz'
  'ingestion for large zoo upper 100m'
  'mol.m-2.s-1'
  'ingest_n_100_lgz, scalar, series'
  'ocean_time'
  'iDbio2(iingest_n_100_lgz)'
  'r2dvar'
  1.0d0

'zloss_n_100_lgz'
  'loss to zooplankton for large zoo upper 100m'
  'mol.m-2.s-1'
  'zloss_n_100_lgz, scalar, series'
  'ocean_time'
  'iDbio2(izloss_n_100_lgz)'
  'r2dvar'
  1.0d0

'hploss_n_100_lgz'
  'loss to higher predators for large zoo upper 100m'
  'mol.m-2.s-1'
  'hploss_n_100_lgz, scalar, series'
  'ocean_time'
  'iDbio2(ihploss_n_100_lgz)'
  'r2dvar'
  1.0d0

'prod_ndet_100_lgz'
  'detrital production for large zoo upper 100m'
  'mol.m-2.s-1'
  'prod_ndet_100_lgz, scalar, series'
  'ocean_time'
  'iDbio2(iprod_ndet_100_lgz)'
  'r2dvar'
  1.0d0

'prod_n_100_bact'
  'nitrogen production for bacteria upper 100m'
  'mol.m-2.s-1'
  'prod_n_100_bact, scalar, series'
  'ocean_time'
  'iDbio2(iprod_n_100_bact)'
  'r2dvar'
  1.0d0

'zloss_n_100_bact'
  'loss to zooplankton for bacteria upper 100m'
  'mol.m-2.s-1'
  'zloss_n_100_bact, scalar, series'
  'ocean_time'
  'iDbio2(izloss_n_100_bact)'
  'r2dvar'
  1.0d0

'prod_n_100_md'
  'nitrogen production for medium phyto upper 100m'
  'mol.m-2.s-1'
  'prod_n_100_md, scalar, series'
  'ocean_time'
  'iDbio2(iprod_n_100_md)'
  'r2dvar'
  1.0d0

'aggloss_n_100_md'
  'loss to aggregation for medium phyto upper 100m'
  'mol.m-2.s-1'
  'aggloss_n_100_md, scalar, series'
  'ocean_time'
  'iDbio2(iaggloss_n_100_md)'
  'r2dvar'
  1.0d0

'zloss_n_100_md'
  'loss to zooplankton for medium phyto upper 100m'
  'mol.m-2.s-1'
  'zloss_n_100_md, scalar, series'
  'ocean_time'
  'iDbio2(izloss_n_100_md)'
  'r2dvar'
  1.0d0

'mesozooprod_200'
  'mesozoo production over top 200m'
  'mgC m-2 day-1'
  'mesozooprod_200, scalar, series'
  'ocean_time'
  'iDbio2(imesozooprod_200)'
  'r2dvar'
  1.0d0

'uptake_din_100'
  'uptake DIN over top 100m'
  'mgC m-2 day-1'
  'uptake_din_100, scalar, series'
  'ocean_time'
  'iDbio2(iuptake_din_100)'
  'r2dvar'
  1.0d0

'uptake_no3_n2_100'
  'uptake NO3+N2 over top 100m'
  'mgC m-2 day-1'
  'uptake_no3_n2_100, scalar, series'
  'ocean_time'
  'iDbio2(iuptake_no3_n2_100)'
  'r2dvar'
  1.0d0

'prod_mesozoo_100'
  'mesozoo production over top 100m'
  'mgC m-2 day-1'
  'prod_mesozoo_100, scalar, series'
  'ocean_time'
  'iDbio2(iprod_mesozoo_100)'
  'r2dvar'
  1.0d0

'z_ratio_100'
  'Z-ratio'
  ''
  'z_ratio_100, scalar, series'
  'ocean_time'
  'iDbio2(iz_ratio_100)'
  'r2dvar'
  1.0d0

'pe_ratio_100'
  'pe-ratio'
  ''
  'pe_ratio_100, scalar, series'
  'ocean_time'
  'iDbio2(ipe_ratio_100)'
  'r2dvar'
  1.0d0

'f_ratio_100'
  'f-ratio'
  ''
  'f_ratio_100, scalar, series'
  'ocean_time'
  'iDbio2(if_ratio_100)'
  'r2dvar'
  1.0d0

'prod_don_100_smz'
  'production DON by small zoo over top 100m'
  'mol.m-2.s-1'
  'prod_don_100_smz, scalar, series'
  'ocean_time'
  'iDbio2(iprod_don_100_smz)'
  'r2dvar'
  1.0d0

'prod_don_100_mdz'
  'production DON by medium zoo over top 100m'
  'mol.m-2.s-1'
  'prod_don_100_mdz, scalar, series'
  'ocean_time'
  'iDbio2(iprod_don_100_mdz)'
  'r2dvar'
  1.0d0

'prod_don_100_lgz'
  'production DON by large zoo over top 100m'
  'mol.m-2.s-1'
  'prod_don_100_lgz, scalar, series'
  'ocean_time'
  'iDbio2(iprod_don_100_lgz)'
  'r2dvar'
  1.0d0

'ndet_btm'
  'N delta to Sediments'
  'mol m-2 s-1'
  'ndet_btm, scalar, series'
  'ocean_time'
  'iDbio2(indet_btm)'
  'r2dvar'
  1.0d0

'chl'
  'Chlorophyll'
  'ug kg-1'
  'chl, scalar, series'
  'ocean_time'
  'iDbio3(ichl)'
  'r3dvar'
  1.0d0

'co3_ion'
  'Carbonate ion'
  'mol/kg'
  'co3_ion, scalar, series'
  'ocean_time'
  'iDbio3(ico3_ion)'
  'r3dvar'
  1.0d0

'htotal'
  'H+ ion concentration'
  'mol/kg'
  'htotal, scalar, series'
  'ocean_time'
  'iDbio3(ihtotal)'
  'r3dvar'
  1.0d0

'irr_mem'
  'Irradiance memory'
  'Watts/m^2'
  'irr_mem, scalar, series'
  'ocean_time'
  'iDbio3(iirr_mem)'
  'r3dvar'
  1.0d0

'irr_mix'
  'Irradiance in mixed layer'
  'Watts/m^2'
  'irr_mix, scalar, series'
  'ocean_time'
  'iDbio3(iirr_mix)'
  'r3dvar'
  1.0d0

'irr_inst'
  'Irradiance instantaneous'
  'Watts/m^2'
  'irr_inst, scalar, series'
  'ocean_time'
  'iDbio3(iirr_inst)'
  'r3dvar'
  1.0d0

'jprod_cadet_arag'
  'production term for cadet_arag'
  'N/A'
  'jprod_cadet_arag, scalar, series'
  'ocean_time'
  'iDbio3(ijprod_cadet_arag)'
  'r3dvar'
  1.0d0

'jprod_cadet_calc'
  'production term for cadet_calc'
  'N/A'
  'jprod_cadet_calc, scalar, series'
  'ocean_time'
  'iDbio3(ijprod_cadet_calc)'
  'r3dvar'
  1.0d0

'jprod_fedet'
  'production term for fedet'
  'N/A'
  'jprod_fedet, scalar, series'
  'ocean_time'
  'iDbio3(ijprod_fedet)'
  'r3dvar'
  1.0d0

'jprod_lithdet'
  'production term for lithdet'
  'N/A'
  'jprod_lithdet, scalar, series'
  'ocean_time'
  'iDbio3(ijprod_lithdet)'
  'r3dvar'
  1.0d0

'jprod_ndet'
  'production term for ndet'
  'N/A'
  'jprod_ndet, scalar, series'
  'ocean_time'
  'iDbio3(ijprod_ndet)'
  'r3dvar'
  1.0d0

'jprod_pdet'
  'production term for pdet'
  'N/A'
  'jprod_pdet, scalar, series'
  'ocean_time'
  'iDbio3(ijprod_pdet)'
  'r3dvar'
  1.0d0

'jprod_sidet'
  'production term for sidet'
  'N/A'
  'jprod_sidet, scalar, series'
  'ocean_time'
  'iDbio3(ijprod_sidet)'
  'r3dvar'
  1.0d0

'jdiss_cadet_arag'
  'dissolution term for cadet_arag'
  'N/A'
  'jdiss_cadet_arag, scalar, series'
  'ocean_time'
  'iDbio3(ijdiss_cadet_arag)'
  'r3dvar'
  1.0d0

'jdiss_cadet_calc'
  'dissolution term for cadet_calc'
  'N/A'
  'jdiss_cadet_calc, scalar, series'
  'ocean_time'
  'iDbio3(ijdiss_cadet_calc)'
  'r3dvar'
  1.0d0

'jremin_ndet'
  'remineralization term for ndet'
  'N/A'
  'jremin_ndet, scalar, series'
  'ocean_time'
  'iDbio3(ijremin_ndet)'
  'r3dvar'
  1.0d0

'jremin_pdet'
  'remineralization term for pdet'
  'N/A'
  'jremin_pdet, scalar, series'
  'ocean_time'
  'iDbio3(ijremin_pdet)'
  'r3dvar'
  1.0d0

'jremin_fedet'
  'remineralization term for fedet'
  'N/A'
  'jremin_fedet, scalar, series'
  'ocean_time'
  'iDbio3(ijremin_fedet)'
  'r3dvar'
  1.0d0

'det_jzloss_n'
  'loss of detritus to zoo term for nitrogen'
  'N/A'
  'det_jzloss_n, scalar, series'
  'ocean_time'
  'iDbio3(idet_jzloss_n)'
  'r3dvar'
  1.0d0

'det_jzloss_p'
  'loss of detritus to zoo term for phosphorus'
  'N/A'
  'det_jzloss_p, scalar, series'
  'ocean_time'
  'iDbio3(idet_jzloss_p)'
  'r3dvar'
  1.0d0

'det_jzloss_fe'
  'loss of detritus to zoo term for iron'
  'N/A'
  'det_jzloss_fe, scalar, series'
  'ocean_time'
  'iDbio3(idet_jzloss_fe)'
  'r3dvar'
  1.0d0

'det_jhploss_n'
  'loss of detritus to higher predator term for nitrogen'
  'N/A'
  'det_jhploss_n, scalar, series'
  'ocean_time'
  'iDbio3(idet_jhploss_n)'
  'r3dvar'
  1.0d0

'det_jhploss_p'
  'loss of detritus to higher predator term for phosphorus'
  'N/A'
  'det_jhploss_p, scalar, series'
  'ocean_time'
  'iDbio3(idet_jhploss_p)'
  'r3dvar'
  1.0d0

'det_jhploss_fe'
  'loss of detritus to higher predator term for iron'
  'N/A'
  'det_jhploss_fe, scalar, series'
  'ocean_time'
  'iDbio3(idet_jhploss_fe)'
  'r3dvar'
  1.0d0

'ndet_b4sink'
  'ndet before sinking is applied'
  'N/A'
  'ndet_b4sink, scalar, series'
  'ocean_time'
  'iDbio3(indet_b4sink)'
  'r3dvar'
  1.0d0

'pdet_b4sink'
  'pdet before sinking is applied'
  'N/A'
  'pdet_b4sink, scalar, series'
  'ocean_time'
  'iDbio3(ipdet_b4sink)'
  'r3dvar'
  1.0d0

'fedet_b4sink'
  'fedet before sinking is applied'
  'N/A'
  'fedet_b4sink, scalar, series'
  'ocean_time'
  'iDbio3(ifedet_b4sink)'
  'r3dvar'
  1.0d0

'ndet_afsink'
  'ndet before sinking is applied'
  'N/A'
  'ndet_afsink, scalar, series'
  'ocean_time'
  'iDbio3(indet_afsink)'
  'r3dvar'
  1.0d0

'pdet_afsink'
  'pdet before sinking is applied'
  'N/A'
  'pdet_afsink, scalar, series'
  'ocean_time'
  'iDbio3(ipdet_afsink)'
  'r3dvar'
  1.0d0

'fedet_afsink'
  'fedet before sinking is applied'
  'N/A'
  'fedet_afsink, scalar, series'
  'ocean_time'
  'iDbio3(ifedet_afsink)'
  'r3dvar'
  1.0d0

'ndet_flx'
  'N delta between layers'
  'mol m-2 s-1'
  'ndet_flx, scalar, series'
  'ocean_time'
  'iDbio3(indet_flx)'
  'r3dvar'
  1.0d0

'co3_sol_arag'
  'solubility of aragonite'
  'N/A'
  'co3_sol_arag, scalar, series'
  'ocean_time'
  'iDbio3(ico3_sol_arag)'
  'r3dvar'
  1.0d0

'co3_sol_calc'
  'solubility of calcite'
  'N/A'
  'co3_sol_calc, scalar, series'
  'ocean_time'
  'iDbio3(ico3_sol_calc)'
  'r3dvar'
  1.0d0

'fe_bulk_flx'
  'iron flux from surface and pen'
  'N/A'
  'fe_bulk_flx, scalar, series'
  'ocean_time'
  'iDbio3(ife_bulk_flx)'
  'r3dvar'
  1.0d0

'omega_cadet_calc'
  'omega calcite'
  '[]'
  'omega_cadet_calc, scalar, series'
  'ocean_time'
  'iDbio3(iomega_cadet_calc)'
  'r3dvar'
  1.0d0

'omega_cadet_arag'
  'omega aragonite'
  '[]'
  'omega_cadet_arag, scalar, series'
  'ocean_time'
  'iDbio3(iomega_cadet_arag)'
  'r3dvar'
  1.0d0

'swdk'
  'shortwave decay function for visible (red+blue)'
  '[]'
  'swdk, scalar, series'
  'ocean_time'
  'iDbio3(iswdk)'
  'r3dvar'
  1.0d0

'mu_mem_sm'
  'aggregation memory small phyto'
  '[]'
  'mu_mem_sm, scalar, series'
  'ocean_time'
  'iDbio3(imu_mem_sm)'
  'r3dvar'
  1.0d0

'mu_mem_di'
  'aggregation memory diatoms '
  '[]'
  'mu_mem_di, scalar, series'
  'ocean_time'
  'iDbio3(imu_mem_di)'
  'r3dvar'
  1.0d0

'mu_mem_lg'
  'aggregation memory large phyto '
  '[]'
  'mu_mem_lg, scalar, series'
  'ocean_time'
  'iDbio3(imu_mem_lg)'
  'r3dvar'
  1.0d0

'mu_mem_md'
  'aggregation memory medium phyto '
  '[]'
  'mu_mem_md, scalar, series'
  'ocean_time'
  'iDbio3(imu_mem_md)'
  'r3dvar'
  1.0d0

'agg_lim_sm'
  'aggregation limitation small phyto'
  '[]'
  'agg_lim_sm, scalar, series'
  'ocean_time'
  'iDbio3(iagg_lim_sm)'
  'r3dvar'
  1.0d0

'agg_lim_di'
  'aggregation limitation diatoms '
  '[]'
  'agg_lim_di, scalar, series'
  'ocean_time'
  'iDbio3(iagg_lim_di)'
  'r3dvar'
  1.0d0

'agg_lim_lg'
  'aggregation limitation large phyto '
  '[]'
  'agg_lim_lg, scalar, series'
  'ocean_time'
  'iDbio3(iagg_lim_lg)'
  'r3dvar'
  1.0d0

'agg_lim_md'
  'aggregation limitation medium phyto '
  'mol/kg'
  'agg_lim_md, scalar, series'
  'ocean_time'
  'iDbio3(iagg_lim_md)'
  'r3dvar'
  1.0d0

'aggloss_sm'
  'losses to aggregation small phyto '
  'mol/kg'
  'aggloss_sm, scalar, series'
  'ocean_time'
  'iDbio3(iaggloss_sm)'
  'r3dvar'
  1.0d0

'aggloss_di'
  'losses to aggregation diazo phyto '
  'mol/kg'
  'aggloss_di, scalar, series'
  'ocean_time'
  'iDbio3(iaggloss_di)'
  'r3dvar'
  1.0d0

'aggloss_lg'
  'losses to aggregation large phyto '
  'mol/kg'
  'aggloss_lg, scalar, series'
  'ocean_time'
  'iDbio3(iaggloss_lg)'
  'r3dvar'
  1.0d0

'aggloss_md'
  'losses to aggregation medium phyto '
  'mol/kg'
  'aggloss_md, scalar, series'
  'ocean_time'
  'iDbio3(iaggloss_md)'
  'r3dvar'
  1.0d0

'virloss_lg'
  'losses to viruses large phyto '
  'mol/kg'
  'virloss_lg, scalar, series'
  'ocean_time'
  'iDbio3(ivirloss_lg)'
  'r3dvar'
  1.0d0

'virloss_sm'
  'losses to viruses small phyto '
  'mol/kg'
  'virloss_sm, scalar, series'
  'ocean_time'
  'iDbio3(ivirloss_sm)'
  'r3dvar'
  1.0d0

'virloss_di'
  'losses to viruses diazo phyto '
  'mol/kg'
  'virloss_di, scalar, series'
  'ocean_time'
  'iDbio3(ivirloss_di)'
  'r3dvar'
  1.0d0

'virloss_md'
  'losses to viruses medium phyto '
  'mol/kg'
  'virloss_md, scalar, series'
  'ocean_time'
  'iDbio3(ivirloss_md)'
  'r3dvar'
  1.0d0

'zloss_lg'
  'losses to zoo large phyto '
  'mol/kg'
  'zloss_lg, scalar, series'
  'ocean_time'
  'iDbio3(izloss_lg)'
  'r3dvar'
  1.0d0

'zloss_sm'
  'losses to zoo small phyto '
  'mol/kg'
  'zloss_sm, scalar, series'
  'ocean_time'
  'iDbio3(izloss_sm)'
  'r3dvar'
  1.0d0

'zloss_di'
  'losses to zoo diazo phyto '
  'mol/kg'
  'zloss_di, scalar, series'
  'ocean_time'
  'iDbio3(izloss_di)'
  'r3dvar'
  1.0d0

'zloss_md'
  'losses to zoo medium phyto '
  'mol/kg'
  'zloss_md, scalar, series'
  'ocean_time'
  'iDbio3(izloss_md)'
  'r3dvar'
  1.0d0

'agg_lim_md'
  'aggregation limitation medium phyto '
  '[]'
  'agg_lim_md, scalar, series'
  'ocean_time'
  'iDbio3(iagg_lim_md)'
  'r3dvar'
  1.0d0

'def_fe_sm'
  ''
  '[]'
  'def_fe_sm, scalar, series'
  'ocean_time'
  'iDbio3(idef_fe_sm)'
  'r3dvar'
  1.0d0

'def_fe_di'
  ''
  '[]'
  'def_fe_di, scalar, series'
  'ocean_time'
  'iDbio3(idef_fe_di)'
  'r3dvar'
  1.0d0

'def_fe_lg'
  ''
  '[]'
  'def_fe_lg, scalar, series'
  'ocean_time'
  'iDbio3(idef_fe_lg)'
  'r3dvar'
  1.0d0

'def_fe_md'
  ''
  '[]'
  'def_fe_md, scalar, series'
  'ocean_time'
  'iDbio3(idef_fe_md)'
  'r3dvar'
  1.0d0

'felim_sm'
  ''
  '[]'
  'felim_sm, scalar, series'
  'ocean_time'
  'iDbio3(ifelim_sm)'
  'r3dvar'
  1.0d0

'felim_di'
  ''
  '[]'
  'felim_di, scalar, series'
  'ocean_time'
  'iDbio3(ifelim_di)'
  'r3dvar'
  1.0d0

'felim_lg'
  ''
  '[]'
  'felim_lg, scalar, series'
  'ocean_time'
  'iDbio3(ifelim_lg)'
  'r3dvar'
  1.0d0

'felim_md'
  ''
  '[]'
  'felim_md, scalar, series'
  'ocean_time'
  'iDbio3(ifelim_md)'
  'r3dvar'
  1.0d0

'no3lim_sm'
  ''
  '[]'
  'no3lim_sm, scalar, series'
  'ocean_time'
  'iDbio3(ino3lim_sm)'
  'r3dvar'
  1.0d0

'no3lim_di'
  ''
  '[]'
  'no3lim_di, scalar, series'
  'ocean_time'
  'iDbio3(ino3lim_di)'
  'r3dvar'
  1.0d0

'no3lim_lg'
  ''
  '[]'
  'no3lim_lg, scalar, series'
  'ocean_time'
  'iDbio3(ino3lim_lg)'
  'r3dvar'
  1.0d0

'no3lim_md'
  ''
  '[]'
  'no3lim_md, scalar, series'
  'ocean_time'
  'iDbio3(ino3lim_md)'
  'r3dvar'
  1.0d0

'nh4lim_sm'
  ''
  '[]'
  'nh4lim_sm, scalar, series'
  'ocean_time'
  'iDbio3(inh4lim_sm)'
  'r3dvar'
  1.0d0

'nh4lim_di'
  ''
  '[]'
  'nh4lim_di, scalar, series'
  'ocean_time'
  'iDbio3(inh4lim_di)'
  'r3dvar'
  1.0d0

'nh4lim_lg'
  ''
  '[]'
  'nh4lim_lg, scalar, series'
  'ocean_time'
  'iDbio3(inh4lim_lg)'
  'r3dvar'
  1.0d0

'nh4lim_md'
  ''
  '[]'
  'nh4lim_md, scalar, series'
  'ocean_time'
  'iDbio3(inh4lim_md)'
  'r3dvar'
  1.0d0

'po4lim_sm'
  ''
  '[]'
  'po4lim_sm, scalar, series'
  'ocean_time'
  'iDbio3(ipo4lim_sm)'
  'r3dvar'
  1.0d0

'po4lim_di'
  ''
  '[]'
  'po4lim_di, scalar, series'
  'ocean_time'
  'iDbio3(ipo4lim_di)'
  'r3dvar'
  1.0d0

'po4lim_lg'
  ''
  '[]'
  'po4lim_lg, scalar, series'
  'ocean_time'
  'iDbio3(ipo4lim_lg)'
  'r3dvar'
  1.0d0

'po4lim_md'
  ''
  '[]'
  'po4lim_md, scalar, series'
  'ocean_time'
  'iDbio3(ipo4lim_md)'
  'r3dvar'
  1.0d0

'chl_di'
  ''
  '[ug/kg]'
  'chl_di, scalar, series'
  'ocean_time'
  'iDbio3(ichl_di)'
  'r3dvar'
  1.0d0

'chl_sm'
  ''
  '[ug/kg]'
  'chl_sm, scalar, series'
  'ocean_time'
  'iDbio3(ichl_sm)'
  'r3dvar'
  1.0d0

'chl_lg'
  ''
  '[ug/kg]'
  'chl_lg, scalar, series'
  'ocean_time'
  'iDbio3(ichl_lg)'
  'r3dvar'
  1.0d0

'chl_md'
  ''
  '[ug/kg]'
  'chl_md, scalar, series'
  'ocean_time'
  'iDbio3(ichl_md)'
  'r3dvar'
  1.0d0

'C_2_chl_di'
  ''
  '[ug/kg]'
  'C_2_chl_di, scalar, series'
  'ocean_time'
  'iDbio3(iC_2_chl_di)'
  'r3dvar'
  1.0d0

'C_2_chl_sm'
  ''
  '[ug/kg]'
  'C_2_chl_sm, scalar, series'
  'ocean_time'
  'iDbio3(iC_2_chl_sm)'
  'r3dvar'
  1.0d0

'C_2_chl_lg'
  ''
  '[ug/kg]'
  'C_2_chl_lg, scalar, series'
  'ocean_time'
  'iDbio3(iC_2_chl_lg)'
  'r3dvar'
  1.0d0

'C_2_chl_md'
  ''
  '[ug/kg]'
  'C_2_chl_md, scalar, series'
  'ocean_time'
  'iDbio3(iC_2_chl_md)'
  'r3dvar'
  1.0d0

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


'alk_east'
  'Alkalinity eastern boundary condition'
  'mol/kg'
  'alk_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ialk)'
  'nulvar'
  1.0d0

'alk_west'
  'Alkalinity western boundary condition'
  'mol/kg'
  'alk_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ialk)'
  'nulvar'
  1.0d0

'alk_north'
  'Alkalinity northern boundary condition'
  'mol/kg'
  'alk_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ialk)'
  'nulvar'
  1.0d0

'alk_south'
  'Alkalinity southern boundary condition'
  'mol/kg'
  'alk_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ialk)'
  'nulvar'
  1.0d0


'cadet_arag_east'
  'Detrital CaCO3 eastern boundary condition'
  'mol/kg'
  'cadet_arag_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,icadet_arag)'
  'nulvar'
  1.0d0

'cadet_arag_west'
  'Detrital CaCO3 western boundary condition'
  'mol/kg'
  'cadet_arag_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,icadet_arag)'
  'nulvar'
  1.0d0

'cadet_arag_north'
  'Detrital CaCO3 northern boundary condition'
  'mol/kg'
  'cadet_arag_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,icadet_arag)'
  'nulvar'
  1.0d0

'cadet_arag_south'
  'Detrital CaCO3 southern boundary condition'
  'mol/kg'
  'cadet_arag_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,icadet_arag)'
  'nulvar'
  1.0d0


'cadet_calc_east'
  'Detrital CaCO3 eastern boundary condition'
  'mol/kg'
  'cadet_calc_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,icadet_calc)'
  'nulvar'
  1.0d0

'cadet_calc_west'
  'Detrital CaCO3 western boundary condition'
  'mol/kg'
  'cadet_calc_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,icadet_calc)'
  'nulvar'
  1.0d0

'cadet_calc_north'
  'Detrital CaCO3 northern boundary condition'
  'mol/kg'
  'cadet_calc_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,icadet_calc)'
  'nulvar'
  1.0d0

'cadet_calc_south'
  'Detrital CaCO3 southern boundary condition'
  'mol/kg'
  'cadet_calc_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,icadet_calc)'
  'nulvar'
  1.0d0


'dic_east'
  'Dissolved Inorganic Carbon eastern boundary condition'
  'mol/kg'
  'dic_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,idic)'
  'nulvar'
  1.0d0

'dic_west'
  'Dissolved Inorganic Carbon western boundary condition'
  'mol/kg'
  'dic_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,idic)'
  'nulvar'
  1.0d0

'dic_north'
  'Dissolved Inorganic Carbon northern boundary condition'
  'mol/kg'
  'dic_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,idic)'
  'nulvar'
  1.0d0

'dic_south'
  'Dissolved Inorganic Carbon southern boundary condition'
  'mol/kg'
  'dic_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,idic)'
  'nulvar'
  1.0d0


'fed_east'
  'Dissolved Iron eastern boundary condition'
  'mol/kg'
  'fed_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ifed)'
  'nulvar'
  1.0d0

'fed_west'
  'Dissolved Iron western boundary condition'
  'mol/kg'
  'fed_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ifed)'
  'nulvar'
  1.0d0

'fed_north'
  'Dissolved Iron northern boundary condition'
  'mol/kg'
  'fed_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ifed)'
  'nulvar'
  1.0d0

'fed_south'
  'Dissolved Iron southern boundary condition'
  'mol/kg'
  'fed_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ifed)'
  'nulvar'
  1.0d0


'fedet_east'
  'Detrital Iron eastern boundary condition'
  'mol/kg'
  'fedet_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ifedet)'
  'nulvar'
  1.0d0

'fedet_west'
  'Detrital Iron western boundary condition'
  'mol/kg'
  'fedet_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ifedet)'
  'nulvar'
  1.0d0

'fedet_north'
  'Detrital Iron northern boundary condition'
  'mol/kg'
  'fedet_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ifedet)'
  'nulvar'
  1.0d0

'fedet_south'
  'Detrital Iron southern boundary condition'
  'mol/kg'
  'fedet_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ifedet)'
  'nulvar'
  1.0d0


'fedi_east'
  'Diazotroph Iron eastern boundary condition'
  'mol/kg'
  'fedi_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ifedi)'
  'nulvar'
  1.0d0

'fedi_west'
  'Diazotroph Iron western boundary condition'
  'mol/kg'
  'fedi_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ifedi)'
  'nulvar'
  1.0d0

'fedi_north'
  'Diazotroph Iron northern boundary condition'
  'mol/kg'
  'fedi_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ifedi)'
  'nulvar'
  1.0d0

'fedi_south'
  'Diazotroph Iron southern boundary condition'
  'mol/kg'
  'fedi_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ifedi)'
  'nulvar'
  1.0d0


'felg_east'
  'Large Phytoplankton Iron eastern boundary condition'
  'mol/kg'
  'felg_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ifelg)'
  'nulvar'
  1.0d0

'felg_west'
  'Large Phytoplankton Iron western boundary condition'
  'mol/kg'
  'felg_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ifelg)'
  'nulvar'
  1.0d0

'felg_north'
  'Large Phytoplankton Iron northern boundary condition'
  'mol/kg'
  'felg_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ifelg)'
  'nulvar'
  1.0d0

'felg_south'
  'Large Phytoplankton Iron southern boundary condition'
  'mol/kg'
  'felg_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ifelg)'
  'nulvar'
  1.0d0


'fesm_east'
  'Small Phytoplankton Iron eastern boundary condition'
  'mol/kg'
  'fesm_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ifesm)'
  'nulvar'
  1.0d0

'fesm_west'
  'Small Phytoplankton Iron western boundary condition'
  'mol/kg'
  'fesm_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ifesm)'
  'nulvar'
  1.0d0

'fesm_north'
  'Small Phytoplankton Iron northern boundary condition'
  'mol/kg'
  'fesm_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ifesm)'
  'nulvar'
  1.0d0

'fesm_south'
  'Small Phytoplankton Iron southern boundary condition'
  'mol/kg'
  'fesm_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ifesm)'
  'nulvar'
  1.0d0


'ldon_east'
  'labile DON eastern boundary condition'
  'mol/kg'
  'ldon_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ildon)'
  'nulvar'
  1.0d0

'ldon_west'
  'labile DON western boundary condition'
  'mol/kg'
  'ldon_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ildon)'
  'nulvar'
  1.0d0

'ldon_north'
  'labile DON northern boundary condition'
  'mol/kg'
  'ldon_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ildon)'
  'nulvar'
  1.0d0

'ldon_south'
  'labile DON southern boundary condition'
  'mol/kg'
  'ldon_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ildon)'
  'nulvar'
  1.0d0


'ldop_east'
  'labile DOP eastern boundary condition'
  'mol/kg'
  'ldop_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ildop)'
  'nulvar'
  1.0d0

'ldop_west'
  'labile DOP western boundary condition'
  'mol/kg'
  'ldop_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ildop)'
  'nulvar'
  1.0d0

'ldop_north'
  'labile DOP northern boundary condition'
  'mol/kg'
  'ldop_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ildop)'
  'nulvar'
  1.0d0

'ldop_south'
  'labile DOP southern boundary condition'
  'mol/kg'
  'ldop_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ildop)'
  'nulvar'
  1.0d0


'lith_east'
  'Lithogenic Aluminosilicate eastern boundary condition'
  'g/kg'
  'lith_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ilith)'
  'nulvar'
  1.0d0

'lith_west'
  'Lithogenic Aluminosilicate western boundary condition'
  'g/kg'
  'lith_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ilith)'
  'nulvar'
  1.0d0

'lith_north'
  'Lithogenic Aluminosilicate northern boundary condition'
  'g/kg'
  'lith_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ilith)'
  'nulvar'
  1.0d0

'lith_south'
  'Lithogenic Aluminosilicate southern boundary condition'
  'g/kg'
  'lith_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ilith)'
  'nulvar'
  1.0d0


'lithdet_east'
  'lithdet eastern boundary condition'
  'g/kg'
  'lithdet_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ilithdet)'
  'nulvar'
  1.0d0

'lithdet_west'
  'lithdet western boundary condition'
  'g/kg'
  'lithdet_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ilithdet)'
  'nulvar'
  1.0d0

'lithdet_north'
  'lithdet northern boundary condition'
  'g/kg'
  'lithdet_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ilithdet)'
  'nulvar'
  1.0d0

'lithdet_south'
  'lithdet southern boundary condition'
  'g/kg'
  'lithdet_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ilithdet)'
  'nulvar'
  1.0d0


'nbact_east'
  'bacterial eastern boundary condition'
  'mol/kg'
  'nbact_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,inbact)'
  'nulvar'
  1.0d0

'nbact_west'
  'bacterial western boundary condition'
  'mol/kg'
  'nbact_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,inbact)'
  'nulvar'
  1.0d0

'nbact_north'
  'bacterial northern boundary condition'
  'mol/kg'
  'nbact_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,inbact)'
  'nulvar'
  1.0d0

'nbact_south'
  'bacterial southern boundary condition'
  'mol/kg'
  'nbact_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,inbact)'
  'nulvar'
  1.0d0


'ndet_east'
  'ndet eastern boundary condition'
  'mol/kg'
  'ndet_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,indet)'
  'nulvar'
  1.0d0

'ndet_west'
  'ndet western boundary condition'
  'mol/kg'
  'ndet_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,indet)'
  'nulvar'
  1.0d0

'ndet_north'
  'ndet northern boundary condition'
  'mol/kg'
  'ndet_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,indet)'
  'nulvar'
  1.0d0

'ndet_south'
  'ndet southern boundary condition'
  'mol/kg'
  'ndet_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,indet)'
  'nulvar'
  1.0d0


'ndi_east'
  'Diazotroph Nitrogen eastern boundary condition'
  'mol/kg'
  'ndi_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,indi)'
  'nulvar'
  1.0d0

'ndi_west'
  'Diazotroph Nitrogen western boundary condition'
  'mol/kg'
  'ndi_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,indi)'
  'nulvar'
  1.0d0

'ndi_north'
  'Diazotroph Nitrogen northern boundary condition'
  'mol/kg'
  'ndi_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,indi)'
  'nulvar'
  1.0d0

'ndi_south'
  'Diazotroph Nitrogen southern boundary condition'
  'mol/kg'
  'ndi_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,indi)'
  'nulvar'
  1.0d0


'nlg_east'
  'Large Phytoplankton Nitrogen eastern boundary condition'
  'mol/kg'
  'nlg_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,inlg)'
  'nulvar'
  1.0d0

'nlg_west'
  'Large Phytoplankton Nitrogen western boundary condition'
  'mol/kg'
  'nlg_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,inlg)'
  'nulvar'
  1.0d0

'nlg_north'
  'Large Phytoplankton Nitrogen northern boundary condition'
  'mol/kg'
  'nlg_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,inlg)'
  'nulvar'
  1.0d0

'nlg_south'
  'Large Phytoplankton Nitrogen southern boundary condition'
  'mol/kg'
  'nlg_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,inlg)'
  'nulvar'
  1.0d0


'nsm_east'
  'Small Phytoplankton Nitrogen eastern boundary condition'
  'mol/kg'
  'nsm_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,insm)'
  'nulvar'
  1.0d0

'nsm_west'
  'Small Phytoplankton Nitrogen western boundary condition'
  'mol/kg'
  'nsm_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,insm)'
  'nulvar'
  1.0d0

'nsm_north'
  'Small Phytoplankton Nitrogen northern boundary condition'
  'mol/kg'
  'nsm_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,insm)'
  'nulvar'
  1.0d0

'nsm_south'
  'Small Phytoplankton Nitrogen southern boundary condition'
  'mol/kg'
  'nsm_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,insm)'
  'nulvar'
  1.0d0


'nh4_east'
  'Ammonia eastern boundary condition'
  'mol/kg'
  'nh4_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,inh4)'
  'nulvar'
  1.0d0

'nh4_west'
  'Ammonia western boundary condition'
  'mol/kg'
  'nh4_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,inh4)'
  'nulvar'
  1.0d0

'nh4_north'
  'Ammonia northern boundary condition'
  'mol/kg'
  'nh4_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,inh4)'
  'nulvar'
  1.0d0

'nh4_south'
  'Ammonia southern boundary condition'
  'mol/kg'
  'nh4_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,inh4)'
  'nulvar'
  1.0d0


'no3_east'
  'Nitrate eastern boundary condition'
  'mol/kg'
  'no3_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ino3)'
  'nulvar'
  1.0d0

'no3_west'
  'Nitrate western boundary condition'
  'mol/kg'
  'no3_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ino3)'
  'nulvar'
  1.0d0

'no3_north'
  'Nitrate northern boundary condition'
  'mol/kg'
  'no3_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ino3)'
  'nulvar'
  1.0d0

'no3_south'
  'Nitrate southern boundary condition'
  'mol/kg'
  'no3_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ino3)'
  'nulvar'
  1.0d0


'o2_east'
  'Oxygen eastern boundary condition'
  'mol/kg'
  'o2_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,io2)'
  'nulvar'
  1.0d0

'o2_west'
  'Oxygen western boundary condition'
  'mol/kg'
  'o2_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,io2)'
  'nulvar'
  1.0d0

'o2_north'
  'Oxygen northern boundary condition'
  'mol/kg'
  'o2_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,io2)'
  'nulvar'
  1.0d0

'o2_south'
  'Oxygen southern boundary condition'
  'mol/kg'
  'o2_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,io2)'
  'nulvar'
  1.0d0


'pdet_east'
  'Detrital Phosphorus eastern boundary condition'
  'mol/kg'
  'pdet_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ipdet)'
  'nulvar'
  1.0d0

'pdet_west'
  'Detrital Phosphorus western boundary condition'
  'mol/kg'
  'pdet_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ipdet)'
  'nulvar'
  1.0d0

'pdet_north'
  'Detrital Phosphorus northern boundary condition'
  'mol/kg'
  'pdet_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ipdet)'
  'nulvar'
  1.0d0

'pdet_south'
  'Detrital Phosphorus southern boundary condition'
  'mol/kg'
  'pdet_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ipdet)'
  'nulvar'
  1.0d0


'po4_east'
  'Phosphate eastern boundary condition'
  'mol/kg'
  'po4_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ipo4)'
  'nulvar'
  1.0d0

'po4_west'
  'Phosphate western boundary condition'
  'mol/kg'
  'po4_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ipo4)'
  'nulvar'
  1.0d0

'po4_north'
  'Phosphate northern boundary condition'
  'mol/kg'
  'po4_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ipo4)'
  'nulvar'
  1.0d0

'po4_south'
  'Phosphate southern boundary condition'
  'mol/kg'
  'po4_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ipo4)'
  'nulvar'
  1.0d0


'srdon_east'
  'Semi-Refractory DON eastern boundary condition'
  'mol/kg'
  'srdon_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,isrdon)'
  'nulvar'
  1.0d0

'srdon_west'
  'Semi-Refractory DON western boundary condition'
  'mol/kg'
  'srdon_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,isrdon)'
  'nulvar'
  1.0d0

'srdon_north'
  'Semi-Refractory DON northern boundary condition'
  'mol/kg'
  'srdon_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,isrdon)'
  'nulvar'
  1.0d0

'srdon_south'
  'Semi-Refractory DON southern boundary condition'
  'mol/kg'
  'srdon_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,isrdon)'
  'nulvar'
  1.0d0


'srdop_east'
  'Semi-Refractory DOP eastern boundary condition'
  'mol/kg'
  'srdop_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,isrdop)'
  'nulvar'
  1.0d0

'srdop_west'
  'Semi-Refractory DOP western boundary condition'
  'mol/kg'
  'srdop_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,isrdop)'
  'nulvar'
  1.0d0

'srdop_north'
  'Semi-Refractory DOP northern boundary condition'
  'mol/kg'
  'srdop_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,isrdop)'
  'nulvar'
  1.0d0

'srdop_south'
  'Semi-Refractory DOP southern boundary condition'
  'mol/kg'
  'srdop_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,isrdop)'
  'nulvar'
  1.0d0


'sldon_east'
  'Semilabile DON eastern boundary condition'
  'mol/kg'
  'sldon_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,isldon)'
  'nulvar'
  1.0d0

'sldon_west'
  'Semilabile DON western boundary condition'
  'mol/kg'
  'sldon_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,isldon)'
  'nulvar'
  1.0d0

'sldon_north'
  'Semilabile DON northern boundary condition'
  'mol/kg'
  'sldon_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,isldon)'
  'nulvar'
  1.0d0

'sldon_south'
  'Semilabile DON southern boundary condition'
  'mol/kg'
  'sldon_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,isldon)'
  'nulvar'
  1.0d0


'sldop_east'
  'Semilabile DOP eastern boundary condition'
  'mol/kg'
  'sldop_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,isldop)'
  'nulvar'
  1.0d0

'sldop_west'
  'Semilabile DOP western boundary condition'
  'mol/kg'
  'sldop_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,isldop)'
  'nulvar'
  1.0d0

'sldop_north'
  'Semilabile DOP northern boundary condition'
  'mol/kg'
  'sldop_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,isldop)'
  'nulvar'
  1.0d0

'sldop_south'
  'Semilabile DOP southern boundary condition'
  'mol/kg'
  'sldop_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,isldop)'
  'nulvar'
  1.0d0


'sidet_east'
  'Detrital Silicon eastern boundary condition'
  'mol/kg'
  'sidet_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,isidet)'
  'nulvar'
  1.0d0

'sidet_west'
  'Detrital Silicon western boundary condition'
  'mol/kg'
  'sidet_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,isidet)'
  'nulvar'
  1.0d0

'sidet_north'
  'Detrital Silicon northern boundary condition'
  'mol/kg'
  'sidet_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,isidet)'
  'nulvar'
  1.0d0

'sidet_south'
  'Detrital Silicon southern boundary condition'
  'mol/kg'
  'sidet_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,isidet)'
  'nulvar'
  1.0d0


'silg_east'
  'Large Phytoplankton Silicon eastern boundary condition'
  'mol/kg'
  'silg_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,isilg)'
  'nulvar'
  1.0d0

'silg_west'
  'Large Phytoplankton Silicon western boundary condition'
  'mol/kg'
  'silg_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,isilg)'
  'nulvar'
  1.0d0

'silg_north'
  'Large Phytoplankton Silicon northern boundary condition'
  'mol/kg'
  'silg_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,isilg)'
  'nulvar'
  1.0d0

'silg_south'
  'Large Phytoplankton Silicon southern boundary condition'
  'mol/kg'
  'silg_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,isilg)'
  'nulvar'
  1.0d0


'sio4_east'
  'Silicate eastern boundary condition'
  'mol/kg'
  'sio4_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,isio4)'
  'nulvar'
  1.0d0

'sio4_west'
  'Silicate western boundary condition'
  'mol/kg'
  'sio4_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,isio4)'
  'nulvar'
  1.0d0

'sio4_north'
  'Silicate northern boundary condition'
  'mol/kg'
  'sio4_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,isio4)'
  'nulvar'
  1.0d0

'sio4_south'
  'Silicate southern boundary condition'
  'mol/kg'
  'sio4_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,isio4)'
  'nulvar'
  1.0d0


'nsmz_east'
  'Small Zooplankton Nitrogen eastern boundary condition'
  'mol/kg'
  'nsmz_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,insmz)'
  'nulvar'
  1.0d0

'nsmz_west'
  'Small Zooplankton Nitrogen western boundary condition'
  'mol/kg'
  'nsmz_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,insmz)'
  'nulvar'
  1.0d0

'nsmz_north'
  'Small Zooplankton Nitrogen northern boundary condition'
  'mol/kg'
  'nsmz_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,insmz)'
  'nulvar'
  1.0d0

'nsmz_south'
  'Small Zooplankton Nitrogen southern boundary condition'
  'mol/kg'
  'nsmz_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,insmz)'
  'nulvar'
  1.0d0


'nmdz_east'
  'Medium-sized zooplankton Nitrogen eastern boundary condition'
  'mol/kg'
  'nmdz_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,inmdz)'
  'nulvar'
  1.0d0

'nmdz_west'
  'Medium-sized zooplankton Nitrogen western boundary condition'
  'mol/kg'
  'nmdz_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,inmdz)'
  'nulvar'
  1.0d0

'nmdz_north'
  'Medium-sized zooplankton Nitrogen northern boundary condition'
  'mol/kg'
  'nmdz_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,inmdz)'
  'nulvar'
  1.0d0

'nmdz_south'
  'Medium-sized zooplankton Nitrogen southern boundary condition'
  'mol/kg'
  'nmdz_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,inmdz)'
  'nulvar'
  1.0d0


'nlgz_east'
  'large Zooplankton Nitrogen eastern boundary condition'
  'mol/kg'
  'nlgz_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,inlgz)'
  'nulvar'
  1.0d0

'nlgz_west'
  'large Zooplankton Nitrogen western boundary condition'
  'mol/kg'
  'nlgz_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,inlgz)'
  'nulvar'
  1.0d0

'nlgz_north'
  'large Zooplankton Nitrogen northern boundary condition'
  'mol/kg'
  'nlgz_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,inlgz)'
  'nulvar'
  1.0d0

'nlgz_south'
  'large Zooplankton Nitrogen southern boundary condition'
  'mol/kg'
  'nlgz_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,inlgz)'
  'nulvar'
  1.0d0

! Coastal diatoms (NVO)

'nmd_east'
  'Medium Phytoplankton Nitrogen eastern boundary condition'
  'mol/kg'
  'nmd_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,inmd)'
  'nulvar'
  1.0d0

'nmd_west'
  'Medium Phytoplankton Nitrogen western boundary condition'
  'mol/kg'
  'nmd_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,inmd)'
  'nulvar'
  1.0d0

'nmd_north'
  'Medium Phytoplankton Nitrogen northern boundary condition'
  'mol/kg'
  'nmd_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,inmd)'
  'nulvar'
  1.0d0

'nmd_south'
  'Medium Phytoplankton Nitrogen southern boundary condition'
  'mol/kg'
  'nmd_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,inmd)'
  'nulvar'
  1.0d0

'simd_east'
  'Medium Phytoplankton Silicon eastern boundary condition'
  'mol/kg'
  'simd_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,isimd)'
  'nulvar'
  1.0d0

'simd_west'
  'Medium Phytoplankton Silicon western boundary condition'
  'mol/kg'
  'simd_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,isimd)'
  'nulvar'
  1.0d0

'simd_north'
  'Medium Phytoplankton Silicon northern boundary condition'
  'mol/kg'
  'simd_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,isimd)'
  'nulvar'
  1.0d0

'simd_south'
  'Medium Phytoplankton Silicon southern boundary condition'
  'mol/kg'
  'simd_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,isimd)'
  'nulvar'
  1.0d0

'femd_east'
  'Medium Phytoplankton Iron eastern boundary condition'
  'mol/kg'
  'femd_east, scalar, series'
  'ocean_time'
  'idTbry(ieast,ifemd)'
  'nulvar'
  1.0d0

'femd_west'
  'Medium Phytoplankton Iron western boundary condition'
  'mol/kg'
  'femd_west, scalar, series'
  'ocean_time'
  'idTbry(iwest,ifemd)'
  'nulvar'
  1.0d0

'femd_north'
  'Medium Phytoplankton Iron northern boundary condition'
  'mol/kg'
  'femd_north, scalar, series'
  'ocean_time'
  'idTbry(inorth,ifemd)'
  'nulvar'
  1.0d0

'femd_south'
  'Medium Phytoplankton Iron southern boundary condition'
  'mol/kg'
  'femd_south, scalar, series'
  'ocean_time'
  'idTbry(isouth,ifemd)'
  'nulvar'
  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_provenance'                                   ! Input
  'observation origin'
  'nondimensional'                                 ! [nondimensional]
  'nulvar'
  'datum'
  'idOpro'
  '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

!-------------------------------------------------------------
! Cobalt variables here again ?
! BEST_NPZ variables
!-------------------------------------------------------------


!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'
!  'idBvar(iBen)'
!  'r3dvar'
!  1.0d0

!'DetBen'
!  'benthic detritus concentration'
!  'milligram carbon meter-2'
!  'benthic detritus, scalar, series'
!  'bio_time'
!  'idBvar(iBenDet)'
!  'r3dvar'
!  1.0d0

!!RD : sedidemnt var for cobalt
'cased'
  'Sediment CaCO3'
  'mol.m-3'
  'cased, scalar, series'
  'ocean_time'
  'idBeTvar(icased)'
  'r2dvar'
  1.0d0

'cadet_arag_btf'
  'aragonite flux to Sediments'
  'mol m-2 s-1'
  'cadet_arag_btf, scalar, series'
  'ocean_time'
  'idBeTvar(icadet_arag_btf)'
  'r2dvar'
  1.0d0

'cadet_calc_btf'
  'calcite flux to Sediments'
  'mol m-2 s-1'
  'cadet_calc_btf, scalar, series'
  'ocean_time'
  'idBeTvar(icadet_calc_btf)'
  'r2dvar'
  1.0d0

'ndet_btf'
  'N flux to Sediments'
  'mol m-2 s-1'
  'ndet_btf, scalar, series'
  'ocean_time'
  'idBeTvar(indet_btf)'
  'r2dvar'
  1.0d0

'pdet_btf'
  'P flux to Sediments'
  'mol m-2 s-1'
  'pdet_btf, scalar, series'
  'ocean_time'
  'idBeTvar(ipdet_btf)'
  'r2dvar'
  1.0d0

'sidet_btf'
  'SiO2 flux to Sediments'
  'mol m-2 s-1'
  'sidet_btf, scalar, series'
  'ocean_time'
  'idBeTvar(isidet_btf)'
  'r2dvar'
  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'
  'snow_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

!
!------------------------------------------------------------------------------
!  Nudging to climatology time scales. We can have a generic variable for all
!  tracer nudging time scales or individual variables for each active and
!  passive tracer. The metadata for the individual variables is set internally
!  in "mod_ncparam".
!------------------------------------------------------------------------------
!

'M2_NudgeCoef'                                     ! Input
  '2D momentum inverse nudging coefficients'
  'day-1'                                          ! [1/s]
  'M2_NudgeCoef, scalar'
  'nulvar'
  'idM2nc'
  'r2dvar'
  1.0

'M3_NudgeCoef'                                     ! Input
  '3D momentum inverse nudging coefficients'
  'day-1'                                          ! [1/s]
  'M3_NudgeCoef, scalar'
  'nulvar'
  'idM3nc'
  'r3dvar'
  1.0

'tracer_NudgeCoef'                                 ! Input
  'generic tracer inverse nudging coefficients'
  'day-1'                                          ! [1/s]
  'tracer_NudgeCoef, scalar'
  'nulvar'
  'idgTnc'
  'r3dvar'
  1.0