1 |
C $Header: /u/gcmpack/models/MITgcmUV/pkg/obcs/Attic/orlanski_init.F,v 1.1.2.1 2001/01/30 21:03:00 adcroft Exp $ |
2 |
C $Name: branch-atmos-merge-freeze $ |
3 |
|
4 |
#include "OBCS_OPTIONS.h" |
5 |
|
6 |
SUBROUTINE ORLANSKI_INIT(bi, bj, myThid ) |
7 |
C /==========================================================\ |
8 |
C | SUBROUTINE ORLANSKI_INIT | |
9 |
C | o Initialise storage and phase speed arrays | |
10 |
C | for Orlanski radiation condition | |
11 |
C |==========================================================| |
12 |
C | | |
13 |
C \==========================================================/ |
14 |
IMPLICIT NONE |
15 |
|
16 |
C === Global variables === |
17 |
#include "SIZE.h" |
18 |
#include "EEPARAMS.h" |
19 |
#include "PARAMS.h" |
20 |
#ifdef ALLOW_ORLANSKI |
21 |
#include "ORLANSKI.h" |
22 |
#endif /* ALLOW_ORLANSKI */ |
23 |
|
24 |
C == Routine arguments == |
25 |
INTEGER bi, bj, myThid |
26 |
|
27 |
#ifdef ALLOW_ORLANSKI |
28 |
|
29 |
C == Local variables == |
30 |
INTEGER I, J, K |
31 |
|
32 |
DO K=1,Nr |
33 |
DO I=1-Olx,sNx+Olx |
34 |
UN_STORE_1(I,K,bi,bj)=0. |
35 |
VN_STORE_1(I,K,bi,bj)=0. |
36 |
TN_STORE_1(I,K,bi,bj)=0. |
37 |
SN_STORE_1(I,K,bi,bj)=0. |
38 |
UN_STORE_2(I,K,bi,bj)=0. |
39 |
VN_STORE_2(I,K,bi,bj)=0. |
40 |
TN_STORE_2(I,K,bi,bj)=0. |
41 |
SN_STORE_2(I,K,bi,bj)=0. |
42 |
UN_STORE_3(I,K,bi,bj)=0. |
43 |
VN_STORE_3(I,K,bi,bj)=0. |
44 |
TN_STORE_3(I,K,bi,bj)=0. |
45 |
SN_STORE_3(I,K,bi,bj)=0. |
46 |
UN_STORE_4(I,K,bi,bj)=0. |
47 |
VN_STORE_4(I,K,bi,bj)=0. |
48 |
TN_STORE_4(I,K,bi,bj)=0. |
49 |
SN_STORE_4(I,K,bi,bj)=0. |
50 |
US_STORE_1(I,K,bi,bj)=0. |
51 |
VS_STORE_1(I,K,bi,bj)=0. |
52 |
TS_STORE_1(I,K,bi,bj)=0. |
53 |
SS_STORE_1(I,K,bi,bj)=0. |
54 |
US_STORE_2(I,K,bi,bj)=0. |
55 |
VS_STORE_2(I,K,bi,bj)=0. |
56 |
TS_STORE_2(I,K,bi,bj)=0. |
57 |
SS_STORE_2(I,K,bi,bj)=0. |
58 |
US_STORE_3(I,K,bi,bj)=0. |
59 |
VS_STORE_3(I,K,bi,bj)=0. |
60 |
TS_STORE_3(I,K,bi,bj)=0. |
61 |
SS_STORE_3(I,K,bi,bj)=0. |
62 |
US_STORE_4(I,K,bi,bj)=0. |
63 |
VS_STORE_4(I,K,bi,bj)=0. |
64 |
TS_STORE_4(I,K,bi,bj)=0. |
65 |
SS_STORE_4(I,K,bi,bj)=0. |
66 |
CVEL_UN(I,K,bi,bj)=0. |
67 |
CVEL_VN(I,K,bi,bj)=0. |
68 |
CVEL_TN(I,K,bi,bj)=0. |
69 |
CVEL_SN(I,K,bi,bj)=0. |
70 |
CVEL_US(I,K,bi,bj)=0. |
71 |
CVEL_VS(I,K,bi,bj)=0. |
72 |
CVEL_TS(I,K,bi,bj)=0. |
73 |
CVEL_SS(I,K,bi,bj)=0. |
74 |
#ifdef ALLOW_NONHYDROSTATIC |
75 |
WN_STORE_1(I,K,bi,bj)=0. |
76 |
WN_STORE_2(I,K,bi,bj)=0. |
77 |
WN_STORE_3(I,K,bi,bj)=0. |
78 |
WN_STORE_4(I,K,bi,bj)=0. |
79 |
WS_STORE_1(I,K,bi,bj)=0. |
80 |
WS_STORE_2(I,K,bi,bj)=0. |
81 |
WS_STORE_3(I,K,bi,bj)=0. |
82 |
WS_STORE_4(I,K,bi,bj)=0. |
83 |
CVEL_WN(I,K,bi,bj)=0. |
84 |
CVEL_WS(I,K,bi,bj)=0. |
85 |
#endif |
86 |
ENDDO |
87 |
DO J=1-Oly,sNy+Oly |
88 |
UE_STORE_1(J,K,bi,bj)=0. |
89 |
VE_STORE_1(J,K,bi,bj)=0. |
90 |
TE_STORE_1(J,K,bi,bj)=0. |
91 |
SE_STORE_1(J,K,bi,bj)=0. |
92 |
UE_STORE_2(J,K,bi,bj)=0. |
93 |
VE_STORE_2(J,K,bi,bj)=0. |
94 |
TE_STORE_2(J,K,bi,bj)=0. |
95 |
SE_STORE_2(J,K,bi,bj)=0. |
96 |
UE_STORE_3(J,K,bi,bj)=0. |
97 |
VE_STORE_3(J,K,bi,bj)=0. |
98 |
TE_STORE_3(J,K,bi,bj)=0. |
99 |
SE_STORE_3(J,K,bi,bj)=0. |
100 |
UE_STORE_4(J,K,bi,bj)=0. |
101 |
VE_STORE_4(J,K,bi,bj)=0. |
102 |
TE_STORE_4(J,K,bi,bj)=0. |
103 |
SE_STORE_4(J,K,bi,bj)=0. |
104 |
UW_STORE_1(J,K,bi,bj)=0. |
105 |
VW_STORE_1(J,K,bi,bj)=0. |
106 |
TW_STORE_1(J,K,bi,bj)=0. |
107 |
SW_STORE_1(J,K,bi,bj)=0. |
108 |
UW_STORE_2(J,K,bi,bj)=0. |
109 |
VW_STORE_2(J,K,bi,bj)=0. |
110 |
TW_STORE_2(J,K,bi,bj)=0. |
111 |
SW_STORE_2(J,K,bi,bj)=0. |
112 |
UW_STORE_3(J,K,bi,bj)=0. |
113 |
VW_STORE_3(J,K,bi,bj)=0. |
114 |
TW_STORE_3(J,K,bi,bj)=0. |
115 |
SW_STORE_3(J,K,bi,bj)=0. |
116 |
UW_STORE_4(J,K,bi,bj)=0. |
117 |
VW_STORE_4(J,K,bi,bj)=0. |
118 |
TW_STORE_4(J,K,bi,bj)=0. |
119 |
SW_STORE_4(J,K,bi,bj)=0. |
120 |
CVEL_UE(J,K,bi,bj)=0. |
121 |
CVEL_VE(J,K,bi,bj)=0. |
122 |
CVEL_TE(J,K,bi,bj)=0. |
123 |
CVEL_SE(J,K,bi,bj)=0. |
124 |
CVEL_UW(J,K,bi,bj)=0. |
125 |
CVEL_VW(J,K,bi,bj)=0. |
126 |
CVEL_TW(J,K,bi,bj)=0. |
127 |
CVEL_SW(J,K,bi,bj)=0. |
128 |
#ifdef ALLOW_NONHYDROSTATIC |
129 |
WE_STORE_1(J,K,bi,bj)=0. |
130 |
WE_STORE_2(J,K,bi,bj)=0. |
131 |
WE_STORE_3(J,K,bi,bj)=0. |
132 |
WE_STORE_4(J,K,bi,bj)=0. |
133 |
WW_STORE_1(J,K,bi,bj)=0. |
134 |
WW_STORE_2(J,K,bi,bj)=0. |
135 |
WW_STORE_3(J,K,bi,bj)=0. |
136 |
WW_STORE_4(J,K,bi,bj)=0. |
137 |
CVEL_WE(J,K,bi,bj)=0. |
138 |
CVEL_WW(J,K,bi,bj)=0. |
139 |
#endif |
140 |
ENDDO |
141 |
ENDDO |
142 |
|
143 |
#endif /* ALLOW_ORLANSKI */ |
144 |
RETURN |
145 |
END |