/[MITgcm]/MITgcm/pkg/bbl/bbl_tendency_apply.F
ViewVC logotype

Contents of /MITgcm/pkg/bbl/bbl_tendency_apply.F

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

Revision 1.2 - (show annotations) (download)
Sun Aug 7 07:08:15 2011 UTC (12 years, 10 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63b, checkpoint63c, checkpoint64
Changes since 1.1: +1 -18 lines 
o adding package bbl (Bottom Boundary Layer)
  description in MITgcm/pkg/bbl/bbl_description.tex
  example/test experiment in MITgcm_contrib/bbl
1 C $Header: /u/gcmpack/MITgcm/pkg/bbl/bbl_tendency_apply.F,v 1.1 2011/08/06 03:13:22 dimitri Exp $
2 C $Name: $
3
4 #include "BBL_OPTIONS.h"
5 C-- File bbl_tendency_apply.F: Routines to apply BBL tendencies
6 C-- Contents
7 C-- o BBL_TENDENCY_APPLY_T
8 C-- o BBL_TENDENCY_APPLY_S
9
10 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
11 CBOP 0
12 C !ROUTINE: BBL_TENDENCY_APPLY_T
13
14 C !INTERFACE:
15 SUBROUTINE BBL_TENDENCY_APPLY_T(
16 & iMin, iMax, jMin, jMax,
17 & bi, bj, k, myTime, myThid )
18
19 C !DESCRIPTION:
20 C Add bbl tendency terms to T (theta) tendency. Routine works
21 C for one level at a time
22
23 C !USES:
24 IMPLICIT NONE
25 #include "SIZE.h"
26 #include "GRID.h"
27 #include "EEPARAMS.h"
28 #include "DYNVARS.h"
29 #include "BBL.h"
30
31 C !INPUT PARAMETERS:
32 INTEGER iMin, iMax, jMin, jMax, k, bi, bj, myThid
33 _RL myTime
34 CEOP
35
36 C !LOCAL VARIABLES:
37 INTEGER i, j
38
39 DO j=jMin,jMax
40 DO i=iMin,iMax
41 IF ( k .EQ. kLowC(i,j,bi,bj) )
42 & gT(i,j,k,bi,bj) = gT(i,j,k,bi,bj)
43 & + maskC(i,j,k,bi,bj) * bbl_TendTheta(i,j,bi,bj)
44 ENDDO
45 ENDDO
46
47 RETURN
48 END
49
50 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
51 CBOP 0
52 C !ROUTINE: BBL_TENDENCY_APPLY_S
53
54 C !INTERFACE:
55 SUBROUTINE BBL_TENDENCY_APPLY_S(
56 & iMin, iMax, jMin, jMax,
57 & bi, bj, k, myTime, myThid )
58
59 C !DESCRIPTION:
60 C Add bbl tendency terms to S tendency. Routine works for one
61 C level at a time
62
63 C !INPUT PARAMETERS:
64 IMPLICIT NONE
65 #include "SIZE.h"
66 #include "GRID.h"
67 #include "EEPARAMS.h"
68 #include "DYNVARS.h"
69 #include "BBL.h"
70
71 C !INPUT PARAMETERS:
72 INTEGER iMin, iMax, jMin, jMax, k, bi, bj, myThid
73 _RL myTime
74 CEOP
75
76 C !LOCAL VARIABLES:
77 INTEGER i, j
78
79 DO j=jMin,jMax
80 DO i=iMin,iMax
81 IF ( k .EQ. kLowC(i,j,bi,bj) )
82 & gS(i,j,k,bi,bj) = gS(i,j,k,bi,bj)
83 & + maskC(i,j,k,bi,bj) * bbl_TendSalt(i,j,bi,bj)
84 ENDDO
85 ENDDO
86
87 RETURN
88 END
  ViewVC Help
Powered by ViewVC 1.1.22