/[MITgcm]/MITgcm/model/src/calc_ivdc.F
ViewVC logotype

Contents of /MITgcm/model/src/calc_ivdc.F

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

Revision 1.7 - (show annotations) (download)
Wed Jun 22 00:25:32 2005 UTC (19 years ago) by jmc
Branch: MAIN
CVS Tags: checkpoint62v, checkpoint57m_post, checkpoint62u, checkpoint62t, checkpoint57s_post, checkpoint58b_post, checkpoint57y_post, checkpoint62c, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint62s, checkpoint58a_post, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint57z_post, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62w, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint58y_post, checkpoint58t_post, checkpoint58m_post, checkpoint57l_post, checkpoint57t_post, checkpoint63g, checkpoint57v_post, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint58w_post, checkpoint57y_pre, checkpoint58o_post, checkpoint58p_post, checkpoint58q_post, checkpoint58e_post, mitgcm_mapl_00, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint58r_post, checkpoint58n_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, checkpoint58k_post, checkpoint62b, checkpoint58v_post, checkpoint58l_post, checkpoint57j_post, checkpoint61f, checkpoint58g_post, checkpoint58x_post, checkpoint61n, checkpoint59j, checkpoint58h_post, checkpoint58j_post, checkpoint57o_post, checkpoint61q, checkpoint57k_post, checkpoint57w_post, checkpoint61e, checkpoint58i_post, checkpoint57x_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post, checkpoint61g, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.6: +2 -2 lines 
"rkSign" replaces "-rkFac" (<- removed)
1 C $Header: /u/gcmpack/MITgcm/model/src/calc_ivdc.F,v 1.6 2004/07/01 21:45:12 jmc Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5 CBOP
6 C !ROUTINE: CALC_IVDC
7 C !INTERFACE:
8 SUBROUTINE CALC_IVDC(
9 I bi, bj, iMin, iMax, jMin, jMax, K,
10 I rhoKm1, rhoKp1,
11 I myTime,myIter, myThid )
12 C !DESCRIPTION: \bv
13 C *==========================================================*
14 C | SUBROUTINE CALC_IVDC
15 C | o Calculates Implicit Vertical Diffusivity for Convection
16 C \==========================================================*
17 C \ev
18 C !USES:
19 IMPLICIT NONE
20 C == Global data ==
21 #include "SIZE.h"
22 #include "EEPARAMS.h"
23 #include "PARAMS.h"
24 #include "DYNVARS.h"
25 #include "GRID.h"
26
27 C !INPUT/OUTPUT PARAMETERS:
28 C == Routine arguments ==
29 C bi,bj,iMin,iMax,jMin,jMax,K :: Loop counters
30 C rhoKm1 :: rho in layer above
31 C rhoKp1 :: rho in layer below
32 C myTime :: Current time in simulation
33 C myIter :: Current iteration in simulation
34 C myThid :: Thread number of this instance of S/R CONVECT
35 INTEGER bi,bj,iMin,iMax,jMin,jMax,K
36 _RL rhoKm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
37 _RL rhoKp1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
38 _RL myTime
39 INTEGER myIter
40 INTEGER myThid
41
42 C !LOCAL VARIABLES:
43 C == Local variables ==
44 C i,j :: Loop counters
45 INTEGER i,j
46 CEOP
47
48 C-- Where statically unstable, mix the heat and salt
49 DO j=jMin,jmax
50 DO i=iMin,imax
51 IF ( _hFacC(i,j,k,bi,bj) .GT. 0. .AND.
52 & (rhokm1(i,j)-rhokp1(i,j))*rkSign*gravitySign.GT.0. ) THEN
53 IVDConvCount(i,j,k,bi,bj) = 1. _d 0
54 ELSE
55 IVDConvCount(i,j,k,bi,bj) = 0. _d 0
56 ENDIF
57 ENDDO
58 ENDDO
59
60 RETURN
61 END
  ViewVC Help
Powered by ViewVC 1.1.22