pkg/monitor/mon_solution.F

C $Header: /u/gcmpack/MITgcm/pkg/monitor/mon_solution.F,v 1.2 2003/12/15 02:13:13 jmc Exp $
C $Name:  $

#include "MONITOR_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP
C     !ROUTINE: MON_SOLUTION

C     !INTERFACE:
      SUBROUTINE MON_SOLUTION(
     I     myIter, myTime, myThid )

C     !DESCRIPTION:
C     Checks that the solutions is within bounds             

C     !USES:
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "GRID.h"
#include "MONITOR.h"

C     !INPUT PARAMETERS:
      INTEGER myIter
      _RL myTime
      INTEGER myThid
CEOP

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