1 |
dimitri |
1.3 |
C $Header: /u/gcmpack/MITgcm/pkg/bbl/bbl_init_varia.F,v 1.2 2012/04/03 00:22:16 jmc Exp $ |
2 |
dimitri |
1.1 |
C $Name: $ |
3 |
|
|
|
4 |
|
|
#include "BBL_OPTIONS.h" |
5 |
|
|
|
6 |
|
|
CBOP |
7 |
|
|
C !ROUTINE: BBL_INIT_VARIA |
8 |
|
|
|
9 |
|
|
C !INTERFACE: |
10 |
|
|
SUBROUTINE BBL_INIT_VARIA( myThid ) |
11 |
|
|
|
12 |
|
|
C !DESCRIPTION: |
13 |
|
|
C Initialize BBL variables |
14 |
|
|
|
15 |
|
|
C !USES: |
16 |
|
|
IMPLICIT NONE |
17 |
|
|
#include "SIZE.h" |
18 |
|
|
#include "EEPARAMS.h" |
19 |
|
|
#include "PARAMS.h" |
20 |
|
|
#include "GRID.h" |
21 |
|
|
#include "DYNVARS.h" |
22 |
|
|
#include "BBL.h" |
23 |
|
|
|
24 |
|
|
C !INPUT/OUTPUT PARAMETERS: |
25 |
|
|
C myThid :: my Thread Id number |
26 |
|
|
INTEGER myThid |
27 |
|
|
CEOP |
28 |
|
|
|
29 |
|
|
#ifdef ALLOW_BBL |
30 |
|
|
|
31 |
|
|
C !LOCAL VARIABLES: |
32 |
jmc |
1.2 |
C i,j,bi,bj :: loop indices |
33 |
dimitri |
1.1 |
C kBot :: k index of bottommost wet grid |
34 |
jmc |
1.2 |
INTEGER i,j,bi,bj,kBot |
35 |
dimitri |
1.1 |
|
36 |
|
|
C Initialize arrays in common blocks. |
37 |
|
|
DO bj = myByLo(myThid), myByHi(myThid) |
38 |
|
|
DO bi = myBxLo(myThid), myBxHi(myThid) |
39 |
jmc |
1.2 |
DO j=1-OLy,sNy+OLy |
40 |
|
|
DO i=1-OLx,sNx+OLx |
41 |
dimitri |
1.1 |
kBot = max(1,kLowC(i,j,bi,bj)) |
42 |
|
|
bbl_theta(i,j,bi,bj) = theta(i,j,kBot,bi,bj) |
43 |
|
|
bbl_salt (i,j,bi,bj) = salt (i,j,kBot,bi,bj) |
44 |
dimitri |
1.3 |
bbl_eta(i,j,bi,bj) = 0. _d 0 |
45 |
dimitri |
1.1 |
bbl_TendTheta(i,j,bi,bj) = 0. _d 0 |
46 |
|
|
bbl_TendSalt (i,j,bi,bj) = 0. _d 0 |
47 |
|
|
ENDDO |
48 |
|
|
ENDDO |
49 |
|
|
ENDDO |
50 |
|
|
ENDDO |
51 |
|
|
|
52 |
|
|
IF ( startTime.EQ.baseTime .AND. nIter0.EQ.0 |
53 |
|
|
& .AND. pickupSuff.EQ.' ' ) THEN |
54 |
|
|
|
55 |
|
|
C Read initialization files, if available. |
56 |
|
|
IF ( bbl_thetaFile .NE. ' ' ) THEN |
57 |
|
|
CALL READ_FLD_XY_RL( bbl_thetaFile, ' ', |
58 |
|
|
& bbl_theta, 0, myThid ) |
59 |
|
|
CALL EXCH_XY_RL( bbl_theta, myThid ) |
60 |
|
|
ENDIF |
61 |
|
|
IF ( bbl_saltFile .NE. ' ' ) THEN |
62 |
|
|
CALL READ_FLD_XY_RL( bbl_saltFile, ' ', |
63 |
|
|
& bbl_salt, 0, myThid ) |
64 |
|
|
CALL EXCH_XY_RL( bbl_salt, myThid ) |
65 |
|
|
ENDIF |
66 |
|
|
IF ( bbl_etaFile .NE. ' ' ) THEN |
67 |
|
|
CALL READ_FLD_XY_RL( bbl_etaFile, ' ', |
68 |
|
|
& bbl_eta, 0, myThid ) |
69 |
|
|
CALL EXCH_XY_RL( bbl_eta, myThid ) |
70 |
|
|
ENDIF |
71 |
|
|
|
72 |
|
|
ELSE |
73 |
|
|
|
74 |
|
|
C Restart from a pickup. |
75 |
|
|
CALL BBL_READ_PICKUP( nIter0, myThid ) |
76 |
|
|
|
77 |
|
|
C End start-from-iter-zero if/else block. |
78 |
|
|
ENDIF |
79 |
|
|
|
80 |
|
|
#endif /* ALLOW_BBL */ |
81 |
|
|
|
82 |
|
|
RETURN |
83 |
|
|
END |