/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_check.F
ViewVC logotype

Contents of /MITgcm/pkg/diagnostics/diagnostics_check.F

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


Revision 1.5 - (show annotations) (download)
Sun Jun 26 16:51:49 2005 UTC (18 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint57o_post, checkpoint57m_post, checkpoint57s_post, checkpoint57k_post, checkpoint57v_post, checkpoint57r_post, checkpoint57n_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, checkpoint57j_post, checkpoint57l_post
Changes since 1.4: +2 -63 lines
change pointers so that 1 diag. can be used several times (with # freq.)

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

  ViewVC Help
Powered by ViewVC 1.1.22