pkg/icefront/icefront_check.F

C $Header: /u/gcmpack/MITgcm/pkg/icefront/icefront_check.F,v 1.2 2010/01/25 22:37:19 dimitri Exp $
C $Name:  $
#include "ICEFRONT_OPTIONS.h"

      SUBROUTINE ICEFRONT_CHECK( myThid )
C     /==========================================================\
C     | SUBROUTINE ICEFRONT_CHECK                                |
C     | o Validate basic package setup and inter-package         |
C     | dependencies.                                            |
C     \==========================================================/
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "ICEFRONT.h"

C     === Routine arguments ===
C     myThid -  Number of this instance of ICEFRONT_CHECK
      INTEGER myThid

#ifdef ALLOW_ICEFRONT

C     === Local variables ===
C     msgBuf      - Informational/error meesage buffer
      CHARACTER*(MAX_LEN_MBUF) msgBuf

      _BEGIN_MASTER(myThid)

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

C     So far, ICEFRONT works only with oceanic z-coordinates
      IF ( buoyancyRelation .NE. 'OCEANIC' ) THEN
         WRITE(msgBuf,'(A)')
     &      'ICEFRONT works only with buoyancyRelation = ''OCEANIC'''
         CALL PRINT_ERROR( msgBuf, myThid )
         STOP 'ABNORMAL END: S/R ICEFRONT_CHECK'
      ENDIF

CML#ifdef ALLOW_NONHYDROSTATIC
CML      IF ( nonHydrostatic ) THEN
CML         WRITE(msgBuf,'(A)')
CML     &      'Currently ICEFRONT and nonHydrostatic cannot be turned'
CML         CALL PRINT_ERROR( msgBuf , 1)
CML         WRITE(msgBuf,'(A)') 'on at the same time'
CML         CALL PRINT_ERROR( msgBuf , 1)
CML         STOP 'ABNORMAL END: S/R ICEFRONT_CHECK'
CML      ENDIF
CML#endif /* ALLOW_NONHYDROSTATIC */

      _END_MASTER(myThid)

#endif /* ALLOW_ICEFRONT */

      RETURN
      END
1	dimitri	1.3	C $Header: /u/gcmpack/MITgcm/pkg/icefront/icefront_check.F,v 1.2 2010/01/25 22:37:19 dimitri Exp $
2	dimitri	1.1	C $Name: $
3			#include "ICEFRONT_OPTIONS.h"
4
5			SUBROUTINE ICEFRONT_CHECK( myThid )
6			C /==========================================================\
7			C \| SUBROUTINE ICEFRONT_CHECK \|
8			C \| o Validate basic package setup and inter-package \|
9			C \| dependencies. \|
10			C \==========================================================/
11			IMPLICIT NONE
12
13			C === Global variables ===
14			#include "SIZE.h"
15			#include "EEPARAMS.h"
16			#include "PARAMS.h"
17			#include "ICEFRONT.h"
18
19			C === Routine arguments ===
20			C myThid - Number of this instance of ICEFRONT_CHECK
21			INTEGER myThid
22
23			#ifdef ALLOW_ICEFRONT
24
25			C === Local variables ===
26			C msgBuf - Informational/error meesage buffer
27			CHARACTER*(MAX_LEN_MBUF) msgBuf
28
29			_BEGIN_MASTER(myThid)
30
31			WRITE(msgBuf,'(A)') 'ICEFRONT_CHECK: #define ALLOW_ICEFRONT'
32			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
33			& SQUEEZE_RIGHT, myThid )
34
35			C So far, ICEFRONT works only with oceanic z-coordinates
36			IF ( buoyancyRelation .NE. 'OCEANIC' ) THEN
37			WRITE(msgBuf,'(A)')
38			& 'ICEFRONT works only with buoyancyRelation = ''OCEANIC'''
39			CALL PRINT_ERROR( msgBuf, myThid )
40			STOP 'ABNORMAL END: S/R ICEFRONT_CHECK'
41			ENDIF
42
43			CML#ifdef ALLOW_NONHYDROSTATIC
44			CML IF ( nonHydrostatic ) THEN
45			CML WRITE(msgBuf,'(A)')
46			CML & 'Currently ICEFRONT and nonHydrostatic cannot be turned'
47			CML CALL PRINT_ERROR( msgBuf , 1)
48			CML WRITE(msgBuf,'(A)') 'on at the same time'
49			CML CALL PRINT_ERROR( msgBuf , 1)
50			CML STOP 'ABNORMAL END: S/R ICEFRONT_CHECK'
51			CML ENDIF
52			CML#endif /* ALLOW_NONHYDROSTATIC */
53
54			_END_MASTER(myThid)
55
56			#endif /* ALLOW_ICEFRONT */
57
58			RETURN
59			END