def_avg_mod Module Reference

Functions/Subroutines
subroutine, public	def_avg (ng, ldef)
subroutine, private	def_avg_nf90 (ng, model, ldef)
subroutine, private	def_avg_pio (ng, model, ldef)

Function/Subroutine Documentation

def_avg()

subroutine, public def_avg_mod::def_avg	(	integer, intent(in)	ng,
		logical, intent(in)	ldef )

Definition at line 78 of file def_avg.F.

!***********************************************************************
!
!  Imported variable declarations.
!
      logical, intent(in) :: ldef
!
      integer, intent(in) :: ng
!
!  Local variable declarations.
!
      integer :: model
!
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__
!
!-----------------------------------------------------------------------
!  Create a new time-averaged file according to IO type.
!-----------------------------------------------------------------------
!
# ifdef ADJOINT
      model=iadm
# else
      model=inlm
# endif
!
      SELECT CASE (avg(ng)%IOtype)
        CASE (io_nf90)
          CALL def_avg_nf90 (ng, model, ldef)
 
# if defined PIO_LIB && defined DISTRIBUTE
        CASE (io_pio)
          CALL def_avg_pio (ng, model, ldef)
# endif
        CASE DEFAULT
          IF (master) WRITE (stdout,10) avg(ng)%IOtype
          exit_flag=3
      END SELECT
      IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
  10  FORMAT (' DEF_AVG - Illegal output file type, io_type = ',i0,     &
     &        /,11x,'Check KeyWord ''OUT_LIB'' in ''roms.in''.')
!
      RETURN

References mod_iounits::avg, def_avg_nf90(), def_avg_pio(), mod_scalars::exit_flag, strings_mod::founderror(), mod_param::iadm, mod_param::inlm, mod_ncparam::io_nf90, mod_ncparam::io_pio, mod_parallel::master, mod_scalars::noerror, and mod_iounits::stdout.

Referenced by ad_output(), output(), rp_output(), and tl_output().

Here is the call graph for this function:

Here is the caller graph for this function:

def_avg_nf90()

subroutine, private def_avg_mod::def_avg_nf90 ( integer, intent(in) ng, integer, intent(in) model, logical, intent(in) ldef )

private

Definition at line 125 of file def_avg.F.

!***********************************************************************
!
      USE mod_netcdf
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, model
!
      logical, intent(in) :: ldef
!
!  Local variable declarations.
!
      logical :: got_var(NV)
!
      integer, parameter :: Natt = 25
 
      integer :: i, ifield, itrc, j, nvd3, nvd4
# if defined WRITE_WATER && defined MASKING
      integer :: xy_pdim, xyz_pdim
# endif
      integer :: recdim, status
 
      integer :: DimIDs(nDimID)
      integer :: p2dgrd(3), t2dgrd(3), u2dgrd(3), v2dgrd(3)
 
# ifdef SOLVE3D
#  ifdef SEDIMENT
      integer :: b3dgrd(4)
#  endif
      integer :: p3dgrd(4), t3dgrd(4), u3dgrd(4), v3dgrd(4), w3dgrd(4)
# endif
!
      real(r8) :: Aval(6)
!
# ifdef ADJOINT
      character (len=21)     :: Prefix
# else
      character (len=13)     :: Prefix
# endif
      character (len=256)    :: ncname
      character (len=MaxLen) :: Vinfo(Natt)
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", def_avg_nf90"
!
      sourcefile=myfile
!
!-----------------------------------------------------------------------
!  Set and report file name.
!-----------------------------------------------------------------------
!
      IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ncname=avg(ng)%name
!
      IF (master) THEN
        IF (ldef) THEN
          WRITE (stdout,10) ng, trim(ncname)
        ELSE
          WRITE (stdout,20) ng, trim(ncname)
        END IF
      END IF
!
!=======================================================================
!  Create a new averages file.
!=======================================================================
!
      define : IF (ldef) THEN
        CALL netcdf_create (ng, model, trim(ncname), avg(ng)%ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
          IF (master) WRITE (stdout,30) trim(ncname)
          RETURN
        END IF
!
!-----------------------------------------------------------------------
!  Define file dimensions.
!-----------------------------------------------------------------------
!
        dimids=0
!
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xi_rho',       &
     &                 iobounds(ng)%xi_rho, dimids( 1))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xi_u',         &
     &                 iobounds(ng)%xi_u, dimids( 2))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xi_v',         &
     &                 iobounds(ng)%xi_v, dimids( 3))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xi_psi',       &
     &                 iobounds(ng)%xi_psi, dimids( 4))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'eta_rho',      &
     &                 iobounds(ng)%eta_rho, dimids( 5))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'eta_u',        &
     &                 iobounds(ng)%eta_u, dimids( 6))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'eta_v',        &
     &                 iobounds(ng)%eta_v, dimids( 7))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'eta_psi',      &
     &                 iobounds(ng)%eta_psi, dimids( 8))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
# if defined WRITE_WATER && defined MASKING
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xy_rho',       &
     &                 iobounds(ng)%xy_rho, dimids(17))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xy_u',         &
     &                 iobounds(ng)%xy_u, dimids(18))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xy_v',         &
     &                 iobounds(ng)%xy_v, dimids(19))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xy_psi',       &
     &                 iobounds(ng)%xy_psi, xy_pdim)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
 
# ifdef SOLVE3D
#  if defined WRITE_WATER && defined MASKING
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xyz_rho',      &
     &                 iobounds(ng)%xy_rho*n(ng), dimids(20))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xyz_u',        &
     &                 iobounds(ng)%xy_u*n(ng), dimids(21))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xyz_v',        &
     &                 iobounds(ng)%xy_v*n(ng), dimids(22))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xyz_w',        &
     &                 iobounds(ng)%xy_rho*(n(ng)+1), dimids(23))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xyz_psi',      &
     &                 iobounds(ng)%xy_psi*n(ng), xyz_pdim)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 's_rho',        &
     &                 n(ng), dimids( 9))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 's_w',          &
     &                 n(ng)+1, dimids(10))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'tracer',       &
     &                 nt(ng), dimids(11))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
#  ifdef SEDIMENT
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'NST',          &
     &                 nst, dimids(32))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'Nbed',         &
     &                 nbed, dimids(16))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
#   if defined WRITE_WATER && defined MASKING
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'xybed',        &
     &                 iobounds(ng)%xy_rho*nbed, dimids(24))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#   endif
#  endif
 
#  ifdef ECOSIM
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'Nbands',       &
     &                 nbands, dimids(33))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'Nphy',         &
     &                 nphy, dimids(25))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'Nbac',         &
     &                 nbac, dimids(26))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'Ndom',         &
     &                 ndom, dimids(27))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'Nfec',         &
     &                 nfec, dimids(28))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
# endif
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'boundary',     &
     &                 4, dimids(14))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
# ifdef FOUR_DVAR
        status=def_dim(ng, model, avg(ng)%ncid, ncname, 'Nstate',       &
     &                 nstatevar(ng), dimids(29))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
 
        status=def_dim(ng, model, avg(ng)%ncid, ncname,                 &
     &                 trim(adjustl(vname(5,idtime))),                  &
     &                 nf90_unlimited, dimids(12))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        recdim=dimids(12)
!
!  Set number of dimensions for output variables.
!
# if defined WRITE_WATER && defined MASKING
        nvd3=2
        nvd4=2
# else
        nvd3=3
        nvd4=4
# endif
!
!  Define dimension vectors for staggered tracer type variables.
!
# if defined WRITE_WATER && defined MASKING
        t2dgrd(1)=dimids(17)
        t2dgrd(2)=dimids(12)
#  ifdef SOLVE3D
        t3dgrd(1)=dimids(20)
        t3dgrd(2)=dimids(12)
#  endif
# else
        t2dgrd(1)=dimids( 1)
        t2dgrd(2)=dimids( 5)
        t2dgrd(3)=dimids(12)
#  ifdef SOLVE3D
        t3dgrd(1)=dimids( 1)
        t3dgrd(2)=dimids( 5)
        t3dgrd(3)=dimids( 9)
        t3dgrd(4)=dimids(12)
#  endif
# endif
!
!  Define dimension vectors for staggered u-momentum type variables.
!
# if defined WRITE_WATER && defined MASKING
        u2dgrd(1)=dimids(18)
        u2dgrd(2)=dimids(12)
#  ifdef SOLVE3D
        u3dgrd(1)=dimids(21)
        u3dgrd(2)=dimids(12)
#  endif
# else
        u2dgrd(1)=dimids( 2)
        u2dgrd(2)=dimids( 6)
        u2dgrd(3)=dimids(12)
#  ifdef SOLVE3D
        u3dgrd(1)=dimids( 2)
        u3dgrd(2)=dimids( 6)
        u3dgrd(3)=dimids( 9)
        u3dgrd(4)=dimids(12)
#  endif
# endif
!
!  Define dimension vectors for staggered v-momentum type variables.
!
# if defined WRITE_WATER && defined MASKING
        v2dgrd(1)=dimids(19)
        v2dgrd(2)=dimids(12)
#  ifdef SOLVE3D
        v3dgrd(1)=dimids(22)
        v3dgrd(2)=dimids(12)
#  endif
# else
        v2dgrd(1)=dimids( 3)
        v2dgrd(2)=dimids( 7)
        v2dgrd(3)=dimids(12)
#  ifdef SOLVE3D
        v3dgrd(1)=dimids( 3)
        v3dgrd(2)=dimids( 7)
        v3dgrd(3)=dimids( 9)
        v3dgrd(4)=dimids(12)
#  endif
# endif
!
!  Define dimension vectors for staggered variables at PSI-points.
!
# if defined WRITE_WATER && defined MASKING
        p2dgrd(1)=xy_pdim
        p2dgrd(2)=dimids(12)
#  ifdef SOLVE3D
        p3dgrd(1)=xyz_pdim
        p3dgrd(2)=dimids(12)
#  endif
# else
        p2dgrd(1)=dimids( 4)
        p2dgrd(2)=dimids( 8)
        p2dgrd(3)=dimids(12)
#  ifdef SOLVE3D
        p3dgrd(1)=dimids( 4)
        p3dgrd(2)=dimids( 8)
        p3dgrd(3)=dimids( 9)
        p3dgrd(4)=dimids(12)
#  endif
# endif
# ifdef SOLVE3D
!
!  Define dimension vector for staggered w-momentum type variables.
!
#  if defined WRITE_WATER && defined MASKING
        w3dgrd(1)=dimids(23)
        w3dgrd(2)=dimids(12)
#  else
        w3dgrd(1)=dimids( 1)
        w3dgrd(2)=dimids( 5)
        w3dgrd(3)=dimids(10)
        w3dgrd(4)=dimids(12)
#  endif
#  ifdef SEDIMENT
!
!  Define dimension vector for sediment bed layer type variables.
!
#   if defined WRITE_WATER && defined MASKING
        b3dgrd(1)=dimids(24)
        b3dgrd(2)=dimids(12)
#   else
        b3dgrd(1)=dimids( 1)
        b3dgrd(2)=dimids( 5)
        b3dgrd(3)=dimids(16)
        b3dgrd(4)=dimids(12)
#   endif
#  endif
# endif
!
!  Initialize unlimited time record dimension.
!
        avg(ng)%Rindex=0
!
!  Initialize local information variable arrays.
!
        DO i=1,natt
          DO j=1,len(vinfo(1))
            vinfo(i)(j:j)=' '
          END DO
        END DO
        DO i=1,6
          aval(i)=0.0_r8
        END DO
!
!  Set long name prefix string.
!
# ifdef ADJOINT
!!      Prefix='time-averaged adjoint'
        prefix='adjoint'
# else
!!      Prefix='time-averaged'
        prefix=char(32)                                   ! blank space
# endif
!
!-----------------------------------------------------------------------
!  Define time-recordless information variables.
!-----------------------------------------------------------------------
!
        CALL def_info (ng, model, avg(ng)%ncid, ncname, dimids)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!-----------------------------------------------------------------------
!  Define time-varying variables.
!-----------------------------------------------------------------------
!
!  Define model time.
!
        vinfo( 1)=vname(1,idtime)
        WRITE (vinfo( 2),'(a,a)') 'averaged ', trim(vname(2,idtime))
        WRITE (vinfo( 3),'(a,a)') 'seconds since ', trim(rclock%string)
        vinfo( 4)=trim(rclock%calendar)
        vinfo(14)=vname(4,idtime)
        vinfo(21)=vname(6,idtime)
        status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idtime),    &
     &                 nf_tout, 1, (/recdim/), aval, vinfo, ncname,     &
     &                 setparaccess = .true.)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!  Define free-surface.
!
        IF (aout(idfsur,ng)) THEN
          vinfo( 1)=vname(1,idfsur)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idfsur))
          vinfo( 3)=vname(3,idfsur)
          vinfo(14)=vname(4,idfsur)
          vinfo(16)=vname(1,idtime)
# if !defined WET_DRY && (defined WRITE_WATER && defined MASKING)
          vinfo(20)='mask_rho'
# endif
          vinfo(21)=vname(6,idfsur)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idfsur,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idfsur),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
# if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided free-surface.
!
        IF (aout(idfsud,ng)) THEN
          vinfo( 1)=vname(1,idfsud)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idfsud))
          vinfo( 3)=vname(3,idfsud)
          vinfo(14)=vname(4,idfsud)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idfsud)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idfsud,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idfsud),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
# endif
!
!  Define 2D momentum in the XI-direction.
!
        IF (aout(idubar,ng)) THEN
          vinfo( 1)=vname(1,idubar)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idubar))
          vinfo( 3)=vname(3,idubar)
          vinfo(14)=vname(4,idubar)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
# endif
          vinfo(21)=vname(6,idubar)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idubar,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idubar),  &
     &                   nf_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
# if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided 2D momentum in the XI-direction.
!
        IF (aout(idu2dd,ng)) THEN
          vinfo( 1)=vname(1,idu2dd)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idu2dd))
          vinfo( 3)=vname(3,idu2dd)
          vinfo(14)=vname(4,idu2dd)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#  endif
          vinfo(21)=vname(6,idu2dd)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idu2dd,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idu2dd),  &
     &                   nf_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
# endif
!
!  Define 2D momentum in the ETA-direction.
!
        IF (aout(idvbar,ng)) THEN
          vinfo( 1)=vname(1,idvbar)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvbar))
          vinfo( 3)=vname(3,idvbar)
          vinfo(14)=vname(4,idvbar)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
# endif
          vinfo(21)=vname(6,idvbar)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvbar,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvbar),  &
     &                   nf_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
# if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided 2D momentum in the ETA-direction.
!
        IF (aout(idv2dd,ng)) THEN
          vinfo( 1)=vname(1,idv2dd)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idv2dd))
          vinfo( 3)=vname(3,idv2dd)
          vinfo(14)=vname(4,idv2dd)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#  endif
          vinfo(21)=vname(6,idv2dd)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idv2dd,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idv2dd),  &
     &                 nf_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
# endif
!
!  Define 2D Eastward momentum component at RHO-points.
!
        IF (aout(idu2de,ng)) THEN
          vinfo( 1)=vname(1,idu2de)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idu2de))
          vinfo( 3)=vname(3,idu2de)
          vinfo(14)=vname(4,idu2de)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
# endif
          vinfo(21)=vname(6,idu2de)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idu2de,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idu2de),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define 2D Northward momentum component at RHO-points.
!
        IF (aout(idv2dn,ng)) THEN
          vinfo( 1)=vname(1,idv2dn)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idv2dn))
          vinfo( 3)=vname(3,idv2dn)
          vinfo(14)=vname(4,idv2dn)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
# endif
          vinfo(21)=vname(6,idv2dn)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idv2dn,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idv2dn),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
# ifdef SOLVE3D
#  ifdef FORWARD_WRITE
!
!  Define time-averaged mass fluxes for 3D momentum coupling.
!
        IF (aout(idufx1,ng)) THEN
          vinfo( 1)=vname(1,idufx1)
          vinfo( 2)=vname(2,idufx1)
          vinfo( 3)=vname(3,idufx1)
          vinfo(14)=vname(4,idufx1)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#   endif
          vinfo(21)=vname(6,idufx1)
          vinfo(22)='coordinates'
          aval(5)=real(u2dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idufx1),  &
     &                   nf_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idufx2,ng)) THEN
          vinfo( 1)=vname(1,idufx2)
          vinfo( 2)=vname(2,idufx2)
          vinfo( 3)=vname(3,idufx2)
          vinfo(14)=vname(4,idufx2)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#   endif
          vinfo(21)=vname(6,idufx2)
          vinfo(22)='coordinates'
          aval(5)=real(u2dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idufx2),  &
     &                   nf_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idvfx1,ng)) THEN
          vinfo( 1)=vname(1,idvfx1)
          vinfo( 2)=vname(2,idvfx1)
          vinfo( 3)=vname(3,idvfx1)
          vinfo(14)=vname(4,idvfx1)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#   endif
          vinfo(21)=vname(6,idvfx1)
          vinfo(22)='coordinates'
          aval(5)=real(v2dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvfx1),  &
     &                   nf_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idvfx2,ng)) THEN
          vinfo( 1)=vname(1,idvfx2)
          vinfo( 2)=vname(2,idvfx2)
          vinfo( 3)=vname(3,idvfx2)
          vinfo(14)=vname(4,idvfx2)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#   endif
          vinfo(21)=vname(6,idvfx2)
          vinfo(22)='coordinates'
          aval(5)=real(v2dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvfx2),  &
     &                   nf_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
!
!  Define 3D momentum component in the XI-direction.
!
        IF (aout(iduvel,ng)) THEN
          vinfo( 1)=vname(1,iduvel)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,iduvel))
          vinfo( 3)=vname(3,iduvel)
          vinfo(14)=vname(4,iduvel)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#  endif
          vinfo(21)=vname(6,iduvel)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,iduvel,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(iduvel),  &
     &                   nf_fout, nvd4, u3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided 3D momentum component in the XI-direction.
!
        IF (aout(idu3dd,ng)) THEN
          vinfo( 1)=vname(1,idu3dd)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idu3dd))
          vinfo( 3)=vname(3,idu3dd)
          vinfo(14)=vname(4,idu3dd)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#   endif
          vinfo(21)=vname(6,idu3dd)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idu3dd,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idu3dd),  &
     &                   nf_fout, nvd4, u3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
!
!  Define 3D momentum component in the ETA-direction.
!
        IF (aout(idvvel,ng)) THEN
          vinfo( 1)=vname(1,idvvel)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvvel))
          vinfo( 3)=vname(3,idvvel)
          vinfo(14)=vname(4,idvvel)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#  endif
          vinfo(21)=vname(6,idvvel)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvvel,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvvel),  &
     &                   nf_fout, nvd4, v3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided 3D momentum component in the ETA-direction.
!
        IF (aout(idv3dd,ng)) THEN
          vinfo( 1)=vname(1,idv3dd)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idv3dd))
          vinfo( 3)=vname(3,idv3dd)
          vinfo(14)=vname(4,idv3dd)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#   endif
          vinfo(21)=vname(6,idv3dd)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idv3dd,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idv3dd),  &
     &                   nf_fout, nvd4, v3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
!
!  Define 3D Eastward momentum component at RHO-points.
!
        IF (aout(idu3de,ng)) THEN
          vinfo( 1)=vname(1,idu3de)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idu3de))
          vinfo( 3)=vname(3,idu3de)
          vinfo(14)=vname(4,idu3de)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idu3de)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idu3de,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idu3de),  &
     &                   nf_fout, nvd4, t3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define 3D Northward momentum component at RHO-points.
