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	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	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	IF ( (SPbrineSconst .LT. 0. _d 0).OR.
55	& (SPbrineSconst .GT. SPbrineSaltmax) ) then
56	WRITE(msgBuf,'(2A,F6.1,A)') 'SALT_PLUME_CHECK:',
57	& ' SPbrineSconst valid range: [0.,',SPbrineSaltmax,'.]'
58	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