/[MITgcm]/MITgcm/pkg/pp81/pp81_init.F
ViewVC logotype

Contents of /MITgcm/pkg/pp81/pp81_init.F

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


Revision 1.1 - (show annotations) (download)
Thu Sep 2 09:11:54 2004 UTC (19 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint55c_post, checkpoint55d_pre, checkpoint54e_post, checkpoint55b_post, checkpoint55, checkpoint55a_post, checkpoint55g_post, checkpoint55f_post, checkpoint54f_post, checkpoint55e_post, checkpoint55d_post
o add two new packages
  - pp81 (Packanowski and Philander, 1981), Richardson number and
    stratification dependent mixing
  - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme

1 C $Header: $
2 C $Name: $
3
4 #include "PP81_OPTIONS.h"
5
6 SUBROUTINE PP81_INIT( myThid )
7 C /==========================================================\
8 C | SUBROUTINE PP81_INIT |
9 C | o Routine to initialize PP81 parameters and variables. |
10 C |==========================================================|
11 C | Initialize PP81 parameters and variables. |
12 C \==========================================================/
13 IMPLICIT NONE
14
15 C === Global variables ===
16 #include "SIZE.h"
17 #include "EEPARAMS.h"
18 #include "PARAMS.h"
19 #include "GRID.h"
20 #include "PP81.h"
21
22 C === Routine arguments ===
23 C myThid - Number of this instance of PP81_INIT
24 INTEGER myThid
25
26 #ifdef ALLOW_PP81
27
28 C === Local variables ===
29 C i,j,k,bi,bj - Loop counters
30
31 INTEGER i, j, k, bi, bj
32
33 C-----------------------------------------------------------------------
34 C Initialize PP variables PPviscAr, PPdiffKrT
35 C-----------------------------------------------------------------------
36
37 do bj = myByLo(myThid), myByHi(myThid)
38 do bi = myBxLo(myThid), myBxHi(myThid)
39 do k = 1, Nr
40 do j = 1-OLy, sNy+OLy
41 do i = 1-OLx, sNx+OLx
42 PPviscAr (i,j,k,bi,bj) = viscAr
43 PPdiffKr (i,j,k,bi,bj) = diffKrT
44 end do
45 end do
46 end do
47 end do
48 end do
49
50 #endif /* ALLOW_PP81 */
51
52 return
53 end

  ViewVC Help
Powered by ViewVC 1.1.22