!
        IF (aout(idv3dn,ng)) THEN
          vinfo( 1)=vname(1,idv3dn)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idv3dn))
          vinfo( 3)=vname(3,idv3dn)
          vinfo(14)=vname(4,idv3dn)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idv3dn)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idv3dn,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idv3dn),  &
     &                   nf_fout, nvd4, t3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define S-coordinate vertical "omega" momentum component.
!
        IF (aout(idovel,ng)) THEN
          vinfo( 1)=vname(1,idovel)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idovel))
          vinfo( 3)=vname(3,idovel)
          vinfo(14)=vname(4,idovel)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idovel)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idovel,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idovel),  &
     &                   nf_fout, nvd4, w3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define "true" vertical momentum component.
!
        IF (aout(idwvel,ng)) THEN
          vinfo( 1)=vname(1,idwvel)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idwvel))
          vinfo( 3)=vname(3,idwvel)
          vinfo(14)=vname(4,idwvel)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idwvel)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idwvel,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idwvel),  &
     &                   nf_fout, nvd4, w3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define tracer type variables.
!
        DO itrc=1,nt(ng)
          IF (aout(idtvar(itrc),ng)) THEN
            vinfo( 1)=vname(1,idtvar(itrc))
            WRITE (vinfo( 2),'(a,1x,a)') prefix,                        &
     &                                   trim(vname(2,idtvar(itrc)))
            vinfo( 3)=vname(3,idtvar(itrc))
            vinfo(14)=vname(4,idtvar(itrc))
            vinfo(16)=vname(1,idtime)
#  ifdef SEDIMENT
            DO i=1,nst
              IF (itrc.eq.idsed(i)) THEN
                WRITE (vinfo(19),40) 1000.0_r8*sd50(i,ng)
              END IF
            END DO
#  endif
#  if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_rho'
#  endif
            vinfo(21)=vname(6,idtvar(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(r3dvar,r8)
            status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Tid(itrc),  &
     &                     nf_fout, nvd4, t3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
 
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided temperature and salinity.
!
        DO itrc=1,nat
          IF (aout(idtrcd(itrc),ng)) THEN
            ifield=idtrcd(itrc)
            vinfo( 1)=vname(1,ifield)
            WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,ifield))
            vinfo( 3)=vname(3,ifield)
            vinfo(14)=vname(4,ifield)
            vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_rho'
#   endif
            vinfo(21)=vname(6,ifield)
            vinfo(22)='coordinates'
            aval(5)=real(r3dvar,r8)
            status=def_var(ng, model, avg(ng)%ncid,                     &
     &                     avg(ng)%Vid(ifield), nf_fout,                &
     &                     nvd4, t3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
#  endif
!
!  Define density anomaly.
!
        IF (aout(iddano,ng)) THEN
          vinfo( 1)=vname(1,iddano)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,iddano))
          vinfo( 3)=vname(3,iddano)
          vinfo(14)=vname(4,iddano)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,iddano)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,iddano,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(iddano),  &
     &                   nf_fout, nvd4, t3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  ifdef LMD_SKPP
!
!  Define depth of surface boundary layer.
!
        IF (aout(idhsbl,ng)) THEN
          vinfo( 1)=vname(1,idhsbl)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idhsbl))
          vinfo( 3)=vname(3,idhsbl)
          vinfo(14)=vname(4,idhsbl)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idhsbl)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idhsbl,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idhsbl),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
#  ifdef LMD_BKPP
!
!  Define depth of bottom boundary layer.
!
        IF (aout(idhbbl,ng)) THEN
          vinfo( 1)=vname(1,idhbbl)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idhbbl))
          vinfo( 3)=vname(3,idhbbl)
          vinfo(14)=vname(4,idhbbl)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idhbbl)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idhbbl,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idhbbl),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
# endif
!
!  Define 2D potential vorticity.
!
        IF (aout(id2dpv,ng)) THEN
          vinfo( 1)=vname(1,id2dpv)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,id2dpv))
          vinfo( 3)=vname(3,id2dpv)
          vinfo(14)=vname(4,id2dpv)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_psi'
# endif
          vinfo(21)=vname(6,id2dpv)
          vinfo(22)='coordinates'
          aval(5)=real(p2dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(id2dpv),  &
     &                   nf_fout, nvd3, p2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define 2D relative vorticity.
!
        IF (aout(id2drv,ng)) THEN
          vinfo( 1)=vname(1,id2drv)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,id2drv))
          vinfo( 3)=vname(3,id2drv)
          vinfo(14)=vname(4,id2drv)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_psi'
# endif
          vinfo(21)=vname(6,id2drv)
          vinfo(22)='coordinates'
          aval(5)=real(p2dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(id2drv),  &
     &                   nf_fout, nvd3, p2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
# ifdef SOLVE3D
!
!  Define 3D potential vorticity.
!
        IF (aout(id3dpv,ng)) THEN
          vinfo( 1)=vname(1,id3dpv)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,id3dpv))
          vinfo( 3)=vname(3,id3dpv)
          vinfo(14)=vname(4,id3dpv)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_psi'
#  endif
          vinfo(21)=vname(6,id3dpv)
          vinfo(22)='coordinates'
          aval(5)=real(p3dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(id3dpv),  &
     &                   nf_fout, nvd4, p3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define 3D relative vorticity.
!
        IF (aout(id3drv,ng)) THEN
          vinfo( 1)=vname(1,id3drv)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,id3drv))
          vinfo( 3)=vname(3,id3drv)
          vinfo(14)=vname(4,id3drv)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_psi'
#  endif
          vinfo(21)=vname(6,id3drv)
          vinfo(22)='coordinates'
          aval(5)=real(p3dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(id3drv),  &
     &                   nf_fout, nvd4, p3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
# endif
!
!  Define quadratic <zeta*zeta> term.
!
        IF (aout(idzzav,ng)) THEN
          vinfo( 1)=vname(1,idzzav)
          vinfo( 2)=trim(vname(2,idzzav))
          vinfo( 3)=vname(3,idzzav)
          vinfo(14)=vname(4,idzzav)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
# endif
          vinfo(21)=vname(6,idzzav)
          vinfo(22)='coordinates'
          aval(5)=real(r2dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idzzav),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <ubar*ubar> term.
!
        IF (aout(idu2av,ng)) THEN
          vinfo( 1)=vname(1,idu2av)
          vinfo( 2)=trim(vname(2,idu2av))
          vinfo( 3)=vname(3,idu2av)
          vinfo(14)=vname(4,idu2av)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
# endif
          vinfo(21)=vname(6,idu2av)
          vinfo(22)='coordinates'
          aval(5)=real(u2dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idu2av),  &
     &                   nf_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <vbar*vbar> term.
!
        IF (aout(idv2av,ng)) THEN
          vinfo( 1)=vname(1,idv2av)
          vinfo( 2)=trim(vname(2,idv2av))
          vinfo( 3)=vname(3,idv2av)
          vinfo(14)=vname(4,idv2av)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
# endif
          vinfo(21)=vname(6,idv2av)
          vinfo(22)='coordinates'
          aval(5)=real(v2dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idv2av),  &
     &                   nf_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
# ifdef SOLVE3D
!
!  Define u-volume flux.
!
        IF (aout(idhuav,ng)) THEN
          vinfo( 1)=vname(1,idhuav)
          vinfo( 2)=trim(vname(2,idhuav))
          vinfo( 3)=vname(3,idhuav)
          vinfo(14)=vname(4,idhuav)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#  endif
          vinfo(21)=vname(6,idhuav)
          vinfo(22)='coordinates'
          aval(5)=real(u3dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idhuav),  &
     &                   nf_fout, nvd4, u3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define v-volume flux.
!
        IF (aout(idhvav,ng)) THEN
          vinfo( 1)=vname(1,idhvav)
          vinfo( 2)=trim(vname(2,idhvav))
          vinfo( 3)=vname(3,idhvav)
          vinfo(14)=vname(4,idhvav)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#  endif
          vinfo(21)=vname(6,idhvav)
          vinfo(22)='coordinates'
          aval(5)=real(v3dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idhvav),  &
     &                   nf_fout, nvd4, v3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <u*u> term.
!
        IF (aout(iduuav,ng)) THEN
          vinfo( 1)=vname(1,iduuav)
          vinfo( 2)=trim(vname(2,iduuav))
          vinfo( 3)=vname(3,iduuav)
          vinfo(14)=vname(4,iduuav)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#  endif
          vinfo(21)=vname(6,iduuav)
          vinfo(22)='coordinates'
          aval(5)=real(u3dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(iduuav),  &
     &                   nf_fout, nvd4, u3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <u*v> term.
!
        IF (aout(iduvav,ng)) THEN
          vinfo( 1)=vname(1,iduvav)
          vinfo( 2)=trim(vname(2,iduvav))
          vinfo( 3)=vname(3,iduvav)
          vinfo(14)=vname(4,iduvav)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,iduvav)
          vinfo(22)='coordinates'
          aval(5)=real(r3dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(iduvav),  &
     &                   nf_fout, nvd4, t3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <v*v> term.
!
        IF (aout(idvvav,ng)) THEN
          vinfo( 1)=vname(1,idvvav)
          vinfo( 2)=trim(vname(2,idvvav))
          vinfo( 3)=vname(3,idvvav)
          vinfo(14)=vname(4,idvvav)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#  endif
          vinfo(21)=vname(6,idvvav)
          vinfo(22)='coordinates'
          aval(5)=real(v3dvar,r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvvav),  &
     &                   nf_fout, nvd4, v3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <t*t> terms.
!
        DO itrc=1,nt(ng)
          IF (aout(idttav(itrc),ng)) THEN
            vinfo( 1)=vname(1,idttav(itrc))
            vinfo( 2)=trim(vname(2,idttav(itrc)))
            vinfo( 3)=vname(3,idttav(itrc))
            vinfo(14)=vname(4,idttav(itrc))
            vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_rho'
#  endif
            vinfo(21)=vname(6,idttav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(r3dvar,r8)
            status=def_var(ng, model, avg(ng)%ncid,                     &
     &                     avg(ng)%Vid(idttav(itrc)), nf_fout,          &
     &                     nvd4, t3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
!
!  Define active tracers volume fluxes.
!
        DO itrc=1,nt(ng)
          IF (aout(ihutav(itrc),ng)) THEN
            vinfo( 1)=vname(1,ihutav(itrc))
            vinfo( 2)=trim(vname(2,ihutav(itrc)))
            vinfo( 3)=vname(3,ihutav(itrc))
            vinfo(14)=vname(4,ihutav(itrc))
            vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_u'
#  endif
            vinfo(21)=vname(6,ihutav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(u3dvar,r8)
            status=def_var(ng, model, avg(ng)%ncid,                     &
     &                     avg(ng)%Vid(ihutav(itrc)), nf_fout,          &
     &                     nvd4, u3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
!
          IF (aout(ihvtav(itrc),ng)) THEN
            vinfo( 1)=vname(1,ihvtav(itrc))
            vinfo( 2)=trim(vname(2,ihvtav(itrc)))
            vinfo( 3)=vname(3,ihvtav(itrc))
            vinfo(14)=vname(4,ihvtav(itrc))
            vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_v'
#  endif
            vinfo(21)=vname(6,ihvtav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(v3dvar,r8)
            status=def_var(ng, model, avg(ng)%ncid,                     &
     &                     avg(ng)%Vid(ihvtav(itrc)), nf_fout,          &
     &                     nvd4, v3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
!
!  Define quadratic <u*t> and <v*t> terms.
!
        DO itrc=1,nt(ng)
          IF (aout(idutav(itrc),ng)) THEN
            vinfo( 1)=vname(1,idutav(itrc))
            vinfo( 2)=trim(vname(2,idutav(itrc)))
            vinfo( 3)=vname(3,idutav(itrc))
            vinfo(14)=vname(4,idutav(itrc))
            vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_u'
#  endif
            vinfo(21)=vname(6,idutav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(u3dvar,r8)
            status=def_var(ng, model, avg(ng)%ncid,                     &
     &                     avg(ng)%Vid(idutav(itrc)), nf_fout,          &
     &                     nvd4, u3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
!
          IF (aout(idvtav(itrc),ng)) THEN
            vinfo( 1)=vname(1,idvtav(itrc))
            vinfo( 2)=trim(vname(2,idvtav(itrc)))
            vinfo( 3)=vname(3,idvtav(itrc))
            vinfo(14)=vname(4,idvtav(itrc))
            vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_v'
#  endif
            vinfo(21)=vname(6,idvtav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(v3dvar,r8)
            status=def_var(ng, model, avg(ng)%ncid,                     &
     &                     avg(ng)%Vid(idvtav(itrc)), nf_fout,          &
     &                     nvd4, v3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
 
#  if defined LMD_MIXING || defined MY25_MIXING || defined GLS_MIXING
!
!  Define vertical viscosity coefficient.
!
        IF (aout(idvvis,ng)) THEN
          vinfo( 1)=vname(1,idvvis)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvvis))
          vinfo( 3)=vname(3,idvvis)
          vinfo(14)=vname(4,idvvis)
          vinfo(16)=vname(1,idtime)
          vinfo(21)=vname(6,idvvis)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvvis,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvvis),  &
     &                   nf_fout, nvd4, w3dgrd, aval, vinfo, ncname,    &
     &                   setfillval = .false.)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define vertical diffusion coefficient for potential temperature.
!
        IF (aout(idtdif,ng)) THEN
          vinfo( 1)=vname(1,idtdif)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idtdif))
          vinfo( 3)=vname(3,idtdif)
          vinfo(14)=vname(4,idtdif)
          vinfo(16)=vname(1,idtime)
          vinfo(21)=vname(6,idtdif)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idtdif,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idtdif),  &
     &                   nf_fout, nvd4, w3dgrd, aval, vinfo, ncname,    &
     &                   setfillval = .false.)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#   ifdef SALINITY
!
!  Define vertical diffusion coefficient for salinity.
!
        IF (aout(idsdif,ng)) THEN
          vinfo( 1)=vname(1,idsdif)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idsdif))
          vinfo( 3)=vname(3,idsdif)
          vinfo(14)=vname(4,idsdif)
          vinfo(16)=vname(1,idtime)
          vinfo(21)=vname(6,idsdif)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idsdif,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idsdif),  &
     &                   nf_fout, nvd4, w3dgrd, aval, vinfo, ncname,    &
     &                   setfillval = .false.)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
#  endif
#  if defined BULK_FLUXES || defined ECOSIM || defined ATM_PRESS
!
!  Define surface air pressure.
!
        IF (aout(idpair,ng)) THEN
          vinfo( 1)=vname(1,idpair)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idpair))
          vinfo( 3)=vname(3,idpair)
          vinfo(14)=vname(4,idpair)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idpair)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idpair,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idpair),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
#  if defined BULK_FLUXES
!
!  Define surface air temperature.
!
        IF (aout(idtair,ng)) THEN
          vinfo( 1)=vname(1,idtair)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idtair))
          vinfo( 3)=vname(3,idtair)
          vinfo(14)=vname(4,idtair)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idtair)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idtair,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idtair),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
#  if defined BULK_FLUXES || defined ECOSIM
!
!  Define surface winds.
!
        IF (aout(iduair,ng)) THEN
          vinfo( 1)=vname(1,iduair)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,iduair))
          vinfo( 3)=vname(3,iduair)
          vinfo(14)=vname(4,iduair)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,iduair)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,iduair,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(iduair),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idvair,ng)) THEN
          vinfo( 1)=vname(1,idvair)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvair))
          vinfo( 3)=vname(3,idvair)
          vinfo(14)=vname(4,idvair)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idvair)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvair,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvair),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define Eastward/Northward surface wind at RHO-points.
!
        IF (aout(iduaie,ng)) THEN
          vinfo( 1)=vname(1,iduaie)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,iduaie))
          vinfo( 3)=vname(3,iduaie)
          vinfo(14)=vname(4,iduaie)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
# endif
          vinfo(21)=vname(6,iduaie)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,iduaie,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(iduaie),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idvain,ng)) THEN
          vinfo( 1)=vname(1,idvain)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvain))
          vinfo( 3)=vname(3,idvain)
          vinfo(14)=vname(4,idvain)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
# endif
          vinfo(21)=vname(6,idvain)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvain,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvain),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
!
!  Define surface net heat flux.
!
        IF (aout(idtsur(itemp),ng)) THEN
          vinfo( 1)=vname(1,idtsur(itemp))
          WRITE (vinfo( 2),'(a,1x,a)') prefix,                          &
     &                                 trim(vname(2,idtsur(itemp)))
          vinfo( 3)=vname(3,idtsur(itemp))
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idtsur(itemp))
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idtsur(itemp))
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idtsur(itemp),ng),r8)
          status=def_var(ng, model, avg(ng)%ncid,                       &
     &                   avg(ng)%Vid(idtsur(itemp)), nf_fout,           &
     &                   nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  ifdef SALINITY
!
!  Define surface net salt flux.
!
        IF (aout(idtsur(isalt),ng)) THEN
          vinfo( 1)=vname(1,idtsur(isalt))
          WRITE (vinfo( 2),'(a,1x,a)') prefix,                          &
     &                                 trim(vname(2,idtsur(isalt)))
          vinfo( 3)=vname(3,idtsur(isalt))
          vinfo(11)='upward flux, freshening (net precipitation)'
          vinfo(12)='downward flux, salting (net evaporation)'
          vinfo(14)=vname(4,idtsur(isalt))
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idtsur(isalt))
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idtsur(isalt),ng),r8)
          status=def_var(ng, model, avg(ng)%ncid,                       &
     &                   avg(ng)%Vid(idtsur(isalt)), nf_fout,           &
     &                   nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
 
#  if defined BULK_FLUXES || defined FRC_COUPLING
!
!  Define latent heat flux.
!
        IF (aout(idlhea,ng)) THEN
          vinfo( 1)=vname(1,idlhea)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idlhea))
          vinfo( 3)=vname(3,idlhea)
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idlhea)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idlhea)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idlhea,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idlhea),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define sensible heat flux.
!
        IF (aout(idshea,ng)) THEN
          vinfo( 1)=vname(1,idshea)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idshea))
          vinfo( 3)=vname(3,idshea)
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idshea)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idshea)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idshea,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idshea),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define net longwave radiation flux.
!
        IF (aout(idlrad,ng)) THEN
          vinfo( 1)=vname(1,idlrad)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idlrad))
          vinfo( 3)=vname(3,idlrad)
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idlrad)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idlrad)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idlrad,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idlrad),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
 
#  ifdef BULK_FLUXES
#   ifdef EMINUSP
!
!  Define evaporation rate.
!
        IF (aout(idevap,ng)) THEN
          vinfo( 1)=vname(1,idevap)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idevap))
          vinfo( 3)=vname(3,idevap)
          vinfo(11)='downward flux, freshening (condensation)'
          vinfo(12)='upward flux, salting (evaporation)'
          vinfo(14)=vname(4,idevap)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idevap)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idevap,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idevap),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define precipitation rate.
!
        IF (aout(idrain,ng)) THEN
          vinfo( 1)=vname(1,idrain)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idrain))
          vinfo( 3)=vname(3,idrain)
          vinfo(12)='downward flux, freshening (precipitation)'
          vinfo(14)=vname(4,idrain)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idrain)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idrain,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idrain),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
