natl_12/code/mon_advcfl.F

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

#include "MONITOR_OPTIONS.h"

      SUBROUTINE MON_ADVCFL(
     I                label, U, rDx, dT,
     I                myThid )
C     /==========================================================\
C     | SUBROUTINE MON_ADVCFL                                    |
C     | o Calculates maximum CFL number                          |
C     |==========================================================|
C     \==========================================================/
      IMPLICIT NONE

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

C     === Routine arguments ===
      CHARACTER*(*) label
      _RL U(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RS rDx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL dT
      INTEGER myThid

C     === Local variables ====
      INTEGER bi,bj,I,J,K
      INTEGER maxI, maxJ, maxK, IG, JG, maxBi, maxBj
      _RL tmpVal,theMax, myMax

      theMax=0.
      maxI  = 0
      maxJ  = 0
      maxK  = 0
      maxBi = 0
      maxBj = 0

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        DO K=1,Nr
         DO J=1,sNy
          DO I=1,sNx
           tmpVal=abs(U(I,J,K,bi,bj))*rDx(I,J,BI,BJ)*dT
           IF ( tmpVal .GT. theMax ) THEN
            maxI  = I
            maxJ  = J
            maxK  = K
            maxBi = bi
            maxBj = bj
           ENDIF
           theMax=max(theMax,tmpVal)
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDDO
      myMax = theMax
      _GLOBAL_MAX_R8(theMax,myThid)
      IF ( theMax .EQ. myMax ) THEN
       IG = myXGlobalLo+maxI-1+(maxBi-1)*sNx
       JG = myYGlobalLo+maxJ-1+(maxBj-1)*sNy
       CALL MON_OUT_I( '_loci',   IG, mon_foot_max, myThid )
       CALL MON_OUT_I( '_locj',   JG, mon_foot_max, myThid )
       CALL MON_OUT_I( '_lock', maxK, mon_foot_max, myThid )
      ENDIF

      CALL MON_OUT_RL( label, theMax, mon_foot_max ,myThid)

      RETURN
      END
1	C $Header: /u/u0/gcmpack/MITgcm/pkg/monitor/mon_advcfl.F,v 1.5 2003/05/13 18:18:05 adcroft Exp $
2	C $Name: $
3
4	#include "MONITOR_OPTIONS.h"
5
6	SUBROUTINE MON_ADVCFL(
7	I label, U, rDx, dT,
8	I myThid )
9	C /==========================================================\
10	C \| SUBROUTINE MON_ADVCFL \|
11	C \| o Calculates maximum CFL number \|
12	C \|==========================================================\|
13	C \==========================================================/
14	IMPLICIT NONE
15
16	C === Global data ===
17	#include "SIZE.h"
18	#include "EEPARAMS.h"
19	#include "MONITOR.h"
20
21	C === Routine arguments ===
22	CHARACTER() label
23	_RL U(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
24	_RS rDx(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
25	_RL dT
26	INTEGER myThid
27
28	C === Local variables ====
29	INTEGER bi,bj,I,J,K
30	INTEGER maxI, maxJ, maxK, IG, JG, maxBi, maxBj
31	_RL tmpVal,theMax, myMax
32
33	theMax=0.
34	maxI = 0
35	maxJ = 0
36	maxK = 0
37	maxBi = 0
38	maxBj = 0
39
40	DO bj=myByLo(myThid),myByHi(myThid)
41	DO bi=myBxLo(myThid),myBxHi(myThid)
42	DO K=1,Nr
43	DO J=1,sNy
44	DO I=1,sNx
45	tmpVal=abs(U(I,J,K,bi,bj))rDx(I,J,BI,BJ)dT
46	IF ( tmpVal .GT. theMax ) THEN
47	maxI = I
48	maxJ = J
49	maxK = K
50	maxBi = bi
51	maxBj = bj
52	ENDIF
53	theMax=max(theMax,tmpVal)
54	ENDDO
55	ENDDO
56	ENDDO
57	ENDDO
58	ENDDO
59	myMax = theMax
60	_GLOBAL_MAX_R8(theMax,myThid)
61	IF ( theMax .EQ. myMax ) THEN
62	IG = myXGlobalLo+maxI-1+(maxBi-1)*sNx
63	JG = myYGlobalLo+maxJ-1+(maxBj-1)*sNy
64	CALL MON_OUT_I( '_loci', IG, mon_foot_max, myThid )
65	CALL MON_OUT_I( '_locj', JG, mon_foot_max, myThid )
66	CALL MON_OUT_I( '_lock', maxK, mon_foot_max, myThid )
67	ENDIF
68
69	CALL MON_OUT_RL( label, theMax, mon_foot_max ,myThid)
70
71	RETURN
72	END