1 |
C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_check.F,v 1.4 2005/05/20 07:28:49 jmc Exp $ |
2 |
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 |
#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 |
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 |
DO n = 2,diagSt_nbLists |
90 |
DO m = 1,n-1 |
91 |
IF ( diagSt_Fname(n).EQ.diagSt_Fname(m) ) THEN |
92 |
WRITE(msgBuf,'(2A)') 'DIAGNOSTICS_CHECK: ', |
93 |
& 'found 2 identical stat_fname:' |
94 |
CALL PRINT_ERROR( msgBuf , myThid ) |
95 |
WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_CHECK: ', |
96 |
& '1rst (m=', m, ' ): ', diagSt_Fname(m) |
97 |
CALL PRINT_ERROR( msgBuf , myThid ) |
98 |
WRITE(msgBuf,'(2A,I3,2A)') 'DIAGNOSTICS_CHECK: ', |
99 |
& ' 2nd (n=', n, ' ): ', diagSt_Fname(n) |
100 |
CALL PRINT_ERROR( msgBuf , myThid ) |
101 |
STOP 'ABNORMAL END: S/R DIAGNOSTICS_CHECK' |
102 |
ENDIF |
103 |
ENDDO |
104 |
ENDDO |
105 |
|
106 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
107 |
C- Check for field that appears 2 times (or more) with differents frequency: |
108 |
C disable this checking since now diagnostics pkg can handle this case. |
109 |
|
110 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |
111 |
|
112 |
_END_MASTER(myThid) |
113 |
|
114 |
RETURN |
115 |
END |