#  endif
#  ifdef SHORTWAVE
!
!  Define net shortwave radiation flux.
!
        IF (aout(idsrad,ng)) THEN
          vinfo( 1)=vname(1,idsrad)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idsrad))
          vinfo( 3)=vname(3,idsrad)
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idsrad)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idsrad)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idsrad,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idsrad),  &
     &                   nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
# endif
!
!  Define surface u-momentum stress.
!
        IF (aout(idusms,ng)) THEN
          vinfo( 1)=vname(1,idusms)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idusms))
          vinfo( 3)=vname(3,idusms)
          vinfo(14)=vname(4,idusms)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
# endif
          vinfo(21)=vname(6,idusms)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idusms,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idusms),  &
     &                   nf_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define surface v-momentum stress.
!
        IF (aout(idvsms,ng)) THEN
          vinfo( 1)=vname(1,idvsms)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvsms))
          vinfo( 3)=vname(3,idvsms)
          vinfo(14)=vname(4,idvsms)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
# endif
          vinfo(21)=vname(6,idvsms)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvsms,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvsms),  &
     &                   nf_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define bottom u-momentum stress.
!
        IF (aout(idubms,ng)) THEN
          vinfo( 1)=vname(1,idubms)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idubms))
          vinfo( 3)=vname(3,idubms)
          vinfo(14)=vname(4,idubms)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
# endif
          vinfo(21)=vname(6,idubms)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idubms,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idubms),  &
     &                   nf_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define bottom v-momentum stress.
!
        IF (aout(idvbms,ng)) THEN
          vinfo( 1)=vname(1,idvbms)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvbms))
          vinfo( 3)=vname(3,idvbms)
          vinfo(14)=vname(4,idvbms)
          vinfo(16)=vname(1,idtime)
# if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
# endif
          vinfo(21)=vname(6,idvbms)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvbms,ng),r8)
          status=def_var(ng, model, avg(ng)%ncid, avg(ng)%Vid(idvbms),  &
     &                   nf_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
# if (defined BBL_MODEL || defined WAVES_OUTPUT) && defined SOLVE3D
!
!-----------------------------------------------------------------------
!  Define the bottom boundary layer model or waves variables.
!-----------------------------------------------------------------------
!
        CALL bbl_def_nf90 (ng, model, ldef, aout, avg,                  &
     &                     t2dgrd, u2dgrd, v2dgrd)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
 
# if defined ICE_MODEL && defined SOLVE3D
!
!-----------------------------------------------------------------------
!  Define the sea-ice model variables.
!-----------------------------------------------------------------------
!
        CALL ice_def_nf90 (ng, model, ldef, aout, avg,                  &
     &                     t2dgrd, u2dgrd, v2dgrd)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
 
# if defined SEDIMENT && defined SOLVE3D
!
!-----------------------------------------------------------------------
!  Define the sediment model variables.
!-----------------------------------------------------------------------
!
        CALL sediment_def_nf90 (ng, model, ldef, aout, avg,             &
     &                          t2dgrd, u2dgrd, v2dgrd,                 &
     &                          b3dgrd)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
 
# if defined WEC_VF && defined SOLVE3D
!
!-----------------------------------------------------------------------
!  Define the Waves Effect on Currents variables.
!-----------------------------------------------------------------------
!
        CALL wec_def_nf90 (ng, model, ldef, aout, avg,                  &
     &                     t2dgrd, u2dgrd, v2dgrd,                      &
     &                     t3dgrd, u3dgrd, v3dgrd, w3dgrd)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
!
!-----------------------------------------------------------------------
!  Leave definition mode.
!-----------------------------------------------------------------------
!
        CALL netcdf_enddef (ng, model, ncname, avg(ng)%ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!-----------------------------------------------------------------------
!  Write out time-recordless, information variables.
!-----------------------------------------------------------------------
!
        CALL wrt_info (ng, model, avg(ng)%ncid, ncname)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      END IF define
!
!=======================================================================
!  Open an existing averages file, check its contents, and prepare
!  for appending data.
!=======================================================================
!
      query : IF (.not.ldef) THEN
        ncname=avg(ng)%name
!
!  Open averages file for read/write.
!
        CALL netcdf_open (ng, model, ncname, 1, avg(ng)%ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
          WRITE (stdout,50) trim(ncname)
          RETURN
        END IF
!
!  Inquire about the dimensions and check for consistency.
!
        CALL netcdf_check_dim (ng, model, ncname,                       &
     &                         ncid = avg(ng)%ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!  Inquire about the variables.
!
        CALL netcdf_inq_var (ng, model, ncname,                         &
     &                       ncid = avg(ng)%ncid)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!  Initialize logical switches.
!
        DO i=1,nv
          got_var(i)=.false.
        END DO
!
!  Scan variable list from input NetCDF and activate switches for
!  average variables. Get variable IDs.
!
        DO i=1,n_var
          IF (trim(var_name(i)).eq.trim(vname(1,idtime))) THEN
            got_var(idtime)=.true.
            avg(ng)%Vid(idtime)=var_id(i)
# ifdef WET_DRY
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idrwet))) THEN
            got_var(idrwet)=.true.
            avg(ng)%Vid(idrwet)=var_id(i)
# endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idfsur))) THEN
            got_var(idfsur)=.true.
            avg(ng)%Vid(idfsur)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubar))) THEN
            got_var(idubar)=.true.
            avg(ng)%Vid(idubar)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbar))) THEN
            got_var(idvbar)=.true.
            avg(ng)%Vid(idvbar)=var_id(i)
# if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idfsud))) THEN
            got_var(idfsud)=.true.
            avg(ng)%Vid(idfsud)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu2dd))) THEN
            got_var(idu2dd)=.true.
            avg(ng)%Vid(idu2dd)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv2dd))) THEN
            got_var(idv2dd)=.true.
            avg(ng)%Vid(idv2dd)=var_id(i)
# endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu2de))) THEN
            got_var(idu2de)=.true.
            avg(ng)%Vid(idu2de)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv2dn))) THEN
            got_var(idv2dn)=.true.
            avg(ng)%Vid(idv2dn)=var_id(i)
# ifdef SOLVE3D
#  ifdef FORWARD_WRITE
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idufx1))) THEN
            got_var(idufx1)=.true.
            avg(ng)%Vid(idufx1)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idufx2))) THEN
            got_var(idufx2)=.true.
            avg(ng)%Vid(idufx2)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvfx1))) THEN
            got_var(idvfx1)=.true.
            avg(ng)%Vid(idvfx1)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvfx2))) THEN
            got_var(idvfx2)=.true.
            avg(ng)%Vid(idvfx2)=var_id(i)
#  endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduvel))) THEN
            got_var(iduvel)=.true.
            avg(ng)%Vid(iduvel)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvvel))) THEN
            got_var(idvvel)=.true.
            avg(ng)%Vid(idvvel)=var_id(i)
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu3dd))) THEN
            got_var(idu3dd)=.true.
            avg(ng)%Vid(idu3dd)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv3dd))) THEN
            got_var(idv3dd)=.true.
            avg(ng)%Vid(idv3dd)=var_id(i)
#  endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu3de))) THEN
            got_var(idu3de)=.true.
            avg(ng)%Vid(idu3de)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv3dn))) THEN
            got_var(idv3dn)=.true.
            avg(ng)%Vid(idv3dn)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idovel))) THEN
            got_var(idovel)=.true.
            avg(ng)%Vid(idovel)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idwvel))) THEN
            got_var(idwvel)=.true.
            avg(ng)%Vid(idwvel)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iddano))) THEN
            got_var(iddano)=.true.
            avg(ng)%Vid(iddano)=var_id(i)
#  ifdef LMD_SKPP
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idhsbl))) THEN
            got_var(idhsbl)=.true.
            avg(ng)%Vid(idhsbl)=var_id(i)
#  endif
#  ifdef LMD_BKPP
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idhbbl))) THEN
            got_var(idhbbl)=.true.
            avg(ng)%Vid(idhbbl)=var_id(i)
#  endif
# endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,id2dpv))) THEN
            got_var(id2dpv)=.true.
            avg(ng)%Vid(id2dpv)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,id2drv))) THEN
            got_var(id2drv)=.true.
            avg(ng)%Vid(id2drv)=var_id(i)
# ifdef SOLVE3D
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,id3dpv))) THEN
            got_var(id3dpv)=.true.
            avg(ng)%Vid(id3dpv)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,id3drv))) THEN
            got_var(id3drv)=.true.
            avg(ng)%Vid(id3drv)=var_id(i)
# endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idzzav))) THEN
            got_var(idzzav)=.true.
            avg(ng)%Vid(idzzav)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu2av))) THEN
            got_var(idu2av)=.true.
            avg(ng)%Vid(idu2av)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv2av))) THEN
            got_var(idv2av)=.true.
            avg(ng)%Vid(idv2av)=var_id(i)
# ifdef SOLVE3D
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idhuav))) THEN
            got_var(idhuav)=.true.
            avg(ng)%Vid(idhuav)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idhvav))) THEN
            got_var(idhvav)=.true.
            avg(ng)%Vid(idhvav)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduuav))) THEN
            got_var(iduuav)=.true.
            avg(ng)%Vid(iduuav)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduvav))) THEN
            got_var(iduvav)=.true.
            avg(ng)%Vid(iduvav)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvvav))) THEN
            got_var(idvvav)=.true.
            avg(ng)%Vid(idvvav)=var_id(i)
#  if defined LMD_MIXING || defined MY25_MIXING || defined GLS_MIXING
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvvis))) THEN
            got_var(idvvis)=.true.
            avg(ng)%Vid(idvvis)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idtdif))) THEN
            got_var(idtdif)=.true.
            avg(ng)%Vid(idtdif)=var_id(i)
#   ifdef SALINITY
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idsdif))) THEN
            got_var(idsdif)=.true.
            avg(ng)%Vid(idsdif)=var_id(i)
#   endif
#  endif
#  if defined BULK_FLUXES || defined ECOSIM || defined ATM_PRESS
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idpair))) THEN
            got_var(idpair)=.true.
            avg(ng)%Vid(idpair)=var_id(i)
#  endif
#  if defined BULK_FLUXES || defined ECOSIM
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduair))) THEN
            got_var(iduair)=.true.
            avg(ng)%Vid(iduair)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvair))) THEN
            got_var(idvair)=.true.
            avg(ng)%Vid(idvair)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduaie))) THEN
            got_var(iduaie)=.true.
            avg(ng)%Vid(iduaie)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvain))) THEN
            got_var(idvain)=.true.
            avg(ng)%Vid(idvain)=var_id(i)
#  endif
          ELSE IF (trim(var_name(i)).eq.                                &
     &             trim(vname(1,idtsur(itemp)))) THEN
            got_var(idtsur(itemp))=.true.
            avg(ng)%Vid(idtsur(itemp))=var_id(i)
#  ifdef SALINITY
          ELSE IF (trim(var_name(i)).eq.                                &
     &             trim(vname(1,idtsur(isalt)))) THEN
            got_var(idtsur(isalt))=.true.
            avg(ng)%Vid(idtsur(isalt))=var_id(i)
#  endif
#  if defined BULK_FLUXES || defined FRC_COUPLING
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idlhea))) THEN
            got_var(idlhea)=.true.
            avg(ng)%Vid(idlhea)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idshea))) THEN
            got_var(idshea)=.true.
            avg(ng)%Vid(idshea)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idlrad))) THEN
            got_var(idlrad)=.true.
            avg(ng)%Vid(idlrad)=var_id(i)
#  endif
#  ifdef BULK_FLUXES
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idtair))) THEN
            got_var(idtair)=.true.
            avg(ng)%Vid(idtair)=var_id(i)
#   ifdef EMINUSP
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idevap))) THEN
            got_var(idevap)=.true.
            avg(ng)%Vid(idevap)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idrain))) THEN
            got_var(idrain)=.true.
            avg(ng)%Vid(idrain)=var_id(i)
#   endif
#  endif
#  ifdef SHORTWAVE
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idsrad))) THEN
            got_var(idsrad)=.true.
            avg(ng)%Vid(idsrad)=var_id(i)
#  endif
# endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idusms))) THEN
            got_var(idusms)=.true.
            avg(ng)%Vid(idusms)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvsms))) THEN
            got_var(idvsms)=.true.
            avg(ng)%Vid(idvsms)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubms))) THEN
            got_var(idubms)=.true.
            avg(ng)%Vid(idubms)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbms))) THEN
            got_var(idvbms)=.true.
            avg(ng)%Vid(idvbms)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubrs))) THEN
            got_var(idubrs)=.true.
            avg(ng)%Vid(idubrs)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbrs))) THEN
            got_var(idvbrs)=.true.
            avg(ng)%Vid(idvbrs)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubws))) THEN
            got_var(idubws)=.true.
            avg(ng)%Vid(idubws)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbws))) THEN
            got_var(idvbws)=.true.
            avg(ng)%Vid(idvbws)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubcs))) THEN
            got_var(idubcs)=.true.
            avg(ng)%Vid(idubcs)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbcs))) THEN
            got_var(idvbcs)=.true.
            avg(ng)%Vid(idvbcs)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduvwc))) THEN
            got_var(iduvwc)=.true.
            avg(ng)%Vid(iduvwc)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubot))) THEN
            got_var(idubot)=.true.
            avg(ng)%Vid(idubot)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbot))) THEN
            got_var(idvbot)=.true.
            avg(ng)%Vid(idvbot)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubur))) THEN
            got_var(idubur)=.true.
            avg(ng)%Vid(idubur)=var_id(i)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbvr))) THEN
            got_var(idvbvr)=.true.
            avg(ng)%Vid(idvbvr)=var_id(i)
          END IF
# ifdef SOLVE3D
          DO itrc=1,nt(ng)
            IF (trim(var_name(i)).eq.trim(vname(1,idtvar(itrc)))) THEN
             got_var(idtvar(itrc))=.true.
             avg(ng)%Tid(itrc)=var_id(i)
            END IF
          END DO
          DO itrc=1,nat
            IF (trim(var_name(i)).eq.trim(vname(1,ihutav(itrc)))) THEN
              got_var(ihutav(itrc))=.true.
              avg(ng)%Vid(ihutav(itrc))=var_id(i)
            ELSE IF (trim(var_name(i)).eq.                              &
     &               trim(vname(1,ihvtav(itrc)))) THEN
              got_var(ihvtav(itrc))=.true.
              avg(ng)%Vid(ihvtav(itrc))=var_id(i)
            ELSE IF (trim(var_name(i)).eq.                              &
     &               trim(vname(1,idutav(itrc)))) THEN
              got_var(idutav(itrc))=.true.
              avg(ng)%Vid(idutav(itrc))=var_id(i)
            ELSE IF (trim(var_name(i)).eq.                              &
     &               trim(vname(1,idvtav(itrc)))) THEN
              got_var(idvtav(itrc))=.true.
              avg(ng)%Vid(idvtav(itrc))=var_id(i)
            ELSE IF (trim(var_name(i)).eq.                              &
     &               trim(vname(1,idttav(itrc)))) THEN
              got_var(idttav(itrc))=.true.
              avg(ng)%Vid(idttav(itrc))=var_id(i)
            END IF
          END DO
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
          DO itrc=1,nat
            IF (trim(var_name(i)).eq.trim(vname(1,idtrcd(itrc)))) THEN
              got_var(idtrcd(itrc))=.true.
              avg(ng)%Vid(idtrcd(itrc))=var_id(i)
            END IF
          END DO
#  endif
# endif
        END DO
!
!  Check if averages variables are available in input NetCDF file.
!
        IF (.not.got_var(idtime)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtime)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idfsur).and.aout(idfsur,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idfsur)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idubar).and.aout(idubar,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idubar)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvbar).and.aout(idvbar,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvbar)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
# if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
        IF (.not.got_var(idfsud).and.aout(idfsud,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idfsud)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idu2dd).and.aout(idu2dd,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu2dd)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv2dd).and.aout(idv2dd,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv2dd)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
# endif
        IF (.not.got_var(idu2de).and.aout(idu2de,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu2de)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv2dn).and.aout(idv2dn,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv2dn)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
