/[MITgcm]/MITgcm/pkg/obcs/obcs_init_variables.F
ViewVC logotype

Diff of /MITgcm/pkg/obcs/obcs_init_variables.F

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

revision 1.1 by adcroft, Tue Jan 30 21:03:00 2001 UTC revision 1.2 by adcroft, Fri Feb 2 21:36:30 2001 UTC
# Line 0  Line 1 
1    C $Header$
2    C $Name$
3    
4    #include "OBCS_OPTIONS.h"
5    
6          SUBROUTINE OBCS_INIT_VARIABLES( myThid )
7    C     /==========================================================\
8    C     | SUBROUTINE OBCS_INIT_VARIABLES                           |
9    C     | o Initialise OBCs variable data                          |
10    C     |==========================================================|
11    C     |                                                          |
12    C     \==========================================================/
13          IMPLICIT NONE
14    
15    C     === Global variables ===
16    #include "SIZE.h"
17    #include "EEPARAMS.h"
18    #include "PARAMS.h"
19    #include "DYNVARS.h"
20    #include "OBCS.h"
21    
22    C     == Routine arguments ==
23    C     myThid -  Number of this instance of INI_DEPTHS
24          INTEGER myThid
25    
26    #ifdef ALLOW_OBCS
27    
28    C     == Local variables ==
29          INTEGER bi, bj
30          INTEGER I, J, K
31    
32          DO bj = myByLo(myThid), myByHi(myThid)
33           DO bi = myBxLo(myThid), myBxHi(myThid)
34    
35            DO K=1,Nr
36             DO I=1-Olx,sNx+Olx
37              OBNu(I,K,bi,bj)=0.
38              OBNv(I,K,bi,bj)=0.
39              OBNt(I,K,bi,bj)=0.
40              OBSu(I,K,bi,bj)=0.
41              OBSv(I,K,bi,bj)=0.
42              OBSt(I,K,bi,bj)=0.
43    #ifdef ALLOW_NONHYDROSTATIC
44              OBNw(I,K,bi,bj)=0.
45              OBSw(I,K,bi,bj)=0.
46    #endif
47             ENDDO
48             DO J=1-Oly,sNy+Oly
49              OBEu(J,K,bi,bj)=0.
50              OBEv(J,K,bi,bj)=0.
51              OBEt(J,K,bi,bj)=0.
52              OBWu(J,K,bi,bj)=0.
53              OBWv(J,K,bi,bj)=0.
54              OBWt(J,K,bi,bj)=0.
55    #ifdef ALLOW_NONHYDROSTATIC
56              OBEw(J,K,bi,bj)=0.
57              OBWw(J,K,bi,bj)=0.
58    #endif
59             ENDDO
60            ENDDO
61    
62           ENDDO
63          ENDDO
64    
65    #ifdef ALLOW_ORLANSKI
66            IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.
67         &      useOrlanskiEast.OR.useOrlanskiWest) THEN
68              CALL ORLANSKI_INIT(bi, bj, myThid)
69            ENDIF
70    #endif /* ALLOW_ORLANSKI */
71    
72    C--   Apply OBCS values to initial conditions for consistancy
73          DO bj = myByLo(myThid), myByHi(myThid)
74           DO bi = myBxLo(myThid), myBxHi(myThid)
75            CALL OBCS_CALC( bi, bj, startTime,
76         &        uVel, vVel, wVel, theta, salt, myThid )
77            DO K=1,Nr
78             CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )
79             CALL OBCS_APPLY_TS( bi, bj, k, theta, salt, myThid )
80            ENDDO
81           ENDDO
82          ENDDO
83    
84    #endif /* ALLOW_OBCS */
85          RETURN
86          END
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22