13 |
_RL dummy(1-OLx:sNx+Olx,1-Oly:sNy+Oly,Nr,Nsx,Nsy) |
_RL dummy(1-OLx:sNx+Olx,1-Oly:sNy+Oly,Nr,Nsx,Nsy) |
14 |
integer i,j,K,bi,bj |
integer i,j,K,bi,bj |
15 |
|
|
|
#ifdef ALLOW_DIAGNOSTICS |
|
|
if(usediagnostics)then |
|
16 |
call fill_diagnostics(myThid,'ETAN ',0,1,0,1,1,etaN) |
call fill_diagnostics(myThid,'ETAN ',0,1,0,1,1,etaN) |
17 |
|
|
18 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
19 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
20 |
do j = 1-OLy,sNy+Oly |
do j = 1-OLy,sNy+Oly |
21 |
do i = 1-OLx,sNx+Olx |
do i = 1-OLx,sNx+Olx |
22 |
dummy(i,j,K,bi,bj) = etaN(i,j,bi,bj)*etaN(i,j,bi,bj) |
dummy(i,j,1,bi,bj) = etaN(i,j,bi,bj)*etaN(i,j,bi,bj) |
23 |
enddo |
enddo |
24 |
enddo |
enddo |
25 |
enddo |
enddo |
26 |
enddo |
enddo |
27 |
call fill_diagnostics(myThid,'ETANSQ ',0,1,0,1,1,etaN) |
call fill_diagnostics(myThid,'ETANSQ ',0,1,0,1,1,dummy) |
28 |
|
|
29 |
call fill_diagnostics(myThid,'UVEL ',0,Nr,0,1,1,uVel) |
call fill_diagnostics(myThid,'UVEL ',0,Nr,0,1,1,uVel) |
30 |
call fill_diagnostics(myThid,'VVEL ',0,Nr,0,1,1,vVel) |
call fill_diagnostics(myThid,'VVEL ',0,Nr,0,1,1,vVel) |
43 |
enddo |
enddo |
44 |
enddo |
enddo |
45 |
enddo |
enddo |
46 |
call fill_diagnostics(myThid,'UVELSQ ',0,Nr,0,1,1,uVel) |
call fill_diagnostics(myThid,'UVELSQ ',0,Nr,0,1,1,dummy) |
47 |
|
|
48 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
49 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
56 |
enddo |
enddo |
57 |
enddo |
enddo |
58 |
enddo |
enddo |
59 |
call fill_diagnostics(myThid,'VVELSQ ',0,Nr,0,1,1,vVel) |
call fill_diagnostics(myThid,'VVELSQ ',0,Nr,0,1,1,dummy) |
60 |
|
|
61 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
62 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
69 |
enddo |
enddo |
70 |
enddo |
enddo |
71 |
enddo |
enddo |
72 |
call fill_diagnostics(myThid,'WVELSQ ',0,Nr,0,1,1,wVel) |
call fill_diagnostics(myThid,'WVELSQ ',0,Nr,0,1,1,dummy) |
73 |
|
|
74 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
75 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
82 |
enddo |
enddo |
83 |
enddo |
enddo |
84 |
enddo |
enddo |
85 |
call fill_diagnostics(myThid,'THETASQ ',0,Nr,0,1,1,theta) |
call fill_diagnostics(myThid,'THETASQ ',0,Nr,0,1,1,dummy) |
86 |
|
|
87 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
88 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
95 |
enddo |
enddo |
96 |
enddo |
enddo |
97 |
enddo |
enddo |
98 |
call fill_diagnostics(myThid,'SALTSQ ',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'SALTSQ ',0,Nr,0,1,1,dummy) |
99 |
|
|
100 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
101 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
108 |
enddo |
enddo |
109 |
enddo |
enddo |
110 |
enddo |
enddo |
111 |
call fill_diagnostics(myThid,'UVELVVEL',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'UVELVVEL',0,Nr,0,1,1,dummy) |
112 |
|
|
113 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
114 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
121 |
enddo |
enddo |
122 |
enddo |
enddo |
123 |
enddo |
enddo |
124 |
call fill_diagnostics(myThid,'UVELTH ',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'UVELTH ',0,Nr,0,1,1,dummy) |
125 |
|
|
126 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
127 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
134 |
enddo |
enddo |
135 |
enddo |
enddo |
136 |
enddo |
enddo |
137 |
call fill_diagnostics(myThid,'VVELTH ',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'VVELTH ',0,Nr,0,1,1,dummy) |
138 |
|
|
139 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
140 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
147 |
enddo |
enddo |
148 |
enddo |
enddo |
149 |
enddo |
enddo |
150 |
call fill_diagnostics(myThid,'WVELTH ',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'WVELTH ',0,Nr,0,1,1,dummy) |
151 |
|
|
152 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
153 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
160 |
enddo |
enddo |
161 |
enddo |
enddo |
162 |
enddo |
enddo |
163 |
call fill_diagnostics(myThid,'UVELSLT ',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'UVELSLT ',0,Nr,0,1,1,dummy) |
164 |
|
|
165 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
166 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
173 |
enddo |
enddo |
174 |
enddo |
enddo |
175 |
enddo |
enddo |
176 |
call fill_diagnostics(myThid,'VVELSLT ',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'VVELSLT ',0,Nr,0,1,1,dummy) |
177 |
|
|
178 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
179 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
186 |
enddo |
enddo |
187 |
enddo |
enddo |
188 |
enddo |
enddo |
189 |
call fill_diagnostics(myThid,'WVELSLT ',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'WVELSLT ',0,Nr,0,1,1,dummy) |
190 |
|
|
191 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
192 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
212 |
enddo |
enddo |
213 |
enddo |
enddo |
214 |
enddo |
enddo |
215 |
call fill_diagnostics(myThid,'VVELMASS',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'VVELMASS',0,Nr,0,1,1,dummy) |
216 |
|
|
217 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
218 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
226 |
enddo |
enddo |
227 |
enddo |
enddo |
228 |
enddo |
enddo |
229 |
call fill_diagnostics(myThid,'UTHMASS ',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'UTHMASS ',0,Nr,0,1,1,dummy) |
230 |
|
|
231 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
232 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
240 |
enddo |
enddo |
241 |
enddo |
enddo |
242 |
enddo |
enddo |
243 |
call fill_diagnostics(myThid,'VTHMASS ',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'VTHMASS ',0,Nr,0,1,1,dummy) |
244 |
|
|
245 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
246 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
254 |
enddo |
enddo |
255 |
enddo |
enddo |
256 |
enddo |
enddo |
257 |
call fill_diagnostics(myThid,'USLTMASS',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'USLTMASS',0,Nr,0,1,1,dummy) |
258 |
|
|
259 |
do bj = myByLo(myThid), myByHi(myThid) |
do bj = myByLo(myThid), myByHi(myThid) |
260 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
268 |
enddo |
enddo |
269 |
enddo |
enddo |
270 |
enddo |
enddo |
271 |
call fill_diagnostics(myThid,'VSLTMASS',0,Nr,0,1,1,salt) |
call fill_diagnostics(myThid,'VSLTMASS',0,Nr,0,1,1,dummy) |
272 |
|
|
|
endif |
|
|
#endif |
|
273 |
|
|
274 |
return |
return |
275 |
end |
end |