# ifdef SOLVE3D
#  ifdef FORWARD_WRITE
        IF (.not.got_var(idufx1).and.aout(idufx1,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idufx1)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idufx2).and.aout(idufx2,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idufx2)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvfx1).and.aout(idvfx1,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvfx1)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvfx2).and.aout(idvfx2,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvfx2)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
        IF (.not.got_var(iduvel).and.aout(iduvel,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduvel)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvvel).and.aout(idvvel,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvvel)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
        IF (.not.got_var(idu3dd).and.aout(idu3dd,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu3dd)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv3dd).and.aout(idv3dd,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv3dd)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
        IF (.not.got_var(idu3de).and.aout(idu3de,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu3de)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv3dn).and.aout(idv3dn,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv3dn)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idovel).and.aout(idovel,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idovel)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idwvel).and.aout(idwvel,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idwvel)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(iddano).and.aout(iddano,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iddano)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  ifdef LMD_SKPP
        IF (.not.got_var(idhsbl).and.aout(idhsbl,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idhsbl)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
#  ifdef LMD_BKPP
        IF (.not.got_var(idhbbl).and.aout(idhbbl,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idhbbl)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
# endif
        IF (.not.got_var(id2dpv).and.aout(id2dpv,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,id2dpv)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(id2drv).and.aout(id2drv,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,id2drv)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
# ifdef SOLVE3D
        IF (.not.got_var(id3dpv).and.aout(id3dpv,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,id3dpv)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(id3drv).and.aout(id3drv,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,id3drv)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
# endif
        IF (.not.got_var(idzzav).and.aout(idzzav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idzzav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idu2av).and.aout(idu2av,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu2av)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv2av).and.aout(idv2av,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv2av)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
# ifdef SOLVE3D
        IF (.not.got_var(idhuav).and.aout(idhuav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idhuav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idhvav).and.aout(idhvav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idhvav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(iduuav).and.aout(iduuav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduuav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(iduvav).and.aout(iduvav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduvav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvvav).and.aout(idvvav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvvav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  if defined LMD_MIXING || defined MY25_MIXING || defined GLS_MIXING
        IF (.not.got_var(idvvis).and.aout(idvvis,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvvis)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idtdif).and.aout(idtdif,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtdif)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   ifdef SALINITY
        IF (.not.got_var(idsdif).and.aout(idsdif,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idsdif)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
#  endif
#  if defined BULK_FLUXES || defined ECOSIM || defined ATM_PRESS
        IF (.not.got_var(idpair).and.aout(idpair,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idpair)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
#  if defined BULK_FLUXES || defined ECOSIM
        IF (.not.got_var(iduair).and.aout(iduair,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduair)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvair).and.aout(idvair,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvair)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(iduaie).and.aout(iduaie,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduaie)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvain).and.aout(idvain,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvain)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
        IF (.not.got_var(idtsur(itemp)).and.aout(idtsur(itemp),ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtsur(itemp))),   &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  ifdef SALINITY
        IF (.not.got_var(idtsur(isalt)).and.aout(idtsur(isalt),ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtsur(isalt))),   &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
#  if defined BULK_FLUXES || defined FRC_COUPLING
        IF (.not.got_var(idlhea).and.aout(idlhea,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idlhea)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idshea).and.aout(idshea,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idshea)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idlrad).and.aout(idlrad,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idlrad)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
#  ifdef BULK_FLUXES
        IF (.not.got_var(idtair).and.aout(idtair,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtair)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   ifdef EMINUSP
        IF (.not.got_var(idevap).and.aout(idevap,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idevap)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idrain).and.aout(idrain,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idrain)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
#  endif
#  ifdef SHORTWAVE
        IF (.not.got_var(idsrad).and.aout(idsrad,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idsrad)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
# endif
        IF (.not.got_var(idusms).and.aout(idusms,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idusms)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvsms).and.aout(idvsms,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvsms)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idubms).and.aout(idubms,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idubms)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvbms).and.aout(idvbms,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvbms)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
# ifdef SOLVE3D
        DO itrc=1,nt(ng)
          IF (.not.got_var(idtvar(itrc)).and.aout(idtvar(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idtvar(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
        END DO
        DO itrc=1,nat
          IF (.not.got_var(ihutav(itrc)).and.aout(ihutav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,ihutav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
          IF (.not.got_var(ihvtav(itrc)).and.aout(ihvtav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,ihvtav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
          IF (.not.got_var(idutav(itrc)).and.aout(idutav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idutav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
          IF (.not.got_var(idvtav(itrc)).and.aout(idvtav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idvtav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
          IF (.not.got_var(idttav(itrc)).and.aout(idttav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idttav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
        END DO
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
        DO itrc=1,nat
          IF (.not.got_var(idtrcd(itrc)).and.aout(idtrcd(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idtrcd(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
        END DO
#  endif
# endif
 
# if (defined BBL_MODEL || defined WAVES_OUTPUT) && defined SOLVE3D
!
!  Scan bottom boundary layer model and waves variables from input
!  NetCDF and activate switches for time-averaged variables. Get
!  variable IDs.
!
        CALL bbl_def_nf90 (ng, model, ldef, aout, avg)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
 
# if defined ICE_MODEL && defined SOLVE3D
!
!  Scan sea-ice variables from input NetCDF and activate switches for
!  time-averaged variables. Get variable IDs.
!
        CALL ice_def_nf90 (ng, model, ldef, aout, avg)
# endif
 
# if defined SEDIMENT && defined SOLVE3D
!
!  Scan sediment model variables from input NetCDF and activate
!  switches for time-averaged variables. Get variable IDs.
!
        CALL sediment_def_nf90 (ng, model, ldef, aout, avg)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
 
# if defined WEC_VF && defined SOLVE3D
!
!  Scan Waves Effect on Currents variables from input NetCDF and
!  activate switches for time-averaged variables. Get variable IDs.
!
        CALL wec_def_nf90 (ng, model, ldef, aout, avg)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
!
!  Set unlimited time record dimension to the appropriate value.
!
        IF (nrst(ng).eq.navg(ng)) THEN
          IF (ndefavg(ng).gt.0) THEN
            avg(ng)%Rindex=((ntstart(ng)-1)-                            &
     &                      ndefavg(ng)*((ntstart(ng)-1)/ndefavg(ng)))/ &
     &                     navg(ng)
          ELSE
            avg(ng)%Rindex=(ntstart(ng)-1)/navg(ng)
          END IF
        ELSE
          avg(ng)%Rindex=rec_size
        END IF
      END IF query
!
!  Set initial average time. Notice that the value is offset by half
!  nAVG*dt so there is not a special case when computing its value
 
# ifdef ADJOINT
!  in "ad_set_avg".
!
        avgtime(ng)=time(ng)+0.5_r8*real(navg(ng),r8)*dt(ng)
# else
!  in "set_avg".
!
      IF (ntsavg(ng).eq.1) THEN
        avgtime(ng)=time(ng)-0.5_r8*real(navg(ng),r8)*dt(ng)
      ELSE
        avgtime(ng)=time(ng)+real(ntsavg(ng),r8)*dt(ng)-                &
     &              0.5_r8*real(navg(ng),r8)*dt(ng)
      END IF
# endif
!
  10  FORMAT (2x,'DEF_AVG_NF90     - creating average file,',t56,       &
     &        'Grid ',i2.2,': ',a)
  20  FORMAT (2x,'DEF_AVG_NF90     - inquiring average file,',t56,      &
     &        'Grid ',i2.2,': ',a)
  30  FORMAT (/,' DEF_AVG_NF90 - unable to create averages NetCDF',     &
     &        ' file: ',a)
  40  FORMAT (1pe11.4,1x,'millimeter')
  50  FORMAT (/,' DEF_AVG_NF90 - unable to open averages NetCDF',       &
     &        ' file: ',a)
  60  FORMAT (/,' DEF_AVG_NF90 - unable to find variable: ',a,2x,       &
     &        ' in averages NetCDF file: ',a)
!
      RETURN

References mod_ncparam::aout, mod_iounits::avg, mod_scalars::avgtime, bbl_output_mod::bbl_def_nf90(), mod_scalars::dt, mod_scalars::exit_flag, strings_mod::founderror(), mod_ncparam::id2dpv, mod_ncparam::id2drv, mod_ncparam::id3dpv, mod_ncparam::id3drv, mod_ncparam::iddano, mod_ncparam::idevap, mod_ncparam::idfsud, mod_ncparam::idfsur, mod_ncparam::idhbbl, mod_ncparam::idhsbl, mod_ncparam::idhuav, mod_ncparam::idhvav, mod_ncparam::idlhea, mod_ncparam::idlrad, mod_ncparam::idovel, mod_ncparam::idpair, mod_ncparam::idrain, mod_ncparam::idrwet, mod_ncparam::idsdif, mod_sediment::idsed, mod_ncparam::idshea, mod_ncparam::idsrad, mod_ncparam::idtair, mod_ncparam::idtdif, mod_ncparam::idtime, mod_ncparam::idtrcd, mod_ncparam::idtsur, mod_ncparam::idttav, mod_ncparam::idtvar, mod_ncparam::idu2av, mod_ncparam::idu2dd, mod_ncparam::idu2de, mod_ncparam::idu3dd, mod_ncparam::idu3de, mod_ncparam::iduaie, mod_ncparam::iduair, mod_ncparam::idubar, mod_ncparam::idubcs, mod_ncparam::idubms, mod_ncparam::idubot, mod_ncparam::idubrs, mod_ncparam::idubur, mod_ncparam::idubws, mod_ncparam::idufx1, mod_ncparam::idufx2, mod_ncparam::idusms, mod_ncparam::idutav, mod_ncparam::iduuav, mod_ncparam::iduvav, mod_ncparam::iduvel, mod_ncparam::iduvwc, mod_ncparam::idv2av, mod_ncparam::idv2dd, mod_ncparam::idv2dn, mod_ncparam::idv3dd, mod_ncparam::idv3dn, mod_ncparam::idvain, mod_ncparam::idvair, mod_ncparam::idvbar, mod_ncparam::idvbcs, mod_ncparam::idvbms, mod_ncparam::idvbot, mod_ncparam::idvbrs, mod_ncparam::idvbvr, mod_ncparam::idvbws, mod_ncparam::idvfx1, mod_ncparam::idvfx2, mod_ncparam::idvsms, mod_ncparam::idvtav, mod_ncparam::idvvav, mod_ncparam::idvvel, mod_ncparam::idvvis, mod_ncparam::idwvel, mod_ncparam::idzzav, mod_ncparam::ihutav, mod_ncparam::ihvtav, mod_ncparam::iinfo, mod_param::iobounds, mod_scalars::isalt, mod_scalars::itemp, mod_parallel::master, mod_param::n, mod_netcdf::n_var, mod_param::nat, mod_scalars::navg, mod_biology::nbac, mod_biology::nbands, mod_param::nbed, mod_scalars::ndefavg, mod_biology::ndom, mod_netcdf::netcdf_check_dim(), mod_netcdf::netcdf_create(), mod_netcdf::netcdf_enddef(), mod_netcdf::netcdf_inq_var(), mod_netcdf::netcdf_open(), mod_netcdf::nf_fout, mod_netcdf::nf_tout, mod_biology::nfec, mod_scalars::noerror, mod_biology::nphy, mod_scalars::nrst, mod_param::nst, mod_fourdvar::nstatevar, mod_param::nt, mod_scalars::ntsavg, mod_scalars::ntstart, mod_ncparam::nv, mod_param::p2dvar, mod_param::p3dvar, mod_param::r2dvar, mod_param::r3dvar, mod_scalars::rclock, mod_netcdf::rec_size, mod_sediment::sd50, sediment_output_mod::sediment_def_nf90(), mod_iounits::sourcefile, mod_iounits::stdout, mod_scalars::time, mod_param::u2dvar, mod_param::u3dvar, mod_param::v2dvar, mod_param::v3dvar, mod_netcdf::var_id, mod_netcdf::var_name, and mod_ncparam::vname.

Referenced by def_avg().

Here is the call graph for this function:

Here is the caller graph for this function:

def_avg_pio()

subroutine, private def_avg_mod::def_avg_pio ( integer, intent(in) ng, integer, intent(in) model, logical, intent(in) ldef )

private

Definition at line 2700 of file def_avg.F.

!***********************************************************************
!
      USE mod_pio_netcdf
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, model
!
      logical, intent(in) :: ldef
!
!  Local variable declarations.
!
      logical :: got_var(NV)
!
      integer, parameter :: Natt = 25
 
      integer :: i, ifield, itrc, j, nvd3, nvd4
#  if defined WRITE_WATER && defined MASKING
      integer :: xy_pdim, xyz_pdim
#  endif
      integer :: recdim, status
 
      integer :: DimIDs(nDimID)
      integer :: p2dgrd(3), t2dgrd(3), u2dgrd(3), v2dgrd(3)
 
#  ifdef SOLVE3D
#   ifdef SEDIMENT
      integer :: b3dgrd(4)
#   endif
      integer :: p3dgrd(4), t3dgrd(4), u3dgrd(4), v3dgrd(4), w3dgrd(4)
#  endif
!
      real(r8) :: Aval(6)
!
#  ifdef ADJOINT
      character (len=21)     :: Prefix
#  else
      character (len=13)     :: Prefix
#  endif
      character (len=256)    :: ncname
      character (len=MaxLen) :: Vinfo(Natt)
 
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__//", def_avg_pio"
!
      sourcefile=myfile
!
!-----------------------------------------------------------------------
!  Set and report file name.
!-----------------------------------------------------------------------
!
      IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      ncname=avg(ng)%name
!
      IF (master) THEN
        IF (ldef) THEN
          WRITE (stdout,10) ng, trim(ncname)
        ELSE
          WRITE (stdout,20) ng, trim(ncname)
        END IF
      END IF
!
!=======================================================================
!  Create a new averages file.
!=======================================================================
!
      define : IF (ldef) THEN
        CALL pio_netcdf_create (ng, model, trim(ncname),                &
     &                          avg(ng)%pioFile)
        IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
          IF (master) WRITE (stdout,30) trim(ncname)
          RETURN
        END IF
!
!-----------------------------------------------------------------------
!  Define file dimensions.
!-----------------------------------------------------------------------
!
        dimids=0
!
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xi_rho',    &
     &                 iobounds(ng)%xi_rho, dimids( 1))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xi_u',      &
     &                 iobounds(ng)%xi_u, dimids( 2))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xi_v',      &
     &                 iobounds(ng)%xi_v, dimids( 3))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xi_psi',    &
     &                 iobounds(ng)%xi_psi, dimids( 4))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'eta_rho',   &
     &                 iobounds(ng)%eta_rho, dimids( 5))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'eta_u',     &
     &                 iobounds(ng)%eta_u, dimids( 6))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'eta_v',     &
     &                 iobounds(ng)%eta_v, dimids( 7))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'eta_psi',   &
     &                 iobounds(ng)%eta_psi, dimids( 8))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
#  if defined WRITE_WATER && defined MASKING
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xy_rho',    &
     &                 iobounds(ng)%xy_rho, dimids(17))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xy_u',      &
     &                 iobounds(ng)%xy_u, dimids(18))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xy_v',      &
     &                 iobounds(ng)%xy_v, dimids(19))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xy_psi',    &
     &                 iobounds(ng)%xy_psi, xy_pdim)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
 
#  ifdef SOLVE3D
#   if defined WRITE_WATER && defined MASKING
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xyz_rho',   &
     &                 iobounds(ng)%xy_rho*n(ng), dimids(20))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xyz_u',     &
     &                 iobounds(ng)%xy_u*n(ng), dimids(21))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xyz_v',     &
     &                 iobounds(ng)%xy_v*n(ng), dimids(22))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xyz_w',     &
     &                 iobounds(ng)%xy_rho*(n(ng)+1), dimids(23))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xyz_psi',   &
     &                 iobounds(ng)%xy_psi*n(ng), xyz_pdim)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#   endif
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 's_rho',     &
     &                 n(ng), dimids( 9))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 's_w',       &
     &                 n(ng)+1, dimids(10))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'tracer',    &
     &                 nt(ng), dimids(11))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
#   ifdef SEDIMENT
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'NST',       &
     &                 nst, dimids(32))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'Nbed',      &
     &                 nbed, dimids(16))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
#    if defined WRITE_WATER && defined MASKING
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'xybed',     &
     &                 iobounds(ng)%xy_rho*nbed, dimids(24))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#    endif
#   endif
 
#   ifdef ECOSIM
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'Nbands',    &
     &                 nbands, dimids(33))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'Nphy',      &
     &                 nphy, dimids(25))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'Nbac',      &
     &                 nbac, dimids(26))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'Ndom',      &
     &                 ndom, dimids(27))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'Nfec',      &
     &                 nfec, dimids(28))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#   endif
#  endif
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'boundary',  &
     &                 4, dimids(14))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
#  ifdef FOUR_DVAR
        status=def_dim(ng, model, avg(ng)%pioFile, ncname, 'Nstate',    &
     &                 nstatevar(ng), dimids(29))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
 
        status=def_dim(ng, model, avg(ng)%pioFile, ncname,              &
     &                 trim(adjustl(vname(5,idtime))),                  &
     &                 pio_unlimited, dimids(12))
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
 
        recdim=dimids(12)
!
!  Set number of dimensions for output variables.
!
#  if defined WRITE_WATER && defined MASKING
        nvd3=2
        nvd4=2
#  else
        nvd3=3
        nvd4=4
#  endif
!
!  Define dimension vectors for staggered tracer type variables.
!
#  if defined WRITE_WATER && defined MASKING
        t2dgrd(1)=dimids(17)
        t2dgrd(2)=dimids(12)
#   ifdef SOLVE3D
        t3dgrd(1)=dimids(20)
        t3dgrd(2)=dimids(12)
#   endif
#  else
        t2dgrd(1)=dimids( 1)
        t2dgrd(2)=dimids( 5)
        t2dgrd(3)=dimids(12)
#   ifdef SOLVE3D
        t3dgrd(1)=dimids( 1)
        t3dgrd(2)=dimids( 5)
        t3dgrd(3)=dimids( 9)
        t3dgrd(4)=dimids(12)
#   endif
#  endif
!
!  Define dimension vectors for staggered u-momentum type variables.
!
#  if defined WRITE_WATER && defined MASKING
        u2dgrd(1)=dimids(18)
        u2dgrd(2)=dimids(12)
#   ifdef SOLVE3D
        u3dgrd(1)=dimids(21)
        u3dgrd(2)=dimids(12)
#   endif
#  else
        u2dgrd(1)=dimids( 2)
        u2dgrd(2)=dimids( 6)
        u2dgrd(3)=dimids(12)
#   ifdef SOLVE3D
        u3dgrd(1)=dimids( 2)
        u3dgrd(2)=dimids( 6)
        u3dgrd(3)=dimids( 9)
        u3dgrd(4)=dimids(12)
#   endif
#  endif
!
!  Define dimension vectors for staggered v-momentum type variables.
!
#  if defined WRITE_WATER && defined MASKING
        v2dgrd(1)=dimids(19)
        v2dgrd(2)=dimids(12)
#   ifdef SOLVE3D
        v3dgrd(1)=dimids(22)
        v3dgrd(2)=dimids(12)
#   endif
#  else
        v2dgrd(1)=dimids( 3)
        v2dgrd(2)=dimids( 7)
        v2dgrd(3)=dimids(12)
#   ifdef SOLVE3D
        v3dgrd(1)=dimids( 3)
        v3dgrd(2)=dimids( 7)
        v3dgrd(3)=dimids( 9)
        v3dgrd(4)=dimids(12)
#   endif
#  endif
!
!  Define dimension vectors for staggered variables at PSI-points.
!
#  if defined WRITE_WATER && defined MASKING
        p2dgrd(1)=xy_pdim
        p2dgrd(2)=dimids(12)
#   ifdef SOLVE3D
        p3dgrd(1)=xyz_pdim
        p3dgrd(2)=dimids(12)
#   endif
#  else
        p2dgrd(1)=dimids( 4)
        p2dgrd(2)=dimids( 8)
        p2dgrd(3)=dimids(12)
#   ifdef SOLVE3D
        p3dgrd(1)=dimids( 4)
        p3dgrd(2)=dimids( 8)
        p3dgrd(3)=dimids( 9)
        p3dgrd(4)=dimids(12)
#   endif
#  endif
#  ifdef SOLVE3D
!
!  Define dimension vector for staggered w-momentum type variables.
!
#   if defined WRITE_WATER && defined MASKING
        w3dgrd(1)=dimids(23)
        w3dgrd(2)=dimids(12)
#   else
        w3dgrd(1)=dimids( 1)
        w3dgrd(2)=dimids( 5)
        w3dgrd(3)=dimids(10)
        w3dgrd(4)=dimids(12)
#   endif
#   ifdef SEDIMENT
!
!  Define dimension vector for sediment bed layer type variables.
!
#    if defined WRITE_WATER && defined MASKING
        b3dgrd(1)=dimids(24)
        b3dgrd(2)=dimids(12)
#    else
        b3dgrd(1)=dimids( 1)
        b3dgrd(2)=dimids( 5)
        b3dgrd(3)=dimids(16)
        b3dgrd(4)=dimids(12)
#    endif
#   endif
#  endif
!
!  Initialize unlimited time record dimension.
!
        avg(ng)%Rindex=0
!
!  Initialize local information variable arrays.
!
        DO i=1,natt
          DO j=1,len(vinfo(1))
            vinfo(i)(j:j)=' '
          END DO
        END DO
        DO i=1,6
          aval(i)=0.0_r8
        END DO
!
!  Set long name prefix string.
!
#  ifdef ADJOINT
!!      Prefix='time-averaged adjoint'
        prefix='adjoint'
#  else
!!      Prefix='time-averaged'
        prefix=char(32)                                   ! blank space
#  endif
!
!-----------------------------------------------------------------------
!  Define time-recordless information variables.
!-----------------------------------------------------------------------
!
        CALL def_info (ng, model, avg(ng)%pioFile, ncname, dimids)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!-----------------------------------------------------------------------
!  Define time-varying variables.
!-----------------------------------------------------------------------
!
!  Define model time.
!
        vinfo( 1)=vname(1,idtime)
        WRITE (vinfo( 2),'(a,a)') 'averaged ', trim(vname(2,idtime))
        WRITE (vinfo( 3),'(a,a)') 'seconds since ', trim(rclock%string)
        vinfo( 4)=trim(rclock%calendar)
        vinfo(14)=vname(4,idtime)
        vinfo(21)=vname(6,idtime)
        avg(ng)%pioVar(idtime)%dkind=pio_tout
        avg(ng)%pioVar(idtime)%gtype=0
