1 |
jmc |
1.3 |
C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_check.F,v 1.2 2004/12/16 21:39:53 jmc Exp $ |
2 |
jmc |
1.1 |
C $Name: $ |
3 |
|
|
|
4 |
|
|
#include "DIAG_OPTIONS.h" |
5 |
|
|
|
6 |
|
|
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
7 |
|
|
CBOP 0 |
8 |
|
|
C !ROUTINE: DIAGNOSTICS_CHECK |
9 |
|
|
|
10 |
|
|
C !INTERFACE: |
11 |
|
|
SUBROUTINE DIAGNOSTICS_CHECK(myThid) |
12 |
|
|
|
13 |
|
|
C !DESCRIPTION: |
14 |
|
|
C Check option and parameter consistency |
15 |
|
|
|
16 |
|
|
C !USES: |
17 |
|
|
IMPLICIT NONE |
18 |
|
|
#include "SIZE.h" |
19 |
|
|
#include "EEPARAMS.h" |
20 |
|
|
#include "PARAMS.h" |
21 |
|
|
#include "DIAGNOSTICS_SIZE.h" |
22 |
|
|
#include "DIAGNOSTICS.h" |
23 |
|
|
|
24 |
|
|
C !INPUT PARAMETERS: |
25 |
|
|
INTEGER myThid |
26 |
|
|
CEOP |
27 |
|
|
|
28 |
|
|
C !LOCAL VARIABLES: |
29 |
|
|
CHARACTER*(MAX_LEN_MBUF) msgBuf |
30 |
|
|
INTEGER k,l,n,m |
31 |
|
|
|
32 |
|
|
_BEGIN_MASTER(myThid) |
33 |
|
|
|
34 |
|
|
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
35 |
|
|
C Check diagnostics parameter consistency |
36 |
|
|
|
37 |
jmc |
1.3 |
#ifdef DIAGNOSTICS_HAS_PICKUP |
38 |
|
|
IF ( diag_pickup_read ) THEN |
39 |
|
|
WRITE(msgBuf,'(2A)') '**CAUTION** (DIAGNOSTICS_CHECK): ', |
40 |
|
|
& 'reading diagnostics previous state' |
41 |
|
|
CALL PRINT_MESSAGE( msgBuf, errorMessageUnit, |
42 |
|
|
& SQUEEZE_RIGHT , myThid) |
43 |
|
|
WRITE(msgBuf,'(2A)') '**CAUTION** ', |
44 |
|
|
& ' from a pickup file can only work if data.diagnostics' |
45 |
|
|
CALL PRINT_MESSAGE( msgBuf, errorMessageUnit, |
46 |
|
|
& SQUEEZE_RIGHT , myThid) |
47 |
|
|
WRITE(msgBuf,'(2A)') '**CAUTION** ', |
48 |
|
|
& ' is not changed (<= further checking not yet implemented)' |
49 |
|
|
CALL PRINT_MESSAGE( msgBuf, errorMessageUnit, |
50 |
|
|
& SQUEEZE_RIGHT , myThid) |
51 |
|
|
ENDIF |
52 |
|
|
#else /* undef DIAGNOSTICS_HAS_PICKUP */ |
53 |
|
|
C- stop if trying to use part of the code that is not compiled: |
54 |
|
|
IF ( diag_pickup_read ) THEN |
55 |
|
|
WRITE(msgBuf,'(3A)') 'DIAGNOSTICS_CHECK: ', |
56 |
|
|
& 'diag_pickup_read is TRUE ', |
57 |
|
|
& 'but DIAGNOSTICS_HAS_PICKUP is "#undef"' |
58 |
|
|
CALL PRINT_ERROR( msgBuf , myThid) |
59 |
|
|
ENDIF |
60 |
|
|
IF ( diag_pickup_write ) THEN |
61 |
|
|
WRITE(msgBuf,'(3A)') 'DIAGNOSTICS_CHECK: ', |
62 |
|
|
& 'diag_pickup_write is TRUE ', |
63 |
|
|
& 'but DIAGNOSTICS_HAS_PICKUP is "#undef"' |
64 |
|
|
CALL PRINT_ERROR( msgBuf , myThid) |
65 |
|
|
ENDIF |
66 |
|
|
IF ( diag_pickup_read .OR. diag_pickup_write ) THEN |
67 |
|
|
STOP 'ABNORMAL END: S/R DIAGNOSTICS_CHECK' |
68 |
|
|
ENDIF |
69 |
|
|
#endif /* DIAGNOSTICS_HAS_PICKUP */ |
70 |
|
|
|
71 |
jmc |
1.1 |
C- File names: |
72 |
|
|
DO n = 2,nlists |
73 |
|
|
DO m = 1,n-1 |
74 |
|
|
IF ( fnames(n).EQ.fnames(m) ) THEN |
75 |
|
|
WRITE(msgBuf,'(2A)') 'DIAGNOSTICS_CHECK: ', |
76 |
|
|
& 'found 2 identical filenames:' |
77 |
|
|
CALL PRINT_ERROR( msgBuf , myThid ) |
78 |
|
|
WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_CHECK: ', |
79 |
|
|
& '1rst (m=', m, ' ): ', fnames(m) |
80 |
|
|
CALL PRINT_ERROR( msgBuf , myThid ) |
81 |
|
|
WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_CHECK: ', |
82 |
|
|
& ' 2nd (n=', n, ' ): ', fnames(n) |
83 |
|
|
CALL PRINT_ERROR( msgBuf , myThid ) |
84 |
|
|
STOP 'ABNORMAL END: S/R DIAGNOSTICS_CHECK' |
85 |
|
|
ENDIF |
86 |
|
|
ENDDO |
87 |
|
|
ENDDO |
88 |
|
|
|
89 |
|
|
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
90 |
|
|
C- Check for field that appears 2 times (or more) with differents frequency: |
91 |
|
|
|
92 |
|
|
DO n = 2,nlists |
93 |
|
|
DO m = 1,n-1 |
94 |
|
|
IF ( freq(m).NE.freq(n) ) THEN |
95 |
|
|
DO k = 1,nActive(n) |
96 |
|
|
DO l = 1,nActive(m) |
97 |
|
|
IF ( flds(k,n).EQ.flds(l,m) ) THEN |
98 |
|
|
WRITE(msgBuf,'(4A)') 'DIAGNOSTICS_CHECK: ', |
99 |
|
|
& 'field : ',flds(k,n),' use 2 different freq. :' |
100 |
|
|
CALL PRINT_ERROR( msgBuf , myThid ) |
101 |
|
|
WRITE(msgBuf,'(2A,I3,A,I6,2A)') 'DIAGNOSTICS_CHECK: ', |
102 |
|
|
& '1rst (m=', m, ' ), freq=', freq(m),' file:',fnames(m) |
103 |
|
|
CALL PRINT_ERROR( msgBuf , myThid ) |
104 |
|
|
WRITE(msgBuf,'(2A,I3,A,I6,2A)') 'DIAGNOSTICS_CHECK: ', |
105 |
|
|
& ' 2nd (n=', n, ' ), freq=', freq(n),' file:',fnames(n) |
106 |
|
|
CALL PRINT_ERROR( msgBuf , myThid ) |
107 |
jmc |
1.2 |
STOP 'ABNORMAL END: S/R DIAGNOSTICS_CHECK' |
108 |
jmc |
1.1 |
ENDIF |
109 |
|
|
ENDDO |
110 |
|
|
ENDDO |
111 |
|
|
ENDIF |
112 |
|
|
ENDDO |
113 |
|
|
ENDDO |
114 |
|
|
|
115 |
|
|
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
116 |
|
|
|
117 |
|
|
_END_MASTER(myThid) |
118 |
|
|
|
119 |
|
|
RETURN |
120 |
|
|
END |