/[MITgcm]/MITgcm_contrib/natl_12/code/mon_advcfl.F
ViewVC logotype

Annotation of /MITgcm_contrib/natl_12/code/mon_advcfl.F

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


Revision 1.1 - (hide annotations) (download)
Tue Aug 5 21:22:44 2003 UTC (21 years, 11 months ago) by cnh
Branch: MAIN
CVS Tags: HEAD
Adding set of files for 1/12 Atlantic configuration

1 cnh 1.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

  ViewVC Help
Powered by ViewVC 1.1.22