!
        status=def_var(ng, model, avg(ng)%pioFile,                      &
     &                 avg(ng)%pioVar(idtime)%vd,                       &
     &                 pio_tout, 1, (/recdim/), aval, vinfo, ncname,    &
     &                 setparaccess = .true.)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!  Define free-surface.
!
        IF (aout(idfsur,ng)) THEN
          vinfo( 1)=vname(1,idfsur)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idfsur))
          vinfo( 3)=vname(3,idfsur)
          vinfo(14)=vname(4,idfsur)
          vinfo(16)=vname(1,idtime)
#  if !defined WET_DRY && (defined WRITE_WATER && defined MASKING)
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idfsur)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idfsur,ng),r8)
          avg(ng)%pioVar(idfsur)%dkind=pio_fout
          avg(ng)%pioVar(idfsur)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idfsur)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided free-surface.
!
        IF (aout(idfsud,ng)) THEN
          vinfo( 1)=vname(1,idfsud)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idfsud))
          vinfo( 3)=vname(3,idfsud)
          vinfo(14)=vname(4,idfsud)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idfsud)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idfsud,ng),r8)
          avg(ng)%pioVar(idfsud)%dkind=pio_fout
          avg(ng)%pioVar(idfsud)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idfsud)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
!
!  Define 2D momentum in the XI-direction.
!
        IF (aout(idubar,ng)) THEN
          vinfo( 1)=vname(1,idubar)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idubar))
          vinfo( 3)=vname(3,idubar)
          vinfo(14)=vname(4,idubar)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#  endif
          vinfo(21)=vname(6,idubar)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idubar,ng),r8)
          avg(ng)%pioVar(idubar)%dkind=pio_fout
          avg(ng)%pioVar(idubar)%gtype=u2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idubar)%vd,                     &
     &                   pio_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided 2D momentum in the XI-direction.
!
        IF (aout(idu2dd,ng)) THEN
          vinfo( 1)=vname(1,idu2dd)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idu2dd))
          vinfo( 3)=vname(3,idu2dd)
          vinfo(14)=vname(4,idu2dd)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#   endif
          vinfo(21)=vname(6,idu2dd)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idu2dd,ng),r8)
          avg(ng)%pioVar(idu2dd)%dkind=pio_fout
          avg(ng)%pioVar(idu2dd)%gtype=u2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idu2dd)%vd,                     &
     &                   pio_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
!
!  Define 2D momentum in the ETA-direction.
!
        IF (aout(idvbar,ng)) THEN
          vinfo( 1)=vname(1,idvbar)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvbar))
          vinfo( 3)=vname(3,idvbar)
          vinfo(14)=vname(4,idvbar)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#  endif
          vinfo(21)=vname(6,idvbar)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvbar,ng),r8)
          avg(ng)%pioVar(idvbar)%dkind=pio_fout
          avg(ng)%pioVar(idvbar)%gtype=v2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvbar)%vd,                     &
     &                   pio_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided 2D momentum in the ETA-direction.
!
        IF (aout(idv2dd,ng)) THEN
          vinfo( 1)=vname(1,idv2dd)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idv2dd))
          vinfo( 3)=vname(3,idv2dd)
          vinfo(14)=vname(4,idv2dd)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#   endif
          vinfo(21)=vname(6,idv2dd)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idv2dd,ng),r8)
          avg(ng)%pioVar(idv2dd)%dkind=pio_fout
          avg(ng)%pioVar(idv2dd)%gtype=v2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idv2dd)%vd,                     &
     &                   pio_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
!
!  Define 2D Eastward momentum component at RHO-points.
!
        IF (aout(idu2de,ng)) THEN
          vinfo( 1)=vname(1,idu2de)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idu2de))
          vinfo( 3)=vname(3,idu2de)
          vinfo(14)=vname(4,idu2de)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idu2de)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idu2de,ng),r8)
          avg(ng)%pioVar(idu2de)%dkind=pio_fout
          avg(ng)%pioVar(idu2de)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idu2de)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define 2D Northward momentum component at RHO-points.
!
        IF (aout(idv2dn,ng)) THEN
          vinfo( 1)=vname(1,idv2dn)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idv2dn))
          vinfo( 3)=vname(3,idv2dn)
          vinfo(14)=vname(4,idv2dn)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idv2dn)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idv2dn,ng),r8)
          avg(ng)%pioVar(idv2dn)%dkind=pio_fout
          avg(ng)%pioVar(idv2dn)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idv2dn)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  ifdef SOLVE3D
#   ifdef FORWARD_WRITE
!
!  Define time-averaged mass fluxes for 3D momentum coupling.
!
        IF (aout(idufx1,ng)) THEN
          vinfo( 1)=vname(1,idufx1)
          vinfo( 2)=vname(2,idufx1)
          vinfo( 3)=vname(3,idufx1)
          vinfo(14)=vname(4,idufx1)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#    endif
          vinfo(21)=vname(6,idufx1)
          vinfo(22)='coordinates'
          aval(5)=real(u2dvar,r8)
          avg(ng)%pioVar(idufx1)%dkind=pio_fout
          avg(ng)%pioVar(idufx1)%gtype=u2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idufx1)%vd,                     &
     &                   pio_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idufx2,ng)) THEN
          vinfo( 1)=vname(1,idufx2)
          vinfo( 2)=vname(2,idufx2)
          vinfo( 3)=vname(3,idufx2)
          vinfo(14)=vname(4,idufx2)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#    endif
          vinfo(21)=vname(6,idufx2)
          vinfo(22)='coordinates'
          aval(5)=real(u2dvar,r8)
          avg(ng)%pioVar(idufx2)%dkind=pio_fout
          avg(ng)%pioVar(idufx2)%gtype=u2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idufx2)%vd,                     &
     &                   pio_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idvfx1,ng)) THEN
          vinfo( 1)=vname(1,idvfx1)
          vinfo( 2)=vname(2,idvfx1)
          vinfo( 3)=vname(3,idvfx1)
          vinfo(14)=vname(4,idvfx1)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#    endif
          vinfo(21)=vname(6,idvfx1)
          vinfo(22)='coordinates'
          aval(5)=real(v2dvar,r8)
          avg(ng)%pioVar(idvfx1)%dkind=pio_fout
          avg(ng)%pioVar(idvfx1)%gtype=v2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvfx1)%vd,                     &
     &                   pio_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idvfx2,ng)) THEN
          vinfo( 1)=vname(1,idvfx2)
          vinfo( 2)=vname(2,idvfx2)
          vinfo( 3)=vname(3,idvfx2)
          vinfo(14)=vname(4,idvfx2)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#    endif
          vinfo(21)=vname(6,idvfx2)
          vinfo(22)='coordinates'
          aval(5)=real(v2dvar,r8)
          avg(ng)%pioVar(idvfx2)%dkind=pio_fout
          avg(ng)%pioVar(idvfx2)%gtype=v2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvfx2)%vd,                     &
     &                   pio_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
!
!  Define 3D momentum component in the XI-direction.
!
        IF (aout(iduvel,ng)) THEN
          vinfo( 1)=vname(1,iduvel)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,iduvel))
          vinfo( 3)=vname(3,iduvel)
          vinfo(14)=vname(4,iduvel)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#   endif
          vinfo(21)=vname(6,iduvel)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,iduvel,ng),r8)
          avg(ng)%pioVar(iduvel)%dkind=pio_fout
          avg(ng)%pioVar(iduvel)%gtype=u3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(iduvel)%vd,                     &
     &                   pio_fout, nvd4, u3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#   if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided 3D momentum component in the XI-direction.
!
        IF (aout(idu3dd,ng)) THEN
          vinfo( 1)=vname(1,idu3dd)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idu3dd))
          vinfo( 3)=vname(3,idu3dd)
          vinfo(14)=vname(4,idu3dd)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#    endif
          vinfo(21)=vname(6,idu3dd)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idu3dd,ng),r8)
          avg(ng)%pioVar(idu3dd)%dkind=pio_fout
          avg(ng)%pioVar(idu3dd)%gtype=u3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idu3dd)%vd,                     &
     &                   pio_fout, nvd4, u3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
!
!  Define 3D momentum component in the ETA-direction.
!
        IF (aout(idvvel,ng)) THEN
          vinfo( 1)=vname(1,idvvel)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvvel))
          vinfo( 3)=vname(3,idvvel)
          vinfo(14)=vname(4,idvvel)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#   endif
          vinfo(21)=vname(6,idvvel)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvvel,ng),r8)
          avg(ng)%pioVar(idvvel)%dkind=pio_fout
          avg(ng)%pioVar(idvvel)%gtype=v3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvvel)%vd,                     &
     &                   pio_fout, nvd4, v3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#   if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided 3D momentum component in the ETA-direction.
!
        IF (aout(idv3dd,ng)) THEN
          vinfo( 1)=vname(1,idv3dd)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idv3dd))
          vinfo( 3)=vname(3,idv3dd)
          vinfo(14)=vname(4,idv3dd)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#    endif
          vinfo(21)=vname(6,idv3dd)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idv3dd,ng),r8)
          avg(ng)%pioVar(idv3dd)%dkind=pio_fout
          avg(ng)%pioVar(idv3dd)%gtype=v3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idv3dd)%vd,                     &
     &                   pio_fout, nvd4, v3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
!
!  Define 3D Eastward momentum component at RHO-points.
!
        IF (aout(idu3de,ng)) THEN
          vinfo( 1)=vname(1,idu3de)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idu3de))
          vinfo( 3)=vname(3,idu3de)
          vinfo(14)=vname(4,idu3de)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idu3de)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idu3de,ng),r8)
          avg(ng)%pioVar(idu3de)%dkind=pio_fout
          avg(ng)%pioVar(idu3de)%gtype=r3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idu3de)%vd,                     &
     &                   pio_fout, nvd4, t3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define 3D Northward momentum component at RHO-points.
!
        IF (aout(idv3dn,ng)) THEN
          vinfo( 1)=vname(1,idv3dn)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idv3dn))
          vinfo( 3)=vname(3,idv3dn)
          vinfo(14)=vname(4,idv3dn)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idv3dn)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idv3dn,ng),r8)
          avg(ng)%pioVar(idv3dn)%dkind=pio_fout
          avg(ng)%pioVar(idv3dn)%gtype=r3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idv3dn)%vd,                     &
     &                   pio_fout, nvd4, t3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define S-coordinate vertical "omega" momentum component.
!
        IF (aout(idovel,ng)) THEN
          vinfo( 1)=vname(1,idovel)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idovel))
          vinfo( 3)=vname(3,idovel)
          vinfo(14)=vname(4,idovel)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idovel)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idovel,ng),r8)
          avg(ng)%pioVar(idovel)%dkind=pio_fout
          avg(ng)%pioVar(idovel)%gtype=w3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idovel)%vd,                     &
     &                   pio_fout, nvd4, w3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define "true" vertical momentum component.
!
        IF (aout(idwvel,ng)) THEN
          vinfo( 1)=vname(1,idwvel)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idwvel))
          vinfo( 3)=vname(3,idwvel)
          vinfo(14)=vname(4,idwvel)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idwvel)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idwvel,ng),r8)
          avg(ng)%pioVar(idwvel)%dkind=pio_fout
          avg(ng)%pioVar(idwvel)%gtype=w3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idwvel)%vd,                     &
     &                   pio_fout, nvd4, w3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define tracer type variables.
!
        DO itrc=1,nt(ng)
          IF (aout(idtvar(itrc),ng)) THEN
            vinfo( 1)=vname(1,idtvar(itrc))
            WRITE (vinfo( 2),'(a,1x,a)') prefix,                        &
     &                                   trim(vname(2,idtvar(itrc)))
            vinfo( 3)=vname(3,idtvar(itrc))
            vinfo(14)=vname(4,idtvar(itrc))
            vinfo(16)=vname(1,idtime)
#   ifdef SEDIMENT
            DO i=1,nst
              IF (itrc.eq.idsed(i)) THEN
                WRITE (vinfo(19),40) 1000.0_r8*sd50(i,ng)
              END IF
            END DO
#   endif
#   if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_rho'
#   endif
            vinfo(21)=vname(6,idtvar(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(r3dvar,r8)
            avg(ng)%pioTrc(itrc)%dkind=pio_fout
            avg(ng)%pioTrc(itrc)%gtype=r3dvar
!
            status=def_var(ng, model, avg(ng)%pioFile,                  &
     &                     avg(ng)%pioTrc(itrc)%vd,                     &
     &                     pio_fout, nvd4, t3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
 
#   if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
!
!  Define detided temperature and salinity.
!
        DO itrc=1,nat
          IF (aout(idtrcd(itrc),ng)) THEN
            ifield=idtrcd(itrc)
            vinfo( 1)=vname(1,ifield)
            WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,ifield))
            vinfo( 3)=vname(3,ifield)
            vinfo(14)=vname(4,ifield)
            vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_rho'
#    endif
            vinfo(21)=vname(6,ifield)
            vinfo(22)='coordinates'
            aval(5)=real(r3dvar,r8)
            avg(ng)%pioVar(ifield)%dkind=pio_fout
            avg(ng)%pioVar(ifield)%gtype=r3dvar
!
            status=def_var(ng, model, avg(ng)%pioFile,                  &
     &                     avg(ng)%pioVar(ifield)%vd,                   &
     &                     pio_fout, nvd4, t3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
#   endif
!
!  Define density anomaly.
!
        IF (aout(iddano,ng)) THEN
          vinfo( 1)=vname(1,iddano)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,iddano))
          vinfo( 3)=vname(3,iddano)
          vinfo(14)=vname(4,iddano)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,iddano)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,iddano,ng),r8)
          avg(ng)%pioVar(iddano)%dkind=pio_fout
          avg(ng)%pioVar(iddano)%gtype=r3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(iddano)%vd,                     &
     &                   pio_fout, nvd4, t3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#   ifdef LMD_SKPP
!
!  Define depth of surface boundary layer.
!
        IF (aout(idhsbl,ng)) THEN
          vinfo( 1)=vname(1,idhsbl)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idhsbl))
          vinfo( 3)=vname(3,idhsbl)
          vinfo(14)=vname(4,idhsbl)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idhsbl)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idhsbl,ng),r8)
          avg(ng)%pioVar(idhsbl)%dkind=pio_fout
          avg(ng)%pioVar(idhsbl)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idhsbl)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
#   ifdef LMD_BKPP
!
!  Define depth of bottom boundary layer.
!
        IF (aout(idhbbl,ng)) THEN
          vinfo( 1)=vname(1,idhbbl)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idhbbl))
          vinfo( 3)=vname(3,idhbbl)
          vinfo(14)=vname(4,idhbbl)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idhbbl)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idhbbl,ng),r8)
          avg(ng)%pioVar(idhbbl)%dkind=pio_fout
          avg(ng)%pioVar(idhbbl)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idhbbl)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
#  endif
!
!  Define 2D potential vorticity.
!
        IF (aout(id2dpv,ng)) THEN
          vinfo( 1)=vname(1,id2dpv)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,id2dpv))
          vinfo( 3)=vname(3,id2dpv)
          vinfo(14)=vname(4,id2dpv)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_psi'
#  endif
          vinfo(21)=vname(6,id2dpv)
          vinfo(22)='coordinates'
          aval(5)=real(p2dvar,r8)
          avg(ng)%pioVar(id2dpv)%dkind=pio_fout
          avg(ng)%pioVar(id2dpv)%gtype=p2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(id2dpv)%vd,                     &
     &                   pio_fout, nvd3, p2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define 2D relative vorticity.
!
        IF (aout(id2drv,ng)) THEN
          vinfo( 1)=vname(1,id2drv)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,id2drv))
          vinfo( 3)=vname(3,id2drv)
          vinfo(14)=vname(4,id2drv)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_psi'
#  endif
          vinfo(21)=vname(6,id2drv)
          vinfo(22)='coordinates'
          aval(5)=real(p2dvar,r8)
          avg(ng)%pioVar(id2drv)%dkind=pio_fout
          avg(ng)%pioVar(id2drv)%gtype=p2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(id2drv)%vd,                     &
     &                   pio_fout, nvd3, p2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  ifdef SOLVE3D
!
!  Define 3D potential vorticity.
!
        IF (aout(id3dpv,ng)) THEN
          vinfo( 1)=vname(1,id3dpv)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,id3dpv))
          vinfo( 3)=vname(3,id3dpv)
          vinfo(14)=vname(4,id3dpv)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_psi'
#   endif
          vinfo(21)=vname(6,id3dpv)
          vinfo(22)='coordinates'
          aval(5)=real(p3dvar,r8)
          avg(ng)%pioVar(id3dpv)%dkind=pio_fout
          avg(ng)%pioVar(id3dpv)%gtype=p3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(id3dpv)%vd,                     &
     &                   pio_fout, nvd4, p3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define 3D relative vorticity.
!
        IF (aout(id3drv,ng)) THEN
          vinfo( 1)=vname(1,id3drv)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,id3drv))
          vinfo( 3)=vname(3,id3drv)
          vinfo(14)=vname(4,id3drv)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_psi'
#   endif
          vinfo(21)=vname(6,id3drv)
          vinfo(22)='coordinates'
          aval(5)=real(p3dvar,r8)
          avg(ng)%pioVar(id3drv)%dkind=pio_fout
          avg(ng)%pioVar(id3drv)%gtype=p3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(id3drv)%vd,                     &
     &                   pio_fout, nvd4, p3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#  endif
!
!  Define quadratic <zeta*zeta> term.
!
        IF (aout(idzzav,ng)) THEN
          vinfo( 1)=vname(1,idzzav)
          vinfo( 2)=trim(vname(2,idzzav))
          vinfo( 3)=vname(3,idzzav)
          vinfo(14)=vname(4,idzzav)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idzzav)
          vinfo(22)='coordinates'
          aval(5)=real(r2dvar,r8)
          avg(ng)%pioVar(idzzav)%dkind=pio_fout
          avg(ng)%pioVar(idzzav)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idzzav)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <ubar*ubar> term.
!
        IF (aout(idu2av,ng)) THEN
          vinfo( 1)=vname(1,idu2av)
          vinfo( 2)=trim(vname(2,idu2av))
          vinfo( 3)=vname(3,idu2av)
          vinfo(14)=vname(4,idu2av)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#  endif
          vinfo(21)=vname(6,idu2av)
          vinfo(22)='coordinates'
          aval(5)=real(u2dvar,r8)
          avg(ng)%pioVar(idu2av)%dkind=pio_fout
          avg(ng)%pioVar(idu2av)%gtype=u2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idu2av)%vd,                     &
     &                   pio_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <vbar*vbar> term.
!
        IF (aout(idv2av,ng)) THEN
          vinfo( 1)=vname(1,idv2av)
          vinfo( 2)=trim(vname(2,idv2av))
          vinfo( 3)=vname(3,idv2av)
          vinfo(14)=vname(4,idv2av)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#  endif
          vinfo(21)=vname(6,idv2av)
          vinfo(22)='coordinates'
          aval(5)=real(v2dvar,r8)
          avg(ng)%pioVar(idv2av)%dkind=pio_fout
          avg(ng)%pioVar(idv2av)%gtype=v2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idv2av)%vd,                     &
     &                   pio_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  ifdef SOLVE3D
