/[MITgcm]/MITgcm/pkg/monitor/mon_advcflw2.F
ViewVC logotype

Contents of /MITgcm/pkg/monitor/mon_advcflw2.F

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

Revision 1.2 - (show annotations) (download)
Tue May 13 18:18:05 2003 UTC (21 years ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint52d_pre, checkpoint52j_pre, checkpoint51o_pre, checkpoint51l_post, checkpoint52l_post, checkpoint52k_post, checkpoint51, checkpoint52, checkpoint52f_post, checkpoint51f_post, checkpoint51d_post, checkpoint51t_post, checkpoint51n_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint51j_post, checkpoint52e_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint52b_pre, checkpoint51l_pre, checkpoint51q_post, checkpoint51b_pre, checkpoint52b_post, checkpoint52c_post, checkpoint51h_pre, checkpoint50f_post, checkpoint50f_pre, checkpoint52f_pre, branchpoint-genmake2, checkpoint51r_post, checkpoint51i_post, checkpoint51b_post, checkpoint51c_post, checkpoint52d_post, checkpoint50g_post, checkpoint52a_pre, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint52j_post, checkpoint50e_post, branch-netcdf, checkpoint51e_post, checkpoint51o_post, checkpoint51f_pre, checkpoint52a_post, checkpoint51g_post, ecco_c52_e35, checkpoint51m_post, checkpoint51a_post, checkpoint51p_post, checkpoint51u_post
Branch point for: branch-genmake2, branch-nonh, tg2-branch, netcdf-sm0, checkpoint51n_branch
Changes since 1.1: +2 -2 lines 
 o split mon_set.F into mon_set_iounit.F and mon_set_pref.F
 o replaced ref's to CPP_OPTIONS with MONITOR_OPTIONS
 o added new s/r monitor_solution.F that checks that model state
   and if unlikely lets the model die cleanly
   - this is to reduce the number of hanging processes we encounter
     if the model dies due to FPEs
1 C $Header: /u/gcmpack/models/MITgcmUV/pkg/monitor/mon_advcflw2.F,v 1.1 2001/07/03 21:34:52 jmc Exp $
2 C $Name: $
3
4 #include "MONITOR_OPTIONS.h"
5
6 SUBROUTINE MON_ADVCFLW2(
7 I label, W, rHFac, rDrF, dT,
8 I myThid )
9 C /==========================================================\
10 C | SUBROUTINE MON_ADVCFLW2 |
11 C | o Calculates maximum CFL number in vertical |
12 C | relevant for tracer Adv. Pb. with Partial Cell |
13 C |==========================================================|
14 C \==========================================================/
15 IMPLICIT NONE
16
17 C === Global data ===
18 #include "SIZE.h"
19 #include "EEPARAMS.h"
20 #include "MONITOR.h"
21
22 C === Routine arguments ===
23 CHARACTER*(*) label
24 _RL W(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
25 _RS rHFac(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
26 _RS rDrF(Nr)
27 _RL dT
28 INTEGER myThid
29
30 C === Local variables ====
31 INTEGER bi,bj,I,J,K
32 _RL tmpVal,theMax
33
34 theMax=0.
35
36 DO bj=myByLo(myThid),myByHi(myThid)
37 DO bi=myBxLo(myThid),myBxHi(myThid)
38 DO K=2,Nr
39 DO J=1,sNy
40 DO I=1,sNx
41 tmpVal=abs(W(I,J,K,bi,bj))*dT
42 & *max( rDrF(K)*rHfac(I,J,K,bi,bj),
43 & rDrF(K-1)*rHfac(I,J,K-1,bi,bj) )
44 theMax=max(theMax,tmpVal)
45 ENDDO
46 ENDDO
47 ENDDO
48 ENDDO
49 ENDDO
50 _GLOBAL_MAX_R8(theMax,myThid)
51
52 CALL MON_OUT_RL( label, theMax, mon_foot_max ,myThid)
53
54 RETURN
55 END
  ViewVC Help
Powered by ViewVC 1.1.22