C***********************************************************************
C***********************************************************************
      SUBROUTINE OUTPUTQ4(KD,SECKD,ITER,DELT,MESH,FILINQ,FILICQ,
     &                    Cd,AHI,NN,NNV,NE,X,Y,IN,BATHY,
     &                    TOTDEP,UBAR,VBAR,Z,U,V,W,Q2,Q2L,RHO,T,S,
     &                    ENZM,ENZH,ENZQ,
     &                    TOTDEPOLD,UBAROLD,VBAROLD,ZOLD,QP)
C Include file for parameter statements
      INCLUDE 'QUODDY.DIM'
C Fixed global variables
      REAL AHI(NNDIM),X(NNDIM),Y(NNDIM),BATHY(NNDIM)
      INTEGER IN(3,NEDIM)
      REAL TOTDEP(NNDIM),UBAR(NNDIM),VBAR(NNDIM)
      REAL Z(NNDIM,NNVDIM),U(NNDIM,NNVDIM),V(NNDIM,NNVDIM),
     &W(NNDIM,NNVDIM),Q2(NNDIM,NNVDIM),Q2L(NNDIM,NNVDIM),
     &RHO(NNDIM,NNVDIM),T(NNDIM,NNVDIM),S(NNDIM,NNVDIM),
     &ENZM(NNDIM,NNVDIM),ENZH(NNDIM,NNVDIM),ENZQ(NNDIM,NNVDIM)
      REAL TOTDEPOLD(NNDIM),UBAROLD(NNDIM),VBAROLD(NNDIM),
     &ZOLD(NNDIM,NNVDIM),QP(NNDIM)
      CHARACTER*72 MESH,FILINQ,FILICQ
C-----------------------------------------------------------------------
C This user subroutine is supplied to act as an interface for the
C   user to output the 2-D and 3-D results at a particular time in the
C   simulation.
C
C ... ... ...
C ... ... ...
C ... ... ...
C
C     &KD,SECKD,DELT,TOTDEP,UBAR,VBAR,TOTDEPOLD,UBAROLD,VBAROLD,
C     &Z,ZOLD,U,V,W,Q2,Q2L,T,S,NN,NNV,NNDIM,NNVDIM)
C
C================Beginning-of-user-specified-instructions===============

C DROG3DDT PARAMETERS
      CHARACTER DROG*3,FILDR3*72  
      DATA DROG/'YES'/           ! Flag for drogue tracking
      DATA FILDR3/'test.ind'/    ! drogue initilization file name
      INTEGER ISHIFT             ! Integer to shift drog starting time
      DATA ISHIFT/0/             ! no time shift for drogue start
C
C CALL DROG3DDT IF DROG=='YES'; ITER SHIFT FOR DELAYED DROG START
C MESH is the 'gridname' for the grid files (.nod,.ele,...) that
C is used to find the corresponding .gr2 filename.
C
      T1=ITER*DELT/3600.  ! TIME IN HOURS
      T2=T1+DELT/3600
      IF (DROG.EQ.'YES') THEN
         IF(ITER.GE.ISHIFT)THEN
            CALL DROG3DDT(NN,NE,NNV,T2,T1,ITER-ISHIFT,
     +                    Z,U,V,W,MESH,FILDR3)
         END IF
      END IF


...
...  etc.etc.etc.
...