!
!  Define u-volume flux.
!
        IF (aout(idhuav,ng)) THEN
          vinfo( 1)=vname(1,idhuav)
          vinfo( 2)=trim(vname(2,idhuav))
          vinfo( 3)=vname(3,idhuav)
          vinfo(14)=vname(4,idhuav)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#   endif
          vinfo(21)=vname(6,idhuav)
          vinfo(22)='coordinates'
          aval(5)=real(u3dvar,r8)
          avg(ng)%pioVar(idhuav)%dkind=pio_fout
          avg(ng)%pioVar(idhuav)%gtype=u3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idhuav)%vd,                     &
     &                   pio_fout, nvd4, u3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define v-volume flux.
!
        IF (aout(idhvav,ng)) THEN
          vinfo( 1)=vname(1,idhvav)
          vinfo( 2)=trim(vname(2,idhvav))
          vinfo( 3)=vname(3,idhvav)
          vinfo(14)=vname(4,idhvav)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#   endif
          vinfo(21)=vname(6,idhvav)
          vinfo(22)='coordinates'
          aval(5)=real(v3dvar,r8)
          avg(ng)%pioVar(idhvav)%dkind=pio_fout
          avg(ng)%pioVar(idhvav)%gtype=v3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idhvav)%vd,                     &
     &                   pio_fout, nvd4, v3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <u*u> term.
!
        IF (aout(iduuav,ng)) THEN
          vinfo( 1)=vname(1,iduuav)
          vinfo( 2)=trim(vname(2,iduuav))
          vinfo( 3)=vname(3,iduuav)
          vinfo(14)=vname(4,iduuav)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#   endif
          vinfo(21)=vname(6,iduuav)
          vinfo(22)='coordinates'
          aval(5)=real(u3dvar,r8)
          avg(ng)%pioVar(iduuav)%dkind=pio_fout
          avg(ng)%pioVar(iduuav)%gtype=u3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(iduuav)%vd,                     &
     &                   pio_fout, nvd4, u3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <u*v> term.
!
        IF (aout(iduvav,ng)) THEN
          vinfo( 1)=vname(1,iduvav)
          vinfo( 2)=trim(vname(2,iduvav))
          vinfo( 3)=vname(3,iduvav)
          vinfo(14)=vname(4,iduvav)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,iduvav)
          vinfo(22)='coordinates'
          aval(5)=real(r3dvar,r8)
          avg(ng)%pioVar(iduvav)%dkind=pio_fout
          avg(ng)%pioVar(iduvav)%gtype=r3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(iduvav)%vd,                     &
     &                   pio_fout, nvd4, t3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <v*v> term.
!
        IF (aout(idvvav,ng)) THEN
          vinfo( 1)=vname(1,idvvav)
          vinfo( 2)=trim(vname(2,idvvav))
          vinfo( 3)=vname(3,idvvav)
          vinfo(14)=vname(4,idvvav)
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#   endif
          vinfo(21)=vname(6,idvvav)
          vinfo(22)='coordinates'
          aval(5)=real(v3dvar,r8)
          avg(ng)%pioVar(idvvav)%dkind=pio_fout
          avg(ng)%pioVar(idvvav)%gtype=v3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvvav)%vd,                     &
     &                   pio_fout, nvd4, v3dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define quadratic <t*t> terms.
!
        DO itrc=1,nt(ng)
          IF (aout(idttav(itrc),ng)) THEN
            vinfo( 1)=vname(1,idttav(itrc))
            vinfo( 2)=trim(vname(2,idttav(itrc)))
            vinfo( 3)=vname(3,idttav(itrc))
            vinfo(14)=vname(4,idttav(itrc))
            vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_rho'
#   endif
            vinfo(21)=vname(6,idttav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(r3dvar,r8)
            avg(ng)%pioVar(idttav(itrc))%dkind=pio_fout
            avg(ng)%pioVar(idttav(itrc))%gtype=r3dvar
!
            status=def_var(ng, model, avg(ng)%pioFile,                  &
     &                     avg(ng)%pioVar(idttav(itrc))%vd,             &
     &                     pio_fout, nvd4, t3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
!
!  Define active tracers volume fluxes.
!
        DO itrc=1,nt(ng)
          IF (aout(ihutav(itrc),ng)) THEN
            vinfo( 1)=vname(1,ihutav(itrc))
            vinfo( 2)=trim(vname(2,ihutav(itrc)))
            vinfo( 3)=vname(3,ihutav(itrc))
            vinfo(14)=vname(4,ihutav(itrc))
            vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_u'
#   endif
            vinfo(21)=vname(6,ihutav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(u3dvar,r8)
            avg(ng)%pioVar(ihutav(itrc))%dkind=pio_fout
            avg(ng)%pioVar(ihutav(itrc))%gtype=u3dvar
!
            status=def_var(ng, model, avg(ng)%pioFile,                  &
     &                     avg(ng)%pioVar(ihutav(itrc))%vd,             &
     &                     pio_fout, nvd4, u3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
!
          IF (aout(ihvtav(itrc),ng)) THEN
            vinfo( 1)=vname(1,ihvtav(itrc))
            vinfo( 2)=trim(vname(2,ihvtav(itrc)))
            vinfo( 3)=vname(3,ihvtav(itrc))
            vinfo(14)=vname(4,ihvtav(itrc))
            vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_v'
#   endif
            vinfo(21)=vname(6,ihvtav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(v3dvar,r8)
            avg(ng)%pioVar(ihvtav(itrc))%dkind=pio_fout
            avg(ng)%pioVar(ihvtav(itrc))%gtype=v3dvar
!
            status=def_var(ng, model, avg(ng)%pioFile,                  &
     &                     avg(ng)%pioVar(ihvtav(itrc))%vd,             &
     &                     pio_fout, nvd4, v3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
!
!  Define quadratic <u*t> and <v*t> terms.
!
        DO itrc=1,nt(ng)
          IF (aout(idutav(itrc),ng)) THEN
            vinfo( 1)=vname(1,idutav(itrc))
            vinfo( 2)=trim(vname(2,idutav(itrc)))
            vinfo( 3)=vname(3,idutav(itrc))
            vinfo(14)=vname(4,idutav(itrc))
            vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_u'
#   endif
            vinfo(21)=vname(6,idutav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(u3dvar,r8)
            avg(ng)%pioVar(idutav(itrc))%dkind=pio_fout
            avg(ng)%pioVar(idutav(itrc))%gtype=u3dvar
!
            status=def_var(ng, model, avg(ng)%pioFile,                  &
     &                     avg(ng)%pioVar(idutav(itrc))%vd,             &
     &                     pio_fout, nvd4, u3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
!
          IF (aout(idvtav(itrc),ng)) THEN
            vinfo( 1)=vname(1,idvtav(itrc))
            vinfo( 2)=trim(vname(2,idvtav(itrc)))
            vinfo( 3)=vname(3,idvtav(itrc))
            vinfo(14)=vname(4,idvtav(itrc))
            vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
            vinfo(20)='mask_v'
#   endif
            vinfo(21)=vname(6,idvtav(itrc))
            vinfo(22)='coordinates'
            aval(5)=real(v3dvar,r8)
            avg(ng)%pioVar(idvtav(itrc))%dkind=pio_fout
            avg(ng)%pioVar(idvtav(itrc))%gtype=v3dvar
!
            status=def_var(ng, model, avg(ng)%pioFile,                  &
     &                     avg(ng)%pioVar(idvtav(itrc))%vd,             &
     &                     pio_fout, nvd4, v3dgrd, aval, vinfo, ncname)
            IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
          END IF
        END DO
 
#   if defined LMD_MIXING || defined MY25_MIXING || defined GLS_MIXING
!
!  Define vertical viscosity coefficient.
!
        IF (aout(idvvis,ng)) THEN
          vinfo( 1)=vname(1,idvvis)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvvis))
          vinfo( 3)=vname(3,idvvis)
          vinfo(14)=vname(4,idvvis)
          vinfo(16)=vname(1,idtime)
          vinfo(21)=vname(6,idvvis)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvvis,ng),r8)
          avg(ng)%pioVar(idvvis)%dkind=pio_fout
          avg(ng)%pioVar(idvvis)%gtype=w3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvvis)%vd,                     &
     &                   pio_fout, nvd4, w3dgrd, aval, vinfo, ncname,   &
     &                   setfillval = .false.)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define vertical diffusion coefficient for potential temperature.
!
        IF (aout(idtdif,ng)) THEN
          vinfo( 1)=vname(1,idtdif)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idtdif))
          vinfo( 3)=vname(3,idtdif)
          vinfo(14)=vname(4,idtdif)
          vinfo(16)=vname(1,idtime)
          vinfo(21)=vname(6,idtdif)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idtdif,ng),r8)
          avg(ng)%pioVar(idtdif)%dkind=pio_fout
          avg(ng)%pioVar(idtdif)%gtype=w3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idtdif)%vd,                     &
     &                   pio_fout, nvd4, w3dgrd, aval, vinfo, ncname,   &
     &                   setfillval = .false.)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#    ifdef SALINITY
!
!  Define vertical diffusion coefficient for salinity.
!
        IF (aout(idsdif,ng)) THEN
          vinfo( 1)=vname(1,idsdif)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idsdif))
          vinfo( 3)=vname(3,idsdif)
          vinfo(14)=vname(4,idsdif)
          vinfo(16)=vname(1,idtime)
          vinfo(21)=vname(6,idsdif)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idsdif,ng),r8)
          avg(ng)%pioVar(idsdif)%dkind=pio_fout
          avg(ng)%pioVar(idsdif)%gtype=w3dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idsdif)%vd,                     &
     &                   pio_fout, nvd4, w3dgrd, aval, vinfo, ncname,   &
     &                   setfillval = .false.)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#    endif
#   endif
#   if defined BULK_FLUXES || defined ECOSIM || defined ATM_PRESS
!
!  Define surface air pressure.
!
        IF (aout(idpair,ng)) THEN
          vinfo( 1)=vname(1,idpair)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idpair))
          vinfo( 3)=vname(3,idpair)
          vinfo(14)=vname(4,idpair)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idpair)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idpair,ng),r8)
          avg(ng)%pioVar(idpair)%dkind=pio_fout
          avg(ng)%pioVar(idpair)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idpair)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
#   if defined BULK_FLUXES
!
!  Define surface air temperature.
!
        IF (aout(idtair,ng)) THEN
          vinfo( 1)=vname(1,idtair)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idtair))
          vinfo( 3)=vname(3,idtair)
          vinfo(14)=vname(4,idtair)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idtair)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idtair,ng),r8)
          avg(ng)%pioVar(idtair)%dkind=pio_fout
          avg(ng)%pioVar(idtair)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idtair)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
#   if defined BULK_FLUXES || defined ECOSIM
!
!  Define surface winds.
!
        IF (aout(iduair,ng)) THEN
          vinfo( 1)=vname(1,iduair)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,iduair))
          vinfo( 3)=vname(3,iduair)
          vinfo(14)=vname(4,iduair)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,iduair)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,iduair,ng),r8)
          avg(ng)%pioVar(iduair)%dkind=pio_fout
          avg(ng)%pioVar(iduair)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(iduair)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idvair,ng)) THEN
          vinfo( 1)=vname(1,idvair)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvair))
          vinfo( 3)=vname(3,idvair)
          vinfo(14)=vname(4,idvair)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idvair)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvair,ng),r8)
          avg(ng)%pioVar(idvair)%dkind=pio_fout
          avg(ng)%pioVar(idvair)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvair)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define Eastward/Northward surface wind at RHO-points.
!
        IF (aout(iduaie,ng)) THEN
          vinfo( 1)=vname(1,iduaie)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,iduaie))
          vinfo( 3)=vname(3,iduaie)
          vinfo(14)=vname(4,iduaie)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,iduaie)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,iduaie,ng),r8)
          avg(ng)%pioVar(iduaie)%dkind=pio_fout
          avg(ng)%pioVar(iduaie)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(iduaie)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
        IF (aout(idvain,ng)) THEN
          vinfo( 1)=vname(1,idvain)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvain))
          vinfo( 3)=vname(3,idvain)
          vinfo(14)=vname(4,idvain)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#  endif
          vinfo(21)=vname(6,idvain)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvain,ng),r8)
          avg(ng)%pioVar(idvain)%dkind=pio_fout
          avg(ng)%pioVar(idvain)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvain)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
!
!  Define surface net heat flux.
!
        IF (aout(idtsur(itemp),ng)) THEN
          vinfo( 1)=vname(1,idtsur(itemp))
          WRITE (vinfo( 2),'(a,1x,a)') prefix,                          &
     &                                 trim(vname(2,idtsur(itemp)))
          vinfo( 3)=vname(3,idtsur(itemp))
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idtsur(itemp))
          vinfo(16)=vname(1,idtime)
#   if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#   endif
          vinfo(21)=vname(6,idtsur(itemp))
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idtsur(itemp),ng),r8)
          avg(ng)%pioVar(idtsur(itemp))%dkind=pio_fout
          avg(ng)%pioVar(idtsur(itemp))%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idtsur(itemp))%vd,              &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#   ifdef SALINITY
!
!  Define surface net salt flux.
!
        IF (aout(idtsur(isalt),ng)) THEN
          vinfo( 1)=vname(1,idtsur(isalt))
          WRITE (vinfo( 2),'(a,1x,a)') prefix,                          &
     &                                 trim(vname(2,idtsur(isalt)))
          vinfo( 3)=vname(3,idtsur(isalt))
          vinfo(11)='upward flux, freshening (net precipitation)'
          vinfo(12)='downward flux, salting (net evaporation)'
          vinfo(14)=vname(4,idtsur(isalt))
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idtsur(isalt))
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idtsur(isalt),ng),r8)
          avg(ng)%pioVar(idtsur(isalt))%dkind=pio_fout
          avg(ng)%pioVar(idtsur(isalt))%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idtsur(isalt))%vd,              &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
 
#   if defined BULK_FLUXES || defined FRC_COUPLING
!
!  Define latent heat flux.
!
        IF (aout(idlhea,ng)) THEN
          vinfo( 1)=vname(1,idlhea)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idlhea))
          vinfo( 3)=vname(3,idlhea)
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idlhea)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idlhea)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idlhea,ng),r8)
          avg(ng)%pioVar(idlhea)%dkind=pio_fout
          avg(ng)%pioVar(idlhea)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idlhea)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define sensible heat flux.
!
        IF (aout(idshea,ng)) THEN
          vinfo( 1)=vname(1,idshea)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idshea))
          vinfo( 3)=vname(3,idshea)
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idshea)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idshea)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idshea,ng),r8)
          avg(ng)%pioVar(idshea)%dkind=pio_fout
          avg(ng)%pioVar(idshea)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idshea)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define net longwave radiation flux.
!
        IF (aout(idlrad,ng)) THEN
          vinfo( 1)=vname(1,idlrad)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idlrad))
          vinfo( 3)=vname(3,idlrad)
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idlrad)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idlrad)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idlrad,ng),r8)
          avg(ng)%pioVar(idlrad)%dkind=pio_fout
          avg(ng)%pioVar(idlrad)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idlrad)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
 
#   ifdef BULK_FLUXES
#    ifdef EMINUSP
!
!  Define evaporation rate.
!
        IF (aout(idevap,ng)) THEN
          vinfo( 1)=vname(1,idevap)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idevap))
          vinfo( 3)=vname(3,idevap)
          vinfo(11)='downward flux, freshening (condensation)'
          vinfo(12)='upward flux, salting (evaporation)'
          vinfo(14)=vname(4,idevap)
          vinfo(16)=vname(1,idtime)
#     if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#     endif
          vinfo(21)=vname(6,idevap)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idevap,ng),r8)
          avg(ng)%pioVar(idevap)%dkind=pio_fout
          avg(ng)%pioVar(idevap)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idevap)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define precipitation rate.
!
        IF (aout(idrain,ng)) THEN
          vinfo( 1)=vname(1,idrain)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idrain))
          vinfo( 3)=vname(3,idrain)
          vinfo(12)='downward flux, freshening (precipitation)'
          vinfo(14)=vname(4,idrain)
          vinfo(16)=vname(1,idtime)
#     if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#     endif
          vinfo(21)=vname(6,idrain)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idrain,ng),r8)
          avg(ng)%pioVar(idrain)%dkind=pio_fout
          avg(ng)%pioVar(idrain)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idrain)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#    endif
#   endif
#   ifdef SHORTWAVE
!
!  Define net shortwave radiation flux.
!
        IF (aout(idsrad,ng)) THEN
          vinfo( 1)=vname(1,idsrad)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idsrad))
          vinfo( 3)=vname(3,idsrad)
          vinfo(11)='upward flux, cooling'
          vinfo(12)='downward flux, heating'
          vinfo(14)=vname(4,idsrad)
          vinfo(16)=vname(1,idtime)
#    if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_rho'
#    endif
          vinfo(21)=vname(6,idsrad)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idsrad,ng),r8)
          avg(ng)%pioVar(idsrad)%dkind=pio_fout
          avg(ng)%pioVar(idsrad)%gtype=r2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idsrad)%vd,                     &
     &                   pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
#   endif
#  endif
!
!  Define surface u-momentum stress.
!
        IF (aout(idusms,ng)) THEN
          vinfo( 1)=vname(1,idusms)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idusms))
          vinfo( 3)=vname(3,idusms)
          vinfo(14)=vname(4,idusms)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#  endif
          vinfo(21)=vname(6,idusms)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idusms,ng),r8)
          avg(ng)%pioVar(idusms)%dkind=pio_fout
          avg(ng)%pioVar(idusms)%gtype=u2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idusms)%vd,                     &
     &                   pio_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define surface v-momentum stress.
!
        IF (aout(idvsms,ng)) THEN
          vinfo( 1)=vname(1,idvsms)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvsms))
          vinfo( 3)=vname(3,idvsms)
          vinfo(14)=vname(4,idvsms)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#  endif
          vinfo(21)=vname(6,idvsms)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvsms,ng),r8)
          avg(ng)%pioVar(idvsms)%dkind=pio_fout
          avg(ng)%pioVar(idvsms)%gtype=v2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvsms)%vd,                     &
     &                   pio_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define bottom u-momentum stress.
!
        IF (aout(idubms,ng)) THEN
          vinfo( 1)=vname(1,idubms)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idubms))
          vinfo( 3)=vname(3,idubms)
          vinfo(14)=vname(4,idubms)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_u'
#  endif
          vinfo(21)=vname(6,idubms)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idubms,ng),r8)
          avg(ng)%pioVar(idubms)%dkind=pio_fout
          avg(ng)%pioVar(idubms)%gtype=u2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idubms)%vd,                     &
     &                   pio_fout, nvd3, u2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
!
!  Define bottom v-momentum stress.
!
        IF (aout(idvbms,ng)) THEN
          vinfo( 1)=vname(1,idvbms)
          WRITE (vinfo( 2),'(a,1x,a)') prefix, trim(vname(2,idvbms))
          vinfo( 3)=vname(3,idvbms)
          vinfo(14)=vname(4,idvbms)
          vinfo(16)=vname(1,idtime)
#  if defined WRITE_WATER && defined MASKING
          vinfo(20)='mask_v'
#  endif
          vinfo(21)=vname(6,idvbms)
          vinfo(22)='coordinates'
          aval(5)=real(iinfo(1,idvbms,ng),r8)
          avg(ng)%pioVar(idvbms)%dkind=pio_fout
          avg(ng)%pioVar(idvbms)%gtype=v2dvar
