pkg/diagnostics/diagnostics_fill_state.F

C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_fill_state.F,v 1.7 2004/06/14 21:54:47 jmc Exp $
C $Name:  $

#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"

       subroutine diagnostics_fill_state(myThid)
       implicit none
#include "SIZE.h"
#include "EEPARAMS.h"
#include "GRID.h"
#include "DYNVARS.h"

      integer myThid 

#ifdef ALLOW_DIAGNOSTICS
      _RL dummy(1-OLx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)
      _RL dummys(1-OLx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
      integer i,j,K,bi,bj
      integer km1

        call fill_diagnostics(myThid,'ETAN    ',0,1,0,1,1,etaN)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
          do j = 1,sNy
          do i = 1,sNx
           dummys(i,j,bi,bj) = etaN(i,j,bi,bj)*etaN(i,j,bi,bj)
          enddo
          enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'ETANSQ  ',0,1,0,1,1,dummys)

        call fill_diagnostics(myThid,'UVEL    ',0,Nr,0,1,1,uVel)
        call fill_diagnostics(myThid,'VVEL    ',0,Nr,0,1,1,vVel)
        call fill_diagnostics(myThid,'WVEL    ',0,Nr,0,1,1,wVel)
        call fill_diagnostics(myThid,'THETA   ',0,Nr,0,1,1,theta)
        call fill_diagnostics(myThid,'SALT    ',0,Nr,0,1,1,salt)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*uVel(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'UVELSQ  ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*vVel(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'VVELSQ  ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = wVel(i,j,K,bi,bj)*wVel(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'WVELSQ  ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = theta(i,j,K,bi,bj)*theta(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'THETASQ ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = salt(i,j,K,bi,bj)*salt(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'SALTSQ  ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*vVel(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'UVELVVEL',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*
     &         0.5*(theta(i,j,K,bi,bj)+theta(i-1,j,K,bi,bj))
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'UVELTH  ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*
     &         0.5*(theta(i,j,K,bi,bj)+theta(i,j-1,K,bi,bj))
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'VVELTH  ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
         do bi = myBxLo(myThid), myBxHi(myThid)
          do K=1,Nr
           km1 = MAX(k-1,1)
           do j = 1,sNy
            do i = 1,sNx
             dummy(i,j,K,bi,bj) = wVel(i,j,K,bi,bj)*0.5*
     &                 (theta(i,j,K,bi,bj)+theta(i,j,km1,bi,bj))
            enddo
           enddo
          enddo
         enddo
        enddo
        call fill_diagnostics(myThid,'WVELTH  ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*
     &         0.5*(salt(i,j,K,bi,bj)+salt(i-1,j,K,bi,bj))
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'UVELSLT ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*
     &         0.5*(salt(i,j,K,bi,bj)+salt(i,j-1,K,bi,bj))
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'VVELSLT ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
         do bi = myBxLo(myThid), myBxHi(myThid)
          do K=1,Nr
           km1 = MAX(k-1,1)
           do j = 1,sNy
            do i = 1,sNx
             dummy(i,j,K,bi,bj) = wVel(i,j,K,bi,bj)*0.5*
     &                 (salt(i,j,K,bi,bj)+salt(i,j,km1,bi,bj))
            enddo
           enddo
          enddo
         enddo
        enddo
        call fill_diagnostics(myThid,'WVELSLT ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*hFacW(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'UVELMASS',0,Nr,0,1,1,salt)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*hFacS(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'VVELMASS',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*
     .         0.5*(theta(i,j,K,bi,bj)+theta(i-1,j,K,bi,bj))
     .                                * hFacW(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'UTHMASS ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*
     .         0.5*(theta(i,j,K,bi,bj)+theta(i,j-1,K,bi,bj))
     .                                * hFacS(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'VTHMASS ',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*
     .         0.5*(salt(i,j,K,bi,bj)+salt(i-1,j,K,bi,bj))
     .                                * hFacW(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'USLTMASS',0,Nr,0,1,1,dummy)

        do bj = myByLo(myThid), myByHi(myThid)
        do bi = myBxLo(myThid), myBxHi(myThid)
         do K=1,Nr
          do j = 1,sNy
          do i = 1,sNx
           dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*
     .         0.5*(salt(i,j,K,bi,bj)+salt(i,j-1,K,bi,bj))
     .                                * hFacS(i,j,K,bi,bj)
          enddo
          enddo
         enddo
        enddo
        enddo
        call fill_diagnostics(myThid,'VSLTMASS',0,Nr,0,1,1,dummy)


#endif /* ALLOW_DIAGNOSTICS */

      return
      end
1	molod	1.8	C $Header: /u/gcmpack/MITgcm/pkg/diagnostics/diagnostics_fill_state.F,v 1.7 2004/06/14 21:54:47 jmc Exp $
2	edhill	1.6	C $Name: $
3
4	jmc	1.7	#include "PACKAGES_CONFIG.h"
5			#include "CPP_OPTIONS.h"
6
7	molod	1.1	subroutine diagnostics_fill_state(myThid)
8	molod	1.2	implicit none
9	molod	1.1	#include "SIZE.h"
10			#include "EEPARAMS.h"
11			#include "GRID.h"
12			#include "DYNVARS.h"
13
14			integer myThid
15	jmc	1.7
16			#ifdef ALLOW_DIAGNOSTICS
17			_RL dummy(1-OLx:sNx+Olx,1-Oly:sNy+Oly,Nr,nSx,nSy)
18	molod	1.8	_RL dummys(1-OLx:sNx+Olx,1-Oly:sNy+Oly,nSx,nSy)
19	molod	1.1	integer i,j,K,bi,bj
20	jmc	1.7	integer km1
21	molod	1.1
22			call fill_diagnostics(myThid,'ETAN ',0,1,0,1,1,etaN)
23
24			do bj = myByLo(myThid), myByHi(myThid)
25			do bi = myBxLo(myThid), myBxHi(myThid)
26	molod	1.4	do j = 1,sNy
27			do i = 1,sNx
28	molod	1.8	dummys(i,j,bi,bj) = etaN(i,j,bi,bj)*etaN(i,j,bi,bj)
29	molod	1.1	enddo
30			enddo
31			enddo
32			enddo
33	molod	1.8	call fill_diagnostics(myThid,'ETANSQ ',0,1,0,1,1,dummys)
34	molod	1.1
35			call fill_diagnostics(myThid,'UVEL ',0,Nr,0,1,1,uVel)
36			call fill_diagnostics(myThid,'VVEL ',0,Nr,0,1,1,vVel)
37			call fill_diagnostics(myThid,'WVEL ',0,Nr,0,1,1,wVel)
38			call fill_diagnostics(myThid,'THETA ',0,Nr,0,1,1,theta)
39			call fill_diagnostics(myThid,'SALT ',0,Nr,0,1,1,salt)
40
41			do bj = myByLo(myThid), myByHi(myThid)
42			do bi = myBxLo(myThid), myBxHi(myThid)
43			do K=1,Nr
44	molod	1.4	do j = 1,sNy
45			do i = 1,sNx
46	molod	1.1	dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*uVel(i,j,K,bi,bj)
47			enddo
48			enddo
49			enddo
50			enddo
51			enddo
52	molod	1.3	call fill_diagnostics(myThid,'UVELSQ ',0,Nr,0,1,1,dummy)
53	molod	1.1
54			do bj = myByLo(myThid), myByHi(myThid)
55			do bi = myBxLo(myThid), myBxHi(myThid)
56			do K=1,Nr
57	molod	1.4	do j = 1,sNy
58			do i = 1,sNx
59	molod	1.1	dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*vVel(i,j,K,bi,bj)
60			enddo
61			enddo
62			enddo
63			enddo
64			enddo
65	molod	1.3	call fill_diagnostics(myThid,'VVELSQ ',0,Nr,0,1,1,dummy)
66	molod	1.1
67			do bj = myByLo(myThid), myByHi(myThid)
68			do bi = myBxLo(myThid), myBxHi(myThid)
69			do K=1,Nr
70	molod	1.4	do j = 1,sNy
71			do i = 1,sNx
72	molod	1.1	dummy(i,j,K,bi,bj) = wVel(i,j,K,bi,bj)*wVel(i,j,K,bi,bj)
73			enddo
74			enddo
75			enddo
76			enddo
77			enddo
78	molod	1.3	call fill_diagnostics(myThid,'WVELSQ ',0,Nr,0,1,1,dummy)
79	molod	1.1
80			do bj = myByLo(myThid), myByHi(myThid)
81			do bi = myBxLo(myThid), myBxHi(myThid)
82			do K=1,Nr
83	molod	1.4	do j = 1,sNy
84			do i = 1,sNx
85	molod	1.1	dummy(i,j,K,bi,bj) = theta(i,j,K,bi,bj)*theta(i,j,K,bi,bj)
86			enddo
87			enddo
88			enddo
89			enddo
90			enddo
91	molod	1.3	call fill_diagnostics(myThid,'THETASQ ',0,Nr,0,1,1,dummy)
92	molod	1.1
93			do bj = myByLo(myThid), myByHi(myThid)
94			do bi = myBxLo(myThid), myBxHi(myThid)
95			do K=1,Nr
96	molod	1.4	do j = 1,sNy
97			do i = 1,sNx
98	molod	1.1	dummy(i,j,K,bi,bj) = salt(i,j,K,bi,bj)*salt(i,j,K,bi,bj)
99			enddo
100			enddo
101			enddo
102			enddo
103			enddo
104	molod	1.3	call fill_diagnostics(myThid,'SALTSQ ',0,Nr,0,1,1,dummy)
105	molod	1.1
106			do bj = myByLo(myThid), myByHi(myThid)
107			do bi = myBxLo(myThid), myBxHi(myThid)
108			do K=1,Nr
109	molod	1.4	do j = 1,sNy
110			do i = 1,sNx
111	molod	1.1	dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*vVel(i,j,K,bi,bj)
112			enddo
113			enddo
114			enddo
115			enddo
116			enddo
117	molod	1.3	call fill_diagnostics(myThid,'UVELVVEL',0,Nr,0,1,1,dummy)
118	molod	1.1
119			do bj = myByLo(myThid), myByHi(myThid)
120			do bi = myBxLo(myThid), myBxHi(myThid)
121			do K=1,Nr
122	molod	1.4	do j = 1,sNy
123			do i = 1,sNx
124	jmc	1.7	dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*
125			& 0.5*(theta(i,j,K,bi,bj)+theta(i-1,j,K,bi,bj))
126	molod	1.1	enddo
127			enddo
128			enddo
129			enddo
130			enddo
131	molod	1.3	call fill_diagnostics(myThid,'UVELTH ',0,Nr,0,1,1,dummy)
132	molod	1.1
133			do bj = myByLo(myThid), myByHi(myThid)
134			do bi = myBxLo(myThid), myBxHi(myThid)
135			do K=1,Nr
136	molod	1.4	do j = 1,sNy
137			do i = 1,sNx
138	jmc	1.7	dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*
139			& 0.5*(theta(i,j,K,bi,bj)+theta(i,j-1,K,bi,bj))
140	molod	1.1	enddo
141			enddo
142			enddo
143			enddo
144			enddo
145	molod	1.3	call fill_diagnostics(myThid,'VVELTH ',0,Nr,0,1,1,dummy)
146	molod	1.1
147			do bj = myByLo(myThid), myByHi(myThid)
148	jmc	1.7	do bi = myBxLo(myThid), myBxHi(myThid)
149			do K=1,Nr
150			km1 = MAX(k-1,1)
151			do j = 1,sNy
152			do i = 1,sNx
153			dummy(i,j,K,bi,bj) = wVel(i,j,K,bi,bj)0.5
154			& (theta(i,j,K,bi,bj)+theta(i,j,km1,bi,bj))
155			enddo
156			enddo
157	molod	1.4	enddo
158			enddo
159			enddo
160	molod	1.3	call fill_diagnostics(myThid,'WVELTH ',0,Nr,0,1,1,dummy)
161	molod	1.1
162			do bj = myByLo(myThid), myByHi(myThid)
163			do bi = myBxLo(myThid), myBxHi(myThid)
164			do K=1,Nr
165	molod	1.4	do j = 1,sNy
166			do i = 1,sNx
167	jmc	1.7	dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*
168			& 0.5*(salt(i,j,K,bi,bj)+salt(i-1,j,K,bi,bj))
169	molod	1.1	enddo
170			enddo
171			enddo
172			enddo
173			enddo
174	molod	1.3	call fill_diagnostics(myThid,'UVELSLT ',0,Nr,0,1,1,dummy)
175	molod	1.1
176			do bj = myByLo(myThid), myByHi(myThid)
177			do bi = myBxLo(myThid), myBxHi(myThid)
178			do K=1,Nr
179	molod	1.4	do j = 1,sNy
180			do i = 1,sNx
181	jmc	1.7	dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*
182			& 0.5*(salt(i,j,K,bi,bj)+salt(i,j-1,K,bi,bj))
183	molod	1.1	enddo
184			enddo
185			enddo
186			enddo
187			enddo
188	molod	1.3	call fill_diagnostics(myThid,'VVELSLT ',0,Nr,0,1,1,dummy)
189	molod	1.1
190			do bj = myByLo(myThid), myByHi(myThid)
191	jmc	1.7	do bi = myBxLo(myThid), myBxHi(myThid)
192			do K=1,Nr
193			km1 = MAX(k-1,1)
194			do j = 1,sNy
195			do i = 1,sNx
196			dummy(i,j,K,bi,bj) = wVel(i,j,K,bi,bj)0.5
197			& (salt(i,j,K,bi,bj)+salt(i,j,km1,bi,bj))
198			enddo
199			enddo
200	molod	1.1	enddo
201			enddo
202			enddo
203	molod	1.3	call fill_diagnostics(myThid,'WVELSLT ',0,Nr,0,1,1,dummy)
204	molod	1.1
205			do bj = myByLo(myThid), myByHi(myThid)
206			do bi = myBxLo(myThid), myBxHi(myThid)
207			do K=1,Nr
208	molod	1.4	do j = 1,sNy
209			do i = 1,sNx
210	molod	1.1	dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*hFacW(i,j,K,bi,bj)
211			enddo
212			enddo
213			enddo
214			enddo
215			enddo
216			call fill_diagnostics(myThid,'UVELMASS',0,Nr,0,1,1,salt)
217
218			do bj = myByLo(myThid), myByHi(myThid)
219			do bi = myBxLo(myThid), myBxHi(myThid)
220			do K=1,Nr
221	molod	1.4	do j = 1,sNy
222			do i = 1,sNx
223	molod	1.1	dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*hFacS(i,j,K,bi,bj)
224			enddo
225			enddo
226			enddo
227			enddo
228			enddo
229	molod	1.3	call fill_diagnostics(myThid,'VVELMASS',0,Nr,0,1,1,dummy)
230	molod	1.1
231			do bj = myByLo(myThid), myByHi(myThid)
232			do bi = myBxLo(myThid), myBxHi(myThid)
233			do K=1,Nr
234	molod	1.4	do j = 1,sNy
235			do i = 1,sNx
236			dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*
237			. 0.5*(theta(i,j,K,bi,bj)+theta(i-1,j,K,bi,bj))
238	molod	1.1	. * hFacW(i,j,K,bi,bj)
239			enddo
240			enddo
241			enddo
242			enddo
243			enddo
244	molod	1.3	call fill_diagnostics(myThid,'UTHMASS ',0,Nr,0,1,1,dummy)
245	molod	1.1
246			do bj = myByLo(myThid), myByHi(myThid)
247			do bi = myBxLo(myThid), myBxHi(myThid)
248			do K=1,Nr
249	molod	1.4	do j = 1,sNy
250			do i = 1,sNx
251			dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*
252			. 0.5*(theta(i,j,K,bi,bj)+theta(i,j-1,K,bi,bj))
253	jmc	1.7	. * hFacS(i,j,K,bi,bj)
254	molod	1.1	enddo
255			enddo
256			enddo
257			enddo
258			enddo
259	molod	1.3	call fill_diagnostics(myThid,'VTHMASS ',0,Nr,0,1,1,dummy)
260	molod	1.1
261			do bj = myByLo(myThid), myByHi(myThid)
262			do bi = myBxLo(myThid), myBxHi(myThid)
263			do K=1,Nr
264	molod	1.4	do j = 1,sNy
265			do i = 1,sNx
266			dummy(i,j,K,bi,bj) = uVel(i,j,K,bi,bj)*
267			. 0.5*(salt(i,j,K,bi,bj)+salt(i-1,j,K,bi,bj))
268	molod	1.1	. * hFacW(i,j,K,bi,bj)
269			enddo
270			enddo
271			enddo
272			enddo
273			enddo
274	molod	1.3	call fill_diagnostics(myThid,'USLTMASS',0,Nr,0,1,1,dummy)
275	molod	1.1
276			do bj = myByLo(myThid), myByHi(myThid)
277			do bi = myBxLo(myThid), myBxHi(myThid)
278			do K=1,Nr
279	molod	1.4	do j = 1,sNy
280			do i = 1,sNx
281			dummy(i,j,K,bi,bj) = vVel(i,j,K,bi,bj)*
282			. 0.5*(salt(i,j,K,bi,bj)+salt(i,j-1,K,bi,bj))
283	jmc	1.7	. * hFacS(i,j,K,bi,bj)
284	molod	1.1	enddo
285			enddo
286			enddo
287			enddo
288			enddo
289	molod	1.3	call fill_diagnostics(myThid,'VSLTMASS',0,Nr,0,1,1,dummy)
290	molod	1.1
291	jmc	1.7
292			#endif /* ALLOW_DIAGNOSTICS */
293	molod	1.1
294			return
295			end