15 |
|
|
16 |
#ifdef ALLOW_DIAGNOSTICS |
#ifdef ALLOW_DIAGNOSTICS |
17 |
_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) |
18 |
|
_RL dummys(1-OLx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy) |
19 |
integer i,j,K,bi,bj |
integer i,j,K,bi,bj |
20 |
integer km1 |
integer km1 |
21 |
|
|
25 |
do bi = myBxLo(myThid), myBxHi(myThid) |
do bi = myBxLo(myThid), myBxHi(myThid) |
26 |
do j = 1,sNy |
do j = 1,sNy |
27 |
do i = 1,sNx |
do i = 1,sNx |
28 |
dummy(i,j,1,bi,bj) = etaN(i,j,bi,bj)*etaN(i,j,bi,bj) |
dummys(i,j,bi,bj) = etaN(i,j,bi,bj)*etaN(i,j,bi,bj) |
29 |
enddo |
enddo |
30 |
enddo |
enddo |
31 |
enddo |
enddo |
32 |
enddo |
enddo |
33 |
call fill_diagnostics(myThid,'ETANSQ ',0,1,0,1,1,dummy) |
call fill_diagnostics(myThid,'ETANSQ ',0,1,0,1,1,dummys) |
34 |
|
|
35 |
call fill_diagnostics(myThid,'UVEL ',0,Nr,0,1,1,uVel) |
call fill_diagnostics(myThid,'UVEL ',0,Nr,0,1,1,uVel) |
36 |
call fill_diagnostics(myThid,'VVEL ',0,Nr,0,1,1,vVel) |
call fill_diagnostics(myThid,'VVEL ',0,Nr,0,1,1,vVel) |