!
          status=def_var(ng, model, avg(ng)%pioFile,                    &
     &                   avg(ng)%pioVar(idvbms)%vd,                     &
     &                   pio_fout, nvd3, v2dgrd, aval, vinfo, ncname)
          IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
        END IF
 
#  if (defined BBL_MODEL || defined WAVES_OUTPUT) && defined SOLVE3D
!
!-----------------------------------------------------------------------
!  Define the bottom boundary layer model or waves variables.
!-----------------------------------------------------------------------
!
        CALL bbl_def_pio (ng, model, ldef, aout, avg,                   &
     &                    t2dgrd, u2dgrd, v2dgrd)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
 
#  if defined ICE_MODEL && defined SOLVE3D
!
!-----------------------------------------------------------------------
!  Define the sea-ice model variables.
!-----------------------------------------------------------------------
!
        CALL ice_def_pio (ng, model, ldef, aout, avg,                   &
     &                    t2dgrd, u2dgrd, v2dgrd)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
 
#  if defined SEDIMENT && defined SOLVE3D
!
!-----------------------------------------------------------------------
!  Define the sediment model variables.
!-----------------------------------------------------------------------
!
        CALL sediment_def_pio (ng, model, ldef, aout, avg,              &
     &                         t2dgrd, u2dgrd, v2dgrd)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
 
#  if defined WEC_VF && defined SOLVE3D
!
!-----------------------------------------------------------------------
!  Define the Waves Effect on Currents variables.
!-----------------------------------------------------------------------
!
        CALL wec_def_pio (ng, model, ldef, aout, avg,                   &
     &                    t2dgrd, u2dgrd, v2dgrd,                       &
     &                    t3dgrd, u3dgrd, v3dgrd, w3dgrd)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
!
!-----------------------------------------------------------------------
!  Leave definition mode.
!-----------------------------------------------------------------------
!
        CALL pio_netcdf_enddef (ng, model, ncname, avg(ng)%pioFile)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!-----------------------------------------------------------------------
!  Write out time-recordless, information variables.
!-----------------------------------------------------------------------
!
        CALL wrt_info (ng, model, avg(ng)%pioFile, ncname)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
      END IF define
!
!=======================================================================
!  Open an existing averages file, check its contents, and prepare
!  for appending data.
!=======================================================================
!
      query : IF (.not.ldef) THEN
        ncname=avg(ng)%name
!
!  Open averages file for read/write.
!
        CALL pio_netcdf_open (ng, model, ncname, 1, avg(ng)%pioFile)
        IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
          WRITE (stdout,50) trim(ncname)
          RETURN
        END IF
!
!  Inquire about the dimensions and check for consistency.
!
        CALL pio_netcdf_check_dim (ng, model, ncname,                   &
     &                             piofile = avg(ng)%pioFile)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!  Inquire about the variables.
!
        CALL pio_netcdf_inq_var (ng, model, ncname,                     &
     &                           piofile = avg(ng)%pioFile)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
!
!  Initialize logical switches.
!
        DO i=1,nv
          got_var(i)=.false.
        END DO
!
!  Scan variable list from input NetCDF and activate switches for
!  average variables. Get variable IDs.
!
        DO i=1,n_var
          IF (trim(var_name(i)).eq.trim(vname(1,idtime))) THEN
            got_var(idtime)=.true.
            avg(ng)%pioVar(idtime)%vd=var_desc(i)
            avg(ng)%pioVar(idtime)%dkind=pio_tout
            avg(ng)%pioVar(idtime)%gtype=0
#  ifdef WET_DRY
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idrwet))) THEN
            got_var(idrwet)=.true.
            avg(ng)%pioVar(idrwet)%vd=var_desc(i)
            avg(ng)%pioVar(idrwet)%dkind=pio_fout
            avg(ng)%pioVar(idrwet)%gtype=r2dvar
#  endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idfsur))) THEN
            got_var(idfsur)=.true.
            avg(ng)%pioVar(idfsur)%vd=var_desc(i)
            avg(ng)%pioVar(idfsur)%dkind=pio_fout
            avg(ng)%pioVar(idfsur)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubar))) THEN
            got_var(idubar)=.true.
            avg(ng)%pioVar(idubar)%vd=var_desc(i)
            avg(ng)%pioVar(idubar)%dkind=pio_fout
            avg(ng)%pioVar(idubar)%gtype=u2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbar))) THEN
            got_var(idvbar)=.true.
            avg(ng)%pioVar(idvbar)%vd=var_desc(i)
            avg(ng)%pioVar(idvbar)%dkind=pio_fout
            avg(ng)%pioVar(idvbar)%gtype=v2dvar
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idfsud))) THEN
            got_var(idfsud)=.true.
            avg(ng)%pioVar(idfsud)%vd=var_desc(i)
            avg(ng)%pioVar(idfsud)%dkind=pio_fout
            avg(ng)%pioVar(idfsud)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu2dd))) THEN
            got_var(idu2dd)=.true.
            avg(ng)%pioVar(idu2dd)%vd=var_desc(i)
            avg(ng)%pioVar(idu2dd)%dkind=pio_fout
            avg(ng)%pioVar(idu2dd)%gtype=u2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv2dd))) THEN
            got_var(idv2dd)=.true.
            avg(ng)%pioVar(idv2dd)%vd=var_desc(i)
            avg(ng)%pioVar(idv2dd)%dkind=pio_fout
            avg(ng)%pioVar(idv2dd)%gtype=v2dvar
#  endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu2de))) THEN
            got_var(idu2de)=.true.
            avg(ng)%pioVar(idu2de)%vd=var_desc(i)
            avg(ng)%pioVar(idu2de)%dkind=pio_fout
            avg(ng)%pioVar(idu2de)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv2dn))) THEN
            got_var(idv2dn)=.true.
            avg(ng)%pioVar(idv2dn)%vd=var_desc(i)
            avg(ng)%pioVar(idv2dn)%dkind=pio_fout
            avg(ng)%pioVar(idv2dn)%gtype=r2dvar
#  ifdef SOLVE3D
#   ifdef FORWARD_WRITE
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idufx1))) THEN
            got_var(idufx1)=.true.
            avg(ng)%pioVar(idufx1)%vd=var_desc(i)
            avg(ng)%pioVar(idufx1)%dkind=pio_fout
            avg(ng)%pioVar(idufx1)%gtype=u2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idufx2))) THEN
            got_var(idufx2)=.true.
            avg(ng)%pioVar(idufx2)%vd=var_desc(i)
            avg(ng)%pioVar(idufx2)%dkind=pio_fout
            avg(ng)%pioVar(idufx2)%gtype=u2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvfx1))) THEN
            got_var(idvfx1)=.true.
            avg(ng)%pioVar(idvfx1)%vd=var_desc(i)
            avg(ng)%pioVar(idvfx1)%dkind=pio_fout
            avg(ng)%pioVar(idvfx1)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvfx2))) THEN
            got_var(idvfx2)=.true.
            avg(ng)%pioVar(idvfx2)%vd=var_desc(i)
            avg(ng)%pioVar(idvfx2)%dkind=pio_fout
            avg(ng)%pioVar(idvfx2)%gtype=v2dvar
#   endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduvel))) THEN
            got_var(iduvel)=.true.
            avg(ng)%pioVar(iduvel)%vd=var_desc(i)
            avg(ng)%pioVar(iduvel)%dkind=pio_fout
            avg(ng)%pioVar(iduvel)%gtype=u3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvvel))) THEN
            got_var(idvvel)=.true.
            avg(ng)%pioVar(idvvel)%vd=var_desc(i)
            avg(ng)%pioVar(idvvel)%dkind=pio_fout
            avg(ng)%pioVar(idvvel)%gtype=v3dvar
#   if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu3dd))) THEN
            got_var(idu3dd)=.true.
            avg(ng)%pioVar(idu3dd)%vd=var_desc(i)
            avg(ng)%pioVar(idu3dd)%dkind=pio_fout
            avg(ng)%pioVar(idu3dd)%gtype=u3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv3dd))) THEN
            got_var(idv3dd)=.true.
            avg(ng)%pioVar(idv3dd)%vd=var_desc(i)
            avg(ng)%pioVar(idv3dd)%dkind=pio_fout
            avg(ng)%pioVar(idv3dd)%gtype=v3dvar
#   endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu3de))) THEN
            got_var(idu3de)=.true.
            avg(ng)%pioVar(idu3de)%vd=var_desc(i)
            avg(ng)%pioVar(idu3de)%dkind=pio_fout
            avg(ng)%pioVar(idu3de)%gtype=r3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv3dn))) THEN
            got_var(idv3dn)=.true.
            avg(ng)%pioVar(idv3dn)%vd=var_desc(i)
            avg(ng)%pioVar(idv3dn)%dkind=pio_fout
            avg(ng)%pioVar(idv3dn)%gtype=r3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idovel))) THEN
            got_var(idovel)=.true.
            avg(ng)%pioVar(idovel)%vd=var_desc(i)
            avg(ng)%pioVar(idovel)%dkind=pio_fout
            avg(ng)%pioVar(idovel)%gtype=w3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idwvel))) THEN
            got_var(idwvel)=.true.
            avg(ng)%pioVar(idwvel)%vd=var_desc(i)
            avg(ng)%pioVar(idwvel)%dkind=pio_fout
            avg(ng)%pioVar(idwvel)%gtype=w3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iddano))) THEN
            got_var(iddano)=.true.
            avg(ng)%pioVar(iddano)%vd=var_desc(i)
            avg(ng)%pioVar(iddano)%dkind=pio_fout
            avg(ng)%pioVar(iddano)%gtype=r3dvar
#   ifdef LMD_SKPP
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idhsbl))) THEN
            got_var(idhsbl)=.true.
            avg(ng)%pioVar(idhsbl)%vd=var_desc(i)
            avg(ng)%pioVar(idhsbl)%dkind=pio_fout
            avg(ng)%pioVar(idhsbl)%gtype=r2dvar
#   endif
#   ifdef LMD_BKPP
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idhbbl))) THEN
            got_var(idhbbl)=.true.
            avg(ng)%pioVar(idhbbl)%vd=var_desc(i)
            avg(ng)%pioVar(idhbbl)%dkind=pio_fout
            avg(ng)%pioVar(idhbbl)%gtype=r2dvar
#   endif
#  endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,id2dpv))) THEN
            got_var(id2dpv)=.true.
            avg(ng)%pioVar(id2dpv)%vd=var_desc(i)
            avg(ng)%pioVar(id2dpv)%dkind=pio_fout
            avg(ng)%pioVar(id2dpv)%gtype=p2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,id2drv))) THEN
            got_var(id2drv)=.true.
            avg(ng)%pioVar(id2drv)%vd=var_desc(i)
            avg(ng)%pioVar(id2drv)%dkind=pio_fout
            avg(ng)%pioVar(id2drv)%gtype=p2dvar
#  ifdef SOLVE3D
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,id3dpv))) THEN
            got_var(id3dpv)=.true.
            avg(ng)%pioVar(id3dpv)%vd=var_desc(i)
            avg(ng)%pioVar(id3dpv)%dkind=pio_fout
            avg(ng)%pioVar(id3dpv)%gtype=p3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,id3drv))) THEN
            got_var(id3drv)=.true.
            avg(ng)%pioVar(id3drv)%vd=var_desc(i)
            avg(ng)%pioVar(id3drv)%dkind=pio_fout
            avg(ng)%pioVar(id3drv)%gtype=p3dvar
#  endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idzzav))) THEN
            got_var(idzzav)=.true.
            avg(ng)%pioVar(idzzav)%vd=var_desc(i)
            avg(ng)%pioVar(idzzav)%dkind=pio_fout
            avg(ng)%pioVar(idzzav)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idu2av))) THEN
            got_var(idu2av)=.true.
            avg(ng)%pioVar(idu2av)%vd=var_desc(i)
            avg(ng)%pioVar(idu2av)%dkind=pio_fout
            avg(ng)%pioVar(idu2av)%gtype=u2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idv2av))) THEN
            got_var(idv2av)=.true.
            avg(ng)%pioVar(idv2av)%vd=var_desc(i)
            avg(ng)%pioVar(idv2av)%dkind=pio_fout
            avg(ng)%pioVar(idv2av)%gtype=v2dvar
#  ifdef SOLVE3D
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idhuav))) THEN
            got_var(idhuav)=.true.
            avg(ng)%pioVar(idhuav)%vd=var_desc(i)
            avg(ng)%pioVar(idhuav)%dkind=pio_fout
            avg(ng)%pioVar(idhuav)%gtype=u3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idhvav))) THEN
            got_var(idhvav)=.true.
            avg(ng)%pioVar(idhvav)%vd=var_desc(i)
            avg(ng)%pioVar(idhvav)%dkind=pio_fout
            avg(ng)%pioVar(idhvav)%gtype=v3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduuav))) THEN
            got_var(iduuav)=.true.
            avg(ng)%pioVar(iduuav)%vd=var_desc(i)
            avg(ng)%pioVar(idhvav)%dkind=pio_fout
            avg(ng)%pioVar(idhvav)%gtype=v3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduvav))) THEN
            got_var(iduvav)=.true.
            avg(ng)%pioVar(iduvav)%vd=var_desc(i)
            avg(ng)%pioVar(iduvav)%dkind=pio_fout
            avg(ng)%pioVar(iduvav)%gtype=r3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvvav))) THEN
            got_var(idvvav)=.true.
            avg(ng)%pioVar(idvvav)%vd=var_desc(i)
            avg(ng)%pioVar(idvvav)%dkind=pio_fout
            avg(ng)%pioVar(idvvav)%gtype=v3dvar
#   if defined LMD_MIXING || defined MY25_MIXING || defined GLS_MIXING
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvvis))) THEN
            got_var(idvvis)=.true.
            avg(ng)%pioVar(idvvis)%vd=var_desc(i)
            avg(ng)%pioVar(idvvis)%dkind=pio_fout
            avg(ng)%pioVar(idvvis)%gtype=w3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idtdif))) THEN
            got_var(idtdif)=.true.
            avg(ng)%pioVar(idtdif)%vd=var_desc(i)
            avg(ng)%pioVar(idtdif)%dkind=pio_fout
            avg(ng)%pioVar(idtdif)%gtype=w3dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idsdif))) THEN
            got_var(idsdif)=.true.
            avg(ng)%pioVar(idsdif)%vd=var_desc(i)
            avg(ng)%pioVar(idsdif)%dkind=pio_fout
            avg(ng)%pioVar(idsdif)%gtype=w3dvar
#   endif
#   if defined BULK_FLUXES || defined ECOSIM || defined ATM_PRESS
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idpair))) THEN
            got_var(idpair)=.true.
            avg(ng)%pioVar(idpair)%vd=var_desc(i)
            avg(ng)%pioVar(idpair)%dkind=pio_fout
            avg(ng)%pioVar(idpair)%gtype=r2dvar
#   endif
#   if defined BULK_FLUXES || defined ECOSIM
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduair))) THEN
            got_var(iduair)=.true.
            avg(ng)%pioVar(iduair)%vd=var_desc(i)
            avg(ng)%pioVar(iduair)%dkind=pio_fout
            avg(ng)%pioVar(iduair)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvair))) THEN
            got_var(idvair)=.true.
            avg(ng)%pioVar(idvair)%vd=var_desc(i)
            avg(ng)%pioVar(idvair)%dkind=pio_fout
            avg(ng)%pioVar(idvair)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduaie))) THEN
            got_var(iduaie)=.true.
            avg(ng)%pioVar(iduaie)%vd=var_desc(i)
            avg(ng)%pioVar(iduaie)%dkind=pio_fout
            avg(ng)%pioVar(iduaie)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvain))) THEN
            got_var(idvain)=.true.
            avg(ng)%pioVar(idvain)%vd=var_desc(i)
            avg(ng)%pioVar(idvain)%dkind=pio_fout
            avg(ng)%pioVar(idvain)%gtype=r2dvar
#   endif
          ELSE IF (trim(var_name(i)).eq.                                &
     &             trim(vname(1,idtsur(itemp)))) THEN
            got_var(idtsur(itemp))=.true.
            avg(ng)%pioVar(idtsur(itemp))%vd=var_desc(i)
            avg(ng)%pioVar(idtsur(itemp))%dkind=pio_fout
            avg(ng)%pioVar(idtsur(itemp))%gtype=r2dvar
#   ifdef SALINITY
          ELSE IF (trim(var_name(i)).eq.                                &
     &             trim(vname(1,idtsur(isalt)))) THEN
            got_var(idtsur(isalt))=.true.
            avg(ng)%pioVar(idtsur(isalt))%vd=var_desc(i)
            avg(ng)%pioVar(idtsur(isalt))%dkind=pio_fout
            avg(ng)%pioVar(idtsur(isalt))%gtype=r2dvar
#   endif
#   if defined BULK_FLUXES || defined FRC_COUPLING
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idlhea))) THEN
            got_var(idlhea)=.true.
            avg(ng)%pioVar(idlhea)%vd=var_desc(i)
            avg(ng)%pioVar(idlhea)%dkind=pio_fout
            avg(ng)%pioVar(idlhea)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idshea))) THEN
            got_var(idshea)=.true.
            avg(ng)%pioVar(idshea)%vd=var_desc(i)
            avg(ng)%pioVar(idshea)%dkind=pio_fout
            avg(ng)%pioVar(idshea)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idlrad))) THEN
            got_var(idlrad)=.true.
            avg(ng)%pioVar(idlrad)%vd=var_desc(i)
            avg(ng)%pioVar(idlrad)%dkind=pio_fout
            avg(ng)%pioVar(idlrad)%gtype=r2dvar
#   endif
#   ifdef BULK_FLUXES
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idtair))) THEN
            got_var(idtair)=.true.
            avg(ng)%pioVar(idtair)%vd=var_desc(i)
            avg(ng)%pioVar(idtair)%dkind=pio_fout
            avg(ng)%pioVar(idtair)%gtype=r2dvar
#    ifdef EMINUSP
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idevap))) THEN
            got_var(idevap)=.true.
            avg(ng)%pioVar(idevap)%vd=var_desc(i)
            avg(ng)%pioVar(idevap)%dkind=pio_fout
            avg(ng)%pioVar(idevap)%gtype=r2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idrain))) THEN
            got_var(idrain)=.true.
            avg(ng)%pioVar(idrain)%vd=var_desc(i)
            avg(ng)%pioVar(idrain)%dkind=pio_fout
            avg(ng)%pioVar(idrain)%gtype=r2dvar
#    endif
#   endif
#   ifdef SHORTWAVE
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idsrad))) THEN
            got_var(idsrad)=.true.
            avg(ng)%pioVar(idsrad)%vd=var_desc(i)
            avg(ng)%pioVar(idsrad)%dkind=pio_fout
            avg(ng)%pioVar(idsrad)%gtype=r2dvar
#   endif
#  endif
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idusms))) THEN
            got_var(idusms)=.true.
            avg(ng)%pioVar(idusms)%vd=var_desc(i)
            avg(ng)%pioVar(idusms)%dkind=pio_fout
            avg(ng)%pioVar(idusms)%gtype=u2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvsms))) THEN
            got_var(idvsms)=.true.
            avg(ng)%pioVar(idvsms)%vd=var_desc(i)
            avg(ng)%pioVar(idvsms)%dkind=pio_fout
            avg(ng)%pioVar(idvsms)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubms))) THEN
            got_var(idubms)=.true.
            avg(ng)%pioVar(idubms)%vd=var_desc(i)
            avg(ng)%pioVar(idubms)%dkind=pio_fout
            avg(ng)%pioVar(idubms)%gtype=u2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbms))) THEN
            got_var(idvbms)=.true.
            avg(ng)%pioVar(idvbms)%vd=var_desc(i)
            avg(ng)%pioVar(idvbms)%dkind=pio_fout
            avg(ng)%pioVar(idvbms)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubrs))) THEN
            got_var(idubrs)=.true.
            avg(ng)%pioVar(idubrs)%vd=var_desc(i)
            avg(ng)%pioVar(idubrs)%dkind=pio_fout
            avg(ng)%pioVar(idubrs)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbrs))) THEN
            got_var(idvbrs)=.true.
            avg(ng)%pioVar(idvbrs)%vd=var_desc(i)
            avg(ng)%pioVar(idvbrs)%dkind=pio_fout
            avg(ng)%pioVar(idvbrs)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubws))) THEN
            got_var(idubws)=.true.
            avg(ng)%pioVar(idubws)%vd=var_desc(i)
            avg(ng)%pioVar(idubws)%dkind=pio_fout
            avg(ng)%pioVar(idubws)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbws))) THEN
            got_var(idvbws)=.true.
            avg(ng)%pioVar(idvbws)%vd=var_desc(i)
            avg(ng)%pioVar(idvbws)%dkind=pio_fout
            avg(ng)%pioVar(idvbws)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubcs))) THEN
            got_var(idubcs)=.true.
            avg(ng)%pioVar(idubcs)%vd=var_desc(i)
            avg(ng)%pioVar(idubcs)%dkind=pio_fout
            avg(ng)%pioVar(idubcs)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbcs))) THEN
            got_var(idvbcs)=.true.
            avg(ng)%pioVar(idvbcs)%vd=var_desc(i)
            avg(ng)%pioVar(idvbcs)%dkind=pio_fout
            avg(ng)%pioVar(idvbcs)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduvwc))) THEN
            got_var(iduvwc)=.true.
            avg(ng)%pioVar(iduvwc)%vd=var_desc(i)
            avg(ng)%pioVar(iduvwc)%dkind=pio_fout
            avg(ng)%pioVar(iduvwc)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubot))) THEN
            got_var(idubot)=.true.
            avg(ng)%pioVar(idubot)%vd=var_desc(i)
            avg(ng)%pioVar(idubot)%dkind=pio_fout
            avg(ng)%pioVar(idubot)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbot))) THEN
            got_var(idvbot)=.true.
            avg(ng)%pioVar(idvbot)%vd=var_desc(i)
            avg(ng)%pioVar(idvbot)%dkind=pio_fout
            avg(ng)%pioVar(idvbot)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubur))) THEN
            got_var(idubur)=.true.
            avg(ng)%pioVar(idubur)%vd=var_desc(i)
            avg(ng)%pioVar(idubur)%dkind=pio_fout
            avg(ng)%pioVar(idubur)%gtype=v2dvar
          ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbvr))) THEN
            got_var(idvbvr)=.true.
            avg(ng)%pioVar(idvbvr)%vd=var_desc(i)
            avg(ng)%pioVar(idvbvr)%dkind=pio_fout
            avg(ng)%pioVar(idvbvr)%gtype=v2dvar
          END IF
#  ifdef SOLVE3D
          DO itrc=1,nt(ng)
            IF (trim(var_name(i)).eq.trim(vname(1,idtvar(itrc)))) THEN
              got_var(idtvar(itrc))=.true.
              avg(ng)%pioTrc(itrc)%vd=var_desc(i)
              avg(ng)%pioTrc(itrc)%dkind=pio_fout
              avg(ng)%pioTrc(itrc)%gtype=r3dvar
            END IF
          END DO
          DO itrc=1,nat
            IF (trim(var_name(i)).eq.trim(vname(1,ihutav(itrc)))) THEN
              got_var(ihutav(itrc))=.true.
              avg(ng)%pioVar(ihutav(itrc))%vd=var_desc(i)
              avg(ng)%pioVar(ihutav(itrc))%dkind=pio_fout
              avg(ng)%pioVar(ihutav(itrc))%gtype=u3dvar
            ELSE IF (trim(var_name(i)).eq.                              &
     &               trim(vname(1,ihvtav(itrc)))) THEN
              got_var(ihvtav(itrc))=.true.
              avg(ng)%pioVar(ihvtav(itrc))%vd=var_desc(i)
              avg(ng)%pioVar(ihvtav(itrc))%dkind=pio_fout
              avg(ng)%pioVar(ihvtav(itrc))%gtype=v3dvar
            ELSE IF (trim(var_name(i)).eq.                              &
     &               trim(vname(1,idutav(itrc)))) THEN
              got_var(idutav(itrc))=.true.
              avg(ng)%pioVar(idutav(itrc))%vd=var_desc(i)
              avg(ng)%pioVar(idutav(itrc))%dkind=pio_fout
              avg(ng)%pioVar(idutav(itrc))%gtype=u3dvar
            ELSE IF (trim(var_name(i)).eq.                              &
     &               trim(vname(1,idvtav(itrc)))) THEN
              got_var(idvtav(itrc))=.true.
              avg(ng)%pioVar(idvtav(itrc))%vd=var_desc(i)
              avg(ng)%pioVar(idvtav(itrc))%dkind=pio_fout
              avg(ng)%pioVar(idvtav(itrc))%gtype=v3dvar
            ELSE IF (trim(var_name(i)).eq.                              &
     &               trim(vname(1,idttav(itrc)))) THEN
              got_var(idttav(itrc))=.true.
              avg(ng)%pioVar(idttav(itrc))%vd=var_desc(i)
              avg(ng)%pioVar(idttav(itrc))%dkind=pio_fout
              avg(ng)%pioVar(idttav(itrc))%gtype=r3dvar
            END IF
          END DO
#   if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
          DO itrc=1,nat
            IF (trim(var_name(i)).eq.trim(vname(1,idtrcd(itrc)))) THEN
              got_var(idtrcd(itrc))=.true.
              avg(ng)%pioVar(idtrcd(itrc))%vd=var_desc(i)
              avg(ng)%pioVar(idtrcd(itrc))%dkind=pio_fout
              avg(ng)%pioVar(idtrcd(itrc))%gtype=r3dvar
            END IF
          END DO
#   endif
#  endif
        END DO
!
!  Check if averages variables are available in input NetCDF file.
!
        IF (.not.got_var(idtime)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtime)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idfsur).and.aout(idfsur,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idfsur)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idubar).and.aout(idubar,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idubar)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvbar).and.aout(idvbar,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvbar)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
        IF (.not.got_var(idfsud).and.aout(idfsud,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idfsud)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idu2dd).and.aout(idu2dd,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu2dd)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv2dd).and.aout(idv2dd,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv2dd)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
        IF (.not.got_var(idu2de).and.aout(idu2de,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu2de)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv2dn).and.aout(idv2dn,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv2dn)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  ifdef SOLVE3D
#   ifdef FORWARD_WRITE
        IF (.not.got_var(idufx1).and.aout(idufx1,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idufx1)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idufx2).and.aout(idufx2,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idufx2)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvfx1).and.aout(idvfx1,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvfx1)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvfx2).and.aout(idvfx2,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvfx2)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
        IF (.not.got_var(iduvel).and.aout(iduvel,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduvel)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvvel).and.aout(idvvel,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvvel)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
        IF (.not.got_var(idu3dd).and.aout(idu3dd,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu3dd)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv3dd).and.aout(idv3dd,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv3dd)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
        IF (.not.got_var(idu3de).and.aout(idu3de,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu3de)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv3dn).and.aout(idv3dn,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv3dn)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idovel).and.aout(idovel,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idovel)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idwvel).and.aout(idwvel,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idwvel)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(iddano).and.aout(iddano,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iddano)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   ifdef LMD_SKPP
        IF (.not.got_var(idhsbl).and.aout(idhsbl,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idhsbl)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
#   ifdef LMD_BKPP
        IF (.not.got_var(idhbbl).and.aout(idhbbl,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idhbbl)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
#  endif
        IF (.not.got_var(id2dpv).and.aout(id2dpv,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,id2dpv)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(id2drv).and.aout(id2drv,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,id2drv)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  ifdef SOLVE3D
        IF (.not.got_var(id3dpv).and.aout(id3dpv,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,id3dpv)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(id3drv).and.aout(id3drv,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,id3drv)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  endif
        IF (.not.got_var(idzzav).and.aout(idzzav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idzzav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idu2av).and.aout(idu2av,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idu2av)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idv2av).and.aout(idv2av,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idv2av)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  ifdef SOLVE3D
        IF (.not.got_var(idhuav).and.aout(idhuav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idhuav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idhvav).and.aout(idhvav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idhvav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(iduuav).and.aout(iduuav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduuav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(iduvav).and.aout(iduvav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduvav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvvav).and.aout(idvvav,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvvav)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   if defined LMD_MIXING || defined MY25_MIXING || defined GLS_MIXING
        IF (.not.got_var(idvvis).and.aout(idvvis,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvvis)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idtdif).and.aout(idtdif,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtdif)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#    ifdef SALINITY
        IF (.not.got_var(idsdif).and.aout(idsdif,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idsdif)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#    endif
#   endif
#   if defined BULK_FLUXES || defined ECOSIM || defined ATM_PRESS
        IF (.not.got_var(idpair).and.aout(idpair,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idpair)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
#   if defined BULK_FLUXES || defined ECOSIM
        IF (.not.got_var(iduair).and.aout(iduair,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduair)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvair).and.aout(idvair,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvair)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(iduaie).and.aout(iduaie,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,iduaie)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvain).and.aout(idvain,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvain)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
        IF (.not.got_var(idtsur(itemp)).and.aout(idtsur(itemp),ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtsur(itemp))),   &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   ifdef SALINITY
        IF (.not.got_var(idtsur(isalt)).and.aout(idtsur(isalt),ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtsur(isalt))),   &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
#   if defined BULK_FLUXES || defined FRC_COUPLING
        IF (.not.got_var(idlhea).and.aout(idlhea,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idlhea)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idshea).and.aout(idshea,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idshea)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idlrad).and.aout(idlrad,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idlrad)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
#   ifdef BULK_FLUXES
        IF (.not.got_var(idtair).and.aout(idtair,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idtair)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#    ifdef EMINUSP
        IF (.not.got_var(idevap).and.aout(idevap,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idevap)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idrain).and.aout(idrain,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idrain)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#    endif
#   endif
#   ifdef SHORTWAVE
        IF (.not.got_var(idsrad).and.aout(idsrad,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idsrad)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#   endif
#  endif
        IF (.not.got_var(idusms).and.aout(idusms,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idusms)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvsms).and.aout(idvsms,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvsms)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idubms).and.aout(idubms,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idubms)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
        IF (.not.got_var(idvbms).and.aout(idvbms,ng)) THEN
          IF (master) WRITE (stdout,60) trim(vname(1,idvbms)),          &
     &                                  trim(ncname)
          exit_flag=3
          RETURN
        END IF
#  ifdef SOLVE3D
        DO itrc=1,nt(ng)
          IF (.not.got_var(idtvar(itrc)).and.aout(idtvar(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idtvar(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
        END DO
        DO itrc=1,nat
          IF (.not.got_var(ihutav(itrc)).and.aout(ihutav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,ihutav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
          IF (.not.got_var(ihvtav(itrc)).and.aout(ihvtav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,ihvtav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
          IF (.not.got_var(idutav(itrc)).and.aout(idutav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idutav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
          IF (.not.got_var(idvtav(itrc)).and.aout(idvtav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idvtav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
          IF (.not.got_var(idttav(itrc)).and.aout(idttav(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idttav(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
        END DO
#   if defined AVERAGES_DETIDE && (defined SSH_TIDES || defined UV_TIDES)
        DO itrc=1,nat
          IF (.not.got_var(idtrcd(itrc)).and.aout(idtrcd(itrc),ng)) THEN
            IF (master) WRITE (stdout,60) trim(vname(1,idtrcd(itrc))),  &
     &                                    trim(ncname)
            exit_flag=3
            RETURN
          END IF
        END DO
#   endif
#  endif
 
#  if (defined BBL_MODEL || defined WAVES_OUTPUT) && defined SOLVE3D
!
!  Scan bottom boundary layer model and waves variables from input
!  NetCDF and activate switches for time-averaged variables. Get
!  variable IDs.
!
        CALL bbl_def_pio (ng, model, ldef, aout, avg)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
 
#  if defined ICE_MODEL && defined SOLVE3D
!
!  Scan sea-ice variables from input NetCDF and activate switches for
!  time-averaged variables. Get variable IDs.
!
        CALL ice_def_pio (ng, model, ldef, aout, avg)
#  endif
 
# if defined SEDIMENT && defined SOLVE3D
!
!  Scan sediment model variables from input NetCDF and activate
!  switches for time-averaged variables. Get variable IDs.
!
        CALL sediment_def_pio (ng, model, ldef, aout, avg)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
# endif
 
#  if defined WEC_VF && defined SOLVE3D
!
!  Scan Waves Effect on Currents variables from input NetCDF and
!  activate switches for time-averaged variables. Get variable IDs.
!
        CALL wec_def_pio (ng, model, ldef, aout, avg)
        IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
#  endif
!
!  Set unlimited time record dimension to the appropriate value.
!
        IF (nrst(ng).eq.navg(ng)) THEN
          IF (ndefavg(ng).gt.0) THEN
            avg(ng)%Rindex=((ntstart(ng)-1)-                            &
     &                      ndefavg(ng)*((ntstart(ng)-1)/ndefavg(ng)))/ &
     &                     navg(ng)
          ELSE
            avg(ng)%Rindex=(ntstart(ng)-1)/navg(ng)
          END IF
        ELSE
          avg(ng)%Rindex=rec_size
        END IF
      END IF query
!
!  Set initial average time. Notice that the value is offset by half
!  nAVG*dt so there is not a special case when computing its value
 
#  ifdef ADJOINT
!  in "ad_set_avg".
!
        avgtime(ng)=time(ng)+0.5_r8*real(navg(ng),r8)*dt(ng)
#  else
!  in "set_avg".
!
      IF (ntsavg(ng).eq.1) THEN
        avgtime(ng)=time(ng)-0.5_r8*real(navg(ng),r8)*dt(ng)
      ELSE
        avgtime(ng)=time(ng)+real(ntsavg(ng),r8)*dt(ng)-                &
     &              0.5_r8*real(navg(ng),r8)*dt(ng)
      END IF
#  endif
!
  10  FORMAT (2x,'DEF_AVG_PIO      - creating average file,',t56,       &
     &        'Grid ',i2.2,': ',a)
  20  FORMAT (2x,'DEF_AVG_PIO      - inquiring average file',t56,       &
     &        'Grid ',i2.2,': ',a)
  30  FORMAT (/,' DEF_AVG_PIO - unable to create averages NetCDF',      &
     &        ' file: ',a)
  40  FORMAT (1pe11.4,1x,'millimeter')
  50  FORMAT (/,' DEF_AVG_PIO - unable to open averages NetCDF'         &
     &        ' file: ',a)
  60  FORMAT (/,' DEF_AVG_PIO - unable to find variable: ',a,2x,        &
     &        ' in averages NetCDF file: ',a)
!
      RETURN

References mod_ncparam::aout, mod_iounits::avg, mod_scalars::avgtime, bbl_output_mod::bbl_def_pio(), mod_scalars::dt, mod_scalars::exit_flag, strings_mod::founderror(), mod_ncparam::id2dpv, mod_ncparam::id2drv, mod_ncparam::id3dpv, mod_ncparam::id3drv, mod_ncparam::iddano, mod_ncparam::idevap, mod_ncparam::idfsud, mod_ncparam::idfsur, mod_ncparam::idhbbl, mod_ncparam::idhsbl, mod_ncparam::idhuav, mod_ncparam::idhvav, mod_ncparam::idlhea, mod_ncparam::idlrad, mod_ncparam::idovel, mod_ncparam::idpair, mod_ncparam::idrain, mod_ncparam::idrwet, mod_ncparam::idsdif, mod_sediment::idsed, mod_ncparam::idshea, mod_ncparam::idsrad, mod_ncparam::idtair, mod_ncparam::idtdif, mod_ncparam::idtime, mod_ncparam::idtrcd, mod_ncparam::idtsur, mod_ncparam::idttav, mod_ncparam::idtvar, mod_ncparam::idu2av, mod_ncparam::idu2dd, mod_ncparam::idu2de, mod_ncparam::idu3dd, mod_ncparam::idu3de, mod_ncparam::iduaie, mod_ncparam::iduair, mod_ncparam::idubar, mod_ncparam::idubcs, mod_ncparam::idubms, mod_ncparam::idubot, mod_ncparam::idubrs, mod_ncparam::idubur, mod_ncparam::idubws, mod_ncparam::idufx1, mod_ncparam::idufx2, mod_ncparam::idusms, mod_ncparam::idutav, mod_ncparam::iduuav, mod_ncparam::iduvav, mod_ncparam::iduvel, mod_ncparam::iduvwc, mod_ncparam::idv2av, mod_ncparam::idv2dd, mod_ncparam::idv2dn, mod_ncparam::idv3dd, mod_ncparam::idv3dn, mod_ncparam::idvain, mod_ncparam::idvair, mod_ncparam::idvbar, mod_ncparam::idvbcs, mod_ncparam::idvbms, mod_ncparam::idvbot, mod_ncparam::idvbrs, mod_ncparam::idvbvr, mod_ncparam::idvbws, mod_ncparam::idvfx1, mod_ncparam::idvfx2, mod_ncparam::idvsms, mod_ncparam::idvtav, mod_ncparam::idvvav, mod_ncparam::idvvel, mod_ncparam::idvvis, mod_ncparam::idwvel, mod_ncparam::idzzav, mod_ncparam::ihutav, mod_ncparam::ihvtav, mod_ncparam::iinfo, mod_param::iobounds, mod_scalars::isalt, mod_scalars::itemp, mod_parallel::master, mod_param::n, mod_param::nat, mod_scalars::navg, mod_biology::nbac, mod_biology::nbands, mod_param::nbed, mod_scalars::ndefavg, mod_biology::ndom, mod_biology::nfec, mod_scalars::noerror, mod_biology::nphy, mod_scalars::nrst, mod_param::nst, mod_fourdvar::nstatevar, mod_param::nt, mod_scalars::ntsavg, mod_scalars::ntstart, mod_ncparam::nv, mod_param::p2dvar, mod_param::p3dvar, mod_pio_netcdf::pio_fout, mod_pio_netcdf::pio_netcdf_check_dim(), mod_pio_netcdf::pio_netcdf_create(), mod_pio_netcdf::pio_netcdf_enddef(), mod_pio_netcdf::pio_netcdf_inq_var(), mod_pio_netcdf::pio_netcdf_open(), mod_pio_netcdf::pio_tout, mod_param::r2dvar, mod_param::r3dvar, mod_scalars::rclock, mod_sediment::sd50, sediment_output_mod::sediment_def_pio(), mod_iounits::sourcefile, mod_iounits::stdout, mod_scalars::time, mod_param::u2dvar, mod_param::u3dvar, mod_param::v2dvar, mod_param::v3dvar, mod_pio_netcdf::var_desc, mod_ncparam::vname, and mod_param::w3dvar.

Referenced by def_avg().

Here is the call graph for this function:

Here is the caller graph for this function: