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

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

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

Revision 1.4 - (hide annotations) (download)
Tue May 29 14:01:39 2001 UTC (22 years, 11 months ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint43a-release1mods, checkpoint40pre3, checkpoint40pre1, checkpoint40pre7, checkpoint40pre6, checkpoint40pre9, checkpoint40pre8, release1_b1, checkpoint43, checkpoint40pre2, release1-branch_tutorials, checkpoint40pre4, release1-branch-end, checkpoint40pre5, release1_beta1, checkpoint42, checkpoint40, checkpoint41, checkpoint44, release1-branch_branchpoint
Branch point for: release1-branch, release1, ecco-branch, release1_coupled
Changes since 1.3: +5 -1 lines 
Merge from branch pre38:
 o essential mods for cubed sphere
 o debugged atmosphere, dynamcis + physics (aim)
 o new packages (mom_vecinv, mom_fluxform, ...)
1 adcroft 1.4 C $Header: /u/gcmpack/models/MITgcmUV/pkg/obcs/obcs_init_variables.F,v 1.3.2.1 2001/05/01 13:17:04 adcroft Exp $
2 adcroft 1.3 C $Name: $
3 adcroft 1.2
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 adcroft 1.4 OBNs(I,K,bi,bj)=0.
41 adcroft 1.2 OBSu(I,K,bi,bj)=0.
42     OBSv(I,K,bi,bj)=0.
43     OBSt(I,K,bi,bj)=0.
44 adcroft 1.4 OBSs(I,K,bi,bj)=0.
45 adcroft 1.2 #ifdef ALLOW_NONHYDROSTATIC
46     OBNw(I,K,bi,bj)=0.
47     OBSw(I,K,bi,bj)=0.
48     #endif
49     ENDDO
50     DO J=1-Oly,sNy+Oly
51     OBEu(J,K,bi,bj)=0.
52     OBEv(J,K,bi,bj)=0.
53     OBEt(J,K,bi,bj)=0.
54 adcroft 1.4 OBEs(J,K,bi,bj)=0.
55 adcroft 1.2 OBWu(J,K,bi,bj)=0.
56     OBWv(J,K,bi,bj)=0.
57     OBWt(J,K,bi,bj)=0.
58 adcroft 1.4 OBWs(J,K,bi,bj)=0.
59 adcroft 1.2 #ifdef ALLOW_NONHYDROSTATIC
60     OBEw(J,K,bi,bj)=0.
61     OBWw(J,K,bi,bj)=0.
62     #endif
63     ENDDO
64     ENDDO
65    
66     #ifdef ALLOW_ORLANSKI
67     IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.
68     & useOrlanskiEast.OR.useOrlanskiWest) THEN
69     CALL ORLANSKI_INIT(bi, bj, myThid)
70     ENDIF
71     #endif /* ALLOW_ORLANSKI */
72 adcroft 1.3
73     ENDDO
74     ENDDO
75 adcroft 1.2
76     C-- Apply OBCS values to initial conditions for consistancy
77     DO bj = myByLo(myThid), myByHi(myThid)
78     DO bi = myBxLo(myThid), myBxHi(myThid)
79     CALL OBCS_CALC( bi, bj, startTime,
80     & uVel, vVel, wVel, theta, salt, myThid )
81     DO K=1,Nr
82     CALL OBCS_APPLY_UV( bi, bj, k, uVel, vVel, myThid )
83     CALL OBCS_APPLY_TS( bi, bj, k, theta, salt, myThid )
84     ENDDO
85     ENDDO
86     ENDDO
87    
88     #endif /* ALLOW_OBCS */
89     RETURN
90     END
  ViewVC Help
Powered by ViewVC 1.1.22