pkg/monitor/mon_solution.F

C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_solution.F,v 1.1 2003/05/13 18:18:05 adcroft Exp $
C $Name:  $

#include "MONITOR_OPTIONS.h"

      SUBROUTINE MON_SOLUTION(
     I                myIter, myTime, myThid )
C     /==========================================================\
C     | SUBROUTINE MON_SOLUTION                                  |
C     | o Checks that the solutions is within bounds             |
C     |==========================================================|
C     \==========================================================/
      IMPLICIT NONE

C     === Global data ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "GRID.h"
#include "MONITOR.h"

C     === Routine arguments ===
      INTEGER myIter
      _RL myTime
      INTEGER myThid

C     === Functions ====

C     === Local variables ====
      CHARACTER*(MAX_LEN_MBUF) msgBuf   
      _RL tMin,tMax,tMean,tSD,tDel2,tVol

C     Statistics for T
      CALL MON_STATS_RL(
     I                Nr, theta, maskC,hFacC,rA,drF,
     O                tMin,tMax,tMean,tSD,tDel2,tVol,
     I                myThid )
      
      IF (tMax-tMin.GT.1.e3) THEN
        _BEGIN_MASTER(myThid)
         WRITE(msgBuf,'(A,1P2E12.3)')
     &    'SOLUTION IS HEADING OUT OF BOUNDS: tMin,tMax=',tMin,tMax
         CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
         WRITE(msgBuf,'(A)') 'STOPPING CALCULATION'
         CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
        _END_MASTER(myThid)

C jmc: add this "if not ..." to avoid beeing stuck when using coupler; 
       IF ( .NOT. useCoupler )
     &  CALL EEDIE
        STOP 'MON_SOLUTION: STOPPED DUE TO EXTREME VALUES OF SOLUTION'
      ENDIF

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_solution.F,v 1.1 2003/05/13 18:18:05 adcroft Exp $
2	C $Name: $
3
4	#include "MONITOR_OPTIONS.h"
5
6	SUBROUTINE MON_SOLUTION(
7	I myIter, myTime, myThid )
8	C /==========================================================\
9	C \| SUBROUTINE MON_SOLUTION \|
10	C \| o Checks that the solutions is within bounds \|
11	C \|==========================================================\|
12	C \==========================================================/
13	IMPLICIT NONE
14
15	C === Global data ===
16	#include "SIZE.h"
17	#include "EEPARAMS.h"
18	#include "PARAMS.h"
19	#include "DYNVARS.h"
20	#include "GRID.h"
21	#include "MONITOR.h"
22
23	C === Routine arguments ===
24	INTEGER myIter
25	_RL myTime
26	INTEGER myThid
27
28	C === Functions ====
29
30	C === Local variables ====
31	CHARACTER*(MAX_LEN_MBUF) msgBuf
32	_RL tMin,tMax,tMean,tSD,tDel2,tVol
33
34	C Statistics for T
35	CALL MON_STATS_RL(
36	I Nr, theta, maskC,hFacC,rA,drF,
37	O tMin,tMax,tMean,tSD,tDel2,tVol,
38	I myThid )
39
40	IF (tMax-tMin.GT.1.e3) THEN
41	_BEGIN_MASTER(myThid)
42	WRITE(msgBuf,'(A,1P2E12.3)')
43	& 'SOLUTION IS HEADING OUT OF BOUNDS: tMin,tMax=',tMin,tMax
44	CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
45	WRITE(msgBuf,'(A)') 'STOPPING CALCULATION'
46	CALL PRINT_MESSAGE( msgBuf, mon_ioUnit, SQUEEZE_RIGHT , 1)
47	_END_MASTER(myThid)
48
49	C jmc: add this "if not ..." to avoid beeing stuck when using coupler;
50	IF ( .NOT. useCoupler )
51	& CALL EEDIE
52	STOP 'MON_SOLUTION: STOPPED DUE TO EXTREME VALUES OF SOLUTION'
53	ENDIF
54
55	RETURN
56	END