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

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

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


Revision 1.1 - (hide annotations) (download)
Thu Sep 2 09:11:54 2004 UTC (20 years, 10 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 mlosch 1.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