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

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

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

Revision 1.9 - (show annotations) (download)
Sat Jul 13 03:12:30 2002 UTC (21 years, 10 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint51k_post, checkpoint47e_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint46l_post, checkpoint46g_pre, checkpoint47c_post, checkpoint50c_post, checkpoint46f_post, checkpoint52d_pre, checkpoint48e_post, checkpoint50c_pre, checkpoint46b_post, checkpoint52j_pre, checkpoint51o_pre, checkpoint51l_post, checkpoint48i_post, checkpoint46l_pre, checkpoint52l_post, checkpoint52k_post, checkpoint54, checkpoint51, checkpoint50, checkpoint53, checkpoint52, checkpoint50d_post, checkpoint52f_post, checkpoint50b_pre, checkpoint51f_post, checkpoint48b_post, checkpoint51d_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint51t_post, checkpoint51n_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint47a_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, checkpoint52e_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint47d_post, checkpoint53d_post, checkpoint46d_pre, checkpoint48d_post, checkpoint48f_post, checkpoint52b_pre, checkpoint54b_post, checkpoint46j_pre, checkpoint51l_pre, checkpoint52m_post, checkpoint48h_post, checkpoint51q_post, checkpoint51b_pre, checkpoint46a_post, checkpoint47g_post, checkpoint52b_post, checkpoint52c_post, checkpoint46j_post, checkpoint51h_pre, checkpoint46k_post, checkpoint48a_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint52f_pre, checkpoint47j_post, checkpoint54a_pre, checkpoint53c_post, branch-exfmods-tag, branchpoint-genmake2, checkpoint54a_post, checkpoint46e_pre, checkpoint51r_post, checkpoint48c_post, checkpoint46b_pre, checkpoint51i_post, checkpoint51b_post, checkpoint51c_post, checkpoint46c_pre, checkpoint53a_post, checkpoint46, checkpoint47b_post, checkpoint46h_pre, checkpoint52d_post, checkpoint53g_post, checkpoint46m_post, checkpoint46a_pre, checkpoint50g_post, checkpoint46g_post, checkpoint52a_pre, checkpoint50h_post, checkpoint52i_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint52h_pre, checkpoint53f_post, checkpoint52j_post, checkpoint47f_post, checkpoint50e_post, checkpoint46i_post, checkpoint46c_post, branch-netcdf, checkpoint50d_pre, checkpoint52n_post, checkpoint53b_pre, checkpoint46e_post, checkpoint51e_post, checkpoint47, checkpoint48, checkpoint49, checkpoint46h_post, checkpoint51o_post, checkpoint51f_pre, checkpoint48g_post, checkpoint53b_post, checkpoint47h_post, checkpoint52a_post, checkpoint51g_post, ecco_c52_e35, checkpoint46d_post, checkpoint50b_post, checkpoint51m_post, checkpoint53d_pre, checkpoint51a_post, checkpoint51p_post, checkpoint51u_post
Branch point for: branch-exfmods-curt, branch-genmake2, branch-nonh, tg2-branch, netcdf-sm0, checkpoint51n_branch
Changes since 1.8: +1 -1 lines 
Merging from release1_p5
o Adjoint-related bug fixes in kpp:
  - kpp_calc: sore of kpphbl avoids recomputation/call to S/R kppmix
  - kpp_routines: store of Rib avoids partial recomputation bug of TAF.
1 C $Header: /u/gcmpack/MITgcm/pkg/kpp/kpp_transport_s.F,v 1.7.4.1 2002/04/08 20:27:13 heimbach Exp $
2 C $Name: $
3
4 #include "KPP_OPTIONS.h"
5
6 subroutine KPP_TRANSPORT_S (
7 I iMin,iMax,jMin,jMax,bi,bj,k,km1,
8 I KappaRS,
9 U df )
10 C /==========================================================\
11 C | o SUBROUTINE KPP_TRANSPORT_S |
12 C | Add non local KPP transport term (ghat) to diffusive |
13 C | salinity flux. |
14 C |==========================================================|
15 C | The nonlocal transport term is nonzero only for scalars |
16 C | in unstable (convective) forcing conditions. |
17 C | Note: KappaRS(i,j,k) is defined at the top of grid cell |
18 C | k while KPPghat(i,j,k,bi,bj) is defined at the |
19 C | bottom of grid cell k. |
20 C | For K=1, KappaRS(i,j,k) = 0 which insures no flux |
21 C | through the surface. |
22 C \==========================================================/
23 IMPLICIT NONE
24
25 C == GLobal variables ==
26 #include "SIZE.h"
27 #include "DYNVARS.h"
28 #include "EEPARAMS.h"
29 #include "PARAMS.h"
30 #include "GRID.h"
31 #include "FFIELDS.h"
32 #include "KPP.h"
33 #include "KPP_PARAMS.h"
34
35 C == Routine arguments ==
36 C iMin,iMax,jMin, - Range of points for which calculation
37 C jMax,bi,bj,k,km1 results will be set.
38 C KappaRS - Total diffusion in vertical for S.
39 C df - Diffusive flux component work array.
40 INTEGER iMin,iMax,jMin,jMax,bi,bj,k,km1
41 _RL KappaRS(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
42 _RL df (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
43
44 #if defined(ALLOW_KPP) && defined(KPP_GHAT)
45
46 C == Local variables ==
47 C I, J - Loop counters
48 INTEGER I, J
49
50 DO j=jMin,jMax
51 DO i=iMin,iMax
52 df(i,j) = df(i,j) - _rA(i,j,bi,bj) *
53 & KappaRS(i,j,k) * KPPghat(i,j,km1,bi,bj) * drF(1) *
54 & SurfaceTendencyS(i,j,bi,bj)
55 ENDDO
56 ENDDO
57
58 #endif /* ALLOW_KPP and KPP_GHAT */
59
60 RETURN
61 END
  ViewVC Help
Powered by ViewVC 1.1.22