rpa_layers/layers/layers_check.F

#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 meesage 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_RL( dZZ, INDEX_NONE, 'dZZ =',
     &  ' /* fine vertical grid spacing for isopycnal interp */')

       CALL WRITE_0D_I( NZZ, INDEX_NONE, 'NZZ =',
     &  ' /* number of levels in the fine vertical grid */')
       
       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're 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_MYPACKAGE */

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