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

Annotation of /MITgcm/pkg/kpp/kpp_transport_t.F

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

Revision 1.7 - (hide annotations) (download)
Tue May 29 14:01:38 2001 UTC (23 years ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint44e_post, checkpoint44f_post, checkpoint43a-release1mods, checkpoint40pre3, checkpoint40pre1, checkpoint40pre7, checkpoint40pre6, checkpoint40pre9, checkpoint40pre8, chkpt44d_post, checkpoint44e_pre, release1_b1, checkpoint43, release1_chkpt44d_post, checkpoint40pre2, release1-branch_tutorials, chkpt44a_post, checkpoint40pre4, chkpt44c_pre, ecco_c44_e18, ecco_c44_e17, ecco_c44_e16, checkpoint44g_post, release1-branch-end, checkpoint44b_post, checkpoint40pre5, chkpt44a_pre, ecco-branch-mod1, ecco-branch-mod2, ecco-branch-mod3, ecco-branch-mod4, ecco-branch-mod5, release1_beta1, checkpoint44b_pre, checkpoint42, checkpoint40, checkpoint41, checkpoint44, chkpt44c_post, checkpoint44f_pre, release1-branch_branchpoint
Branch point for: release1_final, release1-branch, release1, ecco-branch, release1_coupled
Changes since 1.6: +4 -6 lines 
Merge from branch pre38:
 o essential mods for cubed sphere
 o debugged atmosphere, dynamcis + physics (aim)
 o new packages (mom_vecinv, mom_fluxform, ...)
1 adcroft 1.7 C $Header: /u/gcmpack/models/MITgcmUV/pkg/kpp/kpp_transport_t.F,v 1.6.2.1 2001/03/30 23:28:32 jmc Exp $
2     C $Name: $
3 adcroft 1.1
4     #include "KPP_OPTIONS.h"
5    
6     subroutine KPP_TRANSPORT_T (
7     I iMin,iMax,jMin,jMax,bi,bj,k,km1,
8 adcroft 1.7 I KappaRT,
9 adcroft 1.1 U df )
10     C /==========================================================\
11     C | o SUBROUTINE KPP_TRANSPORT_T |
12     C | Add non local KPP transport term (ghat) to diffusive |
13     C | temperature flux. |
14     C |==========================================================|
15     C | The nonlocal transport term is nonzero only for scalars |
16     C | in unstable (convective) forcing conditions. |
17 heimbach 1.2 C | Note: [SurfaceTendencyT + Qsw * (1-KPPfrac)] |
18 adcroft 1.1 C | is the total heat flux penetrating the mixed layer |
19     C | from the surface in deg C / s referenced to the |
20     C | top layer. |
21     C | Note: KappaRT(i,j,k) is defined at the top of grid cell |
22     C | k while KPPghat(i,j,k,bi,bj) is defined at the |
23     C | bottom of grid cell k. |
24     C | For K=1, KappaRT(i,j,k) = 0 which insures no flux |
25     C | through the surface. |
26     C \==========================================================/
27     IMPLICIT NONE
28    
29     C == GLobal variables ==
30     #include "SIZE.h"
31     #include "DYNVARS.h"
32     #include "EEPARAMS.h"
33     #include "PARAMS.h"
34     #include "GRID.h"
35     #include "FFIELDS.h"
36     #include "KPP.h"
37     #include "KPP_PARAMS.h"
38    
39     C == Routine arguments ==
40     C iMin,iMax,jMin, - Range of points for which calculation
41     C jMax,bi,bj,k,km1 results will be set.
42     C KappaRT - Total diffusion in vertical for T.
43     C df - Diffusive flux component work array.
44     INTEGER iMin,iMax,jMin,jMax,bi,bj,k,km1
45     _RL KappaRT(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
46     _RL df (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
47    
48 heimbach 1.2 #if defined(ALLOW_KPP) && defined(KPP_GHAT)
49 adcroft 1.1
50     C == Local variables ==
51     C I, J - Loop counters
52     INTEGER I, J
53    
54 heimbach 1.2 DO j=jMin,jMax
55     DO i=iMin,iMax
56 adcroft 1.7 df(i,j) = df(i,j) - _rA(i,j,bi,bj) *
57 heimbach 1.2 & KappaRT(i,j,k) * KPPghat(i,j,km1,bi,bj) * delZ(1) *
58 adcroft 1.5 & ( SurfaceTendencyT(i,j,bi,bj) -
59     & Qsw(i,j,bi,bj)
60     & *recip_Cp*recip_rhoNil*recip_dRf(1)
61     & * (1.-KPPfrac(i,j,bi,bj))
62 heimbach 1.2 & )
63 adcroft 1.1 ENDDO
64 heimbach 1.2 ENDDO
65 adcroft 1.1
66 heimbach 1.2 #endif /* ALLOW_KPP and KPP_GHAT */
67 adcroft 1.1
68     RETURN
69     END
  ViewVC Help
Powered by ViewVC 1.1.22