/[MITgcm]/MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/salt_plume_check.F
ViewVC logotype

Annotation of /MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/salt_plume_check.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.2 - (hide annotations) (download)
Tue Apr 29 06:49:40 2014 UTC (11 years, 9 months ago) by atn
Branch: MAIN
Changes since 1.1: +4 -4 lines
in progress:
1. add SALT_PLUME_OPTIONS.h to several files;
2. replace SPalpha (vol) with SPbrineSconst (salinity);
3. move diagnostics outside bi,bj loop.

1 atn 1.2 C $Header: /u/gcmpack/MITgcm_contrib/atnguyen/code_21Dec2012_saltplume/salt_plume_check.F,v 1.1 2014/04/20 04:03:07 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     & (SPbrineSconst .GT. 50. _d 0) ) then
56 atn 1.1 WRITE(msgBuf,'(2A)') 'SALT_PLUME_CHECK:',
57 atn 1.2 & ' SPbrineSconst valid range: [0.,50.]'
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

  ViewVC Help
Powered by ViewVC 1.1.22