/[MITgcm]/MITgcm/pkg/kpp/kpp_calc_diff.F
ViewVC logotype

Contents of /MITgcm/pkg/kpp/kpp_calc_diff.F

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

Revision 1.8 - (show annotations) (download)
Mon Aug 13 18:05:26 2001 UTC (22 years, 9 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint48f_post, checkpoint51k_post, checkpoint53f_post, checkpoint47j_post, checkpoint54a_pre, checkpoint55c_post, checkpoint53b_pre, checkpoint48d_pre, checkpoint44b_post, checkpoint51j_post, branch-exfmods-tag, checkpoint47e_post, checkpoint43a-release1mods, checkpoint44h_pre, release1_p12, checkpoint52l_pre, checkpoint48i_post, checkpoint52e_pre, hrcube4, hrcube5, release1_p10, release1_p11, release1_p16, release1_p17, checkpoint52j_post, release1_p15, ecco_c44_e17, checkpoint47f_post, ecco_c44_e16, checkpoint54b_post, checkpoint46j_post, checkpoint47c_post, checkpoint50e_post, checkpoint52e_post, checkpoint50c_post, checkpoint46i_post, checkpoint51n_pre, checkpoint47d_post, ecco_c44_e21, ecco_c44_e26, ecco_c44_e27, ecco_c44_e24, ecco_c44_e25, icebear5, icebear4, checkpoint44f_pre, checkpoint47a_post, icebear3, icebear2, checkpoint46f_post, checkpoint52d_pre, ecco_c50_e33a, checkpoint53c_post, checkpoint53d_post, checkpoint41, checkpoint55d_pre, checkpoint46d_pre, release1_p13_pre, checkpoint48e_post, checkpoint46e_post, checkpoint48h_post, checkpoint48d_post, checkpoint50c_pre, release1-branch_tutorials, checkpoint46c_post, release1_p14, checkpoint44g_post, branchpoint-genmake2, checkpoint46h_pre, checkpoint44h_post, release1_p12_pre, checkpoint44e_post, checkpoint46k_post, checkpoint52j_pre, checkpoint54a_post, checkpoint46e_pre, ecco-branch-mod4, ecco-branch-mod5, branch-netcdf, checkpoint50d_pre, checkpoint45d_post, checkpoint51r_post, checkpoint47i_post, checkpoint52b_pre, checkpoint52n_post, checkpoint51o_pre, checkpoint46j_pre, checkpoint45b_post, checkpoint46b_pre, checkpoint51i_post, checkpoint54d_post, checkpoint47h_post, checkpoint48c_post, checkpoint46l_post, chkpt44a_pre, release1-branch-end, checkpoint54e_post, release1_final_v1, checkpoint55b_post, ecco_c44_e19, checkpoint51e_post, checkpoint51b_post, checkpoint46, ecco_c50_e28, ecco_c44_e20, checkpoint51l_pre, checkpoint52m_post, checkpoint51c_post, ecco_c50_e32, ecco_c50_e33, checkpoint53a_post, ecco_c50_e31, checkpoint55a_post, checkpoint47d_pre, checkpoint44, release1_p13, ecco_c44_e18, checkpoint48, checkpoint49, checkpoint44f_post, checkpoint47b_post, checkpoint51l_post, checkpoint40pre6, checkpoint55g_post, checkpoint51o_post, checkpoint40pre8, checkpoint46l_pre, checkpoint51f_pre, release1_b1, checkpoint55f_post, chkpt44c_post, checkpoint51q_post, checkpoint52l_post, checkpoint52k_post, ecco_c51_e34, chkpt44d_post, ecco_c50_e29, checkpoint42, release1_p9, checkpoint55, checkpoint54, checkpoint53b_post, checkpoint51, checkpoint50, checkpoint53, checkpoint52, release1_p8, checkpoint50d_post, checkpoint52d_post, checkpoint46g_pre, release1_p2, release1_p3, release1_p4, checkpoint51b_pre, release1_p6, checkpoint52a_post, checkpoint46a_post, checkpoint47g_post, checkpoint52b_post, checkpoint53g_post, chkpt44a_post, checkpoint52f_post, checkpoint44b_pre, checkpoint52c_post, release1_p1, checkpoint46m_post, checkpoint48a_post, checkpoint51h_pre, checkpoint46a_pre, ecco_c51_e34e, ecco-branch-mod1, checkpoint50g_post, checkpoint45c_post, checkpoint50b_pre, release1_p5, checkpoint44e_pre, checkpoint51g_post, ecco_c52_e35, checkpoint40pre9, release1_p7, ecco_ice2, ecco_ice1, checkpoint46b_post, checkpoint54f_post, checkpoint51f_post, checkpoint46d_post, ecco-branch-mod2, checkpoint48b_post, checkpoint50b_post, checkpoint46g_post, checkpoint45a_post, ecco_c51_e34d, ecco_c51_e34f, ecco_c51_e34g, ecco_c51_e34a, ecco_c51_e34b, ecco_c51_e34c, checkpoint50f_post, checkpoint50a_post, checkpoint46c_pre, checkpoint50f_pre, checkpoint52a_pre, ecco-branch-mod3, checkpoint43, checkpoint51d_post, ecco_c50_e30, checkpoint48c_pre, release1-branch_branchpoint, ecco_c44_e22, release1_beta1, ecco_c44_e23, checkpoint51m_post, checkpoint51t_post, checkpoint53d_pre, checkpoint47, checkpoint55e_post, checkpoint40pre7, checkpoint54c_post, checkpoint50h_post, checkpoint52i_post, checkpoint51a_post, checkpoint40, checkpoint45, checkpoint46h_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, release1_chkpt44d_post, checkpoint48g_post, checkpoint51i_pre, chkpt44c_pre, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint55d_post
Branch point for: c24_e25_ice, netcdf-sm0, ecco-branch, branch-genmake2, release1, branch-exfmods-curt, release1_coupled, branch-nonh, icebear, tg2-branch, release1_final, checkpoint51n_branch, release1-branch, release1_50yr
Changes since 1.7: +16 -23 lines 
Modifications related to split into thermodynamics.F, dynamics.F
o missing initialisations in dynamics.F added
o some fields no longer needed in dynamics/thermodynamics deleted
o split of calc_diffusivity.F into calc_viscosity.F
  (plus split of kpp_calc_diff.F into kpp_calc_visc.F)
o Modifications of some store directives for TAF
1 C $Header: /u/gcmpack/models/MITgcmUV/pkg/kpp/kpp_calc_diff.F,v 1.7 2001/02/04 14:38:50 cnh Exp $
2 C $Name: checkpoint40pre4 $
3
4 #include "KPP_OPTIONS.h"
5
6 SUBROUTINE KPP_CALC_DIFF(
7 I bi,bj,iMin,iMax,jMin,jMax,K,
8 I maskUp,
9 U KappaRT,KappaRS,
10 I myThid)
11
12 C /==========================================================\
13 C | SUBROUTINE KPP_CALC_DIFF |
14 C | o Add contrubution to net diffusivity from KPP mixing |
15 C \==========================================================/
16 IMPLICIT NONE
17
18 C == GLobal variables ==
19 #include "SIZE.h"
20 #include "EEPARAMS.h"
21 #include "PARAMS.h"
22 #include "DYNVARS.h"
23 #include "GRID.h"
24
25 C == Routine arguments ==
26 C bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
27 C maskUp - land/water mask for Wvel points (above tracer level)
28 C myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
29 C
30 INTEGER bi,bj,iMin,iMax,jMin,jMax,K
31 _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
32 _RL KappaRT(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
33 _RL KappaRS(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
34 INTEGER myThid
35
36 #ifdef ALLOW_KPP
37 #include "KPP.h"
38
39 C == Local variables ==
40 C I, J, K - Loop counters
41 INTEGER i,j
42
43 DO j=jMin,jMax
44 DO i=iMin,iMax
45 KappaRT(i,j,k) = KappaRT(i,j,k)+maskUp(i,j)*
46 & ( KPPdiffKzT(i,j,k,bi,bj)
47 #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
48 & - diffKr(i,j,k,bi,bj) )
49 #else
50 & - diffKrT )
51 #endif
52 ENDDO
53 ENDDO
54
55 DO j=jMin,jMax
56 DO i=iMin,iMax
57 KappaRS(i,j,k) = KappaRS(i,j,k)+maskUp(i,j)*
58 & ( KPPdiffKzS(i,j,k,bi,bj)
59 #if (defined (ALLOW_AUTODIFF_TAMC) && defined (ALLOW_DIFFKR_CONTROL))
60 & - diffKr(i,j,k,bi,bj) )
61 #else
62 & - diffKrS )
63 #endif
64 ENDDO
65 ENDDO
66
67 #endif /* ALLOW_KPP */
68
69 RETURN
70 END
  ViewVC Help
Powered by ViewVC 1.1.22