atnguyen/code_21Dec2012_saltplume/salt_plume_check.F

C $Header: /u/gcmpack/MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/salt_plume_check.F,v 1.2 2014/04/29 06:49:40 atn Exp $
C $Name:  $

#include "SALT_PLUME_OPTIONS.h"

CBOP 0
C     !ROUTINE: SALT_PLUME_CHECK
C     !INTERFACE:
      SUBROUTINE SALT_PLUME_CHECK( myThid )

C     !DESCRIPTION:
C     Check dependances with other packages

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

C     !INPUT PARAMETERS:
C     myThid :: my Thread Id number
      INTEGER myThid
CEOP

#ifdef ALLOW_SALT_PLUME

C     !LOCAL VARIABLES:
C     msgBuf      - Informational/error meesage buffer
      CHARACTER*(MAX_LEN_MBUF) msgBuf

      _BEGIN_MASTER(myThid)

       WRITE(msgBuf,'(A)') 'SALT_PLUME_CHECK: #define SALT_PLUME'
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                     SQUEEZE_RIGHT , myThid )

      IF ( (SPsalFRAC .LT. 0. _d 0).OR.(SPsalFRAC .GT. 1. _d 0) ) then
       WRITE(msgBuf,'(2A)') 'SALT_PLUME_CHECK:',
     &  ' SPsalFRAC needs to lie between 0. and 1.'
       CALL PRINT_ERROR( msgBuf, myThid )
       STOP 'ABNORMAL END: S/R SALT_PLUME_CHECK'
      ENDIF

#ifdef SALT_PLUME_IN_LEADS
      IF ( (SPinflectionPoint .LT. 0. _d 0).OR.
     &     (SPinflectionPoint .GT. 1. _d 0) ) then
       WRITE(msgBuf,'(2A)') 'SALT_PLUME_CHECK:',
     &  ' SPinflectionPoint needs to lie between 0. and 1.'
       CALL PRINT_ERROR( msgBuf, myThid )
       STOP 'ABNORMAL END: S/R SALT_PLUME_CHECK'
      ENDIF
#endif
#ifdef SALT_PLUME_VOLUME
      IF ( (SPbrineSconst .LT. 0. _d 0).OR.
     &     (SPbrineSconst .GT. SPbrineSaltmax) ) then
       WRITE(msgBuf,'(2A,F6.1,A)') 'SALT_PLUME_CHECK:',
     &  ' SPbrineSconst valid range: [0.,',SPbrineSaltmax,'.]'
       CALL PRINT_ERROR( msgBuf, myThid )
       STOP 'ABNORMAL END: S/R SALT_PLUME_CHECK'
      ENDIF
#endif

      _END_MASTER(myThid)

#endif /* ALLOW_SALT_PLUME */

      RETURN
      END
1	atn	1.3	C $Header: /u/gcmpack/MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/salt_plume_check.F,v 1.2 2014/04/29 06:49:40 atn Exp $
2	atn	1.1	C $Name: $
3
4			#include "SALT_PLUME_OPTIONS.h"
5
6			CBOP 0
7			C !ROUTINE: SALT_PLUME_CHECK
8			C !INTERFACE:
9			SUBROUTINE SALT_PLUME_CHECK( myThid )
10
11			C !DESCRIPTION:
12			C Check dependances with other packages
13
14			C !USES:
15			IMPLICIT NONE
16			#include "EEPARAMS.h"
17			#include "SIZE.h"
18			#include "SALT_PLUME.h"
19
20			C !INPUT PARAMETERS:
21			C myThid :: my Thread Id number
22			INTEGER myThid
23			CEOP
24
25			#ifdef ALLOW_SALT_PLUME
26
27			C !LOCAL VARIABLES:
28			C msgBuf - Informational/error meesage buffer
29			CHARACTER*(MAX_LEN_MBUF) msgBuf
30
31			_BEGIN_MASTER(myThid)
32
33			WRITE(msgBuf,'(A)') 'SALT_PLUME_CHECK: #define SALT_PLUME'
34			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
35			& SQUEEZE_RIGHT , myThid )
36
37			IF ( (SPsalFRAC .LT. 0. _d 0).OR.(SPsalFRAC .GT. 1. _d 0) ) then
38			WRITE(msgBuf,'(2A)') 'SALT_PLUME_CHECK:',
39			& ' SPsalFRAC needs to lie between 0. and 1.'
40			CALL PRINT_ERROR( msgBuf, myThid )
41			STOP 'ABNORMAL END: S/R SALT_PLUME_CHECK'
42			ENDIF
43
44			#ifdef SALT_PLUME_IN_LEADS
45			IF ( (SPinflectionPoint .LT. 0. _d 0).OR.
46			& (SPinflectionPoint .GT. 1. _d 0) ) then
47			WRITE(msgBuf,'(2A)') 'SALT_PLUME_CHECK:',
48			& ' SPinflectionPoint needs to lie between 0. and 1.'
49			CALL PRINT_ERROR( msgBuf, myThid )
50			STOP 'ABNORMAL END: S/R SALT_PLUME_CHECK'
51			ENDIF
52			#endif
53			#ifdef SALT_PLUME_VOLUME
54	atn	1.2	IF ( (SPbrineSconst .LT. 0. _d 0).OR.
55	atn	1.3	& (SPbrineSconst .GT. SPbrineSaltmax) ) then
56			WRITE(msgBuf,'(2A,F6.1,A)') 'SALT_PLUME_CHECK:',
57			& ' SPbrineSconst valid range: [0.,',SPbrineSaltmax,'.]'
58	atn	1.1	CALL PRINT_ERROR( msgBuf, myThid )
59			STOP 'ABNORMAL END: S/R SALT_PLUME_CHECK'
60			ENDIF
61			#endif
62
63			_END_MASTER(myThid)
64
65			#endif /* ALLOW_SALT_PLUME */
66
67			RETURN
68			END