C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/salt_plume_diagnostics_fill.F,v 1.3 2014/05/02 06:10:48 atn Exp $
C $Name:  $

#include "SALT_PLUME_OPTIONS.h"

CBOP 0
C     !ROUTINE: SALT_PLUME_DIAGNOSTICS_FILL
C     !INTERFACE:
      SUBROUTINE SALT_PLUME_DIAGNOSTICS_FILL(
     I             bi, bj, myThid )

C     !DESCRIPTION:
C     Provide the information for the diagnostics

C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "SALT_PLUME.h"

C     !INPUT PARAMETERS:
      INTEGER bi,bj
      INTEGER myThid
CEOP

#ifdef ALLOW_SALT_PLUME
#ifdef ALLOW_DIAGNOSTICS
C     !LOCAL VARIABLES:

      IF ( useDiagnostics ) THEN

C-    Salt plume depth (>0), [m]
       CALL DIAGNOSTICS_FILL( SaltPlumeDepth,
     &                             'oceSPDep',0, 1,0,1,1,myThid )

C-    net Salt plume flux into the ocean during freezing (+=down), [psu.kg/m2/s ~ g/m2/s]
       CALL DIAGNOSTICS_FILL( saltPlumeFlux,
     &                             'oceSPflx',0, 1,0,1,1,myThid )
#ifdef SALT_PLUME_VOLUME
       CALL DIAGNOSTICS_FILL(
     &      SPbrineVolFlux,'oceSPvol',0,1,0,1,1,myThid )
Cretire       CALL DIAGNOSTICS_FILL(
Cretire     &      SPkBottom,'SPkbottm',0,1,0,1,1,myThid )
Cretire       CALL DIAGNOSTICS_FILL(
Cretire     &      SPbrineSalt,'SPbrineS',0,1,0,1,1,myThid )
C       CALL DIAGNOSTICS_FILL(
C     &      SPplumek,'PLUMEKB ',0,Nr,0,1,1,myThid )
       CALL DIAGNOSTICS_FILL(
     &      dSPvolSurf2kLev, 'SPVsuf2k',0,Nr,0,1,1,myThid )
Cretire       CALL DIAGNOSTICS_FILL(
Cretire     &      dSPvolBelow2kLev,'SPVp2k  ',0,Nr,0,1,1,myThid )
       CALL DIAGNOSTICS_FILL(
     &      dSPvolkLev2Above,'SPVk2m  ',0,Nr,0,1,1,myThid )
#endif
      ENDIF

#endif /* ALLOW_DIAGNOSTICS */
#endif /* ALLOW_SALT_PLUME */

      RETURN
      END