pkg/layers/layers_check.F

C $Header: /u/gcmpack/MITgcm/pkg/layers/layers_check.F,v 1.3 2009/12/28 02:33:49 jmc Exp $
C $Name:  $

#include "LAYERS_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

      SUBROUTINE LAYERS_CHECK( myThid )

C     Check dependances with other packages

      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "EOS.h"
#include "LAYERS_SIZE.h"
#include "LAYERS.h"

C     myThid   :: my Thread Id number
      INTEGER myThid

C     LOCAL VARIABLES:
C     msgBuf   :: Informational/error message buffer
      CHARACTER*(MAX_LEN_MBUF) msgBuf

#ifdef ALLOW_LAYERS

      _BEGIN_MASTER(myThid)

       WRITE(msgBuf,'(A)') 'LAYERS_CHECK: #define LAYERS'
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                     SQUEEZE_RIGHT , 1)

C--  Print out some key parameters :
       CALL WRITE_0D_I( NZZ, INDEX_NONE, 'NZZ =',
     &  ' /* number of levels in the fine vertical grid */')

       CALL WRITE_1D_RL( dZZf, NZZ, INDEX_K, 'dZZf =',
     &  ' /* fine vertical grid spacing for isopycnal interp */')

       CALL WRITE_1D_RL(layers_G,Nlayers+1, INDEX_K,'layers_G =',
     &  ' /* boundaries of isopycnal-averaging bins */')

C--  Check parameters:


C     For now the package will only work if density ~ temperature
      IF ( (eosType .EQ. 'LINEAR')
     &  .AND. (sBeta .EQ. 0.0 _d 0) ) THEN
C        we are good
      ELSE
        WRITE(msgBuf,'(2A)') 'eosType must be eosType=''LINEAR''',
     &   ' and sBeta must = 0.0'
        CALL PRINT_ERROR( msgBuf , 1)
       STOP 'ABNORMAL END: S/R LAYERS_CHECK'
      ENDIF

      _END_MASTER(myThid)

#endif /* ALLOW_LAYERS */

      RETURN
      END
1	jmc	1.4	C $Header: /u/gcmpack/MITgcm/pkg/layers/layers_check.F,v 1.3 2009/12/28 02:33:49 jmc Exp $
2	rpa	1.1	C $Name: $
3
4			#include "LAYERS_OPTIONS.h"
5
6			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7
8			SUBROUTINE LAYERS_CHECK( myThid )
9
10			C Check dependances with other packages
11
12			IMPLICIT NONE
13			#include "SIZE.h"
14			#include "EEPARAMS.h"
15			#include "PARAMS.h"
16			#include "EOS.h"
17			#include "LAYERS_SIZE.h"
18			#include "LAYERS.h"
19
20			C myThid :: my Thread Id number
21			INTEGER myThid
22
23			C LOCAL VARIABLES:
24	jmc	1.3	C msgBuf :: Informational/error message buffer
25	rpa	1.1	CHARACTER*(MAX_LEN_MBUF) msgBuf
26
27			#ifdef ALLOW_LAYERS
28
29			_BEGIN_MASTER(myThid)
30
31			WRITE(msgBuf,'(A)') 'LAYERS_CHECK: #define LAYERS'
32			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
33			& SQUEEZE_RIGHT , 1)
34
35			C-- Print out some key parameters :
36			CALL WRITE_0D_I( NZZ, INDEX_NONE, 'NZZ =',
37			& ' /* number of levels in the fine vertical grid */')
38
39	rpa	1.2	CALL WRITE_1D_RL( dZZf, NZZ, INDEX_K, 'dZZf =',
40			& ' /* fine vertical grid spacing for isopycnal interp */')
41
42	rpa	1.1	CALL WRITE_1D_RL(layers_G,Nlayers+1, INDEX_K,'layers_G =',
43			& ' /* boundaries of isopycnal-averaging bins */')
44
45			C-- Check parameters:
46
47
48			C For now the package will only work if density ~ temperature
49	jmc	1.3	IF ( (eosType .EQ. 'LINEAR')
50	rpa	1.1	& .AND. (sBeta .EQ. 0.0 _d 0) ) THEN
51	jmc	1.4	C we are good
52	jmc	1.3	ELSE
53	rpa	1.1	WRITE(msgBuf,'(2A)') 'eosType must be eosType=''LINEAR''',
54			& ' and sBeta must = 0.0'
55			CALL PRINT_ERROR( msgBuf , 1)
56			STOP 'ABNORMAL END: S/R LAYERS_CHECK'
57	jmc	1.3	ENDIF
58	rpa	1.1
59			_END_MASTER(myThid)
60
61	jmc	1.3	#endif /* ALLOW_LAYERS */
62	rpa	1.1
63			RETURN
64			END