/[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.3 - (show annotations) (download)
Wed Sep 26 18:09:14 2001 UTC (22 years, 7 months ago) by cnh
Branch: MAIN
CVS Tags: checkpoint44e_post, release1_p13_pre, checkpoint44f_post, checkpoint46b_post, checkpoint43a-release1mods, release1_p13, chkpt44d_post, release1_p8, release1_p9, release1_p1, release1_p2, release1_p3, release1_p4, release1_p5, release1_p6, release1_p7, checkpoint44e_pre, release1_b1, checkpoint43, release1_chkpt44d_post, release1_p11, icebear5, icebear4, icebear3, icebear2, release1-branch_tutorials, checkpoint45d_post, chkpt44a_post, checkpoint44h_pre, checkpoint46a_post, chkpt44c_pre, checkpoint45a_post, ecco_c44_e19, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, release1_p12, release1_p10, release1_p16, release1_p17, release1_p14, release1_p15, checkpoint44g_post, checkpoint45b_post, checkpoint46b_pre, release1-branch-end, release1_final_v1, checkpoint46c_pre, checkpoint46, checkpoint44b_post, checkpoint46a_pre, checkpoint45c_post, ecco_ice2, ecco_ice1, checkpoint44h_post, release1_p12_pre, ecco_c44_e22, ecco_c44_e25, chkpt44a_pre, ecco_c44_e23, ecco_c44_e20, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, ecco-branch-mod1, ecco-branch-mod2, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5, release1_beta1, checkpoint44b_pre, checkpoint42, checkpoint41, checkpoint44, checkpoint45, chkpt44c_post, checkpoint44f_pre, release1-branch_branchpoint
Branch point for: c24_e25_ice, release1_final, release1-branch, release1, ecco-branch, release1_50yr, icebear, release1_coupled
Changes since 1.2: +25 -15 lines 
Bringing comments up to data and formatting for document extraction.
1 C $Header: /u/gcmpack/models/MITgcmUV/model/src/calc_ivdc.F,v 1.2 2001/02/04 14:38:46 cnh Exp $
2 C $Name: $
3
4 #include "CPP_OPTIONS.h"
5 CBOP
6 C !ROUTINE: CALC_IVDC
7 C !INTERFACE:
8 SUBROUTINE CALC_IVDC( bi, bj, iMin, iMax, jMin, jMax, K,
9 & rhoKm1, rhoKp1,
10 & ConvectCount, KappaRT, KappaRS,
11 & 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 EXTERNAL DIFFERENT_MULTIPLE
27 LOGICAL DIFFERENT_MULTIPLE
28
29 C !INPUT/OUTPUT PARAMETERS:
30 C == Routine arguments ==
31 C bi,bj,iMin,iMax,jMin,jMax,K :: Loop counters
32 C rhoKm1 :: rho in layer above
33 C rhoKp1 :: rho in layer below
34 C myTime :: Current time in simulation
35 C myIter :: Current iteration in simulation
36 C myThid :: Thread number of this instance of S/R CONVECT
37 INTEGER bi,bj,iMin,iMax,jMin,jMax,K
38 C KappaRT :: Vertical diffusion coeff. for theta
39 C KappaRS :: Vertical diffusion coeff. for salt
40 C ConvectCount :: Static instability counter
41 _RL rhoKm1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
42 _RL rhoKp1(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
43 _RL ConvectCount(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
44 _RL KappaRT(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
45 _RL KappaRS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
46 _RL myTime
47 INTEGER myIter
48 INTEGER myThid
49
50 C !LOCAL VARIABLES:
51 C == Local variables ==
52 C i,j :: Loop counters
53 INTEGER i,j
54 CEOP
55
56 C-- Where statically unstable, mix the heat and salt
57 DO j=jMin,jmax
58 DO i=iMin,imax
59 IF ( _hFacC(i,j,k,bi,bj) .GT. 0. .AND.
60 & rhokm1(i,j) .GT. rhokp1(i,j) ) THEN
61 KappaRT(i,j,k) = ivdc_kappa
62 KappaRS(i,j,k) = ivdc_kappa
63 ConvectCount(i,j,k) = ConvectCount(i,j,k) + 1.
64 ENDIF
65 ENDDO
66 ENDDO
67
68 RETURN
69 END
  ViewVC Help
Powered by ViewVC